NRQCD Analysis of Bottomonium Production at the Tevatron

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.