---
title: Sub-Ohmic Spin-Boson Model
url: https://www.emergentmind.com/topics/sub-ohmic-spin-boson-model
type: topic
---

# Sub-Ohmic Spin-Boson Model

The sub-Ohmic spin-boson model describes a quantum two-level system coupled linearly to a bosonic environment whose low-frequency spectral density scales as \(J(\omega)\propto \omega^s\) with \(0<s<1\). Relative to the Ohmic case, the stronger infrared weight produces pronounced bath memory, strong renormalization of tunneling, and at zero temperature a delocalized-localized quantum phase transition; closely related work also studies “extremely sub-Ohmic” extensions such as \(s=0\) and \(s=-1\) [2309.00797, 1808.02916]. The model has therefore become a standard setting for dissipative quantum criticality, non-Markovian relaxation, and the analysis of decoherence by low-frequency environments [1111.6623, 1408.1061].

## 1. Canonical formulation and spectral classification

A standard spin-boson Hamiltonian is
\[
\hat{H}=\frac{\varepsilon}{2}\sigma_z-\frac{\Delta}{2}\sigma_x+\sum_k \omega_k b_k^\dagger b_k+\frac{\sigma_z}{2}\sum_k \lambda_k(b_k^\dagger+b_k),
\]
where \(\varepsilon\) is the bias, \(\Delta\) is the tunneling amplitude, and the bath enters through the mode frequencies \(\omega_k\) and couplings \(\lambda_k\) [2309.00797]. The bath is encoded by a spectral density \(J(\omega)\); one common convention is \(J(\omega)=2\alpha \omega_c^{1-s}\omega^s\), while several numerical-renormalization-group formulations use \(J(\omega)=2\pi\alpha\,\omega_c^{1-s}\omega^s\) with a cutoff \(\omega_c\) [2309.00797, 2010.07471].

The exponent \(s\) classifies the environment. The sub-Ohmic regime is \(0<s<1\), the Ohmic case is \(s=1\), and the super-Ohmic case is \(s>1\) [2309.00797]. Smaller \(s\) means stronger low-frequency weight and therefore stronger long-time memory effects from the bath [1111.6623]. Recent work often separates the sub-Ohmic interval into a deep regime \(0<s\le 0.5\) and a shallow regime \(0.5<s<1\), because the critical behavior and the numerical difficulty change qualitatively across that division [2309.00797].

The model is also studied in effective and extended forms. A mechanically engineered analog simulator based on a color center in a free-standing h-BN membrane realizes \(\mathcal J(\nu)\propto \nu^s\) with \(-1\le s\le 0\), with the cases \(s=-1\) and \(s=0\) corresponding to \(1/f\)-type and white-noise environments [1808.02916]. This suggests that the conventional sub-Ohmic family is best viewed as one part of a broader infrared-dominated dissipative landscape.

## 2. Zero-temperature phases and quantum criticality

At zero bias and zero temperature, the canonical sub-Ohmic model has a transition between a delocalized phase with vanishing order parameter and a localized phase with finite magnetization \(M=\langle \sigma_z\rangle\) [1111.6623, 1102.1483]. In a variational Asymmetrically Displaced Oscillator treatment, the transition is continuous for \(0<s<0.5\), with
\[
M\propto (\alpha-\alpha_c)^{1/2},\qquad
\chi\propto |\alpha-\alpha_c|^{-1},
\]
and
\[
\alpha_c=\frac{\sin(\pi s)e^{-s/2}}{2\pi(1-s)}
\left(\frac{\Delta}{\omega_c}\right)^{1-s}
\]
in the scaling limit [1102.1483]. That same work emphasized that the localized phase retains finite coherence \(\langle \sigma_x\rangle\), and interpreted localization as the emergence of a self-consistent low-frequency bath bias rather than a trivial freezing of all tunneling processes [1102.1483].

The universality class has been controversial. Bosonic NRG work on critical susceptibilities reported an interacting quantum critical point with \(x=y=1-\epsilon=s\) and \(\omega/T\) scaling across \(0<\epsilon<1\), and attributed the failure of the quantum-to-classical mapping to a Berry-phase term in the continuum spin path integral [1111.6623]. A displaced-Fock-states expansion instead reported \(\beta=0.5\) for the whole interval \(0<s<1\) and interpreted this as evidence that the system is always above its upper critical dimension [1410.0999]. A later multiple-polaron numerical variational study found mean-field criticality for \(0<s\le 0.5\), non-mean-field \(s\)-dependent criticality for \(0.5<s<1\), and stated that the quantum-to-classical correspondence is fully confirmed over the entire sub-Ohmic range [2309.00797].

These competing claims define a central interpretive fault line in the literature. A plausible implication is that the sub-Ohmic model is not only a problem in dissipative dynamics, but also a sensitive benchmark for how infrared bath physics, bosonic truncation, and continuum limits are handled numerically.

