alphaFeedback

WorksPhysics and Astronomy

NNLO QCD corrections to direct J/ψJ/ψ production with open charm at B factories

AI contributions
  • Generative AI tools
View details

Abstract

We calculate the next-to-next-to-leading-order (NNLO) QCD corrections to direct colour-singlet J/ψJ/ψ production with open charm at B factories, using nested soft-collinear subtraction for massive partons. At s=10.6\sqrt{s}=10.6 GeV, mc=1.5m_c=1.5 GeV and μR=2mcμ_R=2m_c, we obtain an approximate NNLO cross section of 0.324−0.066+0.0410.324^{+0.041}_{-0.066} pb, with scale uncertainties from mc≤μR≤4mcm_c\leqμ_R\leq4m_c. The total cross section is 18.4%18.4\% above NLO and exhibits good perturbative convergence through NNLO, with reduced scale dependence. An NLO-based mass extrapolation to mc=1.4m_c=1.4 GeV, combined with estimated ψ(2S)ψ(2S) feed-down, raises the rate to about 0.590.59 pb. The NNLO correction brings the absolute momentum spectrum near its peak closer to Belle's prompt open-charm data, but shifts the angular shape further away from Belle's prompt charm-associated data. The difference in angular shape motivates more precise measurements at Belle II.

†These authors contributed equally to this work and share first authorship.

*Corresponding authors.

Introduction

Associated production of J/ψJ/\psi and an open charm pair in electron-positron annihilation provides a test of heavy-quarkonium production mechanisms. In nonrelativistic QCD (NRQCD), the production cross section is factorised into perturbative short-distance coefficients and non-perturbative long-distance matrix elements describing quarkonium formation [1]. Early colour-singlet calculations established the contributions from associated open charm and two-gluon production [2, 3], with subsequent studies examining the charm-associated rate and differential distributions [4, 5].

Measurements of prompt charmonium production by BaBar and Belle [6, 7], followed by Belle’s observation of a large double-charm component [8], exposed discrepancies with leading-order (LO) predictions. Zhang and Chao calculated the next-to-leading-order (NLO) QCD correction to the total rate [9]; Gong and Wang performed an independent calculation and studied differential distributions [10]. For commonly used inputs, these calculations found an enhancement by a factor of approximately 1.7–1.8 over LO [9, 10]. Belle reported a charm-associated prompt cross section of 0.74±0.08−0.08+0.09 pb0.74 \pm0.08 {}^{+0.09}_{-0.08}\,\mathrm{pb}, including open-charm and double-charmonium final states, together with momentum and angular distributions [11]. The NLO production-angle distribution, however, differed from the shape measured by Belle [10]. Further studies have examined relativistic corrections [12] and renormalisation-scale setting using the principle of maximum conformality [13].

The complementary J/ψ+ggJ/\psi+gg channel is important for interpreting Belle’s sample without an associated charm pair [11]. Its NLO QCD corrections were calculated independently in refs. [14, 15], and colour-singlet relativistic corrections were obtained in refs. [16, 17]. Further developments include endpoint resummation [18], NLO colour-octet contributions [19], initial-state radiation [20], and the recent mixed O(αsv2)O(\alpha_s v^2) correction to the colour-singlet rate and angular distribution [21].

NLO QCD predictions for these inclusive production channels have been available for more than a decade [9, 10, 14, 15]. More recently, several exclusive processes have been studied at next-to-next-to-leading order (NNLO). For exclusive e+e−→J/ψ+ηce^+e^- \to J/\psi+\eta_c, Feng et al. [22] first obtained NNLO predictions, followed by Huang et al. [23]; analytic two-loop amplitudes are now available as asymptotic expansions in mc2/sm_c^2/s [24]. NNLO corrections to exclusive double-J/ψJ/\psi production have been studied in refs. [25, 26].

Large higher-order corrections in exclusive quarkonium processes have also motivated alternative perturbative prescriptions. An amplitude-level prescription retains the complete modulus square of the truncated amplitude and has been proposed and studied for J/ψ→3γJ/\psi\to3\gamma and exclusive charmonium production [27, 28, 29]. Zhan et al. have applied this prescription at NLO to γ→J/ψ+γ\gamma\to J/\psi+\gamma in ultraperipheral collisions [30]. Separately, Maxia et al. have proposed a soft-improved matching scheme that includes finite soft contributions in the NNLO renormalisation of colour-singlet S-wave NRQCD matrix elements [31].

