Non-equilibrium cotunneling in interacting Josephson junctions

We investigate non-equilibrium transport through interacting superconducting nanojunctions using a Liouville space approach. The formalism allows us to study finite gap effects, and to account for both quasiparticle and Cooper pair tunneling. With focus on the weak tunneling limit, we study the stationary dc and ac current up to second order (cotunneling) in the hybridization energy. We identify the characteristic virtual processes that yield the Andreev and Josephson current and obtain the dependence on the gate and bias voltage for the dc current, the critical current and the phase-dependent dissipative current. In particular, the critical current is characterized by regions in the stability diagram in which its sign changes from positive to negative, resulting in a multitude of 0- transitions. The latter signal the interplay between strong interactions and tunneling at finite bias.