Search

The $X(6900)$ resonance, originally discovered by the \texttt{LHCb} collaboration and later confirmed by both \texttt{ATLAS} and \texttt{CMS} experiments, has sparked broad interests in the fully-charmed tetraquark states. Relative to the mass spectra and decay properties of fully-heavy tetraquarks, our knowledge on their production mechanism is still rather limited. In this work we investigate the inclusive production of fully-charmed $S$-wave tetraquarks at \texttt{LHC} within the nonrelativistic QCD (NRQCD) factorization framework. The partonic cross sections are computed at lowest order in $\alpha_s$ and velocity, while the long-distance NRQCD matrix elements are estimated from phenomenological potential models. We predict the differential $p_T$ spectra of various fully-charmed $S$-wave tetraquarks at the \texttt{LHC}, and compare with the results predicted from the fragmentation mechanism at large $p_T$ end., Comment: 12 pages, 5 figures, 2 tables. SDCs in .m attachment

The failure of observing the $e^+ e^- \to J/ψ+J/ψ$ events at $B$ factories to date is often attributed to the significant negative order-$α_s$ correction. In this work we compute the ${\cal O}(α^2_s)$ correction to this process for the first time. The magnitude of the next-to-next-to-leading order (NNLO) perturbative correction is substantially negative so that the standard NRQCD prediction would suffer from an unphysical, negative cross section. This dilemma may be traced in the fact that the bulk contribution of the fixed-order radiative corrections stems from the perturbative corrections to the $J/ψ$ decay constant. We thus implement an improved NRQCD factorization framework, by decomposing the amplitude into the photon-fragmentation piece and the non-fragmentation piece. With the measured $J/ψ$ decay constant as input, which amounts to resumming a specific class of radiative and relativistic corrections to all orders, the fragmentation-induced production rate can be predicted accurately and serves a benchmark prediction. The non-fragmentation type of the amplitude is then computed through NNLO in $α_s$ and at lowest order in velocity. Both the ${\cal O}(α_s)$ and ${\cal O}(α^2_s)$ corrections in the interference term become positive and exhibit a decent convergence behavior. Our finest prediction is $σ(e^+ e^- \to J/ψ+J/ψ)= 2.13^{+0.30}_{-0.06}$ fb at $\sqrt{s}=10.58$ GeV. With the projected integrated luminosity of 50 ${\rm ab}^{-1}$, the prospect to observe this exclusive process at \texttt{Belle} 2 experiment appears to be bright., 13 pages, 8 figures, 3 tables

In this work, we study the exclusive decay of $Υ$ into $J/ψ$ in association with $η_c$ ($χ_{c0,1,2}$). The decay widths for different helicity configurations are evaluated up to QCD next-to-leading order within the nonrelativistic QCD framework. We find that the QCD corrections notably mitigate the renormalization scale dependence of the decay widths for all the processes. The branching fraction of $Υ\rightarrow J/ψ+χ_{c1}$ is obtained as $3.73^{+5.10+0.10}_{-2.06-1.19}\times 10^{-6}$, which agrees well with the Belle measurement, i.e., ${\rm Br}(Υ\rightarrow J/ψ+χ_{c1})=(3.90\pm1.21\pm0.23)\times10^{-6}$. For the other processes, our results of the branching fractions are compatible with the upper limits given by the Belle experiments, except for $Υ(2S)\to J/ψ+χ_{c1}$, where some tension exists between theory and experiment. Having the polarized decay widths, we study the $J/ψ$ polarization, which turn out to be independent of any nonperturbative parameters. Further, according to our calculation, it is promising to measure all the processes at Super B factory thanks to the high luminosity., 13 pages, 2 figures, 3 tables, to match the version accepted for publication in PRL

