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 pT distribution matched to NNLOA+NNLL, valid for the whole pT 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. © 2018, The Author(s).