Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1012.5615 · Physical Review B

Scaling analysis in the numerical renormalization group study of the sub-Ohmic spin-boson model

Ning-Hua Tong, Yan-Hua Hou

Published Apr 25, 2012• 13 citations
arXivHTMLPDF

Abstract

The spin-boson model has nontrivial quantum phase transitions in the sub-Ohmic regime. For the bath spectra exponent $0\ensuremath{\leqslant}s<1/2$, the bosonic numerical renormalization group (BNRG) study of the exponents $\ensuremath{\beta}$ and $\ensuremath{\delta}$ are hampered by the boson-state truncation, which leads to artificial interacting exponents instead of the correct Gaussian ones. In this paper, guided by a mean-field calculation, we study the order-parameter function $m(\ensuremath{\tau}=\ensuremath{\alpha}\ensuremath{-}{\ensuremath{\alpha}}_{c},\ensuremath{\epsilon},\ensuremath{\Delta})$ using BNRG. Scaling analysis with respect to the boson-state truncation ${N}_{b}$, the logarithmic discretization parameter $\ensuremath{\Lambda}$, and the tunneling strength $\ensuremath{\Delta}$ are carried out. Truncation-induced multiple-power behaviors are observed close to the critical point, with artificial values of $\ensuremath{\beta}$ and $\ensuremath{\delta}$. They cross over to classical behaviors with exponents $\ensuremath{\beta}=1/2$ and $\ensuremath{\delta}=3$ on the intermediate scales of $\ensuremath{\tau}$ and $\ensuremath{\epsilon}$, respectively. We also find $\ensuremath{\tau}/{\ensuremath{\Delta}}^{1\ensuremath{-}s}$ and $\ensuremath{\epsilon}/\ensuremath{\Delta}$ scalings in the function $m(\ensuremath{\tau},\ensuremath{\epsilon},\ensuremath{\Delta})$. The role of boson-state truncation as a scaling variable in the BNRG result for $0\ensuremath{\leqslant}s<1/2$ is identified and its interplay with the logarithmic discretization revealed. Relevance to the validity of quantum-to-classical mapping in other impurity models is discussed.

DOI
10.1103/physrevb.85.144425

Related papers

Criticality of the Mean-Field Spin-Boson Model: Boson State Truncation and Its Scaling Analysis

Yan-hua Hou, Ning-hua Tong

The spin-boson model has nontrivial quantum phase transitions at zero temperature induced by the spin-boson coupling. The bosonic numerical renormalization group (BNRG) study of the critical exponents β and δ of this model is hampered by the effects of boson Hilbert space truncation. Here we analyze the mean-field spin boson model to figure out the scaling behavior of magnetization under the cutoff of boson states N b . We find that the truncation is a strong relevant operator with respect to the Gaussian fixed point in 0 < s < 1/2 and incurs the deviation of the exponents from the classical values. The magnetization at zero bias near the critical point is described by a generalized homogeneous function (GHF) of two variables τ = α – α c and x = 1/N b . The universal function has a double-power form and the powers are obtained analytically as well as numerically. Similarly, m(α = α c ) is found to be a GHF of ϵ and x. In the regime s > 1/2, the truncation produces no effect. Implications of these findings to the BNRG study are discussed.

•Aug 16, 2016•11 citations

Numerical renormalization group for the sub-Ohmic spin-boson model: A conspiracy of errors

Matthias Vojta C

The application of Wilson's numerical renormalization group (NRG) method to dissipative quantum impurity models, in particular the sub-Ohmic spin-boson model, has led to conclusions regarding the quantum critical behavior which are in disagreement with those from other methods and which are by now recognized as erroneous. The errors of NRG remained initially undetected, because NRG delivered an internally consistent set of critical exponents satisfying hyperscaling. Here we discuss how the conspiracy of two errors---the Hilbert-space truncation error and the mass-flow error---could lead to this consistent set of exponents. Remarkably, both errors, albeit of different origin, force the system to obey naive scaling laws even when the physical model violates naive scaling. In particular, we show that a combination of the Hilbert-space truncation and mass-flow errors induce an artificial nonanalytic term in the Landau expansion of the free energy which dominates the critical behavior for bath exponents $s<1/2$.

arXiv:1201.4922•Physical Review B•Mar 15, 2012•37 citations

Numerical renormalization group for quantum impurities in a bosonic bath

R. Bulla, Hyunjung Lee, Ning-Hua Tong +1 more

We present a detailed description of the recently proposed numerical renormalization group method for models of quantum impurities coupled to a bosonic bath. Specifically, the method is applied to the spin-boson model, both in the Ohmic and sub-Ohmic cases. We present various results for static as well as dynamic quantities and discuss details of the numerical implementation, e.g., the discretization of a bosonic bath with arbitrary continuous spectral density, the suitable choice of a finite basis in the bosonic Hilbert space, and questions of convergence with respect to truncation parameters. The method is shown to provide high-accuracy data over the whole range of model parameters and temperatures, which are in agreement with exact results and other numerical data from the literature.