In this paper, we evaluate, up to QCD next-to-next-to-next-to-leading order, the $B_c^*$ decay constant, which, within the nonrelativistic QCD (NRQCD) framework, is factorized as the product of the short-distance coefficient (SDC) and the long-distance matrix element. For the first time, the renormalization constant and the anomalous dimension for the NRQCD vector current composed of $c\bar{b}$, which are functions of the charm quark mass $m_c$, the bottom quark mass $m_b$ and the factorization scale $μ_Λ$, are obtained analytically at $\mathcal{O}(α_s^2)$ and $\mathcal{O}(α_s^3)$. The SDC is calculated up to $\mathcal{O}(α_s^3)$ in perturbative expansion. Explicitly, the $\mathcal{O}(α_s^2)$ correction to the SDC is analytically calculated in terms of logarithmic and polylogarithmic functions of $r\equiv m_b/m_c$, and the $\mathcal{O}(α_s^3)$ correction to the SDC is numerically calculated at a series of values of $r$, ranging from $2.1$ to $4.0$. Surprisingly, we find that the nontrivial part of the SDC at $\mathcal{O}(α_s^3)$ can be well estimated by a linear function of $r$. In addition, We find that the $\mathcal{O}(α_s^2)$ and $\mathcal{O}(α_s^3)$ corrections to the decay constant and decay width are considerable and very significant, which indicates a very poor convergence for the perturbative expansion., 15 pages, 3 figures, 2 tables, a minor typo corrected

We compute the next-to-next-to-leading-order (NNLO) perturbative corrections to the exclusive channels $e^+e^-\to J/ψ+χ_{cJ}$ ($J=0,1,2$) at $\sqrt{s}=10.58$ GeV within the nonrelativistic QCD (NRQCD) factorization framework. The validity of NRQCD factorization at ${\cal O}(α_s^2)$ has been confirmed for these double-charmonium exclusive production processes. We analyze the impact of the NNLO perturbative corrections on the polarized and unpolarized cross sections, as well as the $J/ψ$ angular distributions, which largely reduce the renormalization scale dependence but lower the NLO NRQCD predictions to some extent for $χ_{c0,1}$. With high numerical accuracy, our NNLO predictions for $σ(J/ψ+χ_{c1,2})$ are compatible with the upper limit of the \texttt{Belle} measurement. Although the finest prediction for $σ(J/ψ+χ_{c0})$ is consistent with both \texttt{Belle} and \texttt{BaBar} measurements within uncertainties, there still exists serious tension between the predicted and the measured profiles for the $J/ψ$ angular distribution. Regarding the bright observation prospect of the $e^+e^-\to J/ψ+χ_{c1,2}$ channels in \texttt{Belle 2} experiment, it is interesting to compare the future measurements with our NNLO predictions. The more accurate measurement of $e^+e^-\to J/ψ+χ_{0}$ at \texttt{Belle 2} will also help to clarify the long-standing puzzle of $J/ψ$ angular distribution., 17 pages, 4 tables, 2 figures. b quark loop contributions explicitly included

Inspired by the recent discovery of the $X(6900)$ meson at {\tt LHCb} experiment, we investigate the inclusive production rate of the $C$-odd fully-charmed tetraquarks associated with light hadrons at the $B$ factory within the nonrelativistic QCD (NRQCD) factorization framework. The short-distance coefficient is computed at lowest order in velocity and $\alpha_s$. Employing the diquark-antidiquark model to roughly estimate the long-distance NRQCD matrix elements, we predict the rate for inclusive production of the $1^{+-}$ $T_{4c}$ state and discuss the observation prospects at {\tt Belle 2} experiment., Comment: 11 pages, 3 figures

Within the nonrelativistic QCD (NRQCD) factorization framework, we have computed the ${\mathcal O}(��_s^2)$ corrections to the exclusive production of $P$-wave spin-triplet charmonia $��_{cJ}(J=0,1,2)$ accompanied by a hard photon at $B$ factories. For the first time, we have explicitly verified the validity of NRQCD factorization for exclusive $P$-wave quarkonium production at two-loop order. Unlike the $��_{cJ}$ electromagnetic decay processes, the $\mathcal{O}(��_s^2)$ corrections are found to be smaller than the $\mathcal{O}(��_s)$ corrections in all three channels $e^+e^-\to ��_{c0,1,2}+��$. In particular, the ${\mathcal O}(��_s^2)$ corrections appear moderate for $��_{c1}$ case, and insignificant for $��_{c0}$. In addition, the next-to-next-to-leading order (NNLO) predictions for the production rates of $��_{c0,1}+��$ are insensitive to the renormalization and factorization scales. All of these features may indicate that perturbative expansion in these two channels exhibits a decent convergence behavior. By contrast, both the ${\mathcal O}(��_s)$ and ${\mathcal O}(��_s^2)$ corrections to the $��_{c2}+��$ production cross section are sizable, which even reduces the Born order cross section by one order of magnitude after including the NNLO perturbative correction. Taking the values of the long-distance NRQCD matrix elements from nonrelativistic potential model, our prediction to $��_{c1}+��$ production rate is consistent with the recent {\tt Belle} measurement. The NNLO predictions for the $��_{c0,2}+��$ production rates are much smaller than that for $��_{c1}+��$, which seems to naturally explain why the $e^+e^-\to ��_{c0,2}+��$ channels have escaped experimental detection to date., typos corrected, some references added

