We develop a method to study quantum impurity models, small interacting quantum systems bilinearly coupled to an environment, in the presence of an additional Markovian quantum bath, with a generic nonlinear coupling to the impurity. We aim at computing the evolution operator of the reduced density matrix of the impurity, obtained after tracing out all the environmental degrees of freedom. First, we derive an exact real-time hybridization expansion for this quantity, which generalizes the result obtained in the absence of the additional Markovian dissipation and which could be amenable to stochastic sampling through diagrammatic Monte Carlo.
View Article and Find Full Text PDF