We report on results we recently obtained in Hebey and Thizy [11, 12] for critical stationary Kirchhoff systems in closed manifolds. Let $(M^n,g)$ be a closed $n$-manifold, $n\ge 3$. The critical Kirchhoff systems we consider are written as$$\left(a+b\sum _{j=1}^p{\int }_M{\left|\nabla u_j\right|}^{2}d{v}_g\right){\Delta }_g{u}_i+\sum _{j=1}^p{A}_{ij}{u}_j={\left|U\right|}^{2^{☆}-2}{u}_i$$for all $i=1,\cdots ,p$, where ${\Delta }_g$ is the Laplace-Beltrami operator, $A$ is a $C^1$-map from $M$ into the space $M_s^p\left(ℝ\right)$ of symmetric $p\times p$ matrices with real entries, the $A_{ij}$’s are the components of $A$, $U=(u_1,\cdots ,u_p)$, $\left|U\right|:M\to ℝ$ is the Euclidean norm of $U$, $2^{☆}=\frac{2n}{n-2}$ is the critical Sobolev exponent, and we require that $u_i\ge 0$ in $M$ for all $i=1,\cdots ,p$. We...