“…In last ten years, various numerical techniques were used for calculation of the QCP in the SBM, such as the numerical renormalization group (NRG) [3,4,5], the quantum Monte Carlo (QMC) [8], the method of sparse polynomial space representation [9], the extended coherent state approach [10], and the variational matrix product state approach [11]. Besides, recently an extension of the Silbey-Harris ground state was proposed by Zhao et al [12] and Chin et al [13] to study the QPT in the s = 1/2 sub-Ohmic SBM.…”