D-wave Charmonium Production in e^+e^- Annihilation at sqrt{s}=10.6 GeV

We calculate the D-wave charmonium production in $e^+e^-$ annihilation at BELLE and BABAR at $\sqrt{s}=10.6$ GeV using the nonrelativistic QCD factorization formalism, including color singlet and color octet contributions. We analyze the contributions of various processes for $\delta_J(J=1,2)$ production. We find that both the color singlet and color octet channels may give substantial contributions, and the production rates are estimated to be $\sigma (\delta_1)\simeq 0.043-0.16pb$ and $\sigma (\delta_2)\simeq 0.094-0.29pb$, which, however, are very sensitive to the choice of the color octet matrix elements of the D-wave charmonium states. The measurement of D-wave charmonium production at BABAR and BELLE in the future will be very helpful to test the color octet mechanism and to determine the color octet matrix elements.