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

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.