Abstract
Albeit occurring at zero temperature, quantum critical phenomena have a huge impact on the finite-temperature phase diagram of strongly correlated systems, giving experimental access to their observation. Indeed, the existence of a gapless, zero-temperature quantum critical point induces the existence of an extended region in parameter spaceâthe quantum critical fan (QCF)âcharacterized by power-law temperature dependences of all observables. Identifying experimentally the QCF and its crossovers to other regimes (renormalized classical, quantum disordered) remains nonetheless challenging. Focusing on paradigmatic models of quantum phase transitions, here we show that quantum correlationsâcaptured by the quantum variance of the order parameterâexhibit the temperature scaling associated with the QCF over a parameter region much broader than that revealed by ordinary correlations. The link existing between the quantum variance and the dynamical susceptibility paves the way to an experimental reconstruction of the QCF using spectroscopic techniques.
Similar content being viewed by others
Introduction
Quantum critical (QC) phenomena1,2,3,4,5 represent possibly the most dramatic manifestation of quantum mechanics at the macroscopic scale. Their typical setting involves an Hamiltonian \({\cal H} = {\cal H}_0 + gV\) in which the competition between the two noncommuting terms \({\cal H}_0\) and V, controlled by the parameter g, induces a macroscopic rearrangement of the ground state at a critical value gc, accompanied by the appearance of critical quantum fluctuations of collective observables at all length scales. This complex behavior of correlation and entanglement properties emerges from extensive theoretical work based on exactly solvable microscopic models4, quantum field theory2, as well as numerical studies6. Experiments generally do not have access to ground-state physics, but it was realized at an early stage7 that zero-T quantum critical points affect a sizable portion of the finite-T phase diagram, by inducing the presence of a so-called QC regime, in which the behavior of thermodynamic potentials and of response functions is controlled by the QC point (QCP). Indeed, observables in the QC regime are expected to exhibit thermal QC scaling, namely a power-law dependence on temperature with exponents descending from the critical exponents at the QCP. Strikingly, the QC regime is expected to be wider in parameter space at higher temperatures: as sketched in Fig. 1, it acts as a magnifying lens for the QCP. Even more strikingly, the finite-T QC regime ignores completely the physics of the Tâ=â0 and low-T phases at gââ âgc8, which are generally an ordered phase with a classical analog (for, say, gâ<âgc); and a gapped quantum disordered phase (for gâ>âgc). This implies that, if the temperature is lowered from a point at gââ âgc in the QC regime, a crossover must occur toward a thermodynamic regime which is instead controlled by the presence of long-range order in the ground stateâthe so-called renormalized classical (RC) regime for gâ<âgc; or by the presence of a gap above a disordered ground stateâthe QD regime for gâ>âgc. This is all the more striking, as it shows that a strictly quantum Tâ=â0 phenomenon (the QCP), governed by divergent quantum fluctuations, can have consequences on the phase diagram at temperatures T which are higher than those necessary to melt long-range order via a classical thermal transition.
Many exciting platforms for the exploration of QC phenomena can be found across the physical spectrum4,5,9,10,11,12. But can one reconstruct the QC regime quantitatively? The special scaling properties of the thermodynamic potentials and the dynamical response functions at gâ=âgc and finite T (along the so-called QC trajectory) have been observed in several systems, including magnetic insulators13,14,15,16,17 and heavy-fermion compounds9,18; but the persistence of the QC regime away from gc, and its crossover into the competing low-T regimes, are almost uniquely observed via transport properties in heavy-fermion materials9âthe so-called strange-metal phase in cuprate superconductors is also interpreted as an extended QC regime19 associated with a putative QCP20,21. Hence it is fair to say that the quantitative extent of the QC regime, and its crossovers toward the RC and QD regimes, remain challenging to observe. Quite remarkably, the same observations can be repeated for theoretical calculations on microscopic models, for which the quantitative extent of the QC regime is rarely investigated22. A general scenario (corroborated by the present work) is that different observables exhibit thermal QC scaling over different regions in the (g, T) parameter space. Therefore it is crucial to identify those observables which manifest such a scaling over the broadest possible range.
Here, we propose a constructive definition of the QC regime based on observables that do not admit any classical analog, namely quantum coherence measures, capturing quantum correlations and fluctuations for generic mixed states. Using quantum Monte-Carlo numerical computations and exact-diagonalization on paradigmatic models of quantum phase transitions, we single out a broad, fan-shaped region in the phase diagram where quantum fluctuations of the order parameter obey the expected power-law scaling of the QC fan. Given the link between quantum fluctuations and the dynamical susceptibility of the order parameter, our results open the way to the experimental reconstruction of the QC region using, e.g., spectroscopy on strongly correlated quantum systems.
Results
QC regime from quantum correlations
At first sight, defining the QC regime using the scaling behavior of quantum fluctuations sounds very natural: the Tâ=â0 QCP is characterized by critical quantum fluctuations, and the QC regime, if regarded as echoing the QCP at finite T, should be characterized by enhanced quantum fluctuations as well. The quantum coherence measures of our interest generally belong to the family identified by Petz23,24 as generalizations of the quantum Fisher information (QFI)25,26. Among this family of quantities, we focus on the recently proposed quantum variance (QV)27 of an observable O, possessing a simple definition at thermal equilibrium at inverse temperature βâ=â(kBT)â1; it is the difference between the (total) variance (TV) Var[O]â=âãO2ãâââãOã2 and the susceptibility
where ÏOâ=â(âãOã/âh)h=0â and h is a field coupling to O in the Hamiltonian as \({\cal H} - hO\). Beyond its transparent physical meaning (difference between fluctuations and response function), the QV (like the QFI) has the fundamental property of being an entanglement witnessâdenying separability of the state of the system into clusters of size k (or smaller) when it exceeds a k-dependent bound27,28,29,30. Moreover, unlike the QFI, it has the remarkable property of being directly accessible to state-of-the-art calculations for equilibrium quantum many-body systems at finite T, as e.g., worldline quantum Monte Carlo (the computation of the QFI requires the precise knowledge of the dynamical response function at real frequencies31, which can only be inferred from quantum Monte-Carlo data through analytical continuation, an operation very sensitive to numerical noise). This makes the QV an observable of choice to explore quantum coherence properties across the phase diagram of QC phenomena.
Both the variance and the susceptibility in Eq. (1) can be expressed as integrals of the imaginary part of the dynamical susceptibility32, resulting in the fundamental relationship:
where \({\cal L}(x) = {\mathrm{coth}}\,x - 1/x\) is the Langevin function. Since \({\cal L}(x) \to x/3\) for xâââ0, one sees that VarQ[O] is dominated by modes with frequency Ï such that \(\beta \hbar \omega \,\gtrsim \,1\), namely modes which are mildly (or not at all) affected by thermal fluctuations. A similar expression to Eq. (2) holds for the QFI, with the replacement \({\cal L}(x) \to 4\,{\mathrm{tanh}}(x)\)31. When O is the order parameter of the quantum phase transition of interest, the dynamical susceptibility in the vicinity of the QCP is expected to obey the scaling form
where η, ν, and z are the correlation function exponent, correlation length exponent and dynamical critical exponent of the QCP, respectively, and ΦO is a universal function up to a prefactor2. This directly translates into a scaling Ansatz for the QV:
where Ïâ=â(2âââη)/zâââ1 and \(G_O^{({\mathrm{Q}})}\sim {\int} d\omega \,{\cal L}(\beta \hbar \omega /2)\,{\mathrm{\Phi }}_O\). The TV of the order parameter Var[O] possesses a similar scaling form to Eq. (4), but with a different scaling function \(G_O^{({\mathrm{Q}})} \to G_O\sim {\int} d\omega \,{\mathrm{coth}}(\beta \hbar \omega /2)\,{\mathrm{\Phi }}_O\).
Eq. (4) forms the basis of our constructive definition of the QC regime as detected by quantum correlations. Being controlled by the QCP alone, such a regime must be nearly insensitive to whether the control parameter g lies above or below gc. This defining condition requires that, in the QC regime, the scaling function \(G_O^{({\mathrm{Q}})}(x)\) depend very weakly on its argument, namely \(G_O^{({\mathrm{Q}})}(x) \approx G_O^{({\mathrm{Q}})}(0)\). This leads us then to the following quantitative definition for the QC regime in the (g, T) plane:
where the â sign implies that the above condition is satisfied within some tolerance. As the QC regime is not a phase of matter which is divided from competing phases by sharp boundaries, the tolerance defines operatively the crossover lines toward the other regimes (RC and QD) in the vicinity of the QCP. The constructive definition of the QC regime offered by Eq. (5) identifies the latter regime with the region in the g-T phase diagram in which the T-dependence of the quantum fluctuations of the order parameter is uniquely controlled by the presence of the QCPânamely it is the same (up to some tolerance) as along the so-called QC trajectory (the line of variable T at gâ=âgc). The very fact that such a regime exists in an extended, fan-shaped region, is a fundamental test of the validity of our definition. A similar condition could obviously be formulated for the more conventional TV in the form \({\mathrm{Var}}[O] \approx T^{1 - \frac{{2 - \eta }}{z}}G_O(0)\): as we shall see shortly, this condition in practice singles out only the QC trajectory.
2d transverse field Ising model
We demonstrate our constructive definition of the QC regime using two paradigmatic examples of quantum phase transitions in quantum spin models. We begin by considering the 2d transverse field Ising (TFI) model4
where the indices ãijã and i run over the nearest-neighbor bonds and the sites of a square lattice, respectively; and \(S_i^\alpha\) (αâ=âx, y, z) are Sâ=â1/2 spin operators. A critical value gcâ=â1.522â¦33 of the transverse field divides a ferromagnetic regime (gâ<âgc) from a quantum paramagnetic one (gâ>âgc). We calculate the equilibrium properties of this model on Nâ=âLâÃâL lattices with periodic boundary conditions numerically using Stochastic Series Expansion quantum Monte Carlo34,35, which gives direct access to the TV and QV of most relevant observables27. The temperature scaling of the variances (total and quantum) of the macroscopic order parameter \(J^{\mathrm{z}} = \mathop {\sum}\nolimits_i S_i^{\mathrm{z}}\) in the vicinity of the QCP are contrasted in Fig. 2. Here we take for simplicity Var(Jz)â=âã(Jz)2ã as ãJzãâ=â0 on finite lattices. We have checked that our conclusions do not change when considering a finite-size estimate of the actual variance, namely \(\langle (J^{\mathrm{z}})^2\rangle - N^2m_L^2\) where \(m_L^2 = \langle S_i^{\mathrm{z}}S_{i + L/2}^z\rangle\). The QCP appears to control the T dependence of Var(Jz) only along the QC trajectory [Fig. 2a], where it exhibits the expected power-law dependence \(\sim T^{ - \psi }\); on the contrary the TV is strongly bent upward by the finite-T transition for gâ<âgc, as well as downward by the opening of a spin gap for gâ>âgc.
This picture is completely changed when one looks at quantum fluctuations. Indeed Fig. 2b shows that a power-law QC scaling of the QV as \(\sim T^{ - \psi }\) is manifested not only along the QC trajectory (down to Tâ=â0), but it can be observed also over sizable segments of the T-dependence both above and below the QCP. For gâ<âgc this is due to the fundamental property of the QV to be nearly insensitive to finite-temperature transitions27âas it will be discussed in a future publication, only weak singularities, in the form of inflection points, can appear at Tc. Therefore, the thermal critical region (in which power-law singularities dominate the behavior of the system) is minute in the T-dependence of the QV. Interestingly, a similar observation also applies to the case gâ>âgc in which a finite-T transition is absent. Indeed, unlike thermal fluctuations, quantum fluctuations are much more moderately suppressed by the opening of a gap. This observation can be understood by considering that, associated to the quantum fluctuations of the order parameter, there is an intrinsic quantum coherence length ξQ36 which is always finite at finite T, and much smaller than the ordinary correlation length ξ. Approaching the QCP, ξâââcTâ1/zâ=âc/T, but the opening of a gap (Î) for gââ âgc cuts off the QC growth, as ξ saturates to its ground-state value ξ(Tâ=â0)âââc/Î. Such saturation occurs for TâââÎ. Similarly, one expects that ξQâââcQ/T, but with cQââ<< c; nonetheless, at Tâ=â0, ξQ saturates to the same value ξQ(Tâ=â0)â=âξ(Tâ=â0)â=âc/Î. Therefore, the saturation occurs at lower temperatures, Tâââ(cQ/c)Î. In Supplementary Fig. 1, we show that cQ/câââ1/6. Hence, one may expect that saturation temperatures for ξ and ξQ are in a similar proportion in the gapped phase. In summary, both above and below gc the QV exhibits a clear crossover from a power-law regime varying as TâÏ to a saturating regime, occurring around a temperature Tâ~âTc or Tâ~âÎ: this behavior reveals then the crossover from the QC regime to the RC and QD regimes. We can deduce that the QCP controls the thermodynamics of quantum fluctuations over a region of sizable width, and which is fan-shaped (namely broader in g, the higher T).
In order to quantitatively reconstruct the region in which the T-dependence of the QV is influenced by the existence of the QCP, we establish the following criteria: (1) the QV exhibits a power-law dependence in T with an exponent (namely its logarithmic derivative) reproducing Ïâ=â0.964â¦37 (within some tolerance ε); (2) the coefficient of the power-law dependence (estimated as VarQ(Jz)TÏ) reproduces the one along the QC trajectory G(Q)(0) (within the same tolerance ε). These two criteria are illustrated in Fig. 3a, b: the region matching both criteria is then identified as the QC regime of quantum fluctuations. We would like to stress that the regions identified by criteria (1) or (2) individually can be quite different, and only their intersection can be reliably considered as representative of the QC regime. For instance criterion (2) is satisfied in a sizable portion for gââ¥âgc at low T, but this is a mere coincidence due to the fact that VarQ(Jz)TÏ turns out to be a non-monotonic function, crossing twice the value G(Q)(0)âat a lower and at a higher temperature. But only the higher-temperature crossing occurs with the QV exhibiting a logarithmic derivative compatible with âÏ, and hence complying with criterion 1).
Obviously the extent of such a regime in the phase diagram depends crucially on ε (taken as 10% in Fig. 3): yet it is important to observe that, regardless of the value of ε, its lower boundaries, marking the onset of the crossovers toward the RC and QD regimes, follow faithfully the temperature scales set by Tc and Îâboth scaling as |gâââgc|νz in the vicinity of the QCP. In contrast, applying similar criteria to the scaling of Var(Jz) essentially reconstructs the QC trajectory only, as already anticipated above (Fig. 3c). In Supplementary Notes 3 and 4, we provide evidence of a similar phenomenology of the QV, as well as of the QFI, for the exactly solvable cases of the TFI in dâ=â1 (Supplementary Note 4) and dâ=ââ (Supplementary Note 3). In the dâ=â2 case, the chosen tolerance ε assigns to the QC regime a temperature extent which ranges up to a temperature Tâ~â0.3J. This appears as a reasonable upper bound, as it remains sizeably smaller than the temperature scale Tâ~âJ: at that temperature short-wavelength modes become excited, making the microscopic details of the lattice model overcome the universal power-law temperature scalings which are distinctive of the QC regime.
Heisenberg bilayer
We conclude by considering another paradigmatic model of quantum phase transitions, namely the Sâ=â1/2 Heisenberg antiferromagnet on a square lattice bilayer, with Hamiltonian
comprising intralayer (ãijã) as well as interlayer (ãlmãâ¥) bonds. A quantum phase transition at gcâ=â2.522â¦38 divides a Néel antiferromagnetic regime (with order parameter given by the staggered magnetization \(J_{{\mathrm{st}}}^{\mathrm{z}} = \mathop {\sum}\nolimits_i ( - 1)^iS_i^{\mathrm{z}}\)) from a nonmagnetic dimer-singlet regime. The same kind of analysis of the order-parameter QV, as the one presented above for the 2d TFI model, leads to Fig. 4. There the QC region, where the thermal behavior of quantum fluctuations is governed by the QC point, is shown to be very broad (with a maximum width Îg which is around 20% of gc). On the contrary, a similar analysis based on the behavior of the TV singles out only a narrow region around the QC trajectory (Fig. 4c), similarly to what we reported for the TFI model. In the case at hand the QC-RC crossover is marked by a crossover in the QV from the TâÏ power-law behavior into another power-law behavior, as the RC regime is gapless and without a finite-T phase transition; while the QC-QD crossover is similar to the one observed in the 2d TFI model (low-T saturation of the QV). Despite its sizable width, the QC regime of the QV remains limited to a finite g range, and it does not come close to the limit gâ=â0 corresponding to the most investigated case of the 2d Heisenberg antiferromagnet (2dHAF)âthis holds true even when restricting uniquely to the criterion of the logarithmic derivative. The 2dHAF Heisenberg model has been the subject of an intense search for signatures of a QC-RC crossover in the past8,39,40,41,42: our data (see Supplementary Fig. 2) exclude that such a crossover is visible in the QV of the order parameter.
Discussion
Based on numerical simulations (quantum Monte Carlo) and exact diagonalization on the paradigmatic transverse-field Ising and bilayer Heisenberg models, we have provided evidence that the existence of a zero-temperature QC point (QCP) fully controls the thermodynamics of the quantum fluctuations of the order parameter (estimated via the QV) in a broad, fan-shaped region above the QCP itself. Such a region can be identified with the elusive QC regime, acting as a finite-T magnifying lens of zero-T quantum criticality. The extent of the QC regime, as revealed by quantum fluctuations, far exceeds that of conventional fluctuations propertiesâthe latter contain a large thermal component subject to thermal criticality on one side of the QCP, and to a large suppression due to the opening of a gap on the other side. Therefore we open the unconventional perspective of using a property which bears no classical analog to unveil a finite-temperature regimeâsomewhat reminiscent of the use of entanglement to characterize QCPs at Tâ=â0.
Our proposed constructive definition of the QC regime, summarized by Eqs. (1) and (5), is completely general, and therefore it is immediately applicable to detect the QC regime in numerical as well as field-theoretical studies that have naturally access to the QV. For instance, it could be applied to unveil the QC regime of unconventional âdeconfined" QC points43,44,45. Most importantly, the quantitative relationship between the dynamical susceptibility and quantum fluctuations offers the possibility to access the latter in spectroscopic experiments on strongly correlated materials. The main experimental requirement is the ability to measure the dynamical response to a probe coupling to the order parameter of the quantum phase transition; and the ability to reconstruct the corresponding dynamical susceptibility over the relevant frequency range where its imaginary part Ïâ²â² has its strongest support. These requirements are clearly met by modern experiments in neutron spectroscopy46 or AC susceptometry47,48,49 on bulk materials, to cite a few relevant examples. Our proposed scheme can serve as an effective tool to unveil the existence of zero-T QCPs via finite-T experiments, especially in situations in which other signatures of the existence of a QCP in the low-T thermodynamics prove elusive19,50.
Methods
Monte Carlo calculations
The calculations of the total and quantum variance of the order parameter shown in the main text were obtained using stochastic series expansion (SSE) quantum Monte Carlo (QMC)34 on Nâ=âLâÃâL lattices with periodic boundary conditions. In the case of the transverse-field Ising model the quantization axis was chosen to lie along the field (x) axis35, so that the correlations among the z components of the spins correspond to off-diagonal observables, which are sampled with very high statistics during the directed-loop updates34. A typical simulation consists of 104 MC steps (each comprising as many directed loops as necessary to attempt on average to update all the spins in the SSE) for thermalization, and 105 MC steps to accumulate the statistics for the observables of interest. The QV is straightforwardly obtained as the difference between the integrated equal-time correlation function and the integrated (imaginary-)time-averaged correlations51.
Data availability
The data supporting the findings presented in this paper, as well as the computer codes used to generate them, are available from the authors upon reasonable request.
References
Sondhi, S. L., Girvin, S. M., Carini, J. P. & Shahar, D. Continuous quantum phase transitions. Rev. Mod. Phys. 69, 315â333 (1997).
Sachdev, S. Quantum Phase Transitions (Cambridge University Press, Cambridge, 2011).
Continentino, M. Quantum Scaling in Many-Body Systems (Cambridge University Press, Cambridge, 2017).
Dutta, A. et al. Quantum Phase Transitions in Transverse Field Models (Cambridge University Press, Cambridge, 2015).
Carr, L. (ed.) Understanding Quantum Phase Transitions (CRC Press, Boca Raton FL, 2010).
Kaul, R. K., Melko, R. G. & Sandvik, A. W. Bridging lattice-scale physics and continuum field theory with quantum monte carlo simulations. Annu. Rev. Condens. Matter Phys. 4, 179â215 (2013).
Chakravarty, S., Halperin, B. I. & Nelson, D. R. Two-dimensional quantum Heisenberg antiferromagnet at low temperatures. Phys. Rev. B 39, 2344â2371 (1989).
Chubukov, A. V., Sachdev, S. & Ye, J. Theory of two-dimensional quantum heisenberg antiferromagnets with a nearly critical ground state. Phys. Rev. B 49, 11919â11961 (1994).
Gegenwart, P., Si, Q. & Steglich, F. Quantum criticality in heavy-fermion metals. Nat. Phys. 4, 186â197 (2008).
Trotzky, S. et al. Suppression of the critical temperature for superfluidity near the Mott transition. Nat. Phys. 6, 998â1004 (2010).
Endres, M. et al. The âHiggsâ amplitude mode at the two-dimensional superfluid/Mott insulator transition. Nature 487, 454â458 (2012).
Zhang, X., Hung, C.-L., Tung, S.-K. & Chin, C. Observation of quantum criticality with ultracold atoms in optical lattices. Science 335, 1070â1072 (2012).
Lake, B., Tennant, D. A., Frost, C. D. & Nagler, S. E. Quantum criticality and universal scaling of a quantum antiferromagnet. Nat. Mater. 4, 329â334 (2005).
Mukhopadhyay, S. et al. Quantum-critical spin dynamics in quasi-one-dimensional antiferromagnets. Phys. Rev. Lett. 109, 177206 (2012).
Kinross, A. W. et al. Evolution of quantum fluctuations near the quantum critical point of the transverse field ising chain system CoNb2O6. Phys. Rev. X 4, 031008 (2014).
Povarov, K. Y., Schmidiger, D., Reynolds, N., Bewley, R. & Zheludev, A. Scaling of temporal correlations in an attractive TomonagaâLuttinger spin liquid. Phys. Rev. B 91, 020406 (2015).
Hälg, M., Hüvonen, D., Butch, N. P., Demmel, F. & Zheludev, A. Finite-temperature scaling of spin correlations in a partially magnetized Heisenberg sâ=â1/2 chain. Phys. Rev. B 92, 104416 (2015).
Schröder, A. et al. Onset of antiferromagnetism in heavy-fermion metals. Nature 407, 351â355 (2000).
Sachdev, S. Where is the quantum critical point in the cuprate superconductors? Phys. Status Solidi B 247, 537â543 (2010).
Ramshaw, B. J. et al. Quasiparticle mass enhancement approaching optimal doping in a high-tc superconductor. Science 348, 317â320 (2015).
Badoux, S. et al. Change of carrier density at the pseudogap critical point of a cuprate superconductor. Nature 531, 210â214 (2016).
Kopp, A. & Chakravarty, S. Criticality in correlated quantum matter. Nat. Phys. 1, 53â56 (2005).
Petz, D. Monotone metrics on matrix spaces. Linear Algebra Appl. 244, 81â96 (1996).
Streltsov, A., Adesso, G. & Plenio, M. B. Colloquium: quantum coherence as a resource. Rev. Mod. Phys. 89, 041003 (2017).
Braunstein, S. L. & Caves, C. M. Statistical distance and the geometry of quantum states. Phys. Rev. Lett. 72, 3439â3443 (1994).
Pezzé, L. & Smerzi, A. Quantum theory of phase estimation. In Tino, G. M. & Kasevich, M. A. (eds.) Atom Interferometry, Proceedings of the International School of Physics âEnrico Fermiâ, Course 188 (IOS Press, Amsterdam, 2014).
Frérot, I. & Roscilde, T. Quantum variance: a measure of quantum coherence and quantum correlations for many-body systems. Phys. Rev. B 94, 075121 (2016).
Hyllus, P. et al. Fisher information and multiparticle entanglement. Phys. Rev. A 85, 022321 (2012).
Tóth, G. Multipartite entanglement and high-precision metrology. Phys. Rev. A 85, 022322 (2012).
Pezzè, L., Smerzi, A., Oberthaler, M. K., Schmied, R. & Treutlein, P. Quantum metrology with nonclassical states of atomic ensembles. Rev. Mod. Phys. 90, 035005 (2018).
Hauke, P., Heyl, M., Tagliacozzo, L. & Zoller, P. Measuring multipartite entanglement via dynamic susceptibilities. Nat. Phys. 12, 778â782 (2016).
Foerster, D. Hydrodynamic Fluctuations, Broken Symmetry, and Correlation Functions (CRC Press, Boca Raton FL, 1995).
Blöte, H. W. J. & Deng, Y. Cluster monte carlo simulation of the transverse Ising model. Phys. Rev. E 66, 066110 (2002).
Sandvik, A. W. Computational studies of quantum spin systems. AIP Conf. Proc. 1297, 135â338 (2010).
Frérot, I. & Roscilde, T. Quantum critical metrology. Phys. Rev. Lett. 121, 020402 (2018).
Malpetti, D. & Roscilde, T. Quantum correlations, separability, and quantum coherence length in equilibrium many-body systems. Phys. Rev. Lett. 117, 130401 (2016).
Pelissetto, A. & Vicari, E. Critical phenomena and renormalization-group theory. Phys. Rep. 368, 549â727 (2002).
Sandvik, A. W. & Scalapino, D. J. Order-disorder transition in a two-layer quantum antiferromagnet. Phys. Rev. Lett. 72, 2777â2780 (1994).
Imai, T., Slichter, C. P., Yoshimura, K., Katoh, M. & Kosuge, K. Spinâspin correlation in the quantum critical regime of La2CuO4. Phys. Rev. Lett. 71, 1254â1257 (1993).
Greven, M. et al. Neutron scattering study of the two-dimensional spinsâ=â1/2 square-lattice heisenberg antiferromagnet Sr2CuO2Cl2. Z. für Phys. B Condens. Matter 96, 465â477 (1995).
Rønnow, H. M., McMorrow, D. F. & Harrison, A. High-temperature magnetic correlations in the 2d Sâ=â1/2 antiferromagnet copper formate tetradeuterate. Phys. Rev. Lett. 82, 3152â3155 (1999).
Kim, J.-K. & Troyer, M. Low temperature behavior and crossovers of the square lattice quantum heisenberg antiferromagnet. Phys. Rev. Lett. 80, 2705â2708 (1998).
Senthil, T., Vishwanath, A., Balents, L., Sachdev, S. & Fisher, M. P. A. Deconfined quantum critical points. Science 303, 1490â1494 (2004).
Sandvik, A. W. Evidence for deconfined quantum criticality in a two-dimensional heisenberg model with four-spin interactions. Phys. Rev. Lett. 98, 227202 (2007).
Melko, R. G. & Kaul, R. K. Scaling in the fan of an unconventional quantum critical point. Phys. Rev. Lett. 100, 017203 (2008).
Haravifard, S., Yamani, Z. & Gaulin, B. D. in Neutron ScatteringâMagnetic and Quantum Phenomena, Vol. 48 (eds Fernandez-Alonso, F. & Price, D. L.) Ch. 2 (Academic Press, Elsevier, London UK, 2015).
Yin, L., Xia, J. S., Zapf, V. S., Sullivan, N. S. & Paduan-Filho, A. Direct measurement of the BoseâEinstein condensation universality class in NiCl2â4SC(NH2)2 at ultralow temperatures. Phys. Rev. Lett. 101, 187205 (2008).
Quilliam, J. A., Yaraskavitch, L. R., Dabkowska, H. A., Gaulin, B. D. & Kycia, J. B. Dynamics of the magnetic susceptibility deep in the coulomb phase of the dipolar spin ice material Ho2Ti2O7. Phys. Rev. B 83, 094424 (2011).
Bovo, L., Bloxsom, J. A., Prabhakaran, D., Aeppli, G. & Bramwell, S. T. Brownian motion and quantum dynamics of magnetic monopoles in spin ice. Nat. Commun. 4, 1535 (2013).
Shibauchi, T., Carrington, A. & Matsuda, Y. A quantum critical point lying beneath the superconducting dome in iron pnictides. Annu. Rev. Condens. Matter Phys. 5, 113â135 (2014).
Sandvik, A. W. A generalization of Handscombâs quantum Monte Carlo scheme-application to the 1D Hubbard model. J. Phys. A. Math. Gen. 25, 3667â3682 (1992).
Gabbrielli, M., Smerzi, A. & Pezzè, L. Multipartite entanglement at finite temperature. Sci. Rep. 8, 15663 (2018).
Hesselmann, S. & Wessel, S. Thermal Ising transitions in the vicinity of two-dimensional quantum critical points. Phys. Rev. B 93, 155157 (2016).
Acknowledgments
We thank L. Pezzé and A. Smerzi for fruitful discussions, and for informing us of their recent preprint on related work about the quantum Fisher information around quantum critical points52. This work was supported by ANR (âArtiQâ project). Numerical simulations were run on the PSMN cluster (ENS de Lyon). I.F. acknowledges support from the Spanish Ministry MINECO (National Plan 15 Grant: FISICATEAMO No. FIS2016-79508-P, SEVERO OCHOA No. SEV-2015-0522), Fundació Cellex, Generalitat de Catalunya (AGAUR Grant no. 2017 SGR 1341 and CERCA/Program), ERC AdG OSYRIS, EU FETPRO QUIC, and the National Science Centre, Poland-Symfonia Grant no. 2016/20/W/ST4/00314.
Author information
Authors and Affiliations
Contributions
T.R. performed the quantum Monte Carlo calculations for the Ising and Heisenberg models. I.F. performed the exact-diagonalization calculations for the dâ=â1 and dâ=ââ quantum Ising models. Both authors contributed to the writing of the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Journal peer review information: Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisherâs note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Frérot, I., Roscilde, T. Reconstructing the quantum critical fan of strongly correlated systems using quantum correlations. Nat Commun 10, 577 (2019). https://doi.org/10.1038/s41467-019-08324-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-019-08324-9