Abstract
Superconducting devices subject to strong charging energy interactions and Coulomb blockade are one of the key elements for the development of nanoelectronics and constitute common building blocks of quantum computation platforms and topological superconducting setups. The study of their transport properties is non-trivial and some of their non-perturbative aspects are hard to capture with the most ordinary techniques. Here we present a matrix product state approach to simulate the real-time dynamics of these systems. We propose a study of their transport based on the analysis of the currents after quantum quenches connecting such devices with external leads. Our method is based on the combination of a Wilson chain construction for the leads and a mean-field BCS description for the superconducting scatterers. In particular, we employ a quasiparticle energy eigenbasis which greatly reduces their entanglement growth and we introduce an auxiliary degree of freedom to encode the device tot