When we consider an oscillator $a$ acting with a bath of oscillators $b_i$ with the interaction Hamiltonian reads $$H_{int}=\sum_{i}g_ia b_i^{\dagger}+g_i^*a^{\dagger}b_i,$$ with the free Hamiltonian: $H_{free}=\omega_0a^{\dagger}a+\sum_{i}\omega_ib^{\dagger}b$
It is a standard text book exercise (See, for example, Chapter 7 of this book) to derive the Born-Markov master equation for the density matrix of oscillator $a$ when the bath oscillator $b_i$ are all in thermal ensemble, i.e. $$\rho_{b_i}\propto exp(-\beta b^{\dagger}b)$$
However, I also tried a different scenario, when the bath are in a pure coherent ensemble instead of thermal ensemble, i.e.
$$\rho_{b_i}\propto |\alpha_i\rangle\langle \alpha_i|$$ $$\rho_{bath}=\otimes_i |\alpha_i\rangle\langle \alpha_i|$$
Where $|\alpha_i\rangle$ is the standard quantum harmonic oscillator coherent state. However, it seems that no similar Born-Markov master equation can be derived in this case due to some ill-defined integral in the Born-Markov master equation. Do anyone know any reference for this problem?
Edit: It seems that this problem of applying the Born-Markov approximation resides in the fact that there is some coherence time scale which is very large for the bath. I would actually be very grateful if some one can make this statement more explicit(e.g. by stating what are the things to calculate so that we can evaluate whether BM is a good approximation.)