We present the differential predictions for the rapidity distribution of pseudoscalar Higgs boson through gluon fusion at the LHC. These results are obtained taking into account the soft-virtual (SV) as well as the next-to-soft virtual (NSV) resummation effects to next-to-next-to-leading-logarithmic (NNLL¯) accuracy and matching them to the approximate fixed order next-to-next-to-leading-order (NNLOA) computation. We perform the resummation in two-dimensional Mellin space by limiting ourselves to the contributions only from gluon-gluon (gg) initiated channels. The NNLOA rapidity distribution of pseudoscalar Higgs is obtained by applying a ratio method on the NNLO rapidity distribution of the scalar Higgs boson. We also present the first analytical results of N3LO rapidity distribution of pseudoscalar Higgs at SV+NSV accuracy. The phenomenological impacts of NNLOA+NNLL¯ predictions for 13 TeV LHC are studied. We observe that for mA=125(700) GeV the SV+NSV resummation at the NNLL¯ level brings about 14.76% (11.48%) corrections to the NNLOA results at the central scale value of μR=μF=mA. Further, we find that the sensitivity to the renormalization scale gets improved substantially by the inclusion of NSV resummed predictions at NNLL¯ accuracy. Published by the American Physical Society 2024