The $X(6900)$ resonance, very recently discovered in the double-$J/\psi$ channel at LHCb experiment, has spurred intensive interest in unravelling the nature of the fully charmed tetraquark state. The aim of this paper is to present a model-independent theoretical framework to study the inclusive production of this novel species of exotic hadrons, the resonances composed of four heavy quark (commonly referred to as $T_{4c}$), at large $p_T$ in hadron collision experiments. Appealing to asymptotic freedom and the fact $m_c\gg \Lambda_{\rm QCD}$, we propose that the nonpertubative yet universal gluon-to-$T_{4c}$ fragmentation function, can be decomposed into the product of the perturbatively calculable short-distance coefficient and the long-distance NRQCD matrix elements. We compute the short-distance coefficient at lowest-order in $\alpha_s$ and velocity expansion. Adopting the diquark ansatz to roughly estimate those not-yet-known NRQCD matrix elements, together with the standard QCD factorization theorem, we predict the differential production rates for the $T_{4c/4b}(0^{++})$ and $T_{4c/4b}(2^{++})$ at large $p_T$ in $pp$ collision, which eagerly awaits the confrontation with the future LHC experiments., Comment: 9 pages, 1 table, 3 figures

Within the framework of nonrelativistic QCD (NRQCD) factorization formalism, we compute the helicity amplitude as well as the decay width of $��_{Q2}$ ($Q=c,b$) electromagnetic decay into two photons up to next-to-next-to-leading order (NNLO) in $��_s$ expansion. For the first time, we verify the validity of NRQCD factorization for the D-wave quarkonium decay at NNLO. We find that the $\mathcal{O}(��_s)$ and $\mathcal{O}(��_s^2)$ corrections to the helicity amplitude are negative and moderate, nevertheless both corrections combine to suppress the leading-order prediction for the decay width significantly. By approximating the total decay width of $��_{Q2}$ as the sum of those for the hadronic decay and the electric $E1$ transition, we obtain the branching ratios ${\rm Br}(��_{c2}\to 2��)\approx 5\times10^{-6}$ and ${\rm Br}(��_{b2}\to 2��)\approx 4\times10^{-7}$. To explore the potential measurement on $��_{Q2}$, we further evaluate the production cross section of $��_{Q2}$ at LHCb at the lowest order in $��_s$ expansion. With the kinematic constraint on the longitudinal rapidity $4.5>y>2$ and transverse momentum $P_T>(2-4)m_Q$ for $��_{Q2}$, we find the cross section can reach $2-50$ nb for $��_{c2}$, and $1-22$ pb for $��_{b2}$. Considering the integrated luminosity $\mathcal{L}=10\, {\rm fb}^{-1}$ at $\sqrt{s}=7$ TeV and $\sqrt{s}=13$ TeV, we estimate that there are several hundreds events of $pp\to ��_{c2}\to 2��$. Since the background is relatively clean, it is promising to reconstruct $��_{c2}$ through its electromagnetic decay. On the contrary, due to small branching ratio and production cross section, it is quite challenging to detect $��_{b2}\to 2��$ at LHCb., 14 pages, 1 figures, 2 tables

Within the nonrelativistic QCD (NRQCD) factorization framework, we compute the ${\mathcal O}(��_s v^2)$ corrections to the hadronic decay rate of vector quarkonia, exemplified by $J/��$ and $��$. Setting both the renormalization and NRQCD factorization scales to be $m_Q$, we obtain $��(J/��\to {\rm LH})= 0.0716\frac{��_s^3}{m_c^2} \langle \mathcal{O}_1({}^3S_1)\rangle_{J/��} [1-1.19��_s+(-5.32+3.03��_s)\langle v^2\rangle_{J/��}]$ and $��(��\to {\rm LH})= 0.0716\frac{��_s^3}{m_b^2}\langle\mathcal{O}_1({}^3S_1)\rangle_��[1-1.56��_s+(-5.32+4.61��_s)\langle v^2\rangle_��]$. We confirm the previous calculation of $\mathcal{O}(��_s)$ corrections on a diagram-by-diagram basis, with the accuracy significantly improved. For $J/��$ hadronic decay, we find that the ${\mathcal O}(��_sv^2)$ corrections are moderate and positive, nevertheless unable to counterbalance the huge negative corrections. On the other hand, the effect of ${\mathcal O}(��_sv^2)$ corrections for $��(nS)$ is sensitive to the $\mathcal{O}(v^2)$ NRQCD matrix elements. With the appropriate choice of the NRQCD matrix elements, our theoretical predictions for the decay rates may be consistent with the experimental data for $��(1S,2S)\to {\rm LH}$. As a byproduct, we also present the theoretical predictions for the branching ratio of $J/��(��)\to 3��$ accurate up to $\mathcal{O}(��_s v^2)$., 21 pages, 5 figures, 3 tables

