  In this article we have studied the transverse momentum distribution of the
pseudo-scalar Higgs boson at the Large Hadron Collider (LHC). The small $\pt$
region which provides the bulk of the cross section is not accessible to fixed
order perturbation theory due to the presence of large logarithms in the
series. Using the universal infrared behaviour of the QCD we resum these large
logarithms up to next-to-next-to-leading logarithmic (NNLL) accuracy. We
observe a significant reduction in theoretical uncertainties due to the
unphysical scales at NNLL level compared to the previous order. We present the
$p_T$ distribution matched to NNLO$_A$+NNLL, valid for the whole $p_T$ region
and provide a detailed phenomenological study in the context of both 14 TeV and
13 TeV LHC using different choices of masses, scales and parton distribution
functions which will be useful for the search of such particle at the LHC in
near future.
