Biased Spin-Boson Model

• The biased spin-boson model is a dissipative two-level system where finite bias competes with bath-induced renormalization to govern quantum dynamics.
• It reveals how variations in the bath spectral density and explicit bias drive crossovers from coherent to localized behavior in quantum systems.
• Advanced variational and nonperturbative methods are used to address critical exponent challenges and numerical artifacts in modeling its equilibrium and nonequilibrium regimes.

The biased spin-boson model is a dissipative two-state system in which a two-level degree of freedom with tunneling amplitude $\Delta$ and static bias $\epsilon$ is longitudinally coupled to a bosonic environment. In its standard form,

$$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$

with bath spectral density $J(\omega)$, it provides a canonical setting for quantum dissipation, decoherence of two-level systems, open-system dynamics, and impurity criticality. Finite bias explicitly favors one $\sigma_z$ state and breaks the parity symmetry present at $\epsilon=0$, so several zero-bias phase transitions are replaced by crossover structures once $\epsilon\neq 0$ (Zheng et al., 2017, Nazir et al., 2012, Liu et al., 2016).

1. Hamiltonian, conventions, and bath structure

The basic biased Hamiltonian is used directly in several works on equilibrium, dynamics, and transport,

$$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$

where $\Delta$ is the bare tunneling amplitude, $\epsilon$ is the bias field, and the bath couples through the same $\epsilon$0 operator that defines the localized-state population (Zheng et al., 2017). This longitudinal coincidence is central: the external bias and the dissipative dressing act in the same channel, so bias and bath renormalization do not simply add; they compete and reinforce each other in a strongly model-dependent way.

The bath is usually encoded by a power-law spectral density. A common sub-Ohmic/Ohmic choice is

$\epsilon$1

or, with an exponential cutoff,

$\epsilon$2

Here $\epsilon$3 is sub-Ohmic, $\epsilon$4 is Ohmic, and $\epsilon$5 is super-Ohmic; $\epsilon$6 is the dimensionless dissipation strength and $\epsilon$7 the high-frequency cutoff (Zheng et al., 2017, Liu et al., 2016, Bera et al., 2014). In sub-Ohmic baths, the enhanced low-frequency spectral weight makes the environment especially slow and history-dependent, which is why biased sub-Ohmic dynamics is typically much more non-Markovian than simple weak-coupling intuition suggests.

Several papers use rotated spin conventions. One field-theoretic formulation writes

$\epsilon$8

where the source field $\epsilon$9 is conjugate to the ordering component $$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$0 and therefore plays the role analogous to a bias field in conventional spin-boson notation; this is explicitly described as a spin-axis relabeling of the standard model (Vasin et al., 21 Jan 2025). A different extension introduces hybrid diagonal and off-diagonal system-bath coupling and, after a rotation, generates an effective bias

$$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$1

showing that “bias” can also emerge from basis rotation and tunneling renormalization rather than from an explicit bare $$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$2 term (Acharyya et al., 2020).

2. Equilibrium structure under finite bias

For the sub-Ohmic model at $$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$3, the standard picture is a continuous quantum phase transition at $$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$4 between delocalized and localized phases. Finite bias $$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$5 explicitly breaks the parity symmetry

$$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$6

so the sharp delocalized-localized transition is replaced by a smooth crossover from weak-coupling delocalized-like behavior to strong-coupling localized-like behavior (Zheng et al., 2017). This is one of the most important structural facts about the biased problem: in most equilibrium observables the nearby zero-bias critical fixed point still organizes the physics, but it no longer appears as a true phase transition once the symmetry-breaking field is finite.

The equilibrium dynamical correlation function

$$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$7

has a particularly simple low-frequency law in the biased sub-Ohmic problem. Bosonic NRG finds

$$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$8

for essentially the entire $$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$9 parameter space, including the biased strong-coupling regime, with the only exception at the exact unbiased critical point $J(\omega)$0, where

$J(\omega)$1

