Very weak solutions to degenerate parabolic double-phase systems

We prove a local self-improving property for the gradient of very weak solutions to degenerate parabolic double-phase systems. The result is based on a reverse Hölder inequality with constants that are independent of the solution. Delicate methods are required to avoid a self-referential argument. In particular, we develop a new phase analysis method.