Models of electroweak symmetry breaking with more than a single doublet of Higgs scalars contain a neutral pseudoscalar boson. The production of such a pseudoscalar in hadron collisions proceeds primarily via gluon fusion through a top-quark loop (except for those models in which the pseudoscalar coupling to bottom quarks is strongly enhanced). We compute the QCD corrections to this process in the heavy-quark limit, using an effective Lagrangian derived from the axial anomaly.

9 pages, (BNL number added, 1 typo corrected, PHYZZX format, 4 figures not included, available on request), BNL-49061