How can I apply the energy dissipation form of the Spin-Boson model to both time-dependent and time-independent cases in my code?

[nelyasi]

I am wondering about the exact approach to use the energy dissipation form of the Spin-Boson model, as described in the tutorial, for both time-dependent and time-independent forms. Specifically, in the interaction Hamiltonian, instead of considering the system operator as ( \sigma_z ), I want to use the raising and lowering operators ( \sigma_+ ) and ( \sigma_- ), such that:

[
H_{I} = \frac{1}{2}\sum_{k}g_{k}(\sigma_{-}a^{\dag} + \text{hc})
]

Can I write the code like this?

sigma_p = oqupy.operators.sigma("+")
bath = oqupy.Bath(0.5 * sigma_p, correlations)

[piperfw]

Hi @nelyasi, apologies for not seeing this Q before now (GitHub doesn't seem to notify me of discussions). Unfortunately OQuPy only supports couplings of the form $S \sum_k \left( g_k a_k + \text{H.c.}\right)$ where $S$ is the system operator. There may be other ways for you to model the Jaynes-Cummings type interaction however; see #141

It does seem like we may want to add an explanation of this restriction to the tutorial.