  We perform a perturbative QCD analysis of the gluonic gravitational form
factors (GFFs) of the proton and pion at large momentum transfer. We derive the
explicit factorization formula of the GFFs in terms of the distribution
amplitudes of hadrons. At the leading power, we find that
$A_g^\pi(t)=C_g^\pi(t)\sim 1/(-t)$ for pion, $A_g^p(t)\sim 1/(-t)^2$ and
$C_g^p(t)\sim \ln^2(-t/\Lambda^2)/(-t)^3$ for proton, respectively, where $t$
is the momentum transfer and $\Lambda$ a non-perturbative scale to regulate the
endpoint singularity in $C_g^p$ calculation. Our results provide a unique
perspective of the momentum dependence of the GFFs and will help to improve our
understanding of the internal pressure distributions of hadrons.