In the inclusive case, NNLO results are available for the hadronic decays of ηc\eta_c and ηb\eta_b, whose leading partonic channel is ηQ→gg\eta_Q \to gg [32]. For J/ψJ/\psi production with an open charm pair, the sizeable NLO correction makes it important to test the perturbative stability of both the total rate and differential distributions at NNLO. This process already has a threebody Born final state containing massive coloured particles. Besides the two-loop virtual terms, the NNLO correction involves two gluons or a massless quark pair in the double-real channels, and the one-loop correction to J/ψ+ccˉ+gJ/\psi+ c\bar{c} + g in the real-virtual contribution. Their infrared singularities must be cancelled while retaining a consistent definition of the measured observables across the different final-state multiplicities.

General subtraction methods organise this cancellation locally in phase space. At NLO, the FKS [33] and dipole [34] methods provide standard constructions, with extensions to massive partons [35] and NRQCD production [36]. NNLO approaches include antenna subtraction [37], sector-improved residue subtraction [38], CoLoRFulNNLO [39] and local analytic sector subtraction [40], as well as nested soft-collinear subtraction [41]. Fully differential massive-pair production from a colourless initial state has been treated within the latter scheme [42] and CoLoRFulNNLO [43].

In this work we present the first NNLO QCD calculation of direct colour-singlet J/ψJ/\psi production with an open charm pair in electron-positron annihilation through a virtual photon. The real-emission contributions are treated with the massive-parton formulation of nested soft-collinear subtraction [42]. The heavy-quark mass removes collinear singularities involving the open charm quarks, but overlapping soft and light-parton collinear limits remain in the double-real channels. Soft limits are subtracted before the remaining collinear limits, and the integrated counterterms are combined with the virtual contributions on the corresponding reduced phase spaces. For the double-soft terms, we use the analytic integrals for back-to-back massive emitters of ref. [44]; their extension to arbitrary relative angles has recently been studied in ref. [45].

We combine the real and virtual corrections into finite phase-space integrals, with analytic cancellation of infrared poles on the single-real and Born phase spaces. The subtraction mappings keep the measured quarkonium momentum fixed, allowing us to retain the differential information. We present approximate predictions at s=10.6 GeV\sqrt{s} = 10.6\ \mathrm{GeV} for the total cross section and the production-angle and momentum distributions, including their renormalisation-scale dependence.

§ 2 defines the process, observables and renormalisation conventions. § 3 explains how the NLO and NNLO real and virtual corrections are separated and recombined into finite contributions. The numerical implementation and the phenomenological results are discussed separately in § 4 and § 5, followed by the conclusions in § 6. Appendix § 8 derives the angular reconstruction, and appendix § 9 gives the integrated kernels and pole cancellations. The local counterterms and their integrals, phase-space mappings and auxiliary analytic functions are given in appendices § 10, § 11 and § 12, respectively.

Theoretical framework

Process definition and perturbative expansion

The process studied here is direct J/ψJ/\psi production with open charm in electron-positron annihilation through a virtual photon,

e−(p−)+e+(p+)⟶γ∗(Q)⟶J/ψ(P)+c(p1)+cˉ(p2)+Xlight,Q=p−+p+,Q2=s.(1)e^{-}(p_{-})+e^{+}(p_{+})\longrightarrow\gamma^{*}(Q)\longrightarrow J/\psi(P)+c(p_{1})+\bar{c}(p_{2})+X_{\mathrm{light}}, \qquad Q=p_{-}+p_{+},\qquad Q^{2}=s. \tag*{(1)}

The quarkonium is produced directly in the colour-singlet ccˉ[3S1[1]]c\bar{c}[{}^{3}S_{1}^{[1]}] state. At leading order in the heavy-quark relative velocity vv, the charm quark and antiquark forming the J/ψJ/\psi each carry momentum P/2P/2, with P2=Mψ2=4mc2P^{2}=M_{\psi}^{2}=4m_{c}^{2}. The momenta p1p_{1} and p2p_{2} denote the open charm quark and antiquark, respectively. The incoming beams are unpolarised and the electron mass is set to zero. The light-parton final states are Xlight=∅,g,gg,qqˉX_{\mathrm{light}}=\varnothing,g,gg,q\bar{q}, with q=u,d,sq=u,d,s.