In the paper, we present QCD predictions for $��_{c} + ��$ production at an electron-position collider up to next-to-next-to-leading order (NNLO) accuracy without renormalization scale ambiguities. The NNLO total cross-section for $e^{+}+e^{-}\to��+��_{c}$ using the conventional scale-setting approach has large renormalization scale ambiguities, usually estimated by choosing the renormalization scale to be the $e^+ e^-$ center-of-mass collision energy $\sqrt{s}$. The Principle of Maximum Conformality (PMC) provides a systematic way to eliminate such renormalization scale ambiguities by summing the nonconformal $��$ contributions into the QCD coupling $��_s(Q^2)$. The renormalization group equation then sets the value of $��_s$ for the process. The PMC renormalization scale reflects the virtuality of the underlying process, and the resulting predictions satisfy all of the requirements of renormalization group invariance, including renormalization scheme invariance. After applying the PMC, we obtain a scale-and-scheme independent prediction, $��|_{\rm NNLO, PMC}\simeq 41.18$ fb for $\sqrt{s}$=10.6 GeV. The resulting pQCD series matches the series for conformal theory and thus has no divergent renormalon contributions. The large $K$ factor which contributes to this process reinforces the importance of uncalculated NNNLO and higher-order terms. Using the PMC scale-and-scheme independent conformal series and the $\rm Pad\acute{e}$ approximation approach, we predict $��|_{\rm NNNLO, PMC+Pade} \simeq 21.36$ fb, which is consistent with the recent BELLE measurement $��^{\rm obs}$=$16.58^{+10.51}_{-9.93}$ fb at $\sqrt{s} \simeq 10.6$ GeV. This procedure also provides a first estimate of the NNNLO contribution., 8 pages, 5 figures. Discussions improved and references updated. Revised version to be published in JHEP

We study the $��_c$ photoproduction in $ep$ collisions in this paper. The short-distance coefficients for $c\bar{c}(^1S_0^{[1]})$, $c\bar{c}(^1S_0^{[8]})$, $c\bar{c}(^3S_1^{[8]})$, and $c\bar{c}(^1P_1^{[8]})$ photoproductions are evaluated at leading order in $��_s$ expansion, where the color-singlet contribution is achieved for the first time. We have carefully analyzed different kinematic distributions of the cross sections and found that the color-singlet contribution is considerably suppressed comparing with the color-octet parts. This feature renders the $��_c$ photoproduction process an ideal laboratory to test the color-octet mechanism in nonrelativistic QCD. By taking different sets of long-distance matrix elements, we have observed some apparently distinguishable predictions, which can be utilized to scrutinize the validity of these matrix elements., 6pages, 2 figures

The rare process of $W$ radiative decay into the $B_c$ meson can serve as an interesting play ground for testing two influential perturbative QCD approaches: non-relativistic quantum chromodynamics (NRQCD) and ligh-cone (LC) factorization. Since the $B_c$ meson is a genuine heavy quarkonium which is composed of two different species of heavy quarks, it is natural to employ NRQCD factorization to tackle this exclusive quarkonium production process; on the other hand, since the mass of the $W$ boson is much greater than that of the $B_c$, the collinear factorization is also a viable approach. To fully disentangle the contributions from several distinct scales $M_W$, $m_Q$, $m_Q v$ ($v$ represents the typical velocity of the $c/\bar{b}$ quarks inside $B_c$), and $��_{\rm QCD}$, we carry out a joint analysis that combines both light-cone and NRQCD factorization approaches, and compute the order-$��_s$ correction to $W \to B_c+��$. With the aid of the celebrated Efremov-Radyushkin-Brodsky-Lepage (ERBL) evolution equation, we resum the leading collinear logarithms $��_s \ln\! {M_W\over m_Q}$ to all orders in $��_s$. We also explore some phenomenological implications of our predictions., 16 pages, 3 figures

