The spin-boson model with quadratic coupling is studied using the bosonic numerical renormalization group method.We focus on the dynamical auto-correlation functions Co(ω),with the operator (O) taken as (σ)x,(σ)z,and (X),respectively.In the weak-coupling regime α ＜ αc,these functions show power law ω-dependence in the small frequency limit,with the powers 1 + 2s,1 + 2s,and s,respectively.At the critical point α =αc of the boson-unstable quantum phase transition,the critical exponents Yo of these correlation functions are obtained as Yσx =Yσz =1-2s and yx =-s,respectively.Here s is the bath index and X is the boson displacement operator.Close to the spin flip point,the high frequency peak of Cσx (ω) is broadened significantly and the line shape changes qualitatively,showing enhanced dephasing at the spin flip point.