  Recent data from the CDF collaboration on the production of spin-triplet
bottomonium states at the Tevatron p \bar p collider are analyzed within the
NRQCD factorization formalism. The color-singlet matrix elements are determined
from electromagnetic decays and from potential models. The color-octet matrix
elements are determined by fitting the CDF data on the cross sections for
Upsilon(1S), Upsilon(2S), and Upsilon(3S) at large p_T and the fractions of
Upsilon(1S) coming from chi_b(1P) and chi_b(2P). We use the resulting matrix
elements to predict the cross sections at the Tevatron for the spin-singlet
states eta_b(nS) and h_b(nP). We argue that eta_b(1S) should be observable in
Run II through the decay eta_b -> J/psi + J/psi.
