We expand the study of the pseudoscalar glueball and its first excited state by constructing an interaction Lagrangian which produces the two- and three-body decays of the pseudoscalar glueball, $J^{PC}=0^{-+}$, into the (pseudo)scalar and the excited (pseudo)scalar mesons as well as by constructing other two different chiral Lagrangians which describe the two- and three-body decays of the first excited pseudoscalar glueball, $J^{PC}=0^{*-+}$, into the (pseudo)scalar and the excited (pseudo)scalar mesons. We compute the decay channels for the ground state of a pseudoscalar glueball with a mass of $2.6$ GeV and for the first excited pseudoscalar glueball with a mass $3.7$ GeV, following predictions from Lattice QCD in the quenched approximation. These states and channels are within reach of PANDA experiment at the upcoming FAIR facility experiment and ongoing BESIII experiment. In our approach, the various branching ratios are a parameter-free prediction.