The associated production of Higgs boson with a muon pair, $e^+e^-\to ��^+��^- H$, is one of the golden channels to pin down the properties of the Higgs boson in the prospective Higgs factories exemplified by \textsf{CEPC}. The projected accuracy of the corresponding cross section measurement is about per cent level at CEPC. In this work, we investigate both ${\mathcal O}(��)$ weak correction and the ${\mathcal O}(����_s)$ mixed electroweak-QCD corrections for this channel, appropriately taking into account the effect of finite $Z^0$ width. The $��^+��^-$ invariant mass spectrum is also predicted. The mixed electroweak-QCD correction turns out to reach 1.5\% of the Born-order result, and thereby must be included in future confrontation with the data. We also observe that, after including higher-order corrections, the simplified prediction for the integrated cross section employing the narrow-width-approximation may deviate from our full result by a few per cents., v2, 15 pages, 5 figures, 3 tables. Typos corrected, references added, the version to appear in Chinese Physics C

The next-to-next-to-leading order (NNLO) pQCD prediction for the $����^* \to ��_c$ form factor was evaluated in 2015 using nonrelativistic QCD (NRQCD). A strong discrepancy between the NRQCD prediction and the BaBar measurements was observed. Until now there has been no solution for this puzzle. In this paper, we present a NNLO analysis by applying the Principle of Maximum Conformality (PMC) to set the renormalization scale. By carefully dealing with the light-by-light diagrams at the NNLO level, the resulting high precision PMC prediction agrees with the BaBar measurements within errors, and the conventional renormalization scale uncertainty is eliminated. The PMC is consistent with all of the requirements of the renormalization group, including scheme-independence. The application of the PMC thus provides a rigorous solution for the $����^* \to ��_c$ form factor puzzle, emphasizing the importance of correct renormalization scale-setting. The results also support the applicability of NRQCD to hard exclusive processes involving charmonium., 5 pages, 4 figures. References updated. To be published in Phys.Rev.D