The same study reports $J(\omega)$2 over a wide parameter range, supporting the generalized Shiba relation even for $J(\omega)$3 and $J(\omega)$4 (Zheng et al., 2017). Finite bias therefore cuts off the zero-bias infrared singularity and restores the $J(\omega)$5 continuum tail at sufficiently low frequency.

Near $J(\omega)$6, the biased crossover is controlled by the competition between the unbiased crossover scale

$J(\omega)$7

and a bias scale

$J(\omega)$8

At the critical coupling,

$J(\omega)$9

with

$\sigma_z$0

The regime $\sigma_z$1 is controlled by the nearby unbiased fixed point, whereas for $\sigma_z$2 the bias dominates, suppresses low-frequency continuum weight, and transfers spectral weight into the $\sigma_z$3 component associated with finite $\sigma_z$4 (Zheng et al., 2017).

A rotated-basis field theory expresses the same bias idea through the equation of state

$\sigma_z$5

at criticality, so that

$\sigma_z$6

In that formulation $\sigma_z$7 is the field conjugate to the order parameter and is explicitly identified as the analog of the conventional bias field after spin-axis relabeling (Vasin et al., 21 Jan 2025). By contrast, a mean-field analysis of boson-state truncation shows that the finite-bias exponent can be strongly distorted for $\sigma_z$8 at finite $\sigma_z$9: the true mean-field $\epsilon=0$0 crosses over to an apparent truncation-dominated

$\epsilon=0$1

so low-bias critical fits in truncated calculations are artifact-sensitive unless the $\epsilon=0$2 limit is controlled explicitly (Hou et al., 2010).

3. Nonequilibrium dynamics and bath-induced bias renormalization

Real-time studies make the biased model qualitatively richer than a simple detuned damped two-level system. For the deep sub-Ohmic case $\epsilon=0$3, TDVMPS simulations of

$\epsilon=0$4

show three distinct nonequilibrium regimes. At weak coupling, the spin exhibits damped oscillations close to the bare precession frequency

$\epsilon=0$5

but the late-time behavior departs from a Markovian master equation because slow bath modes acquire a small net displacement and act back on the spin as an additional $\epsilon=0$6-field. At strong coupling, the environment rapidly forms a localized bosonic cloud, strongly suppresses tunneling, and the spin freezes near its initial state. Between them lies a nonperturbative intermediate regime with “sequential, time-retarded” dynamics: the spin first behaves approximately as in weak coupling and is only later dragged into bath-dominated evolution once the slow modes have sufficiently reorganized (Gonzalez-Ballestero et al., 2017).

That study is unusual in explicitly following both spin and environment. The projected mode displacements

$\epsilon=0$7

distinguish two mechanisms. If $\epsilon=0$8 and $\epsilon=0$9 drift together, the bath develops a net displacement and hence an effective dynamical bias field. If they remain strongly anticorrelated, the bath response is predominantly polaronic and suppresses tunneling instead (Gonzalez-Ballestero et al., 2017). The interpretation is that the interaction term generates a time-dependent environmental contribution to the $\epsilon\neq 0$0-field proportional to the mode displacements; this is the paper’s microscopic explanation of environment-induced bias renormalization.

In the weak-coupling Ohmic problem at arbitrary bias, a consistent expansion beyond Bloch-Redfield gives a complementary picture. The renormalized tunneling and Rabi frequency are

$\epsilon\neq 0$1

and the three nonstationary modes are

$\epsilon\neq 0$2

with

$\epsilon\neq 0$3

The time evolution has the generic structure

$\epsilon\neq 0$4

but the preexponential functions are not constants: they contain logarithmic terms in $\epsilon\neq 0$5, branch-cut contributions, and power-law crossovers that are absent from simple Bloch-Redfield theory. After real-time RG resummation, the oscillating modes cross over from exponent $\epsilon\neq 0$6 at exponentially small times to the bias-dependent exponent

