Mixed Grushin Heat Equations with Riesz-Potential Nonlinearities in Marcinkiewicz Spaces: Subcritical, Critical, and Supercritical Regimes

We study the nonlinear evolution equation ${\partial }_{t}u+(G+G^{\delta })u=I_{\alpha }(|u{|}^{\rho })\text{quadon}{\mathbb{R}}^{N+k},$ where $G$ is the nonnegative self adjoint realization of Grushin operator, $0<\delta <1$, and $I_{\alpha }$ is the potential operator with kernel $|z{|}^{-\alpha }$, $0<\alpha <N+k$. We develop a well-posedness theory for initial data in Marcinkiewicz spaces, allowing singular profiles outside the corresponding Lebesgue spaces. Using spectral calculus, subordination, and interpolation, we establish Lebesgue and Lorentz smoothing estimates for the mixed semigroup $e^{-t(G+G^{\delta })}$, which distinguish the second-order behaviour at short times from the fractional decay at large times. In the subcritical regime, considering $\frac{1}{\beta }=\frac{2}{Q}+1-\frac{\alpha }{d},\frac{1}{{\beta }_{\delta }}=\frac{2\delta }{Q}+1-\frac{\alpha }{d},$ where $Q$ is homogenous dimension, we prove local existence, uniqueness, Lipschitz dependence on the initial data, and a blow-up alternative in $L^{p,\infty }$ when $\rho <1+\frac{p}{\beta }.$ At the critical relation $\rho =1+p/\beta$, a Yamazaki-type estimate and Lorentz duality yield global mild solutions for sufficiently small data in $L^{p,\infty }$. In the supercritical regime $\rho >1+p/\beta$, well-posedness is recovered in higher-integrability Marcinkiewicz spaces $L^{q,\infty }$ satisfying $\beta (\rho -1)<q<{\beta }_{\delta }(\rho -1),$ with global existence for sufficiently small initial data. In all cases, the solutions attain their initial data in the weak-$\ast$ sense. These results extend the Marcinkiewicz-space theory to mixed Grushin diffusion with a spatially nonlocal potential source and reveal the role of the two competing diffusion scales in determining the admissible integrability regimes.