NRQCD factorisation gives [1]

dσ=dσ^[e+e−⟶ccˉ[3S1[1]]+ccˉ+Xlight]⟨O1⟩J/ψ(μΛ).(2)\mathrm{d}\sigma=\mathrm{d}\hat{\sigma}\bigl[e^{+}e^{-}\longrightarrow c\bar{c}[{}^{3}S_{1}^{[1]}]+c\bar{c}+X_{\mathrm{light}}\bigr]\langle\mathcal{O}_{1}\rangle_{J/\psi}(\mu_{\Lambda}). \tag*{(2)}

Here dσ^\mathrm{d}\hat{\sigma} is the perturbative short-distance coefficient, and ⟨O1⟩J/ψ\langle\mathcal{O}_{1}\rangle_{J/\psi} is the non-perturbative colour-singlet matrix element in the single-polarisation convention, related to the wave function by ∣R(0)∣2=(2π/3)⟨O1⟩J/ψ\lvert R(0)\rvert^{2}=(2\pi/3)\langle\mathcal{O}_{1}\rangle_{J/\psi}.

We denote the contribution of order α2αs2+n\alpha^{2}\alpha_{s}^{2+n} to the physical cross section σ\sigma, including this matrix element, by σ(n)\sigma^{(n)}, with n=0,1,2n=0,1,2. The fixed-order cross sections are

σLO=σ(0),\sigma_{\mathrm{LO}}=\sigma^{(0)},
σNLO=σ(0)+σ(1),\sigma_{\mathrm{NLO}}=\sigma^{(0)}+\sigma^{(1)},
σNNLO=σ(0)+σ(1)+σ(2).\sigma_{\mathrm{NNLO}}=\sigma^{(0)}+\sigma^{(1)}+\sigma^{(2)}.

Angular reconstruction from the virtual-photon tensor

We obtain the total cross sections σ\sigma and σL\sigma_{L} with the respective photon-polarisation tensors

Πμν=−gμν+QμQνs,ΠLμν=eLμeLν,eLμ=Pμ−(P⋅Q)Qμ/s(P⋅Q)2/s−Mψ2.(3)\Pi^{\mu\nu}=-g^{\mu\nu}+\frac{Q^{\mu}Q^{\nu}}{s}, \qquad \Pi_{L}^{\mu\nu}=e_{L}^{\mu}e_{L}^{\nu}, \qquad e_{L}^{\mu}=\frac{P^{\mu}-(P\cdot Q)Q^{\mu}/s}{\sqrt{(P\cdot Q)^{2}/s-M_{\psi}^{2}}.} \tag*{(3)}

Let θ\theta be the angle between the J/ψJ/\psi and incoming-electron momenta in the centre-of-mass frame. The angular distribution is then

dσdcos⁡θ=38[(σ+σL)+(σ−3σL)cos⁡2θ],αθ=σ−3σLσ+σL.(4)\frac{\mathrm{d}\sigma}{\mathrm{d}\cos\theta} =\frac{3}{8}\left[(\sigma+\sigma_{L})+(\sigma-3\sigma_{L})\cos^{2}\theta\right], \qquad \alpha_{\theta}=\frac{\sigma-3\sigma_{L}}{\sigma+\sigma_{L}}. \tag*{(4)}

The derivation is given in appendix A.

Renormalisation conventions

We use dimensional regularisation with d=4−2ϵd = 4 - 2\epsilon. The charm-quark mass and field are renormalised in the on-shell scheme, and the strong coupling in the MS‾\overline{\mathrm{MS}} scheme with four active flavours. The three light flavours u,d,su,d,s are treated as massless.