$\epsilon\neq 0$7

at exponentially large times (Lindner et al., 2018). Bias therefore changes not only the oscillation frequency but also the long-time asymptotic structure of the decay.

An important reference point is the symmetric sub-Ohmic problem. Numerically exact TEMPO calculations at $\epsilon\neq 0$8 find a three-region dynamical diagram with coherent, incoherent, and strong-coupling “pseudo-coherent” behavior, where the oscillation-like turnaround is controlled by the bath timescale $\epsilon\neq 0$9 rather than by intrinsic tunneling (Otterpohl et al., 2022). Those results are direct only for the unbiased case. A plausible implication is that similar bath-driven turning points can survive under bias, but their interpretation must be reformulated because the relaxation target is no longer symmetry-related.

4. Variational and nonperturbative descriptions of the biased ground state

The biased Ohmic ground state is poorly captured by symmetric polaron pictures. A generalized biased ansatz introduces asymmetric bath displacements

$$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$0

with

$$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$1

leading to the self-consistent renormalized tunneling

$$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$2

For the Ohmic spectrum this gives

$$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$3

Because $$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$4 remains nonzero when $$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$5, the infrared divergence is cut off and $$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$6 stays finite and continuous. The same construction produces an environment-induced bias renormalization

$$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$7

and a smooth ground-state magnetization

$$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$8

with

$$H=-\frac{\Delta}{2}\sigma_x+\frac{\epsilon}{2}\sigma_z+\sum_i \omega_i a_i^\dagger a_i+\frac{\sigma_z}{2}\sum_i \lambda_i(a_i+a_i^\dagger),$$9

This removes the unphysical discontinuous collapse of $\Delta$0 and the spurious jump in $\Delta$1 produced by Silbey-Harris theory for biased systems, and matches Bethe-Ansatz results much more closely in the scaling limit (Nazir et al., 2012).

A more systematic ground-state description is the multi-polaron coherent-state expansion,

$\Delta$2

Its central message is that finite bias does not eliminate antipolarons. High-frequency modes remain polaronic, intermediate-energy modes still develop tunneling-stabilized antipolaronic displacements, and only the low-frequency asymptote changes qualitatively: in the biased case all displacements saturate to the same finite value rather than vanishing as $\Delta$3. The resulting observables

$\Delta$4

$\Delta$5

reproduce Bethe-Ansatz results very accurately; typically $\Delta$6, and sometimes $\Delta$7, is sufficient in the biased regimes studied (Bera et al., 2014).

The displaced-Fock-state improvement of Silbey-Harris pursues the same correction from another angle. Starting from the asymmetric biased ansatz

$\Delta$8

it adds orthogonal displaced Fock sectors order by order. In the biased case, this removes the discontinuous magnetization jump reported in earlier SH-based work and yields smooth, rapidly convergent magnetization curves for both $\Delta$9 and $\epsilon$0; for $\epsilon$1, convergence is essentially reached already by the third-order ansatz (He et al., 2018).

5. Transport, generated bias, and exact structures

In the nonequilibrium spin-boson model with two bosonic reservoirs at different temperatures,

$\epsilon$2

the bias $\epsilon$3 becomes a control parameter for energy transfer rather than only for local polarization. A polaron-transformed NEGF treatment gives a Landauer-like steady-state heat current

$\epsilon$4

and reveals a bias-induced nonmonotonic energy conductance in the intermediate-coupling regime. The physical explanation given is resonance: increasing bias changes the effective spin splitting relative to the thermal window of the baths, so conductance can first increase and later decrease as $\epsilon$5 is tuned (Liu et al., 2016).

Hybrid diagonal/off-diagonal coupling extends the same theme by showing that bias can be generated rather than imposed. Starting from

$\epsilon$6

a unitary rotation maps the model onto an ordinary spin-boson Hamiltonian with

$\epsilon$7

