Relativistic corrections to exclusive χcJ+γ production from e+e- annihilation

We calculate in the nonrelativistic QCD (NRQCD) factorization framework all leading relativistic corrections to the exclusive production of χcJ+γ in e+e- annihilation. In particular, we compute for the first time contributions induced by octet operators with a chromoelectric field. The matching coefficients multiplying production long distance matrix elements (LDMEs) are determined through perturbative matching between QCD and NRQCD at the amplitude level. Technical challenges encountered in the nonrelativistic expansion of the QCD amplitudes are discussed in detail. The main source of uncertainty comes from the not so well known LDMEs. Accounting for it, we provide the following estimates for the production cross sections at s=10.6 GeV: σ(e+e-→χc0+γ)=(1.4±0.3) fb, σ(e+e-→χc1+γ)=(15.0±3.3) fb, and σ(e+e-→χc2+γ)=(4.5±1.4) fb.