## 3. Nonequilibrium dynamics, transient phase diagrams, and hidden dynamical structure

Real-time nonequilibrium studies have shown that the dynamical organization of the sub-Ohmic model is richer than the equilibrium phase diagram alone. A numerically exact inchworm-QMC study extracted transient phase boundaries from the post-quench polarization
\[
\langle \sigma_z(t)\rangle \approx a\cos(\Omega t+\phi)e^{-\gamma_1 t}+be^{-\gamma_2 t}+c,
\]
using \(c\) as a localization diagnostic and \(\Omega\), \(\gamma_1\) as coherence diagnostics [2402.18561]. In that formulation, localization is associated with \(c\neq 0\), while coherence loss can occur either through a smooth overdamping crossover \(\Omega=\gamma_1\) or through a sharp frequency collapse \(\Omega=0\) [2402.18561].

This dynamical viewpoint does not coincide exactly with the equilibrium one. The same work reported that the transient localization threshold \(\alpha^*(s)\) agrees well with equilibrium expectations only for sufficiently small \(s\), while at larger \(s\) the transient and equilibrium critical couplings and apparent critical exponents diverge [2402.18561]. A later finite-temperature extension found that increasing temperature weakens localization and coherence in different ways: the incoherent region expands toward weaker coupling, whereas the transient localization threshold shifts only weakly and, where resolvable, tends to move to higher coupling at small \(s\) [2509.02345]. That study also emphasized that coherence and localization are distinct dynamical properties, since a state can be delocalized yet incoherent [2509.02345].

A further refinement is the “hidden” or pseudo-coherent phase. Using QUAPI/TEMPO, strong-coupling sub-Ohmic dynamics was shown to contain three regimes—coherent, incoherent, and pseudo-coherent—rather than a simple coherent/incoherent dichotomy [2208.10313]. In the pseudo-coherent regime the polarization shows a single minimum without a subsequent maximum, and the timescale of that feature scales as \(t_{\min}\propto 1/\omega_c\), indicating that the spin is slaved to oscillatory bath dynamics rather than executing intrinsic coherent oscillations [2208.10313].

Initial preparation also matters sharply. A time-dependent Davydov-\(D_1\) study at \(s=0.25\) found that under a polarized bath initial condition coherent oscillations persisted up to \(\alpha=0.3\approx 13\alpha_c\), whereas under a factorized initial condition a coherent-incoherent transition occurred at \(\alpha_{\rm CI}^{(f)}\approx 0.1\) [1302.1682]. This establishes that “strong-coupling coherence” in the deep sub-Ohmic regime is not a preparation-independent statement.

Recent entanglement-based work further complicates the picture. A TTN-TDVP-PS study found that the population-based coherent, incoherent, and pseudo-coherent regimes do not map one-to-one onto distinct stationary entanglement phases; instead, the stationary spin entropy defines a single entropy ridge in the \((s,\alpha)\) plane, which follows the population-based boundary only at small \(s\) and remains single-valued inside the incoherent region at larger \(s\) [2606.20313].

## 4. Decoherence, finite temperature, and equilibrium dynamical response

Sub-Ohmic environments are exceptionally effective at dephasing because the infrared sector couples strongly to the nonoscillatory parts of the reduced dynamics. A consistent perturbative treatment showed that the standard Markov approximation fails for the kernel \(I_4\) in the sub-Ohmic regime, replacing simple exponential dephasing by a nonexponential factor \(e^{-K(t)}\) and removing the spurious instantaneous-dephasing pathology of earlier weak-coupling formulas [1408.1061]. In the \(T=0\), \(n\to 0\) limit this treatment gives
\[
\frac{T_2}{T_1}=\frac{2\sin^2\theta}{2-\sin^2\theta},
\]
so arbitrarily small \(T_2/T_1\) ratios arise as the bath coupling becomes predominantly longitudinal [1408.1061]. For \(0<n<1\), the off-diagonal coherence decays as \(\exp[-C\,t^{1-n}]\), i.e. with nonanalytic stretched-exponential-type behavior rather than a simple exponential [1408.1061].

Equilibrium finite-temperature dynamics adds another layer. Full-density-matrix NRG calculations of the symmetrized correlation function \(C(\omega)\) found a thermal peak at
\[
\omega_T\sim T,
\]
with \(C(\omega,T)\) merging with the \(T=0\) result for \(\omega\gg \omega_T\) and deviating strongly from the zero-temperature \(\omega^{\pm s}\) laws for \(\omega\ll \omega_T\) [2010.07471]. The same work interpreted finite temperature as a strong infrared regulator that can mask the zero-temperature crossover scale near criticality [2010.07471].