The associated production of Higgs boson with a hard photon at lepton collider, i.e., $e^+e^-\to H��$, is known to bear a rather small cross section in Standard Model, and can serve as a sensitive probe for the potential new physics signals. Similar to the loop-induced Higgs decay channels $H\to ����, Z��$, the $e^+e^-\to H��$ process also starts at one-loop order provided that the tiny electron mass is neglected. In this work, we calculate the next-to-leading-order (NLO) QCD corrections to this associated $H+��$ production process, which mainly stem from the gluonic dressing to the top quark loop. The QCD corrections are found to be rather modest at lower center-of-mass energy range ($\sqrt{s}, v4, 11 pages, 6 figures, 2 tables; errors in Appendix are fixed; version accepted for publication at PLB

We calculate the next-to-next-to-leading-order (NNLO) perturbative corrections to $P$-wave quarkonia annihilation decay to two photons, in the framework of nonrelativistic QCD (NRQCD) factorization. The order-$��_s^2$ short-distance coefficients associated with each helicity amplitude are presented in a semi-analytic form, including the "light-by-light" contributions. With substantial NNLO corrections, we find disquieting discrepancy when confronting our state-of-the-art predictions with the latest \textsf{BESIII} measurements, especially fail to account for the measured $��_{c2}\to����$ width. Incorporating the effects of spin-dependent forces would even exacerbate the situation, since it lifts the degeneracy between the nonperturbative NRQCD matrix elements of $��_{c0}$ and $��_{c2}$ toward the wrong direction. We also present the order-$��_s^2$ predictions to $��_{b0,2}\to����$, which await the future experimental test., 6 pages, 2 figures, 1 table; v2, typos fixed, minor changes in text

Within the framework of nonrelativisitic QCD factorization formalism, we present the next-to-leading-order relativistic corrections to $��$ exclusive decay into $��_c$ plus $J/��$. The double charmonia can be produced through several immediate channels, i.e. $��\to g^*g^*g^*\to ��_c+J/��$, $��\to g^*g^*��^*\to ��_c+J/��$, and $��\to ��^*\to ��_c+J/��$. The amplitudes of these three channels are obtained accurate up to ${\cal O}(��_s^3 v^2)$, ${\cal O}(����_s^2 v^2)$, and ${\cal O}(����_s^2 v^2)$ respectively, where $v$ indicates the typical heavy quark velocity in bottomonium and/or charmonium rest frame. The decay rates are also presented. We find the next-to-leading-order relativistic corrections to the short distance coefficients as well as the decay rates are both significant and negative, especially for the corrections from bottomonium. More seriously, the decay rates are even brought into negative by including the relativistic corrections, which indicates the poor convergence for the velocity expansion in this kind of processes. Detailed analysis is given in the paper., 11 pages, 2 figures. v2: added an appendix, corrected some typos

Unlike the bewildering situation in the $����^*\to ��$ form factor, a widespread view is that perturbative QCD can decently account for the recent \textsc{BaBar} measurement of $����^*\to ��_c$ transition form factor. The next-to-next-to-leading order (NNLO) perturbative correction to the $����^*\to ��_{c,b}$ form factor, is investigated in the NRQCD factorization framework for the first time. As a byproduct, we obtain by far the most precise order-$��_s^2$ NRQCD matching coefficient for the $��_{c,b}\to ����$ process. After including the substantial negative order-$��_s^2$ correction, the good agreement between NRQCD prediction and the measured $����^*\to ��_c$ form factor is completely ruined over a wide range of momentum transfer squared. This eminent discrepancy casts some doubts on the applicability of NRQCD approach to hard exclusive reactions involving charmonium., 6 pages, 3 figures and 1 table; adding Eqs.(8) and (9) as well as some references, correcting errors in Table 1, updating Fig.3 to include the "light-by-light" contributions, devoting a paragraph to discuss why our strategy of interpreting the NNLO corrections is justified; Accepted by PRL

With the recent LHCb data on $��_c$ production and based on heavy quark spin symmetry, we obtain the long-distance matrix elements for both $��_c$ and $J/��$ productions, among which, the color-singlet one for $��_c$ is obtained directly by the fit of experiment for the first time. Using our long-distance matrix elements, we can provide good description of the $��_c$ and $J/��$ hadroproduction measurements. Our predictions on $J/��$ polarization are in good agreement with the LHCb data and pass through the two sets of CDF measurements in medium $p_t$ region. Considering all the possible uncertainties carefully, we obtained quite narrow bands of the $J/��$ polarization curves.

It has been a long-standing problem that the rare electromagnetic decay process $J/��\to 3��$ is plagued with both large and negative radiative and relativistic corrections. To date it remains futile to make a definite prediction to confront with the branching fraction of $J/��\to 3��$ recently measured by the \textsf{CLEO-c} and \textsf{BESIII} Collaborations. In this work, we investigate the joint perturbative and relativistic correction (i.e. the ${\mathcal O}(��_s v^2)$ correction, where $v$ denotes the characteristic velocity of the charm quark inside the $J/��$) for this decay process, which turns out to be very significant. After incorporating the contribution from this new ingredient, with the reasonable choice of the input parameters, we are able to account for the measured decay rates in a satisfactory degree., 7 pages, 1 figure, version accepted for publication in PRD RC

We investigate the $O(��_s v^2)$ correction to the process of pseudoscalar quarkonium decay to two photons in nonrelativistic QCD (NRQCD) factorization framework. The short-distance coefficient associated with the relative-order $v^2$ NRQCD matrix element is determined to next-to-leading order in $��_s$ through the perturbative matching procedure. Some technical subtleties encountered in calculating the ${O(��_s)$ QCD amplitude are thoroughly addressed., v2, 28 pages, 2 figures and 2 tables, matching the published version; typos corrected, references added, as well as a "Note added in the proof"

In the framework of the nonrelativistic quantum chromodynamics factorization formalism, we study the processes of $��(nS)$ and $��(nS)$ decay into a lepton pair or a charm pair associated with two jets up to the next-to-leading order in velocity expansion. We present the analytic expressions for the differential decay rate to the invariant mass of the lepton pair or charm pair. We find that the ratio of the next-to-leading order short-distance coefficient to the leading order one is in the range from -5.5 to -12.4. The relativistic corrections are so large that they modify the leading order prediction significantly. Utilizing the analytic expressions, we also investigate the relativistic corrections in different kinematic regions and their dependence on the masses of the initial-state quarkonium and the final-state fermion. In addition, we study the momentum distribution of $D^{*+}$ in the process $��(1S)\to c\bar{c}gg\to D^{*+}X$., 26 pages, 4 figures, 3 tables. arXiv admin note: text overlap with arXiv:0707.4056 by other authors