Strong coupling in thermoelectric nanojunctions: a reaction coordinate framework

We study a model of a thermoelectric nanojunction driven by vibrationally-assisted tunneling. We apply the reaction coordinate formalism to derive a master equation governing its thermoelectric performance beyond the weak electron-vibrational coupling limit. Employing full counting statistics we calculate the current flow, thermopower, associated noise, and efficiency without resorting to the weak vibrational coupling approximation. We demonstrate intricacies of the power-efficiency- precision trade-off at strong coupling, showing that the three cannot be maximised simultaneously in our model. Finally, we emphasise the importance of capturing non-additivity when considering strong coupling and multiple environments, demonstrating that an additive treatment of the environments can violate the upper bound on thermoelectric efficiency imposed by Carnot.