  We present the first results on the production of pseudo-scalar through gluon
fusion at the LHC to N$^3$LO in QCD taking into account only soft gluon
effects. We have used the effective Lagrangian that describes the coupling of
pseudo-scalar with the gluons in the large top quark mass limit. We have used
recently available quantities namely the three loop pseudo-scalar form factor
and the third order universal soft function in QCD to achieve this. Along with
the fixed order results, we also present the process dependent resummation
coefficient for threshold resummation to N$^3$LL in QCD. Phenomenological
impact of these threshold N$^3$LO corrections to pseudo-scalar production at
the LHC is presented and their role to reduce the renormalisation scale
dependence is demonstrated.
