The production cross section for pseudo-scalar Higgs bosons at hadron colliders is computed at next-to-next-to-leading order (NNLO) in QCD. The pseudo-scalar Higgs is assumed to couple only to top quarks. The NNLO effects are evaluated using an effective lagrangian where the top quarks are integrated out. The NNLO corrections are similar in size to those found for scalar Higgs boson production.