Accurate and precise theory predictions for collider observables are instrumental for the discovery of new physics. To improve the quality of these predictions, increasingly complex contributions to scattering amplitudes must be evaluated. In certain kinematic regimes this approach no longer suffices due to the occurrence of large logarithmic contributions at every order, which spoil the perturbative series. One such regime is when the final state particle is produced at threshold, rendering additional radiation soft. In this instance, the threshold expansion can be applied to simplify the structure of the scattering process by keeping only leading-power (LP) terms in the soft momenta, allowing for the factorisation of the amplitude and the i...