The paper distinguishes this rotated-frame effective bias from a separate “dynamically generated bias” produced by the renormalized tunneling amplitude $\epsilon$8. In the overdamped regime the long-time population has the form

$\epsilon$9

so that $\epsilon$00 even though no explicit bare $\epsilon$01 term was present in the original Hamiltonian (Acharyya et al., 2020).

The biased Ohmic model also has an exact-overlap and integrability sector. With Hamiltonian

$\epsilon$02

the bias is the local field $\epsilon$03, the renormalized impurity scale is $\epsilon$04, and the magnetization is

$\epsilon$05

For ground-state overlaps between two parameter sets, the orthogonality exponent is

$\epsilon$06

with generalized form

$\epsilon$07

This connects bias variation directly to Anderson orthogonality, and at the Toulouse point $\epsilon$08 the exact magnetization is

$\epsilon$09

(Lukyanov, 2015).

6. Limits, controversies, and open problems

Several technical and conceptual cautions recur across the literature. First, some apparent critical exponents are method-sensitive. In the mean-field spin-boson model, boson-state truncation is a strong relevant operator with respect to the Gaussian fixed point for $\epsilon$10, and the finite-bias law

$\epsilon$11

crosses over at finite $\epsilon$12 to the truncation-dominated

$\epsilon$13

This implies that the low-bias exponent $\epsilon$14 is especially vulnerable to numerical artifacts in truncated calculations (Hou et al., 2010).

Second, the status of Ohmic criticality remains nonuniform across formulations. A Schwinger-Keldysh/Majorana RG analysis claims a second-order quantum phase transition for both $\epsilon$15 and $\epsilon$16, with

$\epsilon$17

and explicitly states that this is contrary to the conventional quantum-classical mapping expectation of a Kosterlitz-Thouless transition in the Ohmic case (Vasin et al., 21 Jan 2025). By contrast, other Ohmic treatments retain the standard BKT picture, for example the hybrid-coupling study with

$\epsilon$18

and the integrable Ohmic discussion that places the Kosterlitz-Thouless transition at $\epsilon$19 (Acharyya et al., 2020, Lukyanov, 2015). These statements are not framed in identical conventions, but they mark a genuine interpretive tension in the present literature.

Third, several influential dynamical results are direct only in the unbiased limit. The TEMPO study of the sub-Ohmic “pseudo-coherent” phase is explicit that $\epsilon$20 throughout and that any extension to $\epsilon$21 is an inference rather than a reported result (Otterpohl et al., 2022). Likewise, exact SSE studies with added Markovian dephasing treat only the unbiased Ohmic model, so their frequency-renormalization formulas are baseline results rather than direct biased-model statements (Kamar et al., 2023).

Finally, rigorous mathematical results are narrower than the standard biased $\epsilon$22-coupled Hamiltonian. The small-coupling ground-state existence theorem for the massless spin-boson model applies to an off-diagonal coupling structure in the energy basis and is not a direct treatment of the conventional biased model; its main structural lesson is that the absence of diagonal coupling terms is what makes the infrared problem tractable in that setting (Hasler et al., 2010). The massive one-boson scattering analysis, although it already describes an energy-asymmetric two-level system with

$\epsilon$23

is likewise not written in the conventional $\epsilon$24 form, even though its resonance structure depends directly on the level asymmetry $\epsilon$25 (Ballesteros et al., 2018).

Taken together, these results support a stable qualitative picture. Finite bias removes the exact $\epsilon$26 symmetry, rounds zero-bias critical singularities into crossovers, produces a biased stationary or equilibrium polarization, and changes both infrared scaling and long-time dynamics. At the same time, the detailed location and meaning of “delocalized,” “localized,” “coherent,” or “pseudo-coherent” behavior remain strongly dependent on bath exponent, cutoff, initial preparation, and approximation scheme. The biased spin-boson model is therefore less a single solved system than a family of closely related impurity problems whose shared core is the competition among tunneling, explicit asymmetry, and bath-induced renormalization.