“…We are now in a position to derive a qubit-energy hybrid master equation for the partially reduced density ρ(E, t) by generalizing the usual steps of a weak coupling Born-Markov treatment [11]. An alternative route in terms of the Nakajima-Zwanzig projection operator technique [33,34] is presented in appendix A, see also [26,27,30,31] for similar applications and the recent work [35]. We will not perform the Davies weak coupling limit prescription here explicitly, as it is the equivalent to the Born-Markov treatment we present below.…”