At relative order αs2\alpha_s^2, the infrared pole associated with the colour-singlet ccˉ[3S1[1]]c\bar{c}[{}^3S_1^{[1]}] pair is absorbed into the MS‾\overline{\mathrm{MS}}-renormalised NRQCD matrix element [46, 47, 48]. Subtracting the term in (B.6) from the ultraviolet-renormalised QCD amplitude implements this matching and introduces the factorisation scale μΛ\mu_\Lambda. We fix μΛ=mc\mu_\Lambda= m_c, independently of the QCD renormalisation scale μR\mu_R. The short-distance amplitudes used below include this NRQCD matching; the remaining QCD infrared poles cancel in the cross section after adding the real-emission contributions.

Organisation of the NLO and NNLO calculation

Let Mccˉ+X(ℓ)\mathcal{M}^{(\ell)}_{c\bar{c}+X} be the ℓ\ell-loop short-distance amplitude coefficient for γ∗→ccˉ[3S1[1]]+ccˉ+X\gamma^* \to c\bar{c}[{}^3S_1^{[1]}] + c\bar{c} + X, and dΦccˉ+Xd\Phi_{c\bar{c}+X} the phase space for the corresponding J/ψ+ccˉ+XJ/\psi+ c\bar{c} + X final state. The common J/ψJ/\psi is understood in these subscripts. With as=αs/(4π)a_s = \alpha_s/(4\pi), the renormalised amplitude is M=M(0)+asM(1)+as2M(2)+⋯\mathcal{M} = \mathcal{M}^{(0)} + a_s\mathcal{M}^{(1)} + a_s^2\mathcal{M}^{(2)} + \cdots. We define the perturbative contributions to the squared matrix element by