Bias does not remove all universal low-frequency structure. For the biased sub-Ohmic model, bosonic NRG found that
\[
C(\omega)\propto \omega^s
\]
throughout the biased parameter space, including the biased strong-coupling regime, except at the special point \((\alpha=\alpha_c,\epsilon=0)\), where
\[
C(\omega)\propto \omega^{-s}
\]
[1701.05831]. The same study also supported the generalized Shiba relation
\[
C(\omega)\propto \chi^2\omega^s
\]
over a wide range of parameters [1701.05831].

Equilibrium dynamics near criticality has also been analyzed through the full set of spin components. A combined bosonic-NRG and Majorana-diagrammatics study concluded that the bosonic self-energy is essential for the description of critical fluctuations, but that many-body vertex corrections must also be included to obtain quantitative agreement with numerical simulations [1106.2655]. This work suggested that long-time out-of-equilibrium dynamics beyond the Bloch-Redfield regime deserves reconsideration in dissipative critical systems [1106.2655].

## 5. Numerical and analytical methodologies

Because the sub-Ohmic model combines strong infrared singularity, large bosonic displacements, and long memory times, no single method is uniformly reliable across all observables and parameter regimes. The main methodological families are summarized below.

| Method family | Representative papers | Main use |
|---|---|---|
| Bosonic NRG / FDM-NRG | [1111.6623], [2010.07471], [1701.05831] | Critical susceptibilities, equilibrium spectra, biased and unbiased low-frequency laws |
| Inchworm QMC | [2402.18561], [2509.02345] | Transient dynamical phase diagrams and finite-temperature real-time polarization |
| Variational coherent-state methods | [1102.1483], [2309.00797] | Ground-state criticality, critical couplings, mean-field versus non-mean-field behavior |
| MPS/DMRG-type approaches | [1211.3464], [1402.5478] | Localized phase with large occupations, chain representations, multi-bath generalizations |
| Path-integral and tensor-network real-time methods | [2208.10313], [2606.20313] | Hidden dynamical phases and entanglement structure |

Several method-specific developments are particularly tied to the sub-Ohmic problem. An MPS representation without explicit local Hilbert-space truncation was introduced specifically to cope with the very large boson occupations in and beyond the localized phase, and it reproduced the mean-field exponent \(1/2\) while allowing extrapolation of infinite-chain critical couplings [1211.3464]. A single-mode approximation constructed from a rotating-wave transformation and NRG-inspired transformations yielded the classical exponents
\[
\beta=\frac12,\quad \delta=3,\quad \gamma=1,\quad x=\frac12,\quad y_t^*=\frac12
\]
for \(0<s<1/2\), and was used to argue that the original bosonic NRG mishandles the crossover temperature \(T^*\) in that regime [1212.1889].

On the dynamical side, a variational surface-hopping algorithm built on the Davydov \(D_1\) ansatz was designed to treat coherent and incoherent population evolution within a single framework, and its hopping rates follow Marcus-theory-like \(\Delta^2\) scaling more closely than conventional semiclassical surface hopping [1306.2865]. This suggests that the sub-Ohmic model has also become a testbed for hybrid quantum-classical algorithm design, not only for critical phenomenology.

## 6. Generalizations, engineered realizations, and broader significance

Several extensions move beyond the canonical one-bath, one-coupling-direction problem. A two-bath model with one bath coupled diagonally to \(\sigma_z\) and another off-diagonally to \(\sigma_x\) was shown, by DMRG with optimized boson basis, to exhibit a second-order phase transition in the deep sub-Ohmic regime between two different doubly degenerate phases rather than the standard delocalized-localized transition [1402.5478]. In a related generalized model with both diagonal and off-diagonal couplings, a Davydov-\(D_1\) analysis found a discontinuous first-order transition between a zero-magnetization state and a finite-magnetization state, while the special case \(s=\bar s\) admits a continuous crossover from a single localized phase to a doubly degenerate localized phase [1310.1548].

The most explicit physical implementation to date is the analog quantum-simulation proposal based on color centers in free-standing h-BN membranes. There the spin-motion coupling is generated by a magnetic-field gradient, and the membrane geometry and tensile strain control the bath exponent. The proposal accesses \(-1\le s\le 0\), with \(s=-1\) realizing \(1/f\) noise and \(s=0\) realizing white noise; in those regimes the calculated dynamics shows coherence revivals at periods set by the bath characteristic frequency and coherent localization of the spin polarization [1808.02916]. Because \(1/f\) noise and white noise are among the most important decoherence sources in solid-state qubits, this platform was proposed not only as a simulator of dissipative many-body physics but also as a testbed for decoherence mechanisms in solid-state devices [1808.02916].

Taken together, these developments show that the sub-Ohmic spin-boson model is not merely a special case of the spin-boson Hamiltonian. It is the regime in which infrared bath structure becomes decisive: for criticality, for the distinction between localization and incoherence, for the validity or failure of Markovian approximations, and for the design of numerical and analog tools capable of resolving strong non-Markovian quantum dissipation [1111.6623, 2402.18561].

Source: https://www.emergentmind.com/topics/sub-ohmic-spin-boson-model