Multi-point correlation functions in the boundary XXZ chain at finite temperature

We consider multi-point correlation functions in the open XXZ chain with longitudinal boundary fields and in a uniform external magnetic field. We show that, at finite temperature, these correlation functions can be written in the quantum transfer matrix framework as sums over thermal form factors. More precisely, and quite remarkably, each term of the sum is given by a simple product of usual matrix elements of the quantum transfer matrix multiplied by a unique factor containing the whole information about the boundary fields. As an example, we provide a detailed expression for the longitudinal spin one-point functions at distance $m$ from the boundary. This work thus solves the long-standing problem of setting up form factor expansions in integrable models subject to open boundary conditions.