The regularity of the boundary of a multidimensional aggregation patch

We consider solutions to the aggregation equation with Newtonian potential where the initial data are the characteristic function of a domain with boundary of class $C^{1+\gamma}$ ,$0<\gamma<1$. Such initial data are known to yield a solution that, going forward in time, retains a patch-like structure with a constant time-dependent density inside an evolving region, which collapses on itself in a finite time, and which, going backward in time, converges in an $L^1$ sense to a self-similar expanding ball solution. In this work, we prove $C^{1+\gamma}$ regularity of the domain's boundary on the time interval on which the solution exists as an $L^\infty$ patch, up to the collapse time going forward in time and for all finite times going backward in time.