Summary: In this paper we study the existence of a global solution for a diffusion problem of Kirchhoff type driven by a nonlocal integro-differential operator. As a particular case, we consider the following parabolic equation involving the fractional \(p\)-Laplacian: \[ \begin{cases} \partial_t u + [u]_{s,p}^{(\lambda - 1) p}(-\Delta)_p^s u = |u|^{q-2}u,\quad & \text{in}\; \Omega \times \mathbb{R}^+\text{,}\quad \partial_t u = \partial u/\partial t\text{,}\\ u(x,0) = u_0(x)\text{,}\quad & \text{in}\; \Omega \text{,} \\ u(x,t)=0\text{,}\quad & \text{in}\; (\mathbb{R}^N \setminus \Omega) \times \mathbb{R}_0^+\text{,}\end{cases} \] where \([u]_{s,p}\) is the Gagliardo \(p\)-seminorm of \(u\), \(\Omega \subset \mathbb{R}^N\) is a bounded domain with Lipschitz boundary \(\partial \Omega\), \(p<q<Np/(N-sp)\) with \(1<p<N/s\) and \(s\in (0,1)\), \(1 \leq \lambda < N/(N-sp)\), \((-\Delta)_p^s\) is the fractional \(p\)-Laplacian. Under some appropriate assumptions, we obtain the existence of a global solution for the problem above by the Galerkin method and potential well theory. It is worth pointing out that the main result covers the degenerate case, that is the coefficient of \((-\Delta)_p^s\) can vanish at zero.