  This article proves an interior gradient higher integrability result for weak
solutions to the parabolic double phase problems. The prototype equation for
the parabolic double phase problem of $p$-Laplace type reads as $$u_t - {\rm
div} \left(|\nabla u|^{p-2}\nabla u+ a(z)|\nabla u|^{q-2}\nabla u\right)=0$$
where $\frac{2n}{n+2}<p\leq q < \infty$ and the coefficient $a(z)$ is a
non-negative H\"{o}lder continuous function on $\Omega_T=\Omega \times (0, T),$
$\Omega \subset \mathbb{R}^{n}.$ Recently, this problem has been studied by
Kim, Kinnunen and Moring [19] for the degenerate case and Kim and
S\"{a}rki\"{o} [20] for the singular case. We introduce a new intrinsic scaling
that can handle the degenerate and the singular case simultaneously. This
scaling can also be used to obtain similar results in the variable exponent
case.
