Dynamical correlation functions of the quadratic coupling spin-Boson model
The spin-boson model with quadratic coupling is studied using the bosonic numerical renormalization group method. We focus on the dynamical auto-correlation functions $C_{O}(\omega)$, with the operator $\hat{O}$ taken as $\hat{\sigma}_x$, $\hat{\sigma}_z$, and $\hat{X}$, respectively. In the weak-coupling regime $\alpha < \alpha_c$, these functions show power law $\omega$-dependence in the small frequency limit, with the powers $1+2s$, $1+2s$, and $s$, respectively. At the critical point $\alpha = \alpha_c$ of the boson-unstable quantum phase transition, the critical exponents $y_{O}$ of these correlation functions are obtained as $y_{\sigma_x} = y_{\sigma_z} = 1-2s$ and $y_{X}=-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_{\sigma_{x}}(\omega)$ is broadened significantly and the line shape changes qualitatively, showing enhanced dephasing at the spin flip point.