<p>Double parton distributions at small distances between the two partons are dominated by a mechanism in which the two observed partons originate from the splitting of a single parton. This contribution can be computed in terms of single-parton distributions and perturbative splitting kernels. We demonstrate that two-loop corrections to these kernels can have a substantial quantitative impact and considerably improve the stability of predictions for double parton scattering. We also consider the impact of heavy quark masses in the two-loop splitting kernels in an approximate manner.</p>