Extended Hybridization Expansion Solver for Impurity Models with Retarded Interactions

We extend the continuous-time hybridization expansion solver to a general form, where the hybridization function and retarded interaction are treated on equal footing. Correlation functions can now be directly obtained via functional derivatives with respect to the bosonic propagators, similar to the measurement of Green's functions. We devise a combinatorial scheme of measuring the correlation function, whose efficiency partially emulates that of the Green's function measurement. The algorithm and numerical methods are validated through application to an impurity model involving both electron-phonon coupling and exchange interactions, a case where the previous hybridization expansion algorithm is not applicable. Our improvement of the hybridization expansion solver promotes its applicability in studies of electron-phonon coupling, the extended dynamical mean field theory, and the dual boson method.