References

  1. [1]G. T. Bodwin, E. Braaten and G. P. Lepage, Rigorous QCD analysis of inclusive annihilation and production of heavy quarkonium, Phys. Rev. D 51 (1995) 1125; erratum D 55 (1997) 5853, hep-ph/9407339.pubmed.ncbi.nlm.nih.gov/10018572
  2. [2]P. L. Cho and A. K. Leibovich, Color-singlet ψQ production at e+e− colliders, Phys. Rev. D 54 (1996) 6690–6695, hep-ph/9606229.
  3. [3]F. Yuan, C.-F. Qiao and K.-T. Chao, Prompt J/ψ production at e+e− colliders, Phys. Rev. D 56 (1997) 321–328, hep-ph/9703438.
  4. [4]K.-Y. Liu, Z.-G. He and K.-T. Chao, Inclusive charmonium production via double c̄c in e+e− annihilation, Phys. Rev. D 69 (2004) 094027, hep-ph/0301218.
  5. [5]K. Hagiwara, E. Kou, Z.-H. Lin, C.-F. Qiao and G.-H. Zhu, Inclusive J/ψ production at e+e− colliders, Phys. Rev. D 70 (2004) 034013, hep-ph/0401246.
  6. [6]BaBar collaboration, B. Aubert et al., Measurement of J/ψ production in continuum e+e− annihilations near √s = 10.6 GeV, Phys. Rev. Lett. 87 (2001) 162002, hep-ex/0106044.
  7. [7]Belle collaboration, K. Abe et al., Production of prompt charmonia in e+e− annihilation at √s ≈ 10.6 GeV, Phys. Rev. Lett. 88 (2002) 052001, hep-ex/0110012.
  8. [8]Belle collaboration, K. Abe et al., Observation of double c̄c production in e+e− annihilation at √s ≈ 10.6 GeV, Phys. Rev. Lett. 89 (2002) 142001, hep-ex/0205104.
  9. [9]Y.-J. Zhang and K.-T. Chao, Double charm production e+e− → J/ψ + c̄c at B factories with next-to-leading order QCD corrections, Phys. Rev. Lett. 98 (2007) 092003, hep-ph/0611086.
  10. [10]B. Gong and J.-X. Wang, Next-to-leading-order QCD corrections to e+e− → J/ψc̄c at the B factories, Phys. Rev. D 80 (2009) 054015, 0904.1103.
  11. [11]Belle collaboration, P. Pakhlov et al., Measurement of the e+e− → J/ψc̄c cross section at √s ≈ 10.6 GeV, Phys. Rev. D 79 (2009) 071101, 0901.2775.
  12. [12]Z.-G. He, Y. Fan and K.-T. Chao, Relativistic corrections to J/ψ exclusive and inclusive double charm production at B factories, Phys. Rev. D 75 (2007) 074011, hep-ph/0702239.
  13. [13]X.-D. Huang, X.-G. Wu, X.-C. Zheng, B. Gong and J.-X. Wang, QCD corrections of e+e− → J/ψ + c̄c using the principle of maximum conformality, Phys. Rev. D 110 (2024) 114010, 2407.14150.
  14. [14]Y.-Q. Ma, Y.-J. Zhang and K.-T. Chao, QCD corrections to e+e− → J/ψ + gg at B factories, Phys. Rev. Lett. 102 (2009) 162002, 0812.5106.
  15. [15]B. Gong and J.-X. Wang, Next-to-leading-order QCD corrections to e+e− → J/ψgg at the B factories, Phys. Rev. Lett. 102 (2009) 162003, 0901.0117.
  16. [16]Z.-G. He, Y. Fan and K.-T. Chao, Relativistic correction to e+e− → J/ψ + gg at B factories and constraint on color-octet matrix elements, Phys. Rev. D 87 (2013) 054036, 0910.3636.
  17. [17]Y. Jia, Color-singlet relativistic correction to inclusive J/ψ production associated with light hadrons at B factories, Phys. Rev. D 82 (2010) 034017, 0912.5498.
  18. [18]A. K. Leibovich and X. Liu, Color-singlet contribution to e+e− → J/ψ + X at the endpoint, Phys. Rev. D 76 (2007) 034005, 0705.3230.
  19. [19]Y.-J. Zhang, Y.-Q. Ma, K. Wang and K.-T. Chao, QCD radiative correction to color-octet J/ψ inclusive production at B factories, Phys. Rev. D 81 (2010) 034015, 0911.2166.
  20. [20]H.-S. Shao, Initial state radiation effects in inclusive J/ψ production at B factories, JHEP 04 (2014) 182, 1402.5840.
  21. [21]C. Li, H.-Y. Liu, X.-D. Huang and W.-L. Sang, Semi-analytical results for e+e− → J/ψ + X_non c̄c up to O(α_sv2) at B factories, JHEP 09 (2026) 128, 2606.17827.
  22. [22]F. Feng, Y. Jia, Z. Mo, W.-L. Sang and J.-Y. Zhang, Next-to-next-to-leading-order QCD corrections to e+e− → J/ψ + ηc at B factories, Phys. Lett. B 850 (2024) 138506, 1901.08447.
  23. [23]X.-D. Huang, B. Gong and J.-X. Wang, Next-to-next-to-leading-order QCD corrections to J/ψ plus ηc production at the B factories, JHEP 02 (2023) 049, 2212.03631.
  24. [24]X. Chen, X. Guan, C.-Q. He, Y.-Q. Ma, J. Wang and D.-J. Zhang, Analytical two-loop amplitudes of e+e− → J/ψ + ηc at B factories, Phys. Rev. D 113 (2026) 074023, 2508.20777.
  25. [25]W.-L. Sang, F. Feng, Y. Jia, Z. Mo, J. Pan and J.-Y. Zhang, Optimized O(α_s^2) correction to exclusive double J/ψ production at B factories, Phys. Rev. Lett. 131 (2023) 161904, 2306.11538.
  26. [26]X.-D. Huang, B. Gong, R.-C. Niu, H.-M. Yu and J.-X. Wang, Next-to-next-to-leading-order QCD corrections to double J/ψ production at the B factories, JHEP 02 (2024) 055, 2311.04751.
  27. [27]C. Zeng, B. Gong, J.-X. Wang, R. Niu, X.-D. Huang and C. Li, Complete next-to-next-to-leading-order QCD correction to J/ψ → 3γ decay, 2603.26199.
  28. [28]X.-D. Huang and J.-X. Wang, Reanalysis of exclusive charmonium production at the B factories, 2609.00637.
  29. [29]B. Gong, X.-D. Huang, J.-X. Wang and C. Zeng, Violation and restoration of probability positivity in perturbative quantum field theory: from quantum mechanics to quantum chromodynamics, 2610.04557.
  30. [30]X.-J. Zhan, T.-F. Feng and X.-G. Wu, Amplitude-level NLO predictions for exclusive J/ψ + γ production in ultraperipheral collisions, 2610.04293.
  31. [31]L. Maxia, H.-S. Shao, L. Simon and G. Wang, Soft contributions stabilize NNLO QCD corrections to quarkonium production and decay, 2606.28494.
  32. [32]F. Feng, Y. Jia and W.-L. Sang, Next-to-next-to-leading-order QCD corrections to the hadronic width of pseudoscalar quarkonium, Phys. Rev. Lett. 119 (2017) 252001, 1707.05758.
  33. [33]S. Frixione, Z. Kunszt and A. Signer, Three-jet cross sections to next-to-leading order, Nucl. Phys. B 467 (1996) 399–442, hep-ph/9512328.
  34. [34]S. Catani and M. H. Seymour, A general algorithm for calculating jet cross sections in NLO QCD, Nucl. Phys. B 485 (1997) 291–419; erratum B 510 (1998) 503–504, hep-ph/9605323.
  35. [35]S. Catani, S. Dittmaier, M. H. Seymour and Z. Trócsányi, The dipole formalism for next-to-leading order QCD calculations with massive partons, Nucl. Phys. B 627 (2002) 189–265, hep-ph/0201036.
  36. [36]M. Butenschön and B. A. Kniehl, Dipole subtraction at next-to-leading order in nonrelativistic-QCD factorization, Nucl. Phys. B 950 (2020) 114843, 1909.03698.
  37. [37]A. Gehrmann-De Ridder, T. Gehrmann and E. W. N. Glover, Antenna subtraction at NNLO, JHEP 09 (2005) 056, hep-ph/0505111.
  38. [38]M. Czakon, A novel subtraction scheme for double-real radiation at NNLO, Phys. Lett. B 693 (2010) 259–268, 1005.0274.
  39. [39]G. Somogyi, Z. Trócsányi and V. Del Duca, A subtraction scheme for computing QCD jet cross sections at NNLO: regularization of doubly-real emissions, JHEP 01 (2007) 070, hep-ph/0609042.
  40. [40]L. Magnea, E. Maina, G. Pelliccioli, C. Signorile-Signorile, P. Torrielli and S. Uccirati, Local analytic sector subtraction at NNLO, JHEP 12 (2018) 107; erratum JHEP 06 (2019) 013, 1806.09570.
  41. [41]F. Caola, K. Melnikov and R. Röntsch, Nested soft-collinear subtractions in NNLO QCD computations, Eur. Phys. J. C 77 (2017) 248, 1702.01352.
  42. [42]A. Behring and W. Bizón, Higgs decay into massive b-quarks at NNLO QCD in the nested soft-collinear subtraction scheme, JHEP 01 (2020) 189, 1911.11524.
  43. [43]G. Somogyi and F. Tramontano, Fully exclusive heavy quark-antiquark pair production from a colourless initial state at NNLO in QCD, JHEP 11 (2020) 142, 2007.15015.
  44. [44]W. Bizón and M. Delto, Analytic double-soft integrated subtraction terms for two massive emitters in a back-to-back kinematics, JHEP 07 (2020) 011, 2004.01663.
  45. [45]M.-M. Long, K. Melnikov and A. Pikelner, Integral of the double-emission eikonal function for two massive emitters at an arbitrary angle, JHEP 03 (2026) 117, 2510.04778.
  46. [46]A. Czarnecki and K. Melnikov, Two-loop QCD corrections to the heavy quark pair production cross section in e⁺e⁻ annihilation near threshold, Phys. Rev. Lett. 80 (1998) 2531–2534, hep-ph/9712222.DOI
  47. [47]M. Beneke, A. Signer and V. A. Smirnov, Two-loop corrections to the leptonic decays of quarkonium, Phys. Rev. Lett. 80 (1998) 2535–2538, hep-ph/9712302.DOI
  48. [48]B. A. Kniehl, A. Onishchenko, J. H. Piclum and M. Steinhauser, Two-loop matching coefficients for heavy quark currents, Phys. Lett. B 638 (2006) 209–213, hep-ph/0604072.DOI
  49. [49]Particle Data Group, S. Navas et al., Review of Particle Physics, Phys. Rev. D 110 (2024) 030001, and 2025 update, Quantum Chromodynamics.
  50. [50]K. G. Chetyrkin, J. H. Kühn and M. Steinhauser, RunDec: a Mathematica package for running and decoupling of the strong coupling and quark masses, Comput. Phys. Commun. 133 (2000) 43–65, hep-ph/0004189.DOI
  51. [51]G. T. Bodwin, H. S. Chung, D. Kang, J. Lee and C. Yu, Improved determination of color-singlet nonrelativistic QCD matrix elements for S-wave charmonium, Phys. Rev. D 77 (2008) 094017, 0710.0994.DOI
  52. [52]F. Lange, J. Usovitsch and Z. Wu, Kira 3: integral reduction with efficient seeding and optimized equation selection, Comput. Phys. Commun. 322 (2026) 109999, 2505.20197.arxiv.org/abs/2505.20197
  53. [53]X. Liu and Y.-Q. Ma, AMFlow: a Mathematica package for Feynman integrals computation via auxiliary mass flow, Comput. Phys. Commun. 283 (2023) 108565, 2201.11669.DOI
  54. [54]J. Kuipers, T. Ueda, J. A. M. Vermaseren and J. Vollinga, FORM version 4.0, Comput. Phys. Commun. 184 (2013) 1453–1467, 1203.6543.
  55. [55]J. Boehm, M. Wittmann, Z. Wu, Y. Xu and Y. Zhang, IBP reduction coefficients made simple, JHEP 12 (2020) 054, 2008.13194.
  56. [56]W. Decker, G.-M. Greuel, G. Pfister and H. Schö​nemann, Singular: a computer algebra system for polynomial computations, www.singular.uni-kl.de.
  57. [57]J.-X. Wang, Progress in FDC project, Nucl. Instrum. Meth. A 534 (2004) 241–245, hep-ph/0407058.DOI
  58. [58]R.-J. Huang, D.-S. Jian, Y.-Q. Ma, D.-M. Mu and W.-H. Wu, Efficient computation of one-loop Feynman integrals and fixed-branch integrals to high orders in ε, Phys. Rev. D 111 (2025) 094028, 2412.21054.arxiv.org/abs/2412.21054
  59. [59]F. Johansson, Arb: efficient arbitrary-precision midpoint-radius interval arithmetic, IEEE Trans. Comput. 66 (2017) 1281–1292, 1611.02831.arxiv.org/abs/1611.02831
  60. [60]G. P. Lepage, Adaptive multidimensional integration: VEGAS enhanced, J. Comput. Phys. 439 (2021) 110386, 2009.05112.DOI
  61. [61]R. Piessens, E. de Doncker-Kapenga, C. W. Überhuber and D. K. Kahaner, QUADPACK: A Subroutine Package for Automatic Integration, Springer (1983), doi:10.1007/978-3-642-61786-7; DQK15 routine.
  62. [62]T. Becher and M. Neubert, Infrared singularities of QCD amplitudes with massive partons, Phys. Rev. D 79 (2009) 125004; erratum D 80 (2009) 109901, 0904.1021.
  63. [63]A. Ferroglia, M. Neubert, B. D. Pecjak and L. L. Yang, Two-loop divergences of massive scattering amplitudes in non-abelian gauge theories, JHEP 11 (2009) 062, 0908.3676.DOI
  64. [64]A. von Manteuffel, R. M. Schabinger and H. X. Zhu, The two-loop soft function for heavy quark pair production at future linear colliders, Phys. Rev. D 92 (2015) 045034, 1408.5134.DOI
  65. [65]I. Bierenbaum, M. Czakon and A. Mitov, The singular behavior of one-loop massive QCD amplitudes with one external soft gluon, Nucl. Phys. B 856 (2012) 228–246, 1107.4384.
  66. [66]M. Czakon and A. Mitov, A simplified expression for the one-loop soft-gluon current with massive fermions, 1804.02069.arxiv.org/abs/1804.02069
  67. [67]S. Alioli, P. Nason, C. Oleari and E. Re, A general framework for implementing NLO calculations in shower Monte Carlo programs: the POWHEG BOX, JHEP 06 (2010) 043, 1002.2581.DOI

Paper details

vlnc:W_a6fn6trx
v1 · latestHistory (1)

published · Fresh — added to the Hub within the last 30 days

added ·

Usage

views · 0

downloads · 0

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Physics and Astronomy

Generative AI tools · View statement →

Looking for related papers…
Contents