  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.
