Threshold of global behavior of solutions to a degenerate drift-diffusion system in between two critical exponents

We consider large time behavior of weak solutions to degenerate drift-diffusion system related to Keller–Segel system. (Formula Presented),where (Formula Presented) and (Formula Presented). There exist two critical diffusion exponents (Formula Presented) and for those cases, large time behavior of solutions is classified by the invariant norms of initial data. We consider the case of the intermediate exponent (Formula Presented) and classify the global existence and finite time blow up of weak solutions by the combination of invariant norms of initial data. Besides we show that the threshold value which classifies the behavior of solutions is characterized by the best possible constant of the modified Hardy–Littlewood–Sobolev inequality: (Formula Presented),where (Formula Presented) and it is given by the radial stationary solution of the system. Here the result is continuous analogue of the known critical cases (Blanchet et al., Calc Var Partial Diff Equ 35:133–168, 2009 and Ogawa, Disc Contin Dyn Syst Ser S 4:875–886, 2011). Analogous result has been obtained in the theory of nonlinear Schrödinger equations. The global behavior of the weak solution is also given and the solution converges to the self-similar Barenbratt solution as time parameter goes to infinity.