arXiv:cond-mat/0407559•Physical Review B•Jan 25, 2005•174 citations

Numerical Renormalization Group for Bosonic Systems and Application to the Sub-Ohmic Spin-Boson Model

R. Bulla, Ning-Hua Tong, Matthias Vojta C

We describe the generalization of Wilson's numerical renormalization group method to quantum impurity models with a bosonic bath, providing a general nonperturbative approach to bosonic impurity models which can access exponentially small energies and temperatures. As an application, we consider the spin-boson model, describing a two-level system coupled to a bosonic bath with power-law spectral density, J(omega) proportional to omega(s). We find clear evidence for a line of continuous quantum phase transitions for sub-Ohmic bath exponents 0<s<1; the line terminates in the well-known Kosterlitz-Thouless transition at s=1. Contact is made with results from perturbative renormalization group, and various other applications are outlined.

arXiv:cond-mat/0306708•Physical Review Letters•Oct 24, 2003•265 citations

Quantum Phase Transition in the Sub-Ohmic Spin-Boson Model: Quantum Monte Carlo Study with a Continuous Imaginary Time Cluster Algorithm

André Winter, Heiko Rieger, Matthias Vojta C +1 more

A continuous time cluster algorithm for two-level systems coupled to a dissipative bosonic bath is presented and applied to the sub-Ohmic spin-boson model. When the power s of the spectral function Jomega proportional, variant omegas is smaller than 1/2, the critical exponents are found to be classical, mean-field like. Potential sources for the discrepancy with recent renormalization group predictions are traced back to the effect of a dangerously irrelevant variable.

arXiv:0807.4716•Physical Review Letters•Jan 20, 2009•156 citations

Quantum criticality in sub-Ohmic systems with three competing terms: beyond conventional spin-boson physics

Nengji Zhou, Yulong Shen, Zhe Sun

Quantum phase transitions (QPTs) in the spin-boson model with/without the rotating-wave approximation (RWA) are systematically investigated through variational calculations using a sub-Ohmic bath with high spectral density. Four cases involving different system-environment interactions are examined, where transition points and critical exponents are accurately determined across varying tunneling strengths. Contrary to prior work, a rich phase diagram is revealed in the tunneling-coupling plane even at the low spectral exponent $s<1/2$, with a novel U(1)-symmetric phase being identified. As coupling increases, a multi-stage QPT sequence arises for the tunneling $0<Δ< Δ^*=0.074(1)$, whereas a single transition occurs beyond this range. Furthermore, an odd-parity phase is found to emerge even under the positive tunneling, exhibiting distinct characteristics relative to the prototype model.

quant-pharXiv:2603.07031•Physical review. B./Physical review. B•Mar 10, 2026•0 citations
The spin-boson model has nontrivial quantum phase transitions in the sub-Ohmic regime. For the bath spectra exponent 0\ensuremath⩽s<1/20\ensuremath{\leqslant}s<1/20\ensuremath⩽s<1/2, the bosonic numerical renormalization group (BNRG) study of the exponents \ensuremathβ\ensuremath{\beta}\ensuremathβ and \ensuremathδ\ensuremath{\delta}\ensuremathδ are hampered by the boson-state truncation, which leads to artificial interacting exponents instead of the correct Gaussian ones. In this paper, guided by a mean-field calculation, we study the order-parameter function m(\ensuremathτ=\ensuremathα\ensuremath−\ensuremathαc,\ensuremathϵ,\ensuremathΔ)m(\ensuremath{\tau}=\ensuremath{\alpha}\ensuremath{-}{\ensuremath{\alpha}}_{c},\ensuremath{\epsilon},\ensuremath{\Delta})m(\ensuremathτ=\ensuremathα\ensuremath−\ensuremathαc​,\ensuremathϵ,\ensuremathΔ) using BNRG. Scaling analysis with respect to the boson-state truncation Nb{N}_{b}Nb​, the logarithmic discretization parameter \ensuremathΛ\ensuremath{\Lambda}\ensuremathΛ, and the tunneling strength \ensuremathΔ\ensuremath{\Delta}\ensuremathΔ are carried out. Truncation-induced multiple-power behaviors are observed close to the critical point, with artificial values of \ensuremathβ\ensuremath{\beta}\ensuremathβ and \ensuremathδ\ensuremath{\delta}\ensuremathδ. They cross over to classical behaviors with exponents \ensuremathβ=1/2\ensuremath{\beta}=1/2\ensuremathβ=1/2 and \ensuremathδ=3\ensuremath{\delta}=3\ensuremathδ=3 on the intermediate scales of \ensuremathτ\ensuremath{\tau}\ensuremathτ and \ensuremathϵ\ensuremath{\epsilon}\ensuremathϵ, respectively. We also find \ensuremathτ/\ensuremathΔ1\ensuremath−s\ensuremath{\tau}/{\ensuremath{\Delta}}^{1\ensuremath{-}s}\ensuremathτ/\ensuremathΔ1\ensuremath−s and \ensuremathϵ/\ensuremathΔ\ensuremath{\epsilon}/\ensuremath{\Delta}\ensuremathϵ/\ensuremathΔ scalings in the function m(\ensuremathτ,\ensuremathϵ,\ensuremathΔ)m(\ensuremath{\tau},\ensuremath{\epsilon},\ensuremath{\Delta})m(\ensuremathτ,\ensuremathϵ,\ensuremathΔ). The role of boson-state truncation as a scaling variable in the BNRG result for 0\ensuremath⩽s<1/20\ensuremath{\leqslant}s<1/20\ensuremath⩽s<1/2 is identified and its interplay with the logarithmic discretization revealed. Relevance to the validity of quantum-to-classical mapping in other impurity models is discussed.
The spin-boson model has nontrivial quantum phase transitions at zero temperature induced by the spin-boson coupling. The bosonic numerical renormalization group (BNRG) study of the critical exponents β and δ of this model is hampered by the effects of boson Hilbert space truncation. Here we analyze the mean-field spin boson model to figure out the scaling behavior of magnetization under the cutoff of boson states N b . We find that the truncation is a strong relevant operator with respect to the Gaussian fixed point in 0 < s < 1/2 and incurs the deviation of the exponents from the classical values. The magnetization at zero bias near the critical point is described by a generalized homogeneous function (GHF) of two variables τ = α – α c and x = 1/N b . The universal function has a double-power form and the powers are obtained analytically as well as numerically. Similarly, m(α = α c ) is found to be a GHF of ϵ and x. In the regime s > 1/2, the truncation produces no effect. Implications of these findings to the BNRG study are discussed.
The application of Wilson's numerical renormalization group (NRG) method to dissipative quantum impurity models, in particular the sub-Ohmic spin-boson model, has led to conclusions regarding the quantum critical behavior which are in disagreement with those from other methods and which are by now recognized as erroneous. The errors of NRG remained initially undetected, because NRG delivered an internally consistent set of critical exponents satisfying hyperscaling. Here we discuss how the conspiracy of two errors---the Hilbert-space truncation error and the mass-flow error---could lead to this consistent set of exponents. Remarkably, both errors, albeit of different origin, force the system to obey naive scaling laws even when the physical model violates naive scaling. In particular, we show that a combination of the Hilbert-space truncation and mass-flow errors induce an artificial nonanalytic term in the Landau expansion of the free energy which dominates the critical behavior for bath exponents s<1/2s<1/2s<1/2.
We present a detailed description of the recently proposed numerical renormalization group method for models of quantum impurities coupled to a bosonic bath. Specifically, the method is applied to the spin-boson model, both in the Ohmic and sub-Ohmic cases. We present various results for static as well as dynamic quantities and discuss details of the numerical implementation, e.g., the discretization of a bosonic bath with arbitrary continuous spectral density, the suitable choice of a finite basis in the bosonic Hilbert space, and questions of convergence with respect to truncation parameters. The method is shown to provide high-accuracy data over the whole range of model parameters and temperatures, which are in agreement with exact results and other numerical data from the literature.
We describe the generalization of Wilson's numerical renormalization group method to quantum impurity models with a bosonic bath, providing a general nonperturbative approach to bosonic impurity models which can access exponentially small energies and temperatures. As an application, we consider the spin-boson model, describing a two-level system coupled to a bosonic bath with power-law spectral density, J(omega) proportional to omega(s). We find clear evidence for a line of continuous quantum phase transitions for sub-Ohmic bath exponents 0<s<1; the line terminates in the well-known Kosterlitz-Thouless transition at s=1. Contact is made with results from perturbative renormalization group, and various other applications are outlined.
Quantum phase transitions (QPTs) in the spin-boson model with/without the rotating-wave approximation (RWA) are systematically investigated through variational calculations using a sub-Ohmic bath with high spectral density. Four cases involving different system-environment interactions are examined, where transition points and critical exponents are accurately determined across varying tunneling strengths. Contrary to prior work, a rich phase diagram is revealed in the tunneling-coupling plane even at the low spectral exponent s<1/2s<1/2s<1/2, with a novel U(1)-symmetric phase being identified. As coupling increases, a multi-stage QPT sequence arises for the tunneling 0<Δ<Δ∗=0.074(1)0<Δ< Δ^*=0.074(1)0<Δ<Δ∗=0.074(1), whereas a single transition occurs beyond this range. Furthermore, an odd-parity phase is found to emerge even under the positive tunneling, exhibiting distinct characteristics relative to the prototype model.