1 Introduction
Precision in calculations of cross sections for hard-scattering processes requires the inclusion of higher-order QCD corrections. A particular class of these corrections arises from the emission of soft gluons. These soft radiative corrections can be formally resummed [1 , 2 , 3 , 4 , 5 ] to all orders in the strong coupling after taking Mellin moments or Laplace transforms of the cross section and using its factorization properties, and the resummed cross section can be used as a generator of fixed-order corrections. Soft-gluon resummation has been successfully developed and applied to processes with massive final states, such as top-quark pair production [3 , 4 , 5 , 6 , 7 , 8 ] , because the cross section receives large contributions from soft-gluon emission near partonic threshold due to the large mass. In all these processes, the soft-gluon corrections are numerically dominant and account for the overwhelming majority of the complete corrections at higher orders.
For processes with complex color flows, such as top production, one does not yet have all the ingredients necessary to calculate the complete set of soft-gluon corrections at next-to-next-to-next-to-leading order (N3 LO). For processes with trivial color flow, and in particular with colorless (e.g. electroweak) final states, the N3 LO soft-gluon corrections are known in fully inclusive kinematics [9 , 10 , 11 ] . However, complete N3 LO soft-gluon contributions have not yet been calculated for a partonic process in single-particle-inclusive (1PI) kinematics. In this paper we provide explicit analytical results for the complete soft-gluon corrections through N3 LO for 1PI partonic processes with no color in the final state. As a concrete application, we employ these results to compute higher-order QCD corrections for charged-Higgs pair production in the Two-Higgs-Doublet Model (2HDM) using 1PI kinematics.
In this work we present the most precise predictions for quark-induced production of H + H − superscript 𝐻 superscript 𝐻 H^{+}H^{-} italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT at the LHC, that occurs through diagrams of Drell-Yan type, in the 2HDM by computing the total cross section and the charged-Higgs rapidity distribution for this process at approximate next-to-next-to-leading order (aNNLO) and approximate N3 LO (aN3 LO) in QCD. Results at aNNLO are derived by adding the second-order soft-plus-virtual corrections to the full next-to-leading order (NLO) QCD calculation, while aN3 LO results are obtained by further adding the third-order soft-gluon corrections.
The 2HDM represents one of the simplest extensions of the Standard Model Higgs sector, introducing an additional S U ( 2 ) L 𝑆 𝑈 subscript 2 𝐿 SU(2)_{L} italic_S italic_U ( 2 ) start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT Higgs doublet. This model gives rise to a very rich phenomenology due to its extended Higgs sector (for a review see [12 ] and references therein). The model predicts the existence of three neutral scalar bosons (h ℎ h italic_h , H 𝐻 H italic_H , and A 𝐴 A italic_A ) and a pair of charged Higgs bosons (H + superscript 𝐻 H^{+} italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT and H − superscript 𝐻 H^{-} italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT ).
The dominant production channel for the charged Higgs at the Large Hadron Collider (LHC) is in association with a top quark p p → t H − → 𝑝 𝑝 𝑡 superscript 𝐻 p\,p\to t\,H^{-} italic_p italic_p → italic_t italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT , for which soft-gluon corrections were calculated in Refs. [14 , 13 , 15 ] , followed by the associated production with a W 𝑊 W italic_W boson p p → H + W − → 𝑝 𝑝 superscript 𝐻 superscript 𝑊 p\,p\to H^{+}W^{-} italic_p italic_p → italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_W start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT (see [16 ] for soft-gluon corrections), and pair production p p → H + H − → 𝑝 𝑝 superscript 𝐻 superscript 𝐻 p\,p\to H^{+}H^{-} italic_p italic_p → italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT . In this work we focus on the last process, which can occur at leading order (LO) through quark-antiquark annihilation (q q ¯ → H + H − → 𝑞 ¯ 𝑞 superscript 𝐻 superscript 𝐻 q\bar{q}\to H^{+}H^{-} italic_q over¯ start_ARG italic_q end_ARG → italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT ) [17 , 18 , 19 ] and gluon fusion (g g → H + H − → 𝑔 𝑔 superscript 𝐻 superscript 𝐻 gg\to H^{+}H^{-} italic_g italic_g → italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT ) [20 , 21 , 22 , 23 , 24 ] . In the case of light quarks in the initial state, the pair production of charged Higgs bosons is of Drell-Yan type, with the exchange of a photon and a Z 𝑍 Z italic_Z boson in the s 𝑠 s italic_s -channel. NLO QCD corrections for this channel have been computed in the 2HDM and the Minimal Supersymmetric Standard Model (MSSM) in [17 , 19 ] . In the case of bottom quarks in the initial states, there are additional diagrams that account for the exchange of a neutral Higgs boson in the s 𝑠 s italic_s -channel and a top quark in the t 𝑡 t italic_t -channel. The NLO QCD corrections to b b ¯ → H + H − → 𝑏 ¯ 𝑏 superscript 𝐻 superscript 𝐻 b\bar{b}\to H^{+}H^{-} italic_b over¯ start_ARG italic_b end_ARG → italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT have been computed in [19 , 18 ] . The charged-Higgs pair production through gluon fusion is a one-loop process already at LO and it has been computed in the 2HDM and the MSSM in [20 , 21 , 22 , 23 , 24 ] .
In models like the type II 2HDM, the quark-antiquark annihilation process q q ¯ → H + H − → 𝑞 ¯ 𝑞 superscript 𝐻 superscript 𝐻 q\bar{q}\to H^{+}H^{-} italic_q over¯ start_ARG italic_q end_ARG → italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT is the dominant channel for low values of tan β 𝛽 \tan\beta roman_tan italic_β , which is the ratio of the vacuum expectation values for the two doublets. On the other hand, for large values of tan β 𝛽 \tan\beta roman_tan italic_β , the g g 𝑔 𝑔 gg italic_g italic_g and b b ¯ 𝑏 ¯ 𝑏 b\bar{b} italic_b over¯ start_ARG italic_b end_ARG channels have cross sections which are enhanced by a tan 4 β superscript 4 𝛽 \tan^{4}\beta roman_tan start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_β factor, and this can compensate for the loop suppression or the small bottom-quark parton densities and can make these channels dominant, as shown in [19 ] . However, recent studies showed that, for type II 2DHM with softly broken Z 2 subscript 𝑍 2 Z_{2} italic_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT symmetry, the size of cos ( β − α ) 𝛽 𝛼 \cos(\beta-\alpha) roman_cos ( italic_β - italic_α ) has to be lower than 𝒪 ( 10 − 2 ) 𝒪 superscript 10 2 {\cal O}(10^{-2}) caligraphic_O ( 10 start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT ) and that low values of tan β 𝛽 \tan\beta roman_tan italic_β are favored [25 , 26 , 27 , 28 ] . Therefore, the calculations presented in this work are particularly relevant for those scenarios.
The paper is organized as follows. In Section 2, we describe the soft-gluon resummation formalism and provide explicit analytical results for the soft-gluon corrections through N3 LO in 1PI kinematics for processes with colorless final states. In Section 3 we introduce the Two-Higgs-Doublet model and apply the resummation formalism of Section 2 to charged-Higgs pair production. In Section 4, we present numerical results for the total cross section through aN3 LO at LHC energies. In Section 5, we present results for the charged-Higgs rapidity distribution for two representative masses of the charged Higgs boson. We conclude in Section 6.
2 Soft-gluon corrections at N3 LO in 1PI kinematics
In this section, we discuss the resummation formalism that we use for the calculation of soft-gluon corrections in 1PI kinematics for partonic processes with colorless final states, e.g. with electroweak final-state particles, such as H + H − superscript 𝐻 superscript 𝐻 H^{+}H^{-} italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT production. The origin of the soft-gluon corrections is from the emission of soft (i.e. low-energy) gluons from the initial state partons, which result in partial cancellations of infrared divergences between real-emission and virtual diagrams close to partonic threshold. We derive the resummed cross section via Laplace transforms, factorization, and renormalization-group evolution. We then provide explicit analytical results for the soft-gluon corrections through N3 LO.
We consider first leading-order partonic processes of the form
q ( p a ) + q ¯ ( p b ) → A ( p 1 ) + B ( p 2 ) , → 𝑞 subscript 𝑝 𝑎 ¯ 𝑞 subscript 𝑝 𝑏 𝐴 subscript 𝑝 1 𝐵 subscript 𝑝 2 q(p_{a})\,+{\bar{q}}(p_{b})\to A(p_{1})\,+\,B(p_{2})\,, italic_q ( italic_p start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) + over¯ start_ARG italic_q end_ARG ( italic_p start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) → italic_A ( italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) + italic_B ( italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) ,
(2.1)
where q 𝑞 q italic_q and q ¯ ¯ 𝑞 {\bar{q}} over¯ start_ARG italic_q end_ARG are quarks and antiquarks in the protons, and A 𝐴 A italic_A and B 𝐵 B italic_B are a pair of colorless particles, each of mass m 𝑚 m italic_m
(e.g. H + H − superscript 𝐻 superscript 𝐻 H^{+}H^{-} italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT ), with particle A 𝐴 A italic_A the observed particle in 1PI kinematics. We will also discuss later the changes needed in the analytical expressions if we instead have gluons in the initial state, i.e. processes g g → A B → 𝑔 𝑔 𝐴 𝐵 gg\to AB italic_g italic_g → italic_A italic_B , or if the two final-state particles have different masses, m A subscript 𝑚 𝐴 m_{A} italic_m start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT and m B subscript 𝑚 𝐵 m_{B} italic_m start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT . We define s = ( p a + p b ) 2 𝑠 superscript subscript 𝑝 𝑎 subscript 𝑝 𝑏 2 s=(p_{a}+p_{b})^{2} italic_s = ( italic_p start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT + italic_p start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , t 1 = ( p a − p 1 ) 2 − m 2 subscript 𝑡 1 superscript subscript 𝑝 𝑎 subscript 𝑝 1 2 superscript 𝑚 2 t_{1}=(p_{a}-p_{1})^{2}-m^{2} italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = ( italic_p start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT - italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , u 1 = ( p b − p 1 ) 2 − m 2 subscript 𝑢 1 superscript subscript 𝑝 𝑏 subscript 𝑝 1 2 superscript 𝑚 2 u_{1}=(p_{b}-p_{1})^{2}-m^{2} italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = ( italic_p start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT - italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , as well as a partonic threshold variable s 4 = s + t 1 + u 1 = ( p 2 + p g ) 2 − m 2 subscript 𝑠 4 𝑠 subscript 𝑡 1 subscript 𝑢 1 superscript subscript 𝑝 2 subscript 𝑝 𝑔 2 superscript 𝑚 2 s_{4}=s+t_{1}+u_{1}=(p_{2}+p_{g})^{2}-m^{2} italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT = italic_s + italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = ( italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + italic_p start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , with p g subscript 𝑝 𝑔 p_{g} italic_p start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT the momentum of an additional gluon in the final state. Near partonic threshold p g → 0 → subscript 𝑝 𝑔 0 p_{g}\to 0 italic_p start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT → 0 and, thus, s 4 → 0 → subscript 𝑠 4 0 s_{4}\to 0 italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT → 0 .
The soft-gluon corrections appear in the perturbative series as logarithms of s 4 subscript 𝑠 4 s_{4} italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT , i.e. [ ( ln k ( s 4 / m 2 ) ) / s 4 ] + subscript delimited-[] superscript 𝑘 subscript 𝑠 4 superscript 𝑚 2 subscript 𝑠 4 [(\ln^{k}(s_{4}/m^{2}))/s_{4}]_{+} [ ( roman_ln start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT / italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ) / italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ] start_POSTSUBSCRIPT + end_POSTSUBSCRIPT , with 0 ≤ k ≤ 2 n − 1 0 𝑘 2 𝑛 1 0\leq k\leq 2n-1 0 ≤ italic_k ≤ 2 italic_n - 1 at n 𝑛 n italic_n th order in the strong coupling, α s subscript 𝛼 𝑠 \alpha_{s} italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT .
The resummation of soft-gluon corrections is a consequence of the factorization properties of the (in general, differential) cross section under Laplace transforms in 1PI kinematics. We first write the differential hadronic cross section, d σ p p → A B 𝑑 subscript 𝜎 → 𝑝 𝑝 𝐴 𝐵 d\sigma_{pp\to AB} italic_d italic_σ start_POSTSUBSCRIPT italic_p italic_p → italic_A italic_B end_POSTSUBSCRIPT , as a convolution of the differential partonic cross section, d σ ^ q q ¯ → A B 𝑑 subscript ^ 𝜎 → 𝑞 ¯ 𝑞 𝐴 𝐵 d{\hat{\sigma}}_{q{\bar{q}}\to AB} italic_d over^ start_ARG italic_σ end_ARG start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT , with the parton distribution functions (pdf), ϕ q / p subscript italic-ϕ 𝑞 𝑝 \phi_{q/p} italic_ϕ start_POSTSUBSCRIPT italic_q / italic_p end_POSTSUBSCRIPT and ϕ q ¯ / p subscript italic-ϕ ¯ 𝑞 𝑝 \phi_{{\bar{q}}/p} italic_ϕ start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG / italic_p end_POSTSUBSCRIPT , as
d σ p p → A B = ∑ q , q ¯ ∫ 𝑑 x a 𝑑 x b ϕ q / p ( x a , μ F ) ϕ q ¯ / p ( x b , μ F ) 𝑑 σ ^ q q ¯ → A B ( s 4 , μ F ) , 𝑑 subscript 𝜎 → 𝑝 𝑝 𝐴 𝐵 subscript 𝑞 ¯ 𝑞
differential-d subscript 𝑥 𝑎 differential-d subscript 𝑥 𝑏 subscript italic-ϕ 𝑞 𝑝 subscript 𝑥 𝑎 subscript 𝜇 𝐹 subscript italic-ϕ ¯ 𝑞 𝑝 subscript 𝑥 𝑏 subscript 𝜇 𝐹 differential-d subscript ^ 𝜎 → 𝑞 ¯ 𝑞 𝐴 𝐵 subscript 𝑠 4 subscript 𝜇 𝐹 d\sigma_{pp\to AB}=\sum_{q,{\bar{q}}}\;\int dx_{a}\,dx_{b}\,\phi_{q/p}(x_{a},%
\mu_{F})\,\phi_{{\bar{q}}/p}(x_{b},\mu_{F})\,d{\hat{\sigma}}_{q{\bar{q}}\to AB%
}(s_{4},\mu_{F})\,, italic_d italic_σ start_POSTSUBSCRIPT italic_p italic_p → italic_A italic_B end_POSTSUBSCRIPT = ∑ start_POSTSUBSCRIPT italic_q , over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT ∫ italic_d italic_x start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_d italic_x start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_ϕ start_POSTSUBSCRIPT italic_q / italic_p end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ) italic_ϕ start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG / italic_p end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT , italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ) italic_d over^ start_ARG italic_σ end_ARG start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT , italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ) ,
(2.2)
where μ F subscript 𝜇 𝐹 \mu_{F} italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT is the factorization scale, and x a subscript 𝑥 𝑎 x_{a} italic_x start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , x b subscript 𝑥 𝑏 x_{b} italic_x start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT are momentum fractions of partons q 𝑞 q italic_q , q ¯ ¯ 𝑞 {\bar{q}} over¯ start_ARG italic_q end_ARG , respectively, in the colliding protons.
We also define the hadron-level variables S = ( P a + P b ) 2 𝑆 superscript subscript 𝑃 𝑎 subscript 𝑃 𝑏 2 S=(P_{a}+P_{b})^{2} italic_S = ( italic_P start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT + italic_P start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , T 1 = ( P a − p 1 ) 2 − m 2 subscript 𝑇 1 superscript subscript 𝑃 𝑎 subscript 𝑝 1 2 superscript 𝑚 2 T_{1}=(P_{a}-p_{1})^{2}-m^{2} italic_T start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = ( italic_P start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT - italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , U 1 = ( P b − p 1 ) 2 − m 2 subscript 𝑈 1 superscript subscript 𝑃 𝑏 subscript 𝑝 1 2 superscript 𝑚 2 U_{1}=(P_{b}-p_{1})^{2}-m^{2} italic_U start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = ( italic_P start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT - italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , where P a subscript 𝑃 𝑎 P_{a} italic_P start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT and P b subscript 𝑃 𝑏 P_{b} italic_P start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT denote the momenta of the colliding protons, and S 4 = S + T 1 + U 1 subscript 𝑆 4 𝑆 subscript 𝑇 1 subscript 𝑈 1 S_{4}=S+T_{1}+U_{1} italic_S start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT = italic_S + italic_T start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_U start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT . Since p a = x a P a subscript 𝑝 𝑎 subscript 𝑥 𝑎 subscript 𝑃 𝑎 p_{a}=x_{a}P_{a} italic_p start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = italic_x start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_P start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT and p b = x b P b subscript 𝑝 𝑏 subscript 𝑥 𝑏 subscript 𝑃 𝑏 p_{b}=x_{b}P_{b} italic_p start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = italic_x start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_P start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT , we have the relations s = x a x b S 𝑠 subscript 𝑥 𝑎 subscript 𝑥 𝑏 𝑆 s=x_{a}x_{b}S italic_s = italic_x start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_x start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_S , t 1 = x a T 1 subscript 𝑡 1 subscript 𝑥 𝑎 subscript 𝑇 1 t_{1}=x_{a}T_{1} italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = italic_x start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_T start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , u 1 = x b U 1 subscript 𝑢 1 subscript 𝑥 𝑏 subscript 𝑈 1 u_{1}=x_{b}U_{1} italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = italic_x start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_U start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , and
S 4 S = s 4 s − ( 1 − x a ) u 1 s − ( 1 − x b ) t 1 s . subscript 𝑆 4 𝑆 subscript 𝑠 4 𝑠 1 subscript 𝑥 𝑎 subscript 𝑢 1 𝑠 1 subscript 𝑥 𝑏 subscript 𝑡 1 𝑠 \frac{S_{4}}{S}=\frac{s_{4}}{s}-(1-x_{a})\frac{u_{1}}{s}-(1-x_{b})\frac{t_{1}}%
{s}\,. divide start_ARG italic_S start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT end_ARG start_ARG italic_S end_ARG = divide start_ARG italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT end_ARG start_ARG italic_s end_ARG - ( 1 - italic_x start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) divide start_ARG italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_s end_ARG - ( 1 - italic_x start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) divide start_ARG italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_s end_ARG .
(2.3)
We then consider the parton-parton cross section d σ q q ¯ → A B 𝑑 subscript 𝜎 → 𝑞 ¯ 𝑞 𝐴 𝐵 d\sigma_{q{\bar{q}}\to AB} italic_d italic_σ start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT , which is of the same form as Eq. (2.2 ) but with incoming partons instead of hadrons [3 , 4 , 5 ] ,
d σ q q ¯ → A B ( S 4 ) = ∫ 𝑑 x a 𝑑 x b ϕ q / q ( x a ) ϕ q ¯ / q ¯ ( x b ) 𝑑 σ ^ q q ¯ → A B ( s 4 ) , 𝑑 subscript 𝜎 → 𝑞 ¯ 𝑞 𝐴 𝐵 subscript 𝑆 4 differential-d subscript 𝑥 𝑎 differential-d subscript 𝑥 𝑏 subscript italic-ϕ 𝑞 𝑞 subscript 𝑥 𝑎 subscript italic-ϕ ¯ 𝑞 ¯ 𝑞 subscript 𝑥 𝑏 differential-d subscript ^ 𝜎 → 𝑞 ¯ 𝑞 𝐴 𝐵 subscript 𝑠 4 d\sigma_{q{\bar{q}}\to AB}(S_{4})=\int dx_{a}\,dx_{b}\,\phi_{q/q}(x_{a})\,\phi%
_{{\bar{q}}/{\bar{q}}}(x_{b})\,d{\hat{\sigma}}_{q{\bar{q}}\to AB}(s_{4})\,, italic_d italic_σ start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT ( italic_S start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ) = ∫ italic_d italic_x start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_d italic_x start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_ϕ start_POSTSUBSCRIPT italic_q / italic_q end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) italic_ϕ start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG / over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) italic_d over^ start_ARG italic_σ end_ARG start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ) ,
(2.4)
and we define its Laplace transform as
d σ ~ q q ¯ → A B ( N ) = ∫ 0 S d S 4 S e − N S 4 / S 𝑑 σ q q ¯ → A B ( S 4 ) . 𝑑 subscript ~ 𝜎 → 𝑞 ¯ 𝑞 𝐴 𝐵 𝑁 superscript subscript 0 𝑆 𝑑 subscript 𝑆 4 𝑆 superscript 𝑒 𝑁 subscript 𝑆 4 𝑆 differential-d subscript 𝜎 → 𝑞 ¯ 𝑞 𝐴 𝐵 subscript 𝑆 4 d{\tilde{\sigma}}_{q{\bar{q}}\to AB}(N)=\int_{0}^{S}\frac{dS_{4}}{S}\,e^{-NS_{%
4}/S}\,d\sigma_{q{\bar{q}}\to AB}(S_{4})\,. italic_d over~ start_ARG italic_σ end_ARG start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT ( italic_N ) = ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_S end_POSTSUPERSCRIPT divide start_ARG italic_d italic_S start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT end_ARG start_ARG italic_S end_ARG italic_e start_POSTSUPERSCRIPT - italic_N italic_S start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT / italic_S end_POSTSUPERSCRIPT italic_d italic_σ start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT ( italic_S start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ) .
(2.5)
Using the expression for S 4 / S subscript 𝑆 4 𝑆 S_{4}/S italic_S start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT / italic_S in Eq. (2.3 ), we can rewrite the Laplace transform, Eq. (2.5 ), of the parton-parton cross section in Eq. (2.4 ) as
d σ ~ q q ¯ → A B ( N ) 𝑑 subscript ~ 𝜎 → 𝑞 ¯ 𝑞 𝐴 𝐵 𝑁 \displaystyle d{\tilde{\sigma}}_{q{\bar{q}}\to AB}(N) italic_d over~ start_ARG italic_σ end_ARG start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT ( italic_N )
= \displaystyle= =
∫ 0 1 𝑑 x a e − N a ( 1 − x a ) ϕ q / q ( x a ) ∫ 0 1 𝑑 x b e − N b ( 1 − x b ) ϕ q ¯ / q ¯ ( x b ) ∫ 0 s d s 4 s e − N s 4 / s 𝑑 σ ^ q q ¯ → A B ( s 4 ) superscript subscript 0 1 differential-d subscript 𝑥 𝑎 superscript 𝑒 subscript 𝑁 𝑎 1 subscript 𝑥 𝑎 subscript italic-ϕ 𝑞 𝑞 subscript 𝑥 𝑎 superscript subscript 0 1 differential-d subscript 𝑥 𝑏 superscript 𝑒 subscript 𝑁 𝑏 1 subscript 𝑥 𝑏 subscript italic-ϕ ¯ 𝑞 ¯ 𝑞 subscript 𝑥 𝑏 superscript subscript 0 𝑠 𝑑 subscript 𝑠 4 𝑠 superscript 𝑒 𝑁 subscript 𝑠 4 𝑠 differential-d subscript ^ 𝜎 → 𝑞 ¯ 𝑞 𝐴 𝐵 subscript 𝑠 4 \displaystyle\int_{0}^{1}dx_{a}e^{-N_{a}(1-x_{a})}\phi_{q/q}(x_{a})\int_{0}^{1%
}dx_{b}e^{-N_{b}(1-x_{b})}\phi_{{\bar{q}}/{\bar{q}}}(x_{b})\int_{0}^{s}\frac{%
ds_{4}}{s}e^{-Ns_{4}/s}d{\hat{\sigma}}_{q{\bar{q}}\to AB}(s_{4}) ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT italic_d italic_x start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_N start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( 1 - italic_x start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) end_POSTSUPERSCRIPT italic_ϕ start_POSTSUBSCRIPT italic_q / italic_q end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT italic_d italic_x start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_N start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( 1 - italic_x start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) end_POSTSUPERSCRIPT italic_ϕ start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG / over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT ( italic_x start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT divide start_ARG italic_d italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT end_ARG start_ARG italic_s end_ARG italic_e start_POSTSUPERSCRIPT - italic_N italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT / italic_s end_POSTSUPERSCRIPT italic_d over^ start_ARG italic_σ end_ARG start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT )
(2.6)
= \displaystyle= =
ϕ ~ q / q ( N a ) ϕ ~ q ¯ / q ¯ ( N b ) d σ ^ ~ q q ¯ → A B ( N ) , subscript ~ italic-ϕ 𝑞 𝑞 subscript 𝑁 𝑎 subscript ~ italic-ϕ ¯ 𝑞 ¯ 𝑞 subscript 𝑁 𝑏 𝑑 subscript ~ ^ 𝜎 → 𝑞 ¯ 𝑞 𝐴 𝐵 𝑁 \displaystyle{\tilde{\phi}}_{q/q}(N_{a})\,{\tilde{\phi}}_{{\bar{q}}/{\bar{q}}}%
(N_{b})\,d{\tilde{\hat{\sigma}}}_{q{\bar{q}}\to AB}(N)\,, over~ start_ARG italic_ϕ end_ARG start_POSTSUBSCRIPT italic_q / italic_q end_POSTSUBSCRIPT ( italic_N start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) over~ start_ARG italic_ϕ end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG / over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT ( italic_N start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) italic_d over~ start_ARG over^ start_ARG italic_σ end_ARG end_ARG start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT ( italic_N ) ,
where N a = N ( − u 1 / s ) subscript 𝑁 𝑎 𝑁 subscript 𝑢 1 𝑠 N_{a}=N(-u_{1}/s) italic_N start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = italic_N ( - italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT / italic_s ) and N b = N ( − t 1 / s ) subscript 𝑁 𝑏 𝑁 subscript 𝑡 1 𝑠 N_{b}=N(-t_{1}/s) italic_N start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = italic_N ( - italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT / italic_s ) .
Next, we introduce a refactorization of the cross section via new functions H q q ¯ → A B subscript 𝐻 → 𝑞 ¯ 𝑞 𝐴 𝐵 H_{q{\bar{q}}\to AB} italic_H start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT , S q q ¯ → A B subscript 𝑆 → 𝑞 ¯ 𝑞 𝐴 𝐵 S_{q{\bar{q}}\to AB} italic_S start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT , ψ q / q subscript 𝜓 𝑞 𝑞 \psi_{q/q} italic_ψ start_POSTSUBSCRIPT italic_q / italic_q end_POSTSUBSCRIPT , ψ q ¯ / q ¯ subscript 𝜓 ¯ 𝑞 ¯ 𝑞 \psi_{{\bar{q}}/{\bar{q}}} italic_ψ start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG / over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT [3 , 4 , 5 ] . The hard function H q q ¯ → A B subscript 𝐻 → 𝑞 ¯ 𝑞 𝐴 𝐵 H_{q{\bar{q}}\to AB} italic_H start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT is purely short-distance and infrared safe while the soft function S q q ¯ → A B subscript 𝑆 → 𝑞 ¯ 𝑞 𝐴 𝐵 S_{q{\bar{q}}\to AB} italic_S start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT describes the emission of noncollinear soft gluons. The coupling of the soft gluons to the partons in the hard-scattering process is described by eikonal (Wilson) lines, i.e. ordered exponentials of the gauge field. The functions ψ q / q subscript 𝜓 𝑞 𝑞 \psi_{q/q} italic_ψ start_POSTSUBSCRIPT italic_q / italic_q end_POSTSUBSCRIPT and ψ q ¯ / q ¯ subscript 𝜓 ¯ 𝑞 ¯ 𝑞 \psi_{{\bar{q}}/{\bar{q}}} italic_ψ start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG / over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT differ from the pdf ϕ q / q subscript italic-ϕ 𝑞 𝑞 \phi_{q/q} italic_ϕ start_POSTSUBSCRIPT italic_q / italic_q end_POSTSUBSCRIPT and ϕ q ¯ / q ¯ subscript italic-ϕ ¯ 𝑞 ¯ 𝑞 \phi_{{\bar{q}}/{\bar{q}}} italic_ϕ start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG / over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT , and they describe collinear emission from the incoming partons [1 , 3 , 4 , 5 ] . The refactorized form of the cross section [5 , 29 ] is then
d σ q q ¯ → A B 𝑑 subscript 𝜎 → 𝑞 ¯ 𝑞 𝐴 𝐵 \displaystyle d{\sigma}_{q{\bar{q}}\to AB} italic_d italic_σ start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT
= \displaystyle= =
∫ 𝑑 w a 𝑑 w b 𝑑 w S ψ q / q ( w a ) ψ q ¯ / q ¯ ( w b ) differential-d subscript 𝑤 𝑎 differential-d subscript 𝑤 𝑏 differential-d subscript 𝑤 𝑆 subscript 𝜓 𝑞 𝑞 subscript 𝑤 𝑎 subscript 𝜓 ¯ 𝑞 ¯ 𝑞 subscript 𝑤 𝑏 \displaystyle\int dw_{a}\,dw_{b}\,dw_{S}\,\psi_{q/q}(w_{a})\,\psi_{{\bar{q}}/{%
\bar{q}}}(w_{b}) ∫ italic_d italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_d italic_w start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_d italic_w start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT italic_ψ start_POSTSUBSCRIPT italic_q / italic_q end_POSTSUBSCRIPT ( italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) italic_ψ start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG / over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT ( italic_w start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT )
(2.7)
× H q q ¯ → A B S q q ¯ → A B ( w S s μ F ) δ ( S 4 S − w S + w a u 1 s + w b t 1 s ) absent subscript 𝐻 → 𝑞 ¯ 𝑞 𝐴 𝐵 subscript 𝑆 → 𝑞 ¯ 𝑞 𝐴 𝐵 subscript 𝑤 𝑆 𝑠 subscript 𝜇 𝐹 𝛿 subscript 𝑆 4 𝑆 subscript 𝑤 𝑆 subscript 𝑤 𝑎 subscript 𝑢 1 𝑠 subscript 𝑤 𝑏 subscript 𝑡 1 𝑠 \displaystyle\times H_{q{\bar{q}}\to AB}\,\,S_{q{\bar{q}}\to AB}\left(\frac{w_%
{S}\sqrt{s}}{\mu_{F}}\right)\;\delta\left(\frac{S_{4}}{S}-w_{S}+w_{a}\frac{u_{%
1}}{s}+w_{b}\frac{t_{1}}{s}\right) × italic_H start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT italic_S start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT ( divide start_ARG italic_w start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT square-root start_ARG italic_s end_ARG end_ARG start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT end_ARG ) italic_δ ( divide start_ARG italic_S start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT end_ARG start_ARG italic_S end_ARG - italic_w start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT + italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT divide start_ARG italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_s end_ARG + italic_w start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT divide start_ARG italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_s end_ARG )
where the w 𝑤 w italic_w ’s are dimensionless weights, with w a subscript 𝑤 𝑎 w_{a} italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT and w b subscript 𝑤 𝑏 w_{b} italic_w start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT for ψ q / q subscript 𝜓 𝑞 𝑞 \psi_{q/q} italic_ψ start_POSTSUBSCRIPT italic_q / italic_q end_POSTSUBSCRIPT and ψ q ¯ / q ¯ subscript 𝜓 ¯ 𝑞 ¯ 𝑞 \psi_{{\bar{q}}/{\bar{q}}} italic_ψ start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG / over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT , respectively, and w S subscript 𝑤 𝑆 w_{S} italic_w start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT for S q q ¯ → A B subscript 𝑆 → 𝑞 ¯ 𝑞 𝐴 𝐵 S_{q{\bar{q}}\to AB} italic_S start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT . The argument in the delta function of Eq. (2.7 ) arises from rewriting Eq. (2.3 ) in terms of the new weights [5 ] , as
S 4 S = w S − w a u 1 s − w b t 1 s . subscript 𝑆 4 𝑆 subscript 𝑤 𝑆 subscript 𝑤 𝑎 subscript 𝑢 1 𝑠 subscript 𝑤 𝑏 subscript 𝑡 1 𝑠 \frac{S_{4}}{S}=w_{S}-w_{a}\frac{u_{1}}{s}-w_{b}\frac{t_{1}}{s}\,. divide start_ARG italic_S start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT end_ARG start_ARG italic_S end_ARG = italic_w start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT - italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT divide start_ARG italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_s end_ARG - italic_w start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT divide start_ARG italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_s end_ARG .
(2.8)
After taking a Laplace transform of Eq. (2.7 ), we have
d σ ~ q q ¯ → A B ( N ) 𝑑 subscript ~ 𝜎 → 𝑞 ¯ 𝑞 𝐴 𝐵 𝑁 \displaystyle d{\tilde{\sigma}}_{q{\bar{q}}\to AB}(N) italic_d over~ start_ARG italic_σ end_ARG start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT ( italic_N )
= \displaystyle= =
∫ 0 1 𝑑 w a e − N a w a ψ q / q ( w a ) ∫ 0 1 𝑑 w b e − N b w b ψ q ¯ / q ¯ ( w b ) superscript subscript 0 1 differential-d subscript 𝑤 𝑎 superscript 𝑒 subscript 𝑁 𝑎 subscript 𝑤 𝑎 subscript 𝜓 𝑞 𝑞 subscript 𝑤 𝑎 superscript subscript 0 1 differential-d subscript 𝑤 𝑏 superscript 𝑒 subscript 𝑁 𝑏 subscript 𝑤 𝑏 subscript 𝜓 ¯ 𝑞 ¯ 𝑞 subscript 𝑤 𝑏 \displaystyle\int_{0}^{1}dw_{a}\,e^{-N_{a}w_{a}}\,\psi_{q/q}(w_{a})\int_{0}^{1%
}dw_{b}\,e^{-N_{b}w_{b}}\,\psi_{{\bar{q}}/{\bar{q}}}(w_{b}) ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT italic_d italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_N start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_ψ start_POSTSUBSCRIPT italic_q / italic_q end_POSTSUBSCRIPT ( italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT italic_d italic_w start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_N start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_w start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_ψ start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG / over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT ( italic_w start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT )
(2.9)
× H q q ¯ → A B ∫ 0 1 𝑑 w s e − N w s S q q ¯ → A B ( w s s μ F ) absent subscript 𝐻 → 𝑞 ¯ 𝑞 𝐴 𝐵 superscript subscript 0 1 differential-d subscript 𝑤 𝑠 superscript 𝑒 𝑁 subscript 𝑤 𝑠 subscript 𝑆 → 𝑞 ¯ 𝑞 𝐴 𝐵 subscript 𝑤 𝑠 𝑠 subscript 𝜇 𝐹 \displaystyle\times\,\,H_{q{\bar{q}}\to AB}\int_{0}^{1}dw_{s}\,e^{-Nw_{s}}\,S_%
{q{\bar{q}}\to AB}\left(\frac{w_{s}\sqrt{s}}{\mu_{F}}\right) × italic_H start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT italic_d italic_w start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - italic_N italic_w start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_S start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT ( divide start_ARG italic_w start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT square-root start_ARG italic_s end_ARG end_ARG start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT end_ARG )
= \displaystyle= =
ψ ~ q / q ( N a ) ψ ~ q ¯ / q ¯ ( N b ) H q q ¯ → A B S ~ q q ¯ → A B ( s N μ F ) . subscript ~ 𝜓 𝑞 𝑞 subscript 𝑁 𝑎 subscript ~ 𝜓 ¯ 𝑞 ¯ 𝑞 subscript 𝑁 𝑏 subscript 𝐻 → 𝑞 ¯ 𝑞 𝐴 𝐵 subscript ~ 𝑆 → 𝑞 ¯ 𝑞 𝐴 𝐵 𝑠 𝑁 subscript 𝜇 𝐹 \displaystyle{\tilde{\psi}}_{q/q}(N_{a})\,\,{\tilde{\psi}}_{{\bar{q}}/{\bar{q}%
}}(N_{b})\,\,H_{q{\bar{q}}\to AB}\,\,{\tilde{S}}_{q{\bar{q}}\to AB}\left(\frac%
{\sqrt{s}}{N\mu_{F}}\right)\,. over~ start_ARG italic_ψ end_ARG start_POSTSUBSCRIPT italic_q / italic_q end_POSTSUBSCRIPT ( italic_N start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) over~ start_ARG italic_ψ end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG / over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT ( italic_N start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) italic_H start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT over~ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT ( divide start_ARG square-root start_ARG italic_s end_ARG end_ARG start_ARG italic_N italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT end_ARG ) .
All N 𝑁 N italic_N -dependence has now been absorbed into the functions ψ ~ q / q subscript ~ 𝜓 𝑞 𝑞 {\tilde{\psi}}_{q/q} over~ start_ARG italic_ψ end_ARG start_POSTSUBSCRIPT italic_q / italic_q end_POSTSUBSCRIPT , ψ ~ q ¯ / q ¯ subscript ~ 𝜓 ¯ 𝑞 ¯ 𝑞 {\tilde{\psi}}_{{\bar{q}}/{\bar{q}}} over~ start_ARG italic_ψ end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG / over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT , and S ~ ~ 𝑆 {\tilde{S}} over~ start_ARG italic_S end_ARG .
Comparing Eqs. (2.6 ) and (2.9 ), we get an expression for the hard-scattering partonic cross section in Laplace transform space,
d σ ^ ~ q q ¯ → A B ( N ) = ψ ~ q / q ( N a ) ψ ~ q ¯ / q ¯ ( N b ) ϕ ~ q / q ( N a ) ϕ ~ q ¯ / q ¯ ( N b ) H q q ¯ → A B S ~ q q ¯ → A B ( s N μ F ) . 𝑑 subscript ~ ^ 𝜎 → 𝑞 ¯ 𝑞 𝐴 𝐵 𝑁 subscript ~ 𝜓 𝑞 𝑞 subscript 𝑁 𝑎 subscript ~ 𝜓 ¯ 𝑞 ¯ 𝑞 subscript 𝑁 𝑏 subscript ~ italic-ϕ 𝑞 𝑞 subscript 𝑁 𝑎 subscript ~ italic-ϕ ¯ 𝑞 ¯ 𝑞 subscript 𝑁 𝑏 subscript 𝐻 → 𝑞 ¯ 𝑞 𝐴 𝐵 subscript ~ 𝑆 → 𝑞 ¯ 𝑞 𝐴 𝐵 𝑠 𝑁 subscript 𝜇 𝐹 d{\tilde{\hat{\sigma}}}_{q{\bar{q}}\to AB}(N)=\frac{{\tilde{\psi}}_{q/q}(N_{a}%
)\,{\tilde{\psi}}_{{\bar{q}}/{\bar{q}}}(N_{b})}{{\tilde{\phi}}_{q/q}(N_{a})\,{%
\tilde{\phi}}_{{\bar{q}}/{\bar{q}}}(N_{b})}\,\,H_{q{\bar{q}}\to AB}\,\,{\tilde%
{S}}_{q{\bar{q}}\to AB}\left(\frac{\sqrt{s}}{N\mu_{F}}\right)\,. italic_d over~ start_ARG over^ start_ARG italic_σ end_ARG end_ARG start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT ( italic_N ) = divide start_ARG over~ start_ARG italic_ψ end_ARG start_POSTSUBSCRIPT italic_q / italic_q end_POSTSUBSCRIPT ( italic_N start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) over~ start_ARG italic_ψ end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG / over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT ( italic_N start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) end_ARG start_ARG over~ start_ARG italic_ϕ end_ARG start_POSTSUBSCRIPT italic_q / italic_q end_POSTSUBSCRIPT ( italic_N start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) over~ start_ARG italic_ϕ end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG / over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT ( italic_N start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) end_ARG italic_H start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT over~ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT ( divide start_ARG square-root start_ARG italic_s end_ARG end_ARG start_ARG italic_N italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT end_ARG ) .
(2.10)
The resummed differential cross section in Laplace-transform space is derived from the renormalization-group evolution of ψ ~ q / q / ϕ ~ q / q subscript ~ 𝜓 𝑞 𝑞 subscript ~ italic-ϕ 𝑞 𝑞 {\tilde{\psi}}_{q/q}/{\tilde{\phi}}_{q/q} over~ start_ARG italic_ψ end_ARG start_POSTSUBSCRIPT italic_q / italic_q end_POSTSUBSCRIPT / over~ start_ARG italic_ϕ end_ARG start_POSTSUBSCRIPT italic_q / italic_q end_POSTSUBSCRIPT , ψ ~ q ¯ / q ¯ / ϕ ~ q ¯ / q ¯ subscript ~ 𝜓 ¯ 𝑞 ¯ 𝑞 subscript ~ italic-ϕ ¯ 𝑞 ¯ 𝑞 {\tilde{\psi}}_{{\bar{q}}/{\bar{q}}}/{\tilde{\phi}}_{{\bar{q}}/{\bar{q}}} over~ start_ARG italic_ψ end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG / over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT / over~ start_ARG italic_ϕ end_ARG start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG / over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT , and S ~ q q ¯ → A B subscript ~ 𝑆 → 𝑞 ¯ 𝑞 𝐴 𝐵 {\tilde{S}}_{q{\bar{q}}\to AB} over~ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT in Eq. (2.10 ), and it is given by the expression [1 , 2 , 3 , 4 , 5 ]
d σ ^ ~ q q ¯ → A B r e s ( N ) 𝑑 subscript superscript ~ ^ 𝜎 𝑟 𝑒 𝑠 → 𝑞 ¯ 𝑞 𝐴 𝐵 𝑁 \displaystyle d{\tilde{\hat{\sigma}}}^{res}_{q{\bar{q}}\to AB}(N) italic_d over~ start_ARG over^ start_ARG italic_σ end_ARG end_ARG start_POSTSUPERSCRIPT italic_r italic_e italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT ( italic_N )
= \displaystyle= =
exp [ E q ( N a ) + E q ¯ ( N b ) ] exp [ 2 ∫ μ F s d μ μ ( γ q / q ( N a , α s ( μ ) ) + γ q ¯ / q ¯ ( N b , α s ( μ ) ) ) ] subscript 𝐸 𝑞 subscript 𝑁 𝑎 subscript 𝐸 ¯ 𝑞 subscript 𝑁 𝑏 2 superscript subscript subscript 𝜇 𝐹 𝑠 𝑑 𝜇 𝜇 subscript 𝛾 𝑞 𝑞 subscript 𝑁 𝑎 subscript 𝛼 𝑠 𝜇 subscript 𝛾 ¯ 𝑞 ¯ 𝑞 subscript 𝑁 𝑏 subscript 𝛼 𝑠 𝜇 \displaystyle\exp\left[E_{q}(N_{a})+E_{\bar{q}}(N_{b})\right]\;\,\exp\left[2%
\int_{\mu_{F}}^{\sqrt{s}}\frac{d\mu}{\mu}\;\left(\gamma_{q/q}\left(N_{a},%
\alpha_{s}(\mu)\right)+\gamma_{{\bar{q}}/{\bar{q}}}\left(N_{b},\alpha_{s}(\mu)%
\right)\right)\right] roman_exp [ italic_E start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT ( italic_N start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) + italic_E start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT ( italic_N start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) ] roman_exp [ 2 ∫ start_POSTSUBSCRIPT italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT square-root start_ARG italic_s end_ARG end_POSTSUPERSCRIPT divide start_ARG italic_d italic_μ end_ARG start_ARG italic_μ end_ARG ( italic_γ start_POSTSUBSCRIPT italic_q / italic_q end_POSTSUBSCRIPT ( italic_N start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_μ ) ) + italic_γ start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG / over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT ( italic_N start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT , italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_μ ) ) ) ]
(2.11)
× H q q ¯ → A B ( α s ( s ) ) S ~ q q ¯ → A B ( α s ( s / N ) ) . absent subscript 𝐻 → 𝑞 ¯ 𝑞 𝐴 𝐵 subscript 𝛼 𝑠 𝑠 subscript ~ 𝑆 → 𝑞 ¯ 𝑞 𝐴 𝐵 subscript 𝛼 𝑠 𝑠 𝑁 \displaystyle\times\,\,H_{q{\bar{q}}\to AB}\left(\alpha_{s}({\sqrt{s}})\right)%
\,\,{\tilde{S}}_{q{\bar{q}}\to AB}\left(\alpha_{s}({\sqrt{s}/N})\right)\,. × italic_H start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT ( italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( square-root start_ARG italic_s end_ARG ) ) over~ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT ( italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( square-root start_ARG italic_s end_ARG / italic_N ) ) .
In the first exponential of Eq. (2.11 ), we have
E q ( N a ) = ∫ 0 1 𝑑 z z N a − 1 − 1 1 − z { ∫ 1 ( 1 − z ) 2 d λ λ A q ( α s ( λ s ) ) + D q [ α s ( ( 1 − z ) 2 s ) ] } , subscript 𝐸 𝑞 subscript 𝑁 𝑎 subscript superscript 1 0 differential-d 𝑧 superscript 𝑧 subscript 𝑁 𝑎 1 1 1 𝑧 superscript subscript 1 superscript 1 𝑧 2 𝑑 𝜆 𝜆 subscript 𝐴 𝑞 subscript 𝛼 𝑠 𝜆 𝑠 subscript 𝐷 𝑞 delimited-[] subscript 𝛼 𝑠 superscript 1 𝑧 2 𝑠 E_{q}(N_{a})=\int^{1}_{0}dz\frac{z^{N_{a}-1}-1}{1-z}\;\left\{\int_{1}^{(1-z)^{%
2}}\frac{d\lambda}{\lambda}A_{q}\left(\alpha_{s}(\lambda s)\right)+D_{q}\left[%
\alpha_{s}((1-z)^{2}s)\right]\right\}\,, italic_E start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT ( italic_N start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) = ∫ start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_d italic_z divide start_ARG italic_z start_POSTSUPERSCRIPT italic_N start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT - 1 end_POSTSUPERSCRIPT - 1 end_ARG start_ARG 1 - italic_z end_ARG { ∫ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 - italic_z ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT divide start_ARG italic_d italic_λ end_ARG start_ARG italic_λ end_ARG italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT ( italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_λ italic_s ) ) + italic_D start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT [ italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( ( 1 - italic_z ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_s ) ] } ,
(2.12)
with an analogous expression for E q ¯ ( N b ) subscript 𝐸 ¯ 𝑞 subscript 𝑁 𝑏 E_{\bar{q}}(N_{b}) italic_E start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT ( italic_N start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) [1 ] .
The first integrand in Eq. (2.12 ) has the perturbative expansion
A q = ( α s / π ) A q ( 1 ) + ( α s / π ) 2 A q ( 2 ) + ( α s / π ) 3 A q ( 3 ) + ⋯ subscript 𝐴 𝑞 subscript 𝛼 𝑠 𝜋 superscript subscript 𝐴 𝑞 1 superscript subscript 𝛼 𝑠 𝜋 2 superscript subscript 𝐴 𝑞 2 superscript subscript 𝛼 𝑠 𝜋 3 superscript subscript 𝐴 𝑞 3 ⋯ A_{q}=(\alpha_{s}/\pi)A_{q}^{(1)}+(\alpha_{s}/\pi)^{2}A_{q}^{(2)}+(\alpha_{s}/%
\pi)^{3}A_{q}^{(3)}+\cdots italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT = ( italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT / italic_π ) italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT + ( italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT / italic_π ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT + ( italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT / italic_π ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT + ⋯ .
The first term in this expansion is A q ( 1 ) = C F superscript subscript 𝐴 𝑞 1 subscript 𝐶 𝐹 A_{q}^{(1)}=C_{F} italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT = italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT [1 ] , where C F = ( N c 2 − 1 ) / ( 2 N c ) subscript 𝐶 𝐹 superscript subscript 𝑁 𝑐 2 1 2 subscript 𝑁 𝑐 C_{F}=(N_{c}^{2}-1)/(2N_{c}) italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT = ( italic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 1 ) / ( 2 italic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ) and N c = 3 subscript 𝑁 𝑐 3 N_{c}=3 italic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = 3 is the number of colors. The second term, A q ( 2 ) superscript subscript 𝐴 𝑞 2 A_{q}^{(2)} italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT , is given by [2 , 30 ]
A q ( 2 ) = C F C A ( 67 36 − ζ 2 2 ) − 5 18 C F n f , superscript subscript 𝐴 𝑞 2 subscript 𝐶 𝐹 subscript 𝐶 𝐴 67 36 subscript 𝜁 2 2 5 18 subscript 𝐶 𝐹 subscript 𝑛 𝑓 A_{q}^{(2)}=C_{F}C_{A}\left(\frac{67}{36}-\frac{\zeta_{2}}{2}\right)-\frac{5}{%
18}C_{F}n_{f}\,, italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT = italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ( divide start_ARG 67 end_ARG start_ARG 36 end_ARG - divide start_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG ) - divide start_ARG 5 end_ARG start_ARG 18 end_ARG italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT ,
(2.13)
where C A = N c subscript 𝐶 𝐴 subscript 𝑁 𝑐 C_{A}=N_{c} italic_C start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT = italic_N start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT and n f subscript 𝑛 𝑓 n_{f} italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT is the number of light-quark flavors (in the next section we will set n f = 5 subscript 𝑛 𝑓 5 n_{f}=5 italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT = 5 for H + H − superscript 𝐻 superscript 𝐻 H^{+}H^{-} italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT production). Finally, the third term, A q ( 3 ) superscript subscript 𝐴 𝑞 3 A_{q}^{(3)} italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT , is given by [31 ]
A q ( 3 ) superscript subscript 𝐴 𝑞 3 \displaystyle A_{q}^{(3)} italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT
= \displaystyle= =
C F C A 2 ( 245 96 − 67 36 ζ 2 + 11 24 ζ 3 + 11 8 ζ 4 ) + C F 2 n f ( − 55 96 + ζ 3 2 ) subscript 𝐶 𝐹 superscript subscript 𝐶 𝐴 2 245 96 67 36 subscript 𝜁 2 11 24 subscript 𝜁 3 11 8 subscript 𝜁 4 superscript subscript 𝐶 𝐹 2 subscript 𝑛 𝑓 55 96 subscript 𝜁 3 2 \displaystyle C_{F}C_{A}^{2}\left(\frac{245}{96}-\frac{67}{36}\zeta_{2}+\frac{%
11}{24}\zeta_{3}+\frac{11}{8}\zeta_{4}\right)+C_{F}^{2}n_{f}\left(-\frac{55}{9%
6}+\frac{\zeta_{3}}{2}\right) italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG 245 end_ARG start_ARG 96 end_ARG - divide start_ARG 67 end_ARG start_ARG 36 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + divide start_ARG 11 end_ARG start_ARG 24 end_ARG italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT + divide start_ARG 11 end_ARG start_ARG 8 end_ARG italic_ζ start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ) + italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT ( - divide start_ARG 55 end_ARG start_ARG 96 end_ARG + divide start_ARG italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG )
(2.14)
+ C F C A n f ( − 209 432 + 5 18 ζ 2 − 7 12 ζ 3 ) − C F n f 2 108 . subscript 𝐶 𝐹 subscript 𝐶 𝐴 subscript 𝑛 𝑓 209 432 5 18 subscript 𝜁 2 7 12 subscript 𝜁 3 subscript 𝐶 𝐹 superscript subscript 𝑛 𝑓 2 108 \displaystyle{}+C_{F}C_{A}n_{f}\left(-\frac{209}{432}+\frac{5}{18}\zeta_{2}-%
\frac{7}{12}\zeta_{3}\right)-C_{F}\frac{n_{f}^{2}}{108}\,. + italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT ( - divide start_ARG 209 end_ARG start_ARG 432 end_ARG + divide start_ARG 5 end_ARG start_ARG 18 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - divide start_ARG 7 end_ARG start_ARG 12 end_ARG italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ) - italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT divide start_ARG italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 108 end_ARG .
Here and below ζ 2 = π 2 / 6 subscript 𝜁 2 superscript 𝜋 2 6 \zeta_{2}=\pi^{2}/6 italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 6 , ζ 3 = 1.202056903 ⋯ subscript 𝜁 3 1.202056903 ⋯ \zeta_{3}=1.202056903\cdots italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT = 1.202056903 ⋯ ,
ζ 4 = π 4 / 90 subscript 𝜁 4 superscript 𝜋 4 90 \zeta_{4}=\pi^{4}/90 italic_ζ start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT = italic_π start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT / 90 , and ζ 5 = 1.036927755 ⋯ subscript 𝜁 5 1.036927755 ⋯ \zeta_{5}=1.036927755\cdots italic_ζ start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT = 1.036927755 ⋯ .
The corresponding expressions for processes with gluons in the initial state (i.e. g g → A B → 𝑔 𝑔 𝐴 𝐵 gg\to AB italic_g italic_g → italic_A italic_B ) are given at these perturbative orders by
A g ( n ) = ( C A / C F ) A q ( n ) superscript subscript 𝐴 𝑔 𝑛 subscript 𝐶 𝐴 subscript 𝐶 𝐹 superscript subscript 𝐴 𝑞 𝑛 A_{g}^{(n)}=(C_{A}/C_{F})A_{q}^{(n)} italic_A start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT = ( italic_C start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT / italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ) italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT , where n = 1 , 2 , 3 𝑛 1 2 3
n=1,2,3 italic_n = 1 , 2 , 3 .
The second integrand in Eq. (2.12 ) can be expanded as D q = ( α s / π ) D q ( 1 ) + ( α s / π ) 2 D q ( 2 ) + ( α s / π ) 3 D q ( 3 ) + ⋯ subscript 𝐷 𝑞 subscript 𝛼 𝑠 𝜋 superscript subscript 𝐷 𝑞 1 superscript subscript 𝛼 𝑠 𝜋 2 superscript subscript 𝐷 𝑞 2 superscript subscript 𝛼 𝑠 𝜋 3 superscript subscript 𝐷 𝑞 3 ⋯ D_{q}=(\alpha_{s}/\pi)D_{q}^{(1)}+(\alpha_{s}/\pi)^{2}D_{q}^{(2)}+(\alpha_{s}/%
\pi)^{3}D_{q}^{(3)}+\cdots italic_D start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT = ( italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT / italic_π ) italic_D start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT + ( italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT / italic_π ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_D start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT + ( italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT / italic_π ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_D start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT + ⋯ . In Feynman gauge, the one-loop term D q ( 1 ) = 0 superscript subscript 𝐷 𝑞 1 0 D_{q}^{(1)}=0 italic_D start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT = 0 , while
at two loops [32 ]
D q ( 2 ) = C F C A ( − 101 54 + 11 6 ζ 2 + 7 4 ζ 3 ) + C F n f ( 7 27 − ζ 2 3 ) , superscript subscript 𝐷 𝑞 2 subscript 𝐶 𝐹 subscript 𝐶 𝐴 101 54 11 6 subscript 𝜁 2 7 4 subscript 𝜁 3 subscript 𝐶 𝐹 subscript 𝑛 𝑓 7 27 subscript 𝜁 2 3 D_{q}^{(2)}=C_{F}C_{A}\left(-\frac{101}{54}+\frac{11}{6}\zeta_{2}+\frac{7}{4}%
\zeta_{3}\right)+C_{F}n_{f}\left(\frac{7}{27}-\frac{\zeta_{2}}{3}\right)\,, italic_D start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT = italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ( - divide start_ARG 101 end_ARG start_ARG 54 end_ARG + divide start_ARG 11 end_ARG start_ARG 6 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + divide start_ARG 7 end_ARG start_ARG 4 end_ARG italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ) + italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT ( divide start_ARG 7 end_ARG start_ARG 27 end_ARG - divide start_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG 3 end_ARG ) ,
(2.15)
and at three loops [9 ]
D q ( 3 ) superscript subscript 𝐷 𝑞 3 \displaystyle D_{q}^{(3)} italic_D start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT
= \displaystyle= =
C F C A 2 ( − 297029 46656 + 6139 648 ζ 2 + 2509 216 ζ 3 − 187 48 ζ 4 − 11 12 ζ 2 ζ 3 − 3 ζ 5 ) subscript 𝐶 𝐹 superscript subscript 𝐶 𝐴 2 297029 46656 6139 648 subscript 𝜁 2 2509 216 subscript 𝜁 3 187 48 subscript 𝜁 4 11 12 subscript 𝜁 2 subscript 𝜁 3 3 subscript 𝜁 5 \displaystyle C_{F}C_{A}^{2}\left(-\frac{297029}{46656}+\frac{6139}{648}\zeta_%
{2}+\frac{2509}{216}\zeta_{3}-\frac{187}{48}\zeta_{4}-\frac{11}{12}\zeta_{2}%
\zeta_{3}-3\zeta_{5}\right) italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( - divide start_ARG 297029 end_ARG start_ARG 46656 end_ARG + divide start_ARG 6139 end_ARG start_ARG 648 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + divide start_ARG 2509 end_ARG start_ARG 216 end_ARG italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT - divide start_ARG 187 end_ARG start_ARG 48 end_ARG italic_ζ start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT - divide start_ARG 11 end_ARG start_ARG 12 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT - 3 italic_ζ start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT )
(2.16)
+ C F C A n f ( 31313 23328 − 1837 648 ζ 2 − 155 72 ζ 3 + 23 24 ζ 4 ) subscript 𝐶 𝐹 subscript 𝐶 𝐴 subscript 𝑛 𝑓 31313 23328 1837 648 subscript 𝜁 2 155 72 subscript 𝜁 3 23 24 subscript 𝜁 4 \displaystyle{}+C_{F}C_{A}n_{f}\left(\frac{31313}{23328}-\frac{1837}{648}\zeta%
_{2}-\frac{155}{72}\zeta_{3}+\frac{23}{24}\zeta_{4}\right) + italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT ( divide start_ARG 31313 end_ARG start_ARG 23328 end_ARG - divide start_ARG 1837 end_ARG start_ARG 648 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - divide start_ARG 155 end_ARG start_ARG 72 end_ARG italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT + divide start_ARG 23 end_ARG start_ARG 24 end_ARG italic_ζ start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT )
+ C F 2 n f ( 1711 1728 − ζ 2 4 − 19 36 ζ 3 − ζ 4 4 ) + C F n f 2 ( − 29 729 + 5 27 ζ 2 + 5 54 ζ 3 ) . superscript subscript 𝐶 𝐹 2 subscript 𝑛 𝑓 1711 1728 subscript 𝜁 2 4 19 36 subscript 𝜁 3 subscript 𝜁 4 4 subscript 𝐶 𝐹 superscript subscript 𝑛 𝑓 2 29 729 5 27 subscript 𝜁 2 5 54 subscript 𝜁 3 \displaystyle{}+C_{F}^{2}n_{f}\left(\frac{1711}{1728}-\frac{\zeta_{2}}{4}-%
\frac{19}{36}\zeta_{3}-\frac{\zeta_{4}}{4}\right)+C_{F}n_{f}^{2}\left(-\frac{2%
9}{729}+\frac{5}{27}\zeta_{2}+\frac{5}{54}\zeta_{3}\right)\,. + italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT ( divide start_ARG 1711 end_ARG start_ARG 1728 end_ARG - divide start_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG - divide start_ARG 19 end_ARG start_ARG 36 end_ARG italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT - divide start_ARG italic_ζ start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG ) + italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( - divide start_ARG 29 end_ARG start_ARG 729 end_ARG + divide start_ARG 5 end_ARG start_ARG 27 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + divide start_ARG 5 end_ARG start_ARG 54 end_ARG italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ) .
The corresponding expressions for processes with gluons in the initial state are given at these perturbative orders by
D g ( n ) = ( C A / C F ) D q ( n ) superscript subscript 𝐷 𝑔 𝑛 subscript 𝐶 𝐴 subscript 𝐶 𝐹 superscript subscript 𝐷 𝑞 𝑛 D_{g}^{(n)}=(C_{A}/C_{F})D_{q}^{(n)} italic_D start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT = ( italic_C start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT / italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ) italic_D start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT , where n = 1 , 2 , 3 𝑛 1 2 3
n=1,2,3 italic_n = 1 , 2 , 3 .
In the second exponent of Eq. (2.11 ), the quantity γ q / q subscript 𝛾 𝑞 𝑞 \gamma_{q/q} italic_γ start_POSTSUBSCRIPT italic_q / italic_q end_POSTSUBSCRIPT is the moment-space
anomalous dimension of the MS ¯ ¯ MS {\overline{\rm MS}} over¯ start_ARG roman_MS end_ARG density ϕ q / q subscript italic-ϕ 𝑞 𝑞 \phi_{q/q} italic_ϕ start_POSTSUBSCRIPT italic_q / italic_q end_POSTSUBSCRIPT
[33 , 34 , 35 , 36 , 37 ] which can be expressed as γ q / q ( N ) = − A q ln N ~ + γ q subscript 𝛾 𝑞 𝑞 𝑁 subscript 𝐴 𝑞 ~ 𝑁 subscript 𝛾 𝑞 \gamma_{q/q}(N)=-A_{q}\ln{\tilde{N}}+\gamma_{q} italic_γ start_POSTSUBSCRIPT italic_q / italic_q end_POSTSUBSCRIPT ( italic_N ) = - italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT roman_ln over~ start_ARG italic_N end_ARG + italic_γ start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT , and similarly for γ q ¯ / q ¯ subscript 𝛾 ¯ 𝑞 ¯ 𝑞 \gamma_{{\bar{q}}/{\bar{q}}} italic_γ start_POSTSUBSCRIPT over¯ start_ARG italic_q end_ARG / over¯ start_ARG italic_q end_ARG end_POSTSUBSCRIPT . The parton anomalous dimension can be expanded as
γ q = ( α s / π ) γ q ( 1 ) + ( α s / π ) 2 γ q ( 2 ) + ⋯ subscript 𝛾 𝑞 subscript 𝛼 𝑠 𝜋 superscript subscript 𝛾 𝑞 1 superscript subscript 𝛼 𝑠 𝜋 2 superscript subscript 𝛾 𝑞 2 ⋯ \gamma_{q}=(\alpha_{s}/\pi)\gamma_{q}^{(1)}+(\alpha_{s}/\pi)^{2}\gamma_{q}^{(2%
)}+\cdots italic_γ start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT = ( italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT / italic_π ) italic_γ start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT + ( italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT / italic_π ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_γ start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT + ⋯
with γ q ( 1 ) = 3 C F / 4 superscript subscript 𝛾 𝑞 1 3 subscript 𝐶 𝐹 4 \gamma_{q}^{(1)}=3C_{F}/4 italic_γ start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT = 3 italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT / 4 and
γ q ( 2 ) = C F 2 ( 3 32 − 3 4 ζ 2 + 3 2 ζ 3 ) + C F C A ( 17 96 + 11 12 ζ 2 − 3 4 ζ 3 ) − C F n f ( 1 48 + ζ 2 6 ) . superscript subscript 𝛾 𝑞 2 superscript subscript 𝐶 𝐹 2 3 32 3 4 subscript 𝜁 2 3 2 subscript 𝜁 3 subscript 𝐶 𝐹 subscript 𝐶 𝐴 17 96 11 12 subscript 𝜁 2 3 4 subscript 𝜁 3 subscript 𝐶 𝐹 subscript 𝑛 𝑓 1 48 subscript 𝜁 2 6 \gamma_{q}^{(2)}=C_{F}^{2}\left(\frac{3}{32}-\frac{3}{4}\zeta_{2}+\frac{3}{2}%
\zeta_{3}\right)+C_{F}C_{A}\left(\frac{17}{96}+\frac{11}{12}\zeta_{2}-\frac{3}%
{4}\zeta_{3}\right)-C_{F}n_{f}\left(\frac{1}{48}+\frac{\zeta_{2}}{6}\right)\,. italic_γ start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT = italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG 3 end_ARG start_ARG 32 end_ARG - divide start_ARG 3 end_ARG start_ARG 4 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + divide start_ARG 3 end_ARG start_ARG 2 end_ARG italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ) + italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ( divide start_ARG 17 end_ARG start_ARG 96 end_ARG + divide start_ARG 11 end_ARG start_ARG 12 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - divide start_ARG 3 end_ARG start_ARG 4 end_ARG italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ) - italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT ( divide start_ARG 1 end_ARG start_ARG 48 end_ARG + divide start_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG 6 end_ARG ) .
(2.17)
For gluons, we have γ g ( 1 ) = β 0 / 4 superscript subscript 𝛾 𝑔 1 subscript 𝛽 0 4 \gamma_{g}^{(1)}=\beta_{0}/4 italic_γ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT = italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT / 4 , where β 0 = 11 C A / 3 − 2 n f / 3 subscript 𝛽 0 11 subscript 𝐶 𝐴 3 2 subscript 𝑛 𝑓 3 \beta_{0}=11C_{A}/3-2n_{f}/3 italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 11 italic_C start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT / 3 - 2 italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT / 3 is the coefficient of the first term in the perturbative expansion of the β 𝛽 \beta italic_β function [38 , 39 ] , and
γ g ( 2 ) = C A 2 ( 2 3 + 3 4 ζ 3 ) − n f ( C F 8 + C A 6 ) . superscript subscript 𝛾 𝑔 2 superscript subscript 𝐶 𝐴 2 2 3 3 4 subscript 𝜁 3 subscript 𝑛 𝑓 subscript 𝐶 𝐹 8 subscript 𝐶 𝐴 6 \gamma_{g}^{(2)}=C_{A}^{2}\left(\frac{2}{3}+\frac{3}{4}\zeta_{3}\right)-n_{f}%
\left(\frac{C_{F}}{8}+\frac{C_{A}}{6}\right)\,. italic_γ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT = italic_C start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG 2 end_ARG start_ARG 3 end_ARG + divide start_ARG 3 end_ARG start_ARG 4 end_ARG italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ) - italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT ( divide start_ARG italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT end_ARG start_ARG 8 end_ARG + divide start_ARG italic_C start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT end_ARG start_ARG 6 end_ARG ) .
(2.18)
By expanding Eq. (2.11 ) to fixed order and inverting back to momentum space, we can calculate the soft-gluon corrections completely through N3 LO.
The Born differential cross section at partonic level can be written as
d 2 σ ^ q q ¯ → A B ( 0 ) d t 1 d u 1 = F q q ¯ → A B B δ ( s 4 ) , superscript 𝑑 2 subscript superscript ^ 𝜎 0 → 𝑞 ¯ 𝑞 𝐴 𝐵 𝑑 subscript 𝑡 1 𝑑 subscript 𝑢 1 subscript superscript 𝐹 𝐵 → 𝑞 ¯ 𝑞 𝐴 𝐵 𝛿 subscript 𝑠 4 \frac{d^{2}{\hat{\sigma}}^{(0)}_{q{\bar{q}}\to AB}}{dt_{1}\,du_{1}}=F^{B}_{q{%
\bar{q}}\to AB}\,\delta(s_{4})\,, divide start_ARG italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over^ start_ARG italic_σ end_ARG start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT end_ARG start_ARG italic_d italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_d italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG = italic_F start_POSTSUPERSCRIPT italic_B end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT italic_δ ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ) ,
(2.19)
where the specific form of F q q ¯ → A B B subscript superscript 𝐹 𝐵 → 𝑞 ¯ 𝑞 𝐴 𝐵 F^{B}_{q{\bar{q}}\to AB} italic_F start_POSTSUPERSCRIPT italic_B end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT depends on the process (see the next section for H + H − superscript 𝐻 superscript 𝐻 H^{+}H^{-} italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT production). If we perturbatively expand H q q ¯ → A B = H q q ¯ → A B ( 0 ) + ( α s / π ) H q q ¯ → A B ( 1 ) + ( α s / π ) 2 H q q ¯ → A B ( 2 ) + ⋯ subscript 𝐻 → 𝑞 ¯ 𝑞 𝐴 𝐵 superscript subscript 𝐻 → 𝑞 ¯ 𝑞 𝐴 𝐵 0 subscript 𝛼 𝑠 𝜋 superscript subscript 𝐻 → 𝑞 ¯ 𝑞 𝐴 𝐵 1 superscript subscript 𝛼 𝑠 𝜋 2 superscript subscript 𝐻 → 𝑞 ¯ 𝑞 𝐴 𝐵 2 ⋯ H_{q{\bar{q}}\to AB}=H_{q{\bar{q}}\to AB}^{(0)}+(\alpha_{s}/\pi)H_{q{\bar{q}}%
\to AB}^{(1)}+(\alpha_{s}/\pi)^{2}H_{q{\bar{q}}\to AB}^{(2)}+\cdots italic_H start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT = italic_H start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT + ( italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT / italic_π ) italic_H start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT + ( italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT / italic_π ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_H start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT + ⋯ , and S q q ¯ → A B = S q q ¯ → A B ( 0 ) + ( α s / π ) S q q ¯ → A B ( 1 ) + ( α s / π ) 2 S q q ¯ → A B ( 2 ) + ⋯ subscript 𝑆 → 𝑞 ¯ 𝑞 𝐴 𝐵 superscript subscript 𝑆 → 𝑞 ¯ 𝑞 𝐴 𝐵 0 subscript 𝛼 𝑠 𝜋 superscript subscript 𝑆 → 𝑞 ¯ 𝑞 𝐴 𝐵 1 superscript subscript 𝛼 𝑠 𝜋 2 superscript subscript 𝑆 → 𝑞 ¯ 𝑞 𝐴 𝐵 2 ⋯ S_{q{\bar{q}}\to AB}=S_{q{\bar{q}}\to AB}^{(0)}+(\alpha_{s}/\pi)S_{q{\bar{q}}%
\to AB}^{(1)}+(\alpha_{s}/\pi)^{2}S_{q{\bar{q}}\to AB}^{(2)}+\cdots italic_S start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT = italic_S start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT + ( italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT / italic_π ) italic_S start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT + ( italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT / italic_π ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_S start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT + ⋯ , then we have F q q ¯ → A B B = H q q ¯ → A B ( 0 ) S q q ¯ → A B ( 0 ) subscript superscript 𝐹 𝐵 → 𝑞 ¯ 𝑞 𝐴 𝐵 superscript subscript 𝐻 → 𝑞 ¯ 𝑞 𝐴 𝐵 0 superscript subscript 𝑆 → 𝑞 ¯ 𝑞 𝐴 𝐵 0 F^{B}_{q{\bar{q}}\to AB}=H_{q{\bar{q}}\to AB}^{(0)}\,S_{q{\bar{q}}\to AB}^{(0)} italic_F start_POSTSUPERSCRIPT italic_B end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT = italic_H start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT italic_S start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT .
The NLO soft-plus-virtual corrections are given by
d 2 σ ^ q q ¯ → A B ( 1 ) d t 1 d u 1 = F q q ¯ → A B B α s ( μ R ) π { c 3 𝒟 1 ( s 4 ) + c 2 𝒟 0 ( s 4 ) + c 1 δ ( s 4 ) } superscript 𝑑 2 subscript superscript ^ 𝜎 1 → 𝑞 ¯ 𝑞 𝐴 𝐵 𝑑 subscript 𝑡 1 𝑑 subscript 𝑢 1 subscript superscript 𝐹 𝐵 → 𝑞 ¯ 𝑞 𝐴 𝐵 subscript 𝛼 𝑠 subscript 𝜇 𝑅 𝜋 subscript 𝑐 3 subscript 𝒟 1 subscript 𝑠 4 subscript 𝑐 2 subscript 𝒟 0 subscript 𝑠 4 subscript 𝑐 1 𝛿 subscript 𝑠 4 \frac{d^{2}{\hat{\sigma}}^{(1)}_{q{\bar{q}}\to AB}}{dt_{1}\,du_{1}}=F^{B}_{q{%
\bar{q}}\to AB}\frac{\alpha_{s}(\mu_{R})}{\pi}\left\{c_{3}\,{\cal D}_{1}(s_{4}%
)+c_{2}\,{\cal D}_{0}(s_{4})+c_{1}\,\delta(s_{4})\right\} divide start_ARG italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over^ start_ARG italic_σ end_ARG start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT end_ARG start_ARG italic_d italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_d italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG = italic_F start_POSTSUPERSCRIPT italic_B end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT divide start_ARG italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ) end_ARG start_ARG italic_π end_ARG { italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT caligraphic_D start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ) + italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT caligraphic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ) + italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_δ ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ) }
(2.20)
where μ R subscript 𝜇 𝑅 \mu_{R} italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT is the renormalization scale and
𝒟 k ( s 4 ) = [ ln k ( s 4 / m 2 ) s 4 ] + . subscript 𝒟 𝑘 subscript 𝑠 4 subscript delimited-[] superscript 𝑘 subscript 𝑠 4 superscript 𝑚 2 subscript 𝑠 4 {\cal D}_{k}(s_{4})=\left[\frac{\ln^{k}(s_{4}/m^{2})}{s_{4}}\right]_{+}\,. caligraphic_D start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ) = [ divide start_ARG roman_ln start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT / italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT end_ARG ] start_POSTSUBSCRIPT + end_POSTSUBSCRIPT .
(2.21)
The coefficients of the 𝒟 1 ( s 4 ) subscript 𝒟 1 subscript 𝑠 4 {\cal D}_{1}(s_{4}) caligraphic_D start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ) and 𝒟 0 ( s 4 ) subscript 𝒟 0 subscript 𝑠 4 {\cal D}_{0}(s_{4}) caligraphic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ) terms in Eq. (2.20 ) are, respectively,
c 3 = 4 C F and c 2 = − 2 C F ln ( t 1 u 1 m 4 ) − 2 C F ln ( μ F 2 s ) . formulae-sequence subscript 𝑐 3 4 subscript 𝐶 𝐹 and
subscript 𝑐 2 2 subscript 𝐶 𝐹 subscript 𝑡 1 subscript 𝑢 1 superscript 𝑚 4 2 subscript 𝐶 𝐹 superscript subscript 𝜇 𝐹 2 𝑠 c_{3}=4C_{F}\qquad{\rm and}\qquad c_{2}=-2\,C_{F}\,\ln\left(\frac{t_{1}u_{1}}{%
m^{4}}\right)-2\,C_{F}\,\ln\left(\frac{\mu_{F}^{2}}{s}\right)\,. italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT = 4 italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_and italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = - 2 italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln ( divide start_ARG italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG ) - 2 italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) .
(2.22)
The coefficient of the δ ( s 4 ) 𝛿 subscript 𝑠 4 \delta(s_{4}) italic_δ ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ) term is
c 1 = C F ln 2 ( − t 1 m 2 ) + C F ln 2 ( − u 1 m 2 ) + [ C F ln ( t 1 u 1 m 4 ) − 2 γ q ( 1 ) ] ln ( μ F 2 s ) + V 1 , subscript 𝑐 1 subscript 𝐶 𝐹 superscript 2 subscript 𝑡 1 superscript 𝑚 2 subscript 𝐶 𝐹 superscript 2 subscript 𝑢 1 superscript 𝑚 2 delimited-[] subscript 𝐶 𝐹 subscript 𝑡 1 subscript 𝑢 1 superscript 𝑚 4 2 superscript subscript 𝛾 𝑞 1 superscript subscript 𝜇 𝐹 2 𝑠 subscript 𝑉 1 c_{1}=C_{F}\,\ln^{2}\left(\frac{-t_{1}}{m^{2}}\right)+C_{F}\,\ln^{2}\left(%
\frac{-u_{1}}{m^{2}}\right)+\left[C_{F}\,\ln\left(\frac{t_{1}u_{1}}{m^{4}}%
\right)-2\gamma_{q}^{(1)}\right]\ln\left(\frac{\mu_{F}^{2}}{s}\right)+V_{1}\,, italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG - italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) + italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG - italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) + [ italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln ( divide start_ARG italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG ) - 2 italic_γ start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT ] roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) + italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ,
(2.23)
where V 1 = ( H q q ¯ → A B ( 1 ) S q q ¯ → A B ( 0 ) + H q q ¯ → A B ( 0 ) S q q ¯ → A B ( 1 ) ) / F q q ¯ → A B B subscript 𝑉 1 superscript subscript 𝐻 → 𝑞 ¯ 𝑞 𝐴 𝐵 1 superscript subscript 𝑆 → 𝑞 ¯ 𝑞 𝐴 𝐵 0 superscript subscript 𝐻 → 𝑞 ¯ 𝑞 𝐴 𝐵 0 superscript subscript 𝑆 → 𝑞 ¯ 𝑞 𝐴 𝐵 1 subscript superscript 𝐹 𝐵 → 𝑞 ¯ 𝑞 𝐴 𝐵 V_{1}=\left(H_{q{\bar{q}}\to AB}^{(1)}\,S_{q{\bar{q}}\to AB}^{(0)}+H_{q{\bar{q%
}}\to AB}^{(0)}\,S_{q{\bar{q}}\to AB}^{(1)}\right)/F^{B}_{q{\bar{q}}\to AB} italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = ( italic_H start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT italic_S start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT + italic_H start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT italic_S start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT ) / italic_F start_POSTSUPERSCRIPT italic_B end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT is the contribution from the one-loop virtual corrections. For Drell-Yan type processes with quark-antiquark annihilation, such as H + H − superscript 𝐻 superscript 𝐻 H^{+}H^{-} italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT production, V 1 subscript 𝑉 1 V_{1} italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT is given by
the expression V 1 = 2 C F ( − 2 + ζ 2 ) subscript 𝑉 1 2 subscript 𝐶 𝐹 2 subscript 𝜁 2 V_{1}=2C_{F}(-2+\zeta_{2}) italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 2 italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ( - 2 + italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) [40 , 41 , 42 ] . For gluon-initiated processes, we replace C F subscript 𝐶 𝐹 C_{F} italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT by C A subscript 𝐶 𝐴 C_{A} italic_C start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT in Eqs. (2.22 ) and (2.23 ), and we use the corresponding V 1 subscript 𝑉 1 V_{1} italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT (for example, see [43 , 44 ] for the NLO virtual corrections for g g → H → 𝑔 𝑔 𝐻 gg\to H italic_g italic_g → italic_H ).
The NNLO soft-plus-virtual corrections are given by
d 2 σ ^ q q ¯ → A B ( 2 ) d t 1 d u 1 superscript 𝑑 2 subscript superscript ^ 𝜎 2 → 𝑞 ¯ 𝑞 𝐴 𝐵 𝑑 subscript 𝑡 1 𝑑 subscript 𝑢 1 \displaystyle\frac{d^{2}{\hat{\sigma}}^{(2)}_{q{\bar{q}}\to AB}}{dt_{1}\ du_{1}} divide start_ARG italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over^ start_ARG italic_σ end_ARG start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT end_ARG start_ARG italic_d italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_d italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG
= \displaystyle= =
F q q ¯ → A B B α s 2 ( μ R ) π 2 { 1 2 c 3 2 𝒟 3 ( s 4 ) + [ 3 2 c 3 c 2 − β 0 4 c 3 ] 𝒟 2 ( s 4 ) \displaystyle F^{B}_{q{\bar{q}}\to AB}\frac{\alpha_{s}^{2}(\mu_{R})}{\pi^{2}}%
\left\{\frac{1}{2}c_{3}^{2}\,{\cal D}_{3}(s_{4})+\left[\frac{3}{2}c_{3}c_{2}-%
\frac{\beta_{0}}{4}c_{3}\right]{\cal D}_{2}(s_{4})\right. italic_F start_POSTSUPERSCRIPT italic_B end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT divide start_ARG italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ) end_ARG start_ARG italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG { divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT caligraphic_D start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ) + [ divide start_ARG 3 end_ARG start_ARG 2 end_ARG italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ] caligraphic_D start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT )
+ [ c 3 c 1 + c 2 2 − ζ 2 c 3 2 − β 0 2 c 2 − β 0 C F ln ( μ F 2 μ R 2 ) + 4 A q ( 2 ) ] 𝒟 1 ( s 4 ) delimited-[] subscript 𝑐 3 subscript 𝑐 1 superscript subscript 𝑐 2 2 subscript 𝜁 2 superscript subscript 𝑐 3 2 subscript 𝛽 0 2 subscript 𝑐 2 subscript 𝛽 0 subscript 𝐶 𝐹 superscript subscript 𝜇 𝐹 2 superscript subscript 𝜇 𝑅 2 4 superscript subscript 𝐴 𝑞 2 subscript 𝒟 1 subscript 𝑠 4 \displaystyle{}+\left[c_{3}c_{1}+c_{2}^{2}-\zeta_{2}c_{3}^{2}-\frac{\beta_{0}}%
{2}c_{2}-\beta_{0}C_{F}\ln\left(\frac{\mu_{F}^{2}}{\mu_{R}^{2}}\right)+4A_{q}^%
{(2)}\right]{\cal D}_{1}(s_{4}) + [ italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) + 4 italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT ] caligraphic_D start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT )
+ [ c 2 c 1 − ζ 2 c 3 c 2 + ζ 3 c 3 2 + β 0 4 c 2 ln ( μ R 2 s ) \displaystyle{}+\left[c_{2}c_{1}-\zeta_{2}c_{3}c_{2}+\zeta_{3}c_{3}^{2}+\frac{%
\beta_{0}}{4}c_{2}\ln\left(\frac{\mu_{R}^{2}}{s}\right)\right. + [ italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG )
− β 0 2 C F ln 2 ( − t 1 m 2 ) − β 0 2 C F ln 2 ( − u 1 m 2 ) − 2 A q ( 2 ) ln ( t 1 u 1 m 4 ) + 2 D q ( 2 ) subscript 𝛽 0 2 subscript 𝐶 𝐹 superscript 2 subscript 𝑡 1 superscript 𝑚 2 subscript 𝛽 0 2 subscript 𝐶 𝐹 superscript 2 subscript 𝑢 1 superscript 𝑚 2 2 superscript subscript 𝐴 𝑞 2 subscript 𝑡 1 subscript 𝑢 1 superscript 𝑚 4 2 superscript subscript 𝐷 𝑞 2 \displaystyle\;{}-\frac{\beta_{0}}{2}C_{F}\ln^{2}\left(\frac{-t_{1}}{m^{2}}%
\right)-\frac{\beta_{0}}{2}C_{F}\ln^{2}\left(\frac{-u_{1}}{m^{2}}\right)-2A_{q%
}^{(2)}\ln\left(\frac{t_{1}u_{1}}{m^{4}}\right)+2D_{q}^{(2)} - divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG - italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) - divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG - italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) - 2 italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT roman_ln ( divide start_ARG italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG ) + 2 italic_D start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT
+ β 0 4 C F ln 2 ( μ F 2 s ) − 2 A q ( 2 ) ln ( μ F 2 s ) ] 𝒟 0 ( s 4 ) \displaystyle\;\left.{}+\frac{\beta_{0}}{4}C_{F}\ln^{2}\left(\frac{\mu_{F}^{2}%
}{s}\right)-2A_{q}^{(2)}\ln\left(\frac{\mu_{F}^{2}}{s}\right)\right]{\cal D}_{%
0}(s_{4}) + divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) - 2 italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) ] caligraphic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT )
+ [ V 2 + 1 2 ( c 1 2 − V 1 2 ) − ζ 2 2 c 2 2 + ζ 3 c 3 c 2 + β 0 6 C F ( ln 3 ( − t 1 m 2 ) + ln 3 ( − u 1 m 2 ) ) \displaystyle{}+\left[V_{2}+\frac{1}{2}\left(c_{1}^{2}-V_{1}^{2}\right)-\frac{%
\zeta_{2}}{2}c_{2}^{2}+\zeta_{3}c_{3}c_{2}+\frac{\beta_{0}}{6}C_{F}\left(\ln^{%
3}\left(\frac{-t_{1}}{m^{2}}\right)+\ln^{3}\left(\frac{-u_{1}}{m^{2}}\right)%
\right)\right. + [ italic_V start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) - divide start_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 6 end_ARG italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ( roman_ln start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ( divide start_ARG - italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) + roman_ln start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ( divide start_ARG - italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) )
+ ( β 0 4 C F + A q ( 2 ) ) ( ln 2 ( − t 1 m 2 ) + ln 2 ( − u 1 m 2 ) ) + β 0 4 c 1 ln ( μ R 2 s ) − 2 γ q ( 2 ) ln ( μ F 2 s ) subscript 𝛽 0 4 subscript 𝐶 𝐹 superscript subscript 𝐴 𝑞 2 superscript 2 subscript 𝑡 1 superscript 𝑚 2 superscript 2 subscript 𝑢 1 superscript 𝑚 2 subscript 𝛽 0 4 subscript 𝑐 1 superscript subscript 𝜇 𝑅 2 𝑠 2 superscript subscript 𝛾 𝑞 2 superscript subscript 𝜇 𝐹 2 𝑠 \displaystyle\;{}+\left(\frac{\beta_{0}}{4}C_{F}+A_{q}^{(2)}\right)\left(\ln^{%
2}\left(\frac{-t_{1}}{m^{2}}\right)+\ln^{2}\left(\frac{-u_{1}}{m^{2}}\right)%
\right)+\frac{\beta_{0}}{4}c_{1}\ln\left(\frac{\mu_{R}^{2}}{s}\right)-2\gamma_%
{q}^{(2)}\ln\left(\frac{\mu_{F}^{2}}{s}\right) + ( divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT + italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT ) ( roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG - italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) + roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG - italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) ) + divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) - 2 italic_γ start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG )
+ A q ( 2 ) ln ( t 1 u 1 m 4 ) ln ( μ F 2 s ) + β 0 8 ( 2 γ q ( 1 ) − C F ln ( t 1 u 1 m 4 ) ) ln 2 ( μ F 2 s ) ] δ ( s 4 ) } \displaystyle\;\left.\left.{}+A_{q}^{(2)}\ln\left(\frac{t_{1}u_{1}}{m^{4}}%
\right)\ln\left(\frac{\mu_{F}^{2}}{s}\right)+\frac{\beta_{0}}{8}\left(2\gamma_%
{q}^{(1)}-C_{F}\ln\left(\frac{t_{1}u_{1}}{m^{4}}\right)\right)\ln^{2}\left(%
\frac{\mu_{F}^{2}}{s}\right)\right]\delta(s_{4})\right\} + italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT roman_ln ( divide start_ARG italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG ) roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) + divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 8 end_ARG ( 2 italic_γ start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT - italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln ( divide start_ARG italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG ) ) roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) ] italic_δ ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ) }
where V 2 = ( H q q ¯ → A B ( 2 ) S q q ¯ → A B ( 0 ) + H q q ¯ → A B ( 0 ) S q q ¯ → A B ( 2 ) + H q q ¯ → A B ( 1 ) S q q ¯ → A B ( 1 ) ) / F q q ¯ → A B B subscript 𝑉 2 superscript subscript 𝐻 → 𝑞 ¯ 𝑞 𝐴 𝐵 2 superscript subscript 𝑆 → 𝑞 ¯ 𝑞 𝐴 𝐵 0 superscript subscript 𝐻 → 𝑞 ¯ 𝑞 𝐴 𝐵 0 superscript subscript 𝑆 → 𝑞 ¯ 𝑞 𝐴 𝐵 2 superscript subscript 𝐻 → 𝑞 ¯ 𝑞 𝐴 𝐵 1 superscript subscript 𝑆 → 𝑞 ¯ 𝑞 𝐴 𝐵 1 subscript superscript 𝐹 𝐵 → 𝑞 ¯ 𝑞 𝐴 𝐵 V_{2}=\left(H_{q{\bar{q}}\to AB}^{(2)}\,S_{q{\bar{q}}\to AB}^{(0)}+H_{q{\bar{q%
}}\to AB}^{(0)}\,S_{q{\bar{q}}\to AB}^{(2)}+H_{q{\bar{q}}\to AB}^{(1)}\,S_{q{%
\bar{q}}\to AB}^{(1)}\right)/F^{B}_{q{\bar{q}}\to AB} italic_V start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = ( italic_H start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT italic_S start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT + italic_H start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT italic_S start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT + italic_H start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT italic_S start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT ) / italic_F start_POSTSUPERSCRIPT italic_B end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT is the contribution from the virtual two-loop corrections.
For Drell-Yan type processes with quark-antiquark annihilation, such as H + H − superscript 𝐻 superscript 𝐻 H^{+}H^{-} italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT production, V 2 subscript 𝑉 2 V_{2} italic_V start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT is given by the expression [45 ]
V 2 subscript 𝑉 2 \displaystyle V_{2} italic_V start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT
= \displaystyle= =
C F 2 ( 511 64 − 35 8 ζ 2 − 15 4 ζ 3 + ζ 2 2 10 ) + C F C A ( − 1535 192 + 37 9 ζ 2 + 7 4 ζ 3 − 3 20 ζ 2 2 ) superscript subscript 𝐶 𝐹 2 511 64 35 8 subscript 𝜁 2 15 4 subscript 𝜁 3 superscript subscript 𝜁 2 2 10 subscript 𝐶 𝐹 subscript 𝐶 𝐴 1535 192 37 9 subscript 𝜁 2 7 4 subscript 𝜁 3 3 20 superscript subscript 𝜁 2 2 \displaystyle C_{F}^{2}\left(\frac{511}{64}-\frac{35}{8}\zeta_{2}-\frac{15}{4}%
\zeta_{3}+\frac{\zeta_{2}^{2}}{10}\right)+C_{F}C_{A}\left(-\frac{1535}{192}+%
\frac{37}{9}\zeta_{2}+\frac{7}{4}\zeta_{3}-\frac{3}{20}\zeta_{2}^{2}\right) italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG 511 end_ARG start_ARG 64 end_ARG - divide start_ARG 35 end_ARG start_ARG 8 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - divide start_ARG 15 end_ARG start_ARG 4 end_ARG italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT + divide start_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 10 end_ARG ) + italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ( - divide start_ARG 1535 end_ARG start_ARG 192 end_ARG + divide start_ARG 37 end_ARG start_ARG 9 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + divide start_ARG 7 end_ARG start_ARG 4 end_ARG italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT - divide start_ARG 3 end_ARG start_ARG 20 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT )
(2.25)
+ C F n f ( 127 96 − 7 9 ζ 2 + ζ 3 2 ) . subscript 𝐶 𝐹 subscript 𝑛 𝑓 127 96 7 9 subscript 𝜁 2 subscript 𝜁 3 2 \displaystyle{}+C_{F}n_{f}\left(\frac{127}{96}-\frac{7}{9}\zeta_{2}+\frac{%
\zeta_{3}}{2}\right)\,. + italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT ( divide start_ARG 127 end_ARG start_ARG 96 end_ARG - divide start_ARG 7 end_ARG start_ARG 9 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + divide start_ARG italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG ) .
For processes g g → A B → 𝑔 𝑔 𝐴 𝐵 gg\to AB italic_g italic_g → italic_A italic_B , we replace in Eq. (LABEL:NNLO ) all explicit appearances of C F subscript 𝐶 𝐹 C_{F} italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT by C A subscript 𝐶 𝐴 C_{A} italic_C start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT , and we also use the corresponding expressions for A g subscript 𝐴 𝑔 A_{g} italic_A start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , D g subscript 𝐷 𝑔 D_{g} italic_D start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , γ g subscript 𝛾 𝑔 \gamma_{g} italic_γ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , c 3 subscript 𝑐 3 c_{3} italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT , c 2 subscript 𝑐 2 c_{2} italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , and c 1 subscript 𝑐 1 c_{1} italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT as detailed above, as well as the appropriate expressions for V 1 subscript 𝑉 1 V_{1} italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and V 2 subscript 𝑉 2 V_{2} italic_V start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT (for example, see [46 ] for the NNLO virtual corrections for g g → H → 𝑔 𝑔 𝐻 gg\to H italic_g italic_g → italic_H ).
The N3 LO soft-gluon corrections are given by
d 2 σ ^ q q ¯ → A B ( 3 ) d t 1 d u 1 superscript 𝑑 2 subscript superscript ^ 𝜎 3 → 𝑞 ¯ 𝑞 𝐴 𝐵 𝑑 subscript 𝑡 1 𝑑 subscript 𝑢 1 \displaystyle\frac{d^{2}{\hat{\sigma}}^{(3)}_{q{\bar{q}}\to AB}}{dt_{1}\,du_{1}} divide start_ARG italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over^ start_ARG italic_σ end_ARG start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT end_ARG start_ARG italic_d italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_d italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG
= \displaystyle= =
F q q ¯ → A B B α s 3 ( μ R ) π 3 { 1 8 c 3 3 𝒟 5 ( s 4 ) + [ 5 8 c 3 2 c 2 − 5 24 c 3 2 β 0 ] 𝒟 4 ( s 4 ) \displaystyle F^{B}_{q{\bar{q}}\to AB}\frac{\alpha_{s}^{3}(\mu_{R})}{\pi^{3}}%
\left\{\frac{1}{8}c_{3}^{3}\;{\cal D}_{5}(s_{4})+\left[\frac{5}{8}c_{3}^{2}c_{%
2}-\frac{5}{24}c_{3}^{2}\beta_{0}\right]\;{\cal D}_{4}(s_{4})\right. italic_F start_POSTSUPERSCRIPT italic_B end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT divide start_ARG italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ( italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT ) end_ARG start_ARG italic_π start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG { divide start_ARG 1 end_ARG start_ARG 8 end_ARG italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT caligraphic_D start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ) + [ divide start_ARG 5 end_ARG start_ARG 8 end_ARG italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - divide start_ARG 5 end_ARG start_ARG 24 end_ARG italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ] caligraphic_D start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT )
(2.26)
+ [ c 3 c 2 2 + 1 2 c 3 2 c 1 − ζ 2 c 3 3 + β 0 2 12 c 3 − 5 6 β 0 c 3 c 2 − β 0 C F c 3 ln ( μ F 2 μ R 2 ) + 4 c 3 A q ( 2 ) ] 𝒟 3 ( s 4 ) delimited-[] subscript 𝑐 3 superscript subscript 𝑐 2 2 1 2 superscript subscript 𝑐 3 2 subscript 𝑐 1 subscript 𝜁 2 superscript subscript 𝑐 3 3 superscript subscript 𝛽 0 2 12 subscript 𝑐 3 5 6 subscript 𝛽 0 subscript 𝑐 3 subscript 𝑐 2 subscript 𝛽 0 subscript 𝐶 𝐹 subscript 𝑐 3 superscript subscript 𝜇 𝐹 2 superscript subscript 𝜇 𝑅 2 4 subscript 𝑐 3 superscript subscript 𝐴 𝑞 2 subscript 𝒟 3 subscript 𝑠 4 \displaystyle{}+\left[c_{3}c_{2}^{2}+\frac{1}{2}c_{3}^{2}c_{1}-\zeta_{2}c_{3}^%
{3}+\frac{\beta_{0}^{2}}{12}c_{3}-\frac{5}{6}\beta_{0}c_{3}c_{2}-\beta_{0}C_{F%
}c_{3}\ln\left(\frac{\mu_{F}^{2}}{\mu_{R}^{2}}\right)+4c_{3}A_{q}^{(2)}\right]%
\;{\cal D}_{3}(s_{4}) + [ italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 12 end_ARG italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT - divide start_ARG 5 end_ARG start_ARG 6 end_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) + 4 italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT ] caligraphic_D start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT )
+ [ 3 2 c 3 c 2 c 1 + 1 2 c 2 3 − 3 ζ 2 c 3 2 c 2 + 5 2 ζ 3 c 3 3 − β 0 4 c 3 c 1 + 9 8 β 0 ζ 2 c 3 2 \displaystyle{}+\left[\frac{3}{2}c_{3}c_{2}c_{1}+\frac{1}{2}c_{2}^{3}-3\zeta_{%
2}c_{3}^{2}c_{2}+\frac{5}{2}\zeta_{3}c_{3}^{3}-\frac{\beta_{0}}{4}c_{3}c_{1}+%
\frac{9}{8}\beta_{0}\zeta_{2}c_{3}^{2}\right. + [ divide start_ARG 3 end_ARG start_ARG 2 end_ARG italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT - 3 italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + divide start_ARG 5 end_ARG start_ARG 2 end_ARG italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT - divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + divide start_ARG 9 end_ARG start_ARG 8 end_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT
+ ( 3 c 2 − β 0 ) ( − β 0 4 c 2 − β 0 2 C F ln ( μ F 2 μ R 2 ) + 2 A q ( 2 ) ) − C F β 1 4 − 3 2 c 3 X 1 ] 𝒟 2 ( s 4 ) \displaystyle\;\left.{}+\left(3c_{2}-\beta_{0}\right)\left(-\frac{\beta_{0}}{4%
}c_{2}-\frac{\beta_{0}}{2}C_{F}\ln\left(\frac{\mu_{F}^{2}}{\mu_{R}^{2}}\right)%
+2A_{q}^{(2)}\right)-C_{F}\frac{\beta_{1}}{4}-\frac{3}{2}c_{3}X_{1}\right]{%
\cal D}_{2}(s_{4}) + ( 3 italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ( - divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) + 2 italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT ) - italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT divide start_ARG italic_β start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG - divide start_ARG 3 end_ARG start_ARG 2 end_ARG italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_X start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ] caligraphic_D start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT )
+ [ 1 2 c 3 c 1 2 + c 2 2 c 1 − ζ 2 c 3 2 c 1 − 5 2 ζ 2 c 3 c 2 2 + 5 ζ 3 c 3 2 c 2 + 5 4 ζ 2 2 c 3 3 − 15 4 ζ 4 c 3 3 \displaystyle{}+\left[\frac{1}{2}c_{3}c_{1}^{2}+c_{2}^{2}c_{1}-\zeta_{2}c_{3}^%
{2}c_{1}-\frac{5}{2}\zeta_{2}c_{3}c_{2}^{2}+5\zeta_{3}c_{3}^{2}c_{2}+\frac{5}{%
4}\zeta_{2}^{2}c_{3}^{3}-\frac{15}{4}\zeta_{4}c_{3}^{3}\right. + [ divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - divide start_ARG 5 end_ARG start_ARG 2 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 5 italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + divide start_ARG 5 end_ARG start_ARG 4 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT - divide start_ARG 15 end_ARG start_ARG 4 end_ARG italic_ζ start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT
− β 0 2 4 ζ 2 c 3 − 5 3 β 0 ζ 3 c 3 2 + β 0 ζ 2 c 3 c 2 + ( 2 c 1 − 5 ζ 2 c 3 ) ( − β 0 4 c 2 − β 0 2 C F ln ( μ F 2 μ R 2 ) + 2 A q ( 2 ) ) superscript subscript 𝛽 0 2 4 subscript 𝜁 2 subscript 𝑐 3 5 3 subscript 𝛽 0 subscript 𝜁 3 superscript subscript 𝑐 3 2 subscript 𝛽 0 subscript 𝜁 2 subscript 𝑐 3 subscript 𝑐 2 2 subscript 𝑐 1 5 subscript 𝜁 2 subscript 𝑐 3 subscript 𝛽 0 4 subscript 𝑐 2 subscript 𝛽 0 2 subscript 𝐶 𝐹 superscript subscript 𝜇 𝐹 2 superscript subscript 𝜇 𝑅 2 2 superscript subscript 𝐴 𝑞 2 \displaystyle\;{}-\frac{\beta_{0}^{2}}{4}\zeta_{2}c_{3}-\frac{5}{3}\beta_{0}%
\zeta_{3}c_{3}^{2}+\beta_{0}\zeta_{2}c_{3}c_{2}+\left(2c_{1}-5\zeta_{2}c_{3}%
\right)\left(-\frac{\beta_{0}}{4}c_{2}-\frac{\beta_{0}}{2}C_{F}\ln\left(\frac{%
\mu_{F}^{2}}{\mu_{R}^{2}}\right)+2A_{q}^{(2)}\right) - divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 4 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT - divide start_ARG 5 end_ARG start_ARG 3 end_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + ( 2 italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - 5 italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ) ( - divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) + 2 italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT )
+ ( β 0 − 2 c 2 ) X 1 + c 3 X 0 + 4 A q ( 3 ) + C F β 0 2 4 ln 2 ( μ F 2 μ R 2 ) − 2 β 0 A q ( 2 ) ln ( μ F 2 μ R 2 ) subscript 𝛽 0 2 subscript 𝑐 2 subscript 𝑋 1 subscript 𝑐 3 subscript 𝑋 0 4 superscript subscript 𝐴 𝑞 3 subscript 𝐶 𝐹 superscript subscript 𝛽 0 2 4 superscript 2 superscript subscript 𝜇 𝐹 2 superscript subscript 𝜇 𝑅 2 2 subscript 𝛽 0 superscript subscript 𝐴 𝑞 2 superscript subscript 𝜇 𝐹 2 superscript subscript 𝜇 𝑅 2 \displaystyle\;{}+(\beta_{0}-2c_{2})X_{1}+c_{3}X_{0}+4A_{q}^{(3)}+C_{F}\frac{%
\beta_{0}^{2}}{4}\ln^{2}\left(\frac{\mu_{F}^{2}}{\mu_{R}^{2}}\right)-2\beta_{0%
}A_{q}^{(2)}\ln\left(\frac{\mu_{F}^{2}}{\mu_{R}^{2}}\right) + ( italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - 2 italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) italic_X start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_X start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + 4 italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT + italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 4 end_ARG roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) - 2 italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG )
+ C F β 1 4 ln ( μ R 2 s ) + C F β 1 4 ln ( t 1 u 1 m 4 ) ] 𝒟 1 ( s 4 ) \displaystyle\;\left.{}+C_{F}\frac{\beta_{1}}{4}\ln\left(\frac{\mu_{R}^{2}}{s}%
\right)+C_{F}\frac{\beta_{1}}{4}\ln\left(\frac{t_{1}u_{1}}{m^{4}}\right)\right%
]\;{\cal D}_{1}(s_{4}) + italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT divide start_ARG italic_β start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) + italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT divide start_ARG italic_β start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG roman_ln ( divide start_ARG italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG ) ] caligraphic_D start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT )
+ [ 1 2 c 2 c 1 2 + 3 ζ 5 c 3 3 − 15 4 ζ 4 c 3 2 c 2 − 2 ζ 2 ζ 3 c 3 3 + ζ 3 c 3 2 c 1 + 2 ζ 3 c 3 c 2 2 + 5 4 ζ 2 2 c 3 2 c 2 − ζ 2 c 3 c 2 c 1 − ζ 2 2 c 2 3 \displaystyle{}+\left[\frac{1}{2}c_{2}c_{1}^{2}+3\zeta_{5}c_{3}^{3}-\frac{15}{%
4}\zeta_{4}c_{3}^{2}c_{2}-2\zeta_{2}\zeta_{3}c_{3}^{3}+\zeta_{3}c_{3}^{2}c_{1}%
+2\zeta_{3}c_{3}c_{2}^{2}+\frac{5}{4}\zeta_{2}^{2}c_{3}^{2}c_{2}-\zeta_{2}c_{3%
}c_{2}c_{1}-\frac{\zeta_{2}}{2}c_{2}^{3}\right. + [ divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 3 italic_ζ start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT - divide start_ARG 15 end_ARG start_ARG 4 end_ARG italic_ζ start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - 2 italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT + italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 2 italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG 5 end_ARG start_ARG 4 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - divide start_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT
+ β 0 12 c 3 ( 15 ζ 4 c 3 − 8 ζ 3 c 2 − 6 ζ 2 2 c 3 + 3 ζ 2 c 1 ) subscript 𝛽 0 12 subscript 𝑐 3 15 subscript 𝜁 4 subscript 𝑐 3 8 subscript 𝜁 3 subscript 𝑐 2 6 superscript subscript 𝜁 2 2 subscript 𝑐 3 3 subscript 𝜁 2 subscript 𝑐 1 \displaystyle\;{}+\frac{\beta_{0}}{12}c_{3}\left(15\zeta_{4}c_{3}-8\zeta_{3}c_%
{2}-6\zeta_{2}^{2}c_{3}+3\zeta_{2}c_{1}\right) + divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 12 end_ARG italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ( 15 italic_ζ start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT - 8 italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - 6 italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT + 3 italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT )
+ ( 4 ζ 3 c 3 − 3 ζ 2 c 2 ) ( − β 0 4 c 2 − β 0 2 C F ln ( μ F 2 μ R 2 ) + 2 A q ( 2 ) ) 4 subscript 𝜁 3 subscript 𝑐 3 3 subscript 𝜁 2 subscript 𝑐 2 subscript 𝛽 0 4 subscript 𝑐 2 subscript 𝛽 0 2 subscript 𝐶 𝐹 superscript subscript 𝜇 𝐹 2 superscript subscript 𝜇 𝑅 2 2 superscript subscript 𝐴 𝑞 2 \displaystyle\;{}+\left(4\zeta_{3}c_{3}-3\zeta_{2}c_{2}\right)\left(-\frac{%
\beta_{0}}{4}c_{2}-\frac{\beta_{0}}{2}C_{F}\ln\left(\frac{\mu_{F}^{2}}{\mu_{R}%
^{2}}\right)+2A_{q}^{(2)}\right) + ( 4 italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT - 3 italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) ( - divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) + 2 italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT )
+ ( ζ 2 c 3 − c 1 ) X 1 + c 2 X 0 − β 0 2 4 C F ( ln 2 ( − t 1 m 2 ) + ln 2 ( − u 1 m 2 ) ) ln ( μ R 2 s ) subscript 𝜁 2 subscript 𝑐 3 subscript 𝑐 1 subscript 𝑋 1 subscript 𝑐 2 subscript 𝑋 0 superscript subscript 𝛽 0 2 4 subscript 𝐶 𝐹 superscript 2 subscript 𝑡 1 superscript 𝑚 2 superscript 2 subscript 𝑢 1 superscript 𝑚 2 superscript subscript 𝜇 𝑅 2 𝑠 \displaystyle\;{}+\left(\zeta_{2}c_{3}-c_{1}\right)X_{1}+c_{2}X_{0}-\frac{%
\beta_{0}^{2}}{4}C_{F}\left(\ln^{2}\left(\frac{-t_{1}}{m^{2}}\right)+\ln^{2}%
\left(\frac{-u_{1}}{m^{2}}\right)\right)\ln\left(\frac{\mu_{R}^{2}}{s}\right) + ( italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT - italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) italic_X start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_X start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 4 end_ARG italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ( roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG - italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) + roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG - italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) ) roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG )
+ β 0 2 16 c 2 ln 2 ( μ R 2 s ) + β 0 2 8 C F ln ( μ F 2 s ) ln 2 ( μ R 2 s ) superscript subscript 𝛽 0 2 16 subscript 𝑐 2 superscript 2 superscript subscript 𝜇 𝑅 2 𝑠 superscript subscript 𝛽 0 2 8 subscript 𝐶 𝐹 superscript subscript 𝜇 𝐹 2 𝑠 superscript 2 superscript subscript 𝜇 𝑅 2 𝑠 \displaystyle\;{}+\frac{\beta_{0}^{2}}{16}c_{2}\ln^{2}\left(\frac{\mu_{R}^{2}}%
{s}\right)+\frac{\beta_{0}^{2}}{8}C_{F}\ln\left(\frac{\mu_{F}^{2}}{s}\right)%
\ln^{2}\left(\frac{\mu_{R}^{2}}{s}\right) + divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 16 end_ARG italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) + divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 8 end_ARG italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG )
− β 1 8 C F ( ln 2 ( − t 1 m 2 ) + ln 2 ( − u 1 m 2 ) ) + Y 0 ] 𝒟 0 ( s 4 ) } \displaystyle\;\left.\left.{}-\frac{\beta_{1}}{8}C_{F}\left(\ln^{2}\left(\frac%
{-t_{1}}{m^{2}}\right)+\ln^{2}\left(\frac{-u_{1}}{m^{2}}\right)\right)+Y_{0}%
\right]\;{\cal D}_{0}(s_{4})\right\} - divide start_ARG italic_β start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG 8 end_ARG italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ( roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG - italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) + roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG - italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) ) + italic_Y start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ] caligraphic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ) }
where
X 1 subscript 𝑋 1 \displaystyle X_{1} italic_X start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT
= \displaystyle= =
β 0 4 ζ 2 c 3 − β 0 4 c 2 ln ( μ R 2 s ) + β 0 2 C F ln 2 ( − t 1 m 2 ) + β 0 2 C F ln 2 ( − u 1 m 2 ) subscript 𝛽 0 4 subscript 𝜁 2 subscript 𝑐 3 subscript 𝛽 0 4 subscript 𝑐 2 superscript subscript 𝜇 𝑅 2 𝑠 subscript 𝛽 0 2 subscript 𝐶 𝐹 superscript 2 subscript 𝑡 1 superscript 𝑚 2 subscript 𝛽 0 2 subscript 𝐶 𝐹 superscript 2 subscript 𝑢 1 superscript 𝑚 2 \displaystyle\frac{\beta_{0}}{4}\zeta_{2}c_{3}-\frac{\beta_{0}}{4}c_{2}\ln%
\left(\frac{\mu_{R}^{2}}{s}\right)+\frac{\beta_{0}}{2}C_{F}\ln^{2}\left(\frac{%
-t_{1}}{m^{2}}\right)+\frac{\beta_{0}}{2}C_{F}\ln^{2}\left(\frac{-u_{1}}{m^{2}%
}\right) divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT - divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) + divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG - italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) + divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG - italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG )
(2.27)
+ 2 A q ( 2 ) ln ( t 1 u 1 m 4 ) − 2 D q ( 2 ) − β 0 4 C F ln 2 ( μ F 2 s ) + 2 A q ( 2 ) ln ( μ F 2 s ) , 2 superscript subscript 𝐴 𝑞 2 subscript 𝑡 1 subscript 𝑢 1 superscript 𝑚 4 2 superscript subscript 𝐷 𝑞 2 subscript 𝛽 0 4 subscript 𝐶 𝐹 superscript 2 superscript subscript 𝜇 𝐹 2 𝑠 2 superscript subscript 𝐴 𝑞 2 superscript subscript 𝜇 𝐹 2 𝑠 \displaystyle{}+2A_{q}^{(2)}\ln\left(\frac{t_{1}u_{1}}{m^{4}}\right)-2D_{q}^{(%
2)}-\frac{\beta_{0}}{4}C_{F}\ln^{2}\left(\frac{\mu_{F}^{2}}{s}\right)+2A_{q}^{%
(2)}\ln\left(\frac{\mu_{F}^{2}}{s}\right)\,, + 2 italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT roman_ln ( divide start_ARG italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG ) - 2 italic_D start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT - divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) + 2 italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) ,
X 0 subscript 𝑋 0 \displaystyle X_{0} italic_X start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT
= \displaystyle= =
V 2 − 1 2 V 1 2 − 1 4 ζ 2 2 c 3 2 + 3 4 ζ 4 c 3 2 + β 0 4 c 1 ln ( μ R 2 s ) + β 0 6 ζ 3 c 3 − β 0 4 ζ 2 c 2 − β 0 2 ζ 2 C F ln ( μ F 2 μ R 2 ) subscript 𝑉 2 1 2 superscript subscript 𝑉 1 2 1 4 superscript subscript 𝜁 2 2 superscript subscript 𝑐 3 2 3 4 subscript 𝜁 4 superscript subscript 𝑐 3 2 subscript 𝛽 0 4 subscript 𝑐 1 superscript subscript 𝜇 𝑅 2 𝑠 subscript 𝛽 0 6 subscript 𝜁 3 subscript 𝑐 3 subscript 𝛽 0 4 subscript 𝜁 2 subscript 𝑐 2 subscript 𝛽 0 2 subscript 𝜁 2 subscript 𝐶 𝐹 superscript subscript 𝜇 𝐹 2 superscript subscript 𝜇 𝑅 2 \displaystyle V_{2}-\frac{1}{2}V_{1}^{2}-\frac{1}{4}\zeta_{2}^{2}c_{3}^{2}+%
\frac{3}{4}\zeta_{4}c_{3}^{2}+\frac{\beta_{0}}{4}c_{1}\ln\left(\frac{\mu_{R}^{%
2}}{s}\right)+\frac{\beta_{0}}{6}\zeta_{3}c_{3}-\frac{\beta_{0}}{4}\zeta_{2}c_%
{2}-\frac{\beta_{0}}{2}\zeta_{2}C_{F}\ln\left(\frac{\mu_{F}^{2}}{\mu_{R}^{2}}\right) italic_V start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - divide start_ARG 1 end_ARG start_ARG 4 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG 3 end_ARG start_ARG 4 end_ARG italic_ζ start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) + divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 6 end_ARG italic_ζ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT - divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG )
(2.28)
+ 2 A q ( 2 ) ζ 2 − 2 γ q ( 2 ) ln ( μ F 2 s ) + β 0 8 [ 2 γ q ( 1 ) − C F ln ( t 1 u 1 m 4 ) ] ln 2 ( μ F 2 s ) 2 superscript subscript 𝐴 𝑞 2 subscript 𝜁 2 2 superscript subscript 𝛾 𝑞 2 superscript subscript 𝜇 𝐹 2 𝑠 subscript 𝛽 0 8 delimited-[] 2 superscript subscript 𝛾 𝑞 1 subscript 𝐶 𝐹 subscript 𝑡 1 subscript 𝑢 1 superscript 𝑚 4 superscript 2 superscript subscript 𝜇 𝐹 2 𝑠 \displaystyle{}+2A_{q}^{(2)}\zeta_{2}-2\gamma_{q}^{(2)}\ln\left(\frac{\mu_{F}^%
{2}}{s}\right)+\frac{\beta_{0}}{8}\left[2\gamma_{q}^{(1)}-C_{F}\ln\left(\frac{%
t_{1}u_{1}}{m^{4}}\right)\right]\ln^{2}\left(\frac{\mu_{F}^{2}}{s}\right) + 2 italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT italic_ζ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - 2 italic_γ start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) + divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 8 end_ARG [ 2 italic_γ start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT - italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln ( divide start_ARG italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG ) ] roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG )
+ A q ( 2 ) ln ( t 1 u 1 m 4 ) ln ( μ F 2 s ) + β 0 6 C F ln 3 ( − t 1 m 2 ) + β 0 6 C F ln 3 ( − u 1 m 2 ) superscript subscript 𝐴 𝑞 2 subscript 𝑡 1 subscript 𝑢 1 superscript 𝑚 4 superscript subscript 𝜇 𝐹 2 𝑠 subscript 𝛽 0 6 subscript 𝐶 𝐹 superscript 3 subscript 𝑡 1 superscript 𝑚 2 subscript 𝛽 0 6 subscript 𝐶 𝐹 superscript 3 subscript 𝑢 1 superscript 𝑚 2 \displaystyle{}+A_{q}^{(2)}\ln\left(\frac{t_{1}u_{1}}{m^{4}}\right)\ln\left(%
\frac{\mu_{F}^{2}}{s}\right)+\frac{\beta_{0}}{6}C_{F}\ln^{3}\left(\frac{-t_{1}%
}{m^{2}}\right)+\frac{\beta_{0}}{6}C_{F}\ln^{3}\left(\frac{-u_{1}}{m^{2}}\right) + italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT roman_ln ( divide start_ARG italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG ) roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) + divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 6 end_ARG italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ( divide start_ARG - italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) + divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 6 end_ARG italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT roman_ln start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ( divide start_ARG - italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG )
+ ( C F β 0 4 + A q ( 2 ) ) [ ln 2 ( − t 1 m 2 ) + ln 2 ( − u 1 m 2 ) ] , subscript 𝐶 𝐹 subscript 𝛽 0 4 superscript subscript 𝐴 𝑞 2 delimited-[] superscript 2 subscript 𝑡 1 superscript 𝑚 2 superscript 2 subscript 𝑢 1 superscript 𝑚 2 \displaystyle{}+\left(C_{F}\frac{\beta_{0}}{4}+A_{q}^{(2)}\right)\left[\ln^{2}%
\left(\frac{-t_{1}}{m^{2}}\right)+\ln^{2}\left(\frac{-u_{1}}{m^{2}}\right)%
\right]\,, + ( italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG + italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT ) [ roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG - italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) + roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG - italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) ] ,
and
Y 0 subscript 𝑌 0 \displaystyle Y_{0} italic_Y start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT
= \displaystyle= =
− C F β 0 2 24 [ ln 3 ( μ F 2 μ R 2 ) + ln 3 ( μ R 2 s ) ] + 1 16 ( C F β 1 + 8 β 0 A q ( 2 ) ) [ ln 2 ( μ F 2 μ R 2 ) − ln 2 ( μ R 2 s ) ] subscript 𝐶 𝐹 superscript subscript 𝛽 0 2 24 delimited-[] superscript 3 superscript subscript 𝜇 𝐹 2 superscript subscript 𝜇 𝑅 2 superscript 3 superscript subscript 𝜇 𝑅 2 𝑠 1 16 subscript 𝐶 𝐹 subscript 𝛽 1 8 subscript 𝛽 0 superscript subscript 𝐴 𝑞 2 delimited-[] superscript 2 superscript subscript 𝜇 𝐹 2 superscript subscript 𝜇 𝑅 2 superscript 2 superscript subscript 𝜇 𝑅 2 𝑠 \displaystyle-C_{F}\frac{\beta_{0}^{2}}{24}\left[\ln^{3}\left(\frac{\mu_{F}^{2%
}}{\mu_{R}^{2}}\right)+\ln^{3}\left(\frac{\mu_{R}^{2}}{s}\right)\right]+\frac{%
1}{16}\left(C_{F}\beta_{1}+8\beta_{0}A_{q}^{(2)}\right)\left[\ln^{2}\left(%
\frac{\mu_{F}^{2}}{\mu_{R}^{2}}\right)-\ln^{2}\left(\frac{\mu_{R}^{2}}{s}%
\right)\right] - italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 24 end_ARG [ roman_ln start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) + roman_ln start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) ] + divide start_ARG 1 end_ARG start_ARG 16 end_ARG ( italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT italic_β start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 8 italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT ) [ roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) - roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) ]
(2.29)
− ( A q ( 2 ) β 0 + C F β 1 8 ) ln ( t 1 u 1 m 4 ) ln ( μ R 2 s ) + D q ( 2 ) β 0 ln ( μ R 2 s ) − 2 A q ( 3 ) ln ( μ F 2 s ) + 2 D q ( 3 ) superscript subscript 𝐴 𝑞 2 subscript 𝛽 0 subscript 𝐶 𝐹 subscript 𝛽 1 8 subscript 𝑡 1 subscript 𝑢 1 superscript 𝑚 4 superscript subscript 𝜇 𝑅 2 𝑠 superscript subscript 𝐷 𝑞 2 subscript 𝛽 0 superscript subscript 𝜇 𝑅 2 𝑠 2 superscript subscript 𝐴 𝑞 3 superscript subscript 𝜇 𝐹 2 𝑠 2 superscript subscript 𝐷 𝑞 3 \displaystyle{}-\left(A_{q}^{(2)}\beta_{0}+C_{F}\frac{\beta_{1}}{8}\right)\ln%
\left(\frac{t_{1}u_{1}}{m^{4}}\right)\ln\left(\frac{\mu_{R}^{2}}{s}\right)+D_{%
q}^{(2)}\beta_{0}\ln\left(\frac{\mu_{R}^{2}}{s}\right)-2A_{q}^{(3)}\ln\left(%
\frac{\mu_{F}^{2}}{s}\right)+2D_{q}^{(3)} - ( italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT divide start_ARG italic_β start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG 8 end_ARG ) roman_ln ( divide start_ARG italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG ) roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) + italic_D start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) - 2 italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT roman_ln ( divide start_ARG italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_s end_ARG ) + 2 italic_D start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT
− C F 6 β 0 2 [ ln 3 ( − t 1 m 2 ) + ln 3 ( − u 1 m 2 ) ] − β 0 4 ( C F β 0 + 4 A q ( 2 ) ) [ ln 2 ( − t 1 m 2 ) + ln 2 ( − u 1 m 2 ) ] subscript 𝐶 𝐹 6 superscript subscript 𝛽 0 2 delimited-[] superscript 3 subscript 𝑡 1 superscript 𝑚 2 superscript 3 subscript 𝑢 1 superscript 𝑚 2 subscript 𝛽 0 4 subscript 𝐶 𝐹 subscript 𝛽 0 4 superscript subscript 𝐴 𝑞 2 delimited-[] superscript 2 subscript 𝑡 1 superscript 𝑚 2 superscript 2 subscript 𝑢 1 superscript 𝑚 2 \displaystyle{}-\frac{C_{F}}{6}\beta_{0}^{2}\left[\ln^{3}\left(\frac{-t_{1}}{m%
^{2}}\right)+\ln^{3}\left(\frac{-u_{1}}{m^{2}}\right)\right]-\frac{\beta_{0}}{%
4}\left(C_{F}\beta_{0}+4A_{q}^{(2)}\right)\left[\ln^{2}\left(\frac{-t_{1}}{m^{%
2}}\right)+\ln^{2}\left(\frac{-u_{1}}{m^{2}}\right)\right] - divide start_ARG italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT end_ARG start_ARG 6 end_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [ roman_ln start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ( divide start_ARG - italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) + roman_ln start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ( divide start_ARG - italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) ] - divide start_ARG italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG 4 end_ARG ( italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + 4 italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT ) [ roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG - italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) + roman_ln start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG - italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) ]
+ ( β 0 D q ( 2 ) − 2 A q ( 3 ) ) ln ( t 1 u 1 m 4 ) subscript 𝛽 0 superscript subscript 𝐷 𝑞 2 2 superscript subscript 𝐴 𝑞 3 subscript 𝑡 1 subscript 𝑢 1 superscript 𝑚 4 \displaystyle{}+\left(\beta_{0}D_{q}^{(2)}-2A_{q}^{(3)}\right)\ln\left(\frac{t%
_{1}u_{1}}{m^{4}}\right) + ( italic_β start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_D start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT - 2 italic_A start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT ) roman_ln ( divide start_ARG italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_m start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG )
with β 1 = 34 C A 2 / 3 − 2 C F n f − 10 C A n f / 3 subscript 𝛽 1 34 superscript subscript 𝐶 𝐴 2 3 2 subscript 𝐶 𝐹 subscript 𝑛 𝑓 10 subscript 𝐶 𝐴 subscript 𝑛 𝑓 3 \beta_{1}=34C_{A}^{2}/3-2C_{F}n_{f}-10C_{A}n_{f}/3 italic_β start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 34 italic_C start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 3 - 2 italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT - 10 italic_C start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT italic_f end_POSTSUBSCRIPT / 3 [47 , 48 , 49 ] the coefficient of the second term in the perturbative expansion of the β 𝛽 \beta italic_β function.
For processes g g → A B → 𝑔 𝑔 𝐴 𝐵 gg\to AB italic_g italic_g → italic_A italic_B , we replace in Eqs. (2.26 ), (2.27 ), (2.28 ), and (2.29 ) all explicit appearances of C F subscript 𝐶 𝐹 C_{F} italic_C start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT by C A subscript 𝐶 𝐴 C_{A} italic_C start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT , and we use the corresponding expressions for A g subscript 𝐴 𝑔 A_{g} italic_A start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , D g subscript 𝐷 𝑔 D_{g} italic_D start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , γ g subscript 𝛾 𝑔 \gamma_{g} italic_γ start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT , c 3 subscript 𝑐 3 c_{3} italic_c start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT , c 2 subscript 𝑐 2 c_{2} italic_c start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , and c 1 subscript 𝑐 1 c_{1} italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT as detailed above, as well as the appropriate expressions for V 1 subscript 𝑉 1 V_{1} italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and V 2 subscript 𝑉 2 V_{2} italic_V start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT .
Finally, we consider the more general case where the observed particle A 𝐴 A italic_A has mass m A subscript 𝑚 𝐴 m_{A} italic_m start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT , which is different from the mass m B subscript 𝑚 𝐵 m_{B} italic_m start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT of the other final-state particle B 𝐵 B italic_B . In that case, we simply replace m 𝑚 m italic_m by m A subscript 𝑚 𝐴 m_{A} italic_m start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT everywhere it appears explicitly in Eqs. (2.21 ) through (2.29 ), and we also redefine and use t 1 = t − m B 2 subscript 𝑡 1 𝑡 superscript subscript 𝑚 𝐵 2 t_{1}=t-m_{B}^{2} italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = italic_t - italic_m start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT and u 1 = u − m B 2 subscript 𝑢 1 𝑢 superscript subscript 𝑚 𝐵 2 u_{1}=u-m_{B}^{2} italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = italic_u - italic_m start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT everywhere. With these simple changes all the expressions above can be used for the more general case of unequal masses.
It is worth emphasizing that the expressions for the higher-order soft-gluon corrections presented in this section are universal and can be applied to any electroweak production process that occurs through quark-antiquark annihilation or gluon-gluon fusion. The only process-dependent parts are encoded into the virtual correction terms V 1 subscript 𝑉 1 V_{1} italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and V 2 subscript 𝑉 2 V_{2} italic_V start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , which depend on the specific process. For example, explicit expressions of those terms for various electroweak processes can be found in [11 ] .
3 Charged Higgs pair production in the 2HDM
In this section we discuss the production of a charged Higgs pair via quark-antiquark annihilation in the 2HDM that occurs through diagrams of Drell-Yan type. For concreteness, we consider a CP-conserving type II 2HDM with a softly-broken Z 2 subscript 𝑍 2 Z_{2} italic_Z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT symmetry defined by the transformations Φ 1 → Φ 1 → subscript Φ 1 subscript Φ 1 \Phi_{1}\to\Phi_{1} roman_Φ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT → roman_Φ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and Φ 2 → − Φ 2 → subscript Φ 2 subscript Φ 2 \Phi_{2}\to-\Phi_{2} roman_Φ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT → - roman_Φ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , where Φ 1 subscript Φ 1 \Phi_{1} roman_Φ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and Φ 2 subscript Φ 2 \Phi_{2} roman_Φ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT are the two Higgs doublets. In this model, the Higgs doublet Φ 1 subscript Φ 1 \Phi_{1} roman_Φ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT is responsible for the Yukawa interactions that generate masses for the down-type quarks and charged leptons, while Φ 2 subscript Φ 2 \Phi_{2} roman_Φ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT is responsible for the Yukawa interactions that generate masses for the up-type quarks. The Higgs potential is given by:
V 2 H D M subscript 𝑉 2 H D M \displaystyle V_{\rm 2HDM} italic_V start_POSTSUBSCRIPT 2 roman_H roman_D roman_M end_POSTSUBSCRIPT
= \displaystyle= =
m 1 2 | Φ 1 | 2 + m 2 2 | Φ 2 | 2 − m 12 2 ( Φ 1 † Φ 2 + Φ 2 † Φ 1 ) + λ 1 2 | Φ 1 | 4 + λ 2 2 | Φ 2 | 4 superscript subscript 𝑚 1 2 superscript subscript Φ 1 2 superscript subscript 𝑚 2 2 superscript subscript Φ 2 2 superscript subscript 𝑚 12 2 superscript subscript Φ 1 † subscript Φ 2 superscript subscript Φ 2 † subscript Φ 1 subscript 𝜆 1 2 superscript subscript Φ 1 4 subscript 𝜆 2 2 superscript subscript Φ 2 4 \displaystyle m_{1}^{2}|\Phi_{1}|^{2}+m_{2}^{2}|\Phi_{2}|^{2}-m_{12}^{2}(\Phi_%
{1}^{\dagger}\Phi_{2}+\Phi_{2}^{\dagger}\Phi_{1})+\frac{\lambda_{1}}{2}|\Phi_{%
1}|^{4}+\frac{\lambda_{2}}{2}|\Phi_{2}|^{4} italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | roman_Φ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | roman_Φ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( roman_Φ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT roman_Φ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + roman_Φ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT roman_Φ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) + divide start_ARG italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG | roman_Φ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT + divide start_ARG italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG | roman_Φ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT
(3.1)
+ λ 3 | Φ 1 | 2 | Φ 2 | 2 + λ 4 | Φ 1 † Φ 2 | 2 + λ 5 2 [ ( Φ 1 † Φ 2 ) 2 + ( Φ 2 † Φ 1 ) 2 ] , subscript 𝜆 3 superscript subscript Φ 1 2 superscript subscript Φ 2 2 subscript 𝜆 4 superscript superscript subscript Φ 1 † subscript Φ 2 2 subscript 𝜆 5 2 delimited-[] superscript superscript subscript Φ 1 † subscript Φ 2 2 superscript superscript subscript Φ 2 † subscript Φ 1 2 \displaystyle+\lambda_{3}|\Phi_{1}|^{2}|\Phi_{2}|^{2}+\lambda_{4}|\Phi_{1}^{%
\dagger}\Phi_{2}|^{2}+\frac{\lambda_{5}}{2}[(\Phi_{1}^{\dagger}\Phi_{2})^{2}+(%
\Phi_{2}^{\dagger}\Phi_{1})^{2}]\,, + italic_λ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT | roman_Φ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | roman_Φ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_λ start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT | roman_Φ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT roman_Φ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG italic_λ start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG [ ( roman_Φ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT roman_Φ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( roman_Φ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT roman_Φ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] ,
where all the masses (m 1 subscript 𝑚 1 m_{1} italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , m 2 subscript 𝑚 2 m_{2} italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , m 12 subscript 𝑚 12 m_{12} italic_m start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT ) and quartic couplings (λ 1 subscript 𝜆 1 \lambda_{1} italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , λ 2 subscript 𝜆 2 \lambda_{2} italic_λ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , λ 3 subscript 𝜆 3 \lambda_{3} italic_λ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT , λ 4 subscript 𝜆 4 \lambda_{4} italic_λ start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT , λ 5 subscript 𝜆 5 \lambda_{5} italic_λ start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT ) are taken to be real in order to avoid CP violation. After electroweak symmetry breaking, this model predicts the existence of three neutral scalar fields (two CP-even states h ℎ h italic_h and H 𝐻 H italic_H , and a CP-odd state A 𝐴 A italic_A ) and a pair of charged Higgs bosons H ± superscript 𝐻 plus-or-minus H^{\pm} italic_H start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT . The eight parameters of the scalar Higgs potential in Eq. (3.1 ) are usually rewritten in terms of more physical ones: a common choice is to consider the electroweak vacuum expectation value υ = 246 𝜐 246 \upsilon=246 italic_υ = 246 GeV, the four Higgs masses m h subscript 𝑚 ℎ m_{h} italic_m start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT , m H subscript 𝑚 𝐻 m_{H} italic_m start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT , m A subscript 𝑚 𝐴 m_{A} italic_m start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT , m H + subscript 𝑚 superscript 𝐻 m_{H^{+}} italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT , the tangent tan β 𝛽 \tan\beta roman_tan italic_β , the cosine cos ( β − α ) 𝛽 𝛼 \cos(\beta-\alpha) roman_cos ( italic_β - italic_α ) , and the discrete symmetry soft-breaking term m 12 subscript 𝑚 12 m_{12} italic_m start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT (for more details see Section III of [50 ] ). For this particular model, recent studies [25 , 26 , 27 , 28 ] showed that the size of cos ( β − α ) 𝛽 𝛼 \cos(\beta-\alpha) roman_cos ( italic_β - italic_α ) has to be lower than 𝒪 ( 10 − 2 ) 𝒪 superscript 10 2 {\cal O}(10^{-2}) caligraphic_O ( 10 start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT ) and that low values of tan β 𝛽 \tan\beta roman_tan italic_β are favored.
Figure 1: LO Feynman diagrams for quark-induced H + H − superscript 𝐻 superscript 𝐻 H^{+}H^{-} italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT production at the LHC.
As already mentioned, charged-Higgs pair production is a process that can occur through quark-antiquark annihilation and gluon-fusion channels. In this work we consider only quark-induced production which is the dominant production mechanism at low tan β 𝛽 \tan\beta roman_tan italic_β values [19 ] . Leading-order Feynman diagrams for quark-induced production are shown in Fig. 1 . In the case of light quarks in the initial state, we have diagrams of Drell-Yan type where the production of a pair of charged Higgs bosons occurs through the exchange of a photon and a Z 𝑍 Z italic_Z boson in the s 𝑠 s italic_s -channel (left diagram of Fig. 1 ). In the case of bottom quarks, in addition to Drell-Yan type diagrams, the production of charged Higgs bosons also occurs through the exchange of the CP-even neutral Higgs bosons (h , H ℎ 𝐻
h,H italic_h , italic_H ) in the s 𝑠 s italic_s -channel (middle diagram of Fig. 1 ) and the top-quark in the t 𝑡 t italic_t -channel (right diagram of Fig. 1 ). These extra diagrams might interfere with Drell-Yan type ones and introduce some dependence on the scalar potential parameters tan β 𝛽 \tan\beta roman_tan italic_β , cos ( β − α ) 𝛽 𝛼 \cos(\beta-\alpha) roman_cos ( italic_β - italic_α ) and m 12 2 superscript subscript 𝑚 12 2 m_{12}^{2} italic_m start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT .
In this work, we consider only Drell-Yan type production of a pair of charged Higgs bosons in the 5-flavor scheme where we include also a massless bottom quark in the initial state. In this way, our H + H − superscript 𝐻 superscript 𝐻 H^{+}H^{-} italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT production process depends only on the mass of the final-state charged Higgs and is independent of the other 2HDM scalar potential parameters, like tan β 𝛽 \tan\beta roman_tan italic_β , cos ( β − α ) 𝛽 𝛼 \cos(\beta-\alpha) roman_cos ( italic_β - italic_α ) , and m 12 2 superscript subscript 𝑚 12 2 m_{12}^{2} italic_m start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . In considering the total contribution of both light and heavy quarks in the initial state, we are neglecting some bottom-quark contributions (middle and right diagrams of Fig. 1 ), which might give a sizable contribution, depending on the specific values of the Higgs potential parameters. However, for small values of tan β 𝛽 \tan\beta roman_tan italic_β favored by recent fits, i.e. 2 < tan β < 4 2 𝛽 4 2<\tan\beta<4 2 < roman_tan italic_β < 4 , these contributions are of the order of 2-3% at LO and NLO, as also shown in [19 ] , which are smaller than the aNNLO corrections for the Drell-Yan type processes and within the theoretical uncertainty of our calculation obtained from scale and pdf variation. We would like to emphasize that the goal of this section is to provide the most precise prediction for processes that occur through Drell-Yan type diagrams, which for light quarks represent their complete contribution. Providing the total b b ¯ 𝑏 ¯ 𝑏 b\bar{b} italic_b over¯ start_ARG italic_b end_ARG contribution, which depends on the values of several Higgs potential parameters, is beyond the scope of this paper.
Therefore, the quantity F q q ¯ → A B B subscript superscript 𝐹 𝐵 → 𝑞 ¯ 𝑞 𝐴 𝐵 F^{B}_{q{\bar{q}}\to AB} italic_F start_POSTSUPERSCRIPT italic_B end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_q over¯ start_ARG italic_q end_ARG → italic_A italic_B end_POSTSUBSCRIPT entering in the Born differential cross section in Eq. (2.19 ) computed using Drell-Yan type diagrams is given by
F u u ¯ → H + H − B subscript superscript 𝐹 𝐵 → 𝑢 ¯ 𝑢 superscript 𝐻 superscript 𝐻 \displaystyle F^{B}_{u{\bar{u}}\to H^{+}H^{-}} italic_F start_POSTSUPERSCRIPT italic_B end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_u over¯ start_ARG italic_u end_ARG → italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT end_POSTSUBSCRIPT
= \displaystyle= =
1 16 π s 2 α 2 π 2 m Z 4 27 m W 4 t 1 u 1 + m H + 2 ( t 1 + u 1 ) ( m W 2 − m Z 2 ) 2 ( m Z 2 − s ) 2 s 2 1 16 𝜋 superscript 𝑠 2 superscript 𝛼 2 superscript 𝜋 2 superscript subscript 𝑚 𝑍 4 27 superscript subscript 𝑚 𝑊 4 subscript 𝑡 1 subscript 𝑢 1 superscript subscript 𝑚 superscript 𝐻 2 subscript 𝑡 1 subscript 𝑢 1 superscript superscript subscript 𝑚 𝑊 2 superscript subscript 𝑚 𝑍 2 2 superscript superscript subscript 𝑚 𝑍 2 𝑠 2 superscript 𝑠 2 \displaystyle\frac{1}{16\pi s^{2}}\frac{\alpha^{2}\pi^{2}m_{Z}^{4}}{27m_{W}^{4%
}}\frac{t_{1}\,u_{1}+m_{H^{+}}^{2}(t_{1}+u_{1})}{(m_{W}^{2}-m_{Z}^{2})^{2}(m_{%
Z}^{2}-s)^{2}s^{2}} divide start_ARG 1 end_ARG start_ARG 16 italic_π italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG italic_α start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG 27 italic_m start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG divide start_ARG italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) end_ARG start_ARG ( italic_m start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_m start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_s ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG
(3.2)
× [ 128 m W 8 + 17 m Z 4 s 2 − 32 m W 6 ( 8 m Z 2 + s ) − 4 m W 2 m Z 2 s ( 20 m Z 2 + 7 s ) \displaystyle\times\left[128m_{W}^{8}+17m_{Z}^{4}s^{2}-32m_{W}^{6}(8m_{Z}^{2}+%
s)-4m_{W}^{2}m_{Z}^{2}s(20m_{Z}^{2}+7s)\right. × [ 128 italic_m start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT + 17 italic_m start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 32 italic_m start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT ( 8 italic_m start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_s ) - 4 italic_m start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_s ( 20 italic_m start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 7 italic_s )
+ 4 m W 4 ( 32 m Z 4 + 28 m Z 2 s + 5 s 2 ) ] \displaystyle\left.\quad\;{}+4m_{W}^{4}(32m_{Z}^{4}+28m_{Z}^{2}s+5s^{2})\right] + 4 italic_m start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT ( 32 italic_m start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT + 28 italic_m start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_s + 5 italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ]
for up-type quarks, while for down-type quarks we have
F d d ¯ → H + H − B subscript superscript 𝐹 𝐵 → 𝑑 ¯ 𝑑 superscript 𝐻 superscript 𝐻 \displaystyle F^{B}_{d{\bar{d}}\to H^{+}H^{-}} italic_F start_POSTSUPERSCRIPT italic_B end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_d over¯ start_ARG italic_d end_ARG → italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT end_POSTSUBSCRIPT
= \displaystyle= =
1 16 π s 2 α 2 π 2 m Z 4 27 m W 4 t 1 u 1 + m H + 2 ( t 1 + u 1 ) ( m W 2 − m Z 2 ) 2 ( m Z 2 − s ) 2 s 2 1 16 𝜋 superscript 𝑠 2 superscript 𝛼 2 superscript 𝜋 2 superscript subscript 𝑚 𝑍 4 27 superscript subscript 𝑚 𝑊 4 subscript 𝑡 1 subscript 𝑢 1 superscript subscript 𝑚 superscript 𝐻 2 subscript 𝑡 1 subscript 𝑢 1 superscript superscript subscript 𝑚 𝑊 2 superscript subscript 𝑚 𝑍 2 2 superscript superscript subscript 𝑚 𝑍 2 𝑠 2 superscript 𝑠 2 \displaystyle\frac{1}{16\pi s^{2}}\frac{\alpha^{2}\pi^{2}m_{Z}^{4}}{27m_{W}^{4%
}}\frac{t_{1}\,u_{1}+m_{H^{+}}^{2}(t_{1}+u_{1})}{(m_{W}^{2}-m_{Z}^{2})^{2}(m_{%
Z}^{2}-s)^{2}s^{2}} divide start_ARG 1 end_ARG start_ARG 16 italic_π italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG italic_α start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG 27 italic_m start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG divide start_ARG italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_u start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) end_ARG start_ARG ( italic_m start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_m start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_m start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_s ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG
(3.3)
× [ 32 m W 8 + 5 m Z 4 s 2 + 16 m W 6 ( − 4 m Z 2 + s ) − 8 m W 2 m Z 2 s ( m Z 2 + 2 s ) \displaystyle\times\left[32m_{W}^{8}+5m_{Z}^{4}s^{2}+16m_{W}^{6}(-4m_{Z}^{2}+s%
)-8m_{W}^{2}m_{Z}^{2}s(m_{Z}^{2}+2s)\right. × [ 32 italic_m start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT + 5 italic_m start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 16 italic_m start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT ( - 4 italic_m start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_s ) - 8 italic_m start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_s ( italic_m start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 2 italic_s )
+ 4 m W 4 ( 8 m Z 4 − 2 m Z 2 s + 5 s 2 ) ] . \displaystyle\left.\quad\;{}+4m_{W}^{4}(8m_{Z}^{4}-2m_{Z}^{2}s+5s^{2})\right]\,. + 4 italic_m start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT ( 8 italic_m start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT - 2 italic_m start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_s + 5 italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ] .
These expressions are used in Eq. (2.19 ) to compute the Born cross section, and in Eqs. (2.20 ), (LABEL:NNLO ), and (2.26 ) to derive the higher-order soft-gluon corrections through N3 LO.
4 Total cross sections for H + H − superscript 𝐻 superscript 𝐻 H^{+}H^{-} italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT production at the LHC
In this section we present results for the total cross section of charged-Higgs pair production via quark-antiquark annihilation q q ¯ → H + H − → 𝑞 ¯ 𝑞 superscript 𝐻 superscript 𝐻 q\bar{q}\to H^{+}H^{-} italic_q over¯ start_ARG italic_q end_ARG → italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT at the LHC, for two center-of-mass energies, namely 13 TeV and 13.6 TeV. As already discussed in the previous section, this cross section depends only on m H + subscript 𝑚 superscript 𝐻 m_{H^{+}} italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT and is independent of the other Higgs potential parameters. In our calculation we work in the 5-flavor scheme and we include the bottom quark in the initial state, which is taken to be massless. Lower bounds on the mass of the charged Higgs of the order of 570-800 GeV have been obtained in Ref. [51 ] . In light of these constraints, the charged-Higgs mass is taken here to vary between 500 and 1500 GeV. Lower values of m H + subscript 𝑚 superscript 𝐻 m_{H^{+}} italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT are excluded, while higher values, although allowed, are not considered because they give very small cross sections. The complete LO and NLO QCD results are calculated using MadGraph5_aMC@NLO [52 ] and an ad hoc UFO model for the CP-conserving type II 2HDM, which can be found at [53 ] . We take m W = 80.377 subscript 𝑚 𝑊 80.377 m_{W}=80.377 italic_m start_POSTSUBSCRIPT italic_W end_POSTSUBSCRIPT = 80.377 GeV, m Z = 91.1876 subscript 𝑚 𝑍 91.1876 m_{Z}=91.1876 italic_m start_POSTSUBSCRIPT italic_Z end_POSTSUBSCRIPT = 91.1876 GeV, and α − 1 = 127.9 superscript 𝛼 1 127.9 \alpha^{-1}=127.9 italic_α start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = 127.9 . We use the same MSHT20 aN3 LO pdf set [54 ] for computing results at every perturbative order, in order to show how each order in the series contributes to the aN3 LO cross section.
For the total cross section, the central results are obtained by setting the factorization and renormalization scales both equal to m H + subscript 𝑚 superscript 𝐻 m_{H^{+}} italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT . Scale uncertainties are obtained by varying independently μ R subscript 𝜇 𝑅 \mu_{R} italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT and μ F subscript 𝜇 𝐹 \mu_{F} italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT between m H + / 2 subscript 𝑚 superscript 𝐻 2 m_{H^{+}}/2 italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT / 2 and 2 m H + 2 subscript 𝑚 superscript 𝐻 2m_{H^{+}} 2 italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT , and pdf uncertainties are also computed. The cross sections at aNNLO are computed by adding the second-order soft-plus-virtual QCD corrections to the complete NLO result. The third-order soft-gluon corrections are further added to derive the result at aN3 LO.
Figure 2: The total cross sections at LO, NLO, aNNLO, and aN3 LO for q q ¯ → H + H − → 𝑞 ¯ 𝑞 superscript 𝐻 superscript 𝐻 q\bar{q}\to H^{+}H^{-} italic_q over¯ start_ARG italic_q end_ARG → italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT production in p p 𝑝 𝑝 pp italic_p italic_p collisions at 13 TeV energy with μ F = μ R = m H + subscript 𝜇 𝐹 subscript 𝜇 𝑅 subscript 𝑚 superscript 𝐻 \mu_{F}=\mu_{R}=m_{H^{+}} italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT = italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT = italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT using MSHT20 aN3 LO pdf. The inset plot displays the NLO/LO, aNNLO/LO, and aN3 LO/LO K 𝐾 K italic_K -factors.
Figure 3: The total cross sections at LO, NLO, aNNLO, and aN3 LO for q q ¯ → H + H − → 𝑞 ¯ 𝑞 superscript 𝐻 superscript 𝐻 q\bar{q}\to H^{+}H^{-} italic_q over¯ start_ARG italic_q end_ARG → italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT production in p p 𝑝 𝑝 pp italic_p italic_p collisions at 13.6 TeV energy with μ F = μ R = m H + subscript 𝜇 𝐹 subscript 𝜇 𝑅 subscript 𝑚 superscript 𝐻 \mu_{F}=\mu_{R}=m_{H^{+}} italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT = italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT = italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT using MSHT20 aN3 LO pdf. The inset plot displays the NLO/LO, aNNLO/LO, and aN3 LO/LO K 𝐾 K italic_K -factors.
In Fig. 2 and Fig. 3 we show the total cross section for H + H − superscript 𝐻 superscript 𝐻 H^{+}H^{-} italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT production at the LHC with 13 TeV and 13.6 TeV energy, respectively. These plots show the central value, obtained with μ = m H + 𝜇 subscript 𝑚 superscript 𝐻 \mu=m_{H^{+}} italic_μ = italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT (where μ = μ R = μ F 𝜇 subscript 𝜇 𝑅 subscript 𝜇 𝐹 \mu=\mu_{R}=\mu_{F} italic_μ = italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT = italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT ), of the LO, NLO, aNNLO, and aN3 LO cross sections as a function of the mass of the charged Higgs boson. The inset plots display the K 𝐾 K italic_K -factors of the higher-order cross sections relative to the LO result. At each perturbative order, we find that the K 𝐾 K italic_K -factors slowly increase with increasing charged-Higgs mass. For instance, at both 13 TeV and 13.6 TeV, the NLO K 𝐾 K italic_K -factor varies from 1.15 for m H + = 500 subscript 𝑚 superscript 𝐻 500 m_{H^{+}}=500 italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT = 500 GeV to 1.20 for m H + = 1500 subscript 𝑚 superscript 𝐻 1500 m_{H^{+}}=1500 italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT = 1500 GeV. The NNLO soft-plus-virtual corrections provide an additional 5% to 6% increase, and the N3 LO soft-gluon corrections give a further 0.5% to 0.6% increase. Thus, the aN3 LO K 𝐾 K italic_K -factor varies from 1.21 for m H + = 500 subscript 𝑚 superscript 𝐻 500 m_{H^{+}}=500 italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT = 500 GeV to 1.26 for m H + = 1500 subscript 𝑚 superscript 𝐻 1500 m_{H^{+}}=1500 italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT = 1500 GeV.
Table 1: The total cross section for p p → H + H − → 𝑝 𝑝 superscript 𝐻 superscript 𝐻 pp\to H^{+}H^{-} italic_p italic_p → italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT via quark-antiquark annihilation at LO, NLO, aNNLO, and aN3 LO at the LHC with S = 13 𝑆 13 \sqrt{S}=13 square-root start_ARG italic_S end_ARG = 13 TeV and MSHT20 aN3 LO pdf. The central results are for μ F = μ R = m H + subscript 𝜇 𝐹 subscript 𝜇 𝑅 subscript 𝑚 superscript 𝐻 \mu_{F}=\mu_{R}=m_{H^{+}} italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT = italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT = italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT and are shown together with scale and pdf uncertainties.
Table 2: The total cross section for p p → H + H − → 𝑝 𝑝 superscript 𝐻 superscript 𝐻 pp\to H^{+}H^{-} italic_p italic_p → italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT via quark-antiquark annihilation at LO, NLO, aNNLO, and aN3 LO at the LHC with S = 13.6 𝑆 13.6 \sqrt{S}=13.6 square-root start_ARG italic_S end_ARG = 13.6 TeV and MSHT20 aN3 LO pdf. The central results are for μ F = μ R = m H + subscript 𝜇 𝐹 subscript 𝜇 𝑅 subscript 𝑚 superscript 𝐻 \mu_{F}=\mu_{R}=m_{H^{+}} italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT = italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT = italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT and are shown together with scale and pdf uncertainties.
In Table 1 and Table 2 we present some numerical values for the total cross section of H + H − superscript 𝐻 superscript 𝐻 H^{+}H^{-} italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT production in proton-proton collisions, with a center-of-mass energy of 13 TeV and 13.6 TeV, respectively. These cross sections are shown for select charged-Higgs masses, namely 600, 800, 1000, and 1200 GeV. They are calculated at LO, NLO, aNNLO, and aN3 LO and they are shown together with scale and pdf uncertainties. We want to stress that the soft-gluon corrections are numerically dominant. The result derived with only NLO soft-gluon corrections differs from the exact NLO cross section by only around one percent or less for the masses and energies considered here. Furthermore, the aNNLO result would change by only 2 per mille if we did not include the δ ( s 4 ) 𝛿 subscript 𝑠 4 \delta(s_{4}) italic_δ ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ) terms (i.e. the virtual corrections) at NNLO; the soft-gluon corrections are numerically by far the dominant contribution to the higher-order QCD corrections.
As shown in Table 1 for the results at 13 TeV energy, the scale uncertainty at LO is large and is reduced significantly at NLO. The higher-order soft-plus-virtual corrections at NNLO decrease the uncertainty even more. Finally, the addition of the N3 LO soft-gluon corrections further decreases this uncertainty for larger charged-Higgs masses but not for smaller ones. The pdf uncertainty varies from about ± 3 % plus-or-minus percent 3 \pm 3\% ± 3 % for a mass of 600 GeV to ± 7 % plus-or-minus percent 7 \pm 7\% ± 7 % for a 1200 GeV mass, at each perturbative order. The results in Table 2 at 13.6 TeV energy show similar behavior for the scale uncertainty at different perturbative orders. Also in this case, the pdf uncertainty varies at each perturbative order from about ± 3 % plus-or-minus percent 3 \pm 3\% ± 3 % for m H + = 600 subscript 𝑚 superscript 𝐻 600 m_{H^{+}}=600 italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT = 600 GeV to ± 6 % plus-or-minus percent 6 \pm 6\% ± 6 % for m H + = 1200 subscript 𝑚 superscript 𝐻 1200 m_{H^{+}}=1200 italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT = 1200 GeV.
Regarding pdf uncertainties, we can see that, for both center-of-mass energies, they increase with increasing mass of the charged Higgs boson. Moreover, they become bigger than the scale ones already at NLO and therefore, starting at that order, they should be considered the dominant source of theoretical uncertainty. As argued in [54 ] , with the use of aN3 LO pdf, the factorisation scale variation is already contained within the predicted pdf uncertainties.
It is also interesting to study the relative contributions of the different logarithmic powers to the higher-order corrections. As an example, we consider a charged-Higgs mass of 600 GeV and an energy of 13.6 TeV. Setting μ F = μ R = m H + subscript 𝜇 𝐹 subscript 𝜇 𝑅 subscript 𝑚 superscript 𝐻 \mu_{F}=\mu_{R}=m_{H^{+}} italic_μ start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT = italic_μ start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT = italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT and using the notation of Eq. (2.21 ), the NNLO corrections from the 𝒟 3 subscript 𝒟 3 {\cal D}_{3} caligraphic_D start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT terms alone are 0.0134 fb, from the 𝒟 3 subscript 𝒟 3 {\cal D}_{3} caligraphic_D start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT terms plus the 𝒟 2 subscript 𝒟 2 {\cal D}_{2} caligraphic_D start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT terms they are 0.0189 fb, from 𝒟 3 + 𝒟 2 + 𝒟 1 subscript 𝒟 3 subscript 𝒟 2 subscript 𝒟 1 {\cal D}_{3}+{\cal D}_{2}+{\cal D}_{1} caligraphic_D start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT + caligraphic_D start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + caligraphic_D start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT they are 0.0150 fb, and from 𝒟 3 + 𝒟 2 + 𝒟 1 + 𝒟 0 subscript 𝒟 3 subscript 𝒟 2 subscript 𝒟 1 subscript 𝒟 0 {\cal D}_{3}+{\cal D}_{2}+{\cal D}_{1}+{\cal D}_{0} caligraphic_D start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT + caligraphic_D start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + caligraphic_D start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + caligraphic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT they are 0.0118 fb. By further adding the NNLO δ ( s 4 ) 𝛿 subscript 𝑠 4 \delta(s_{4}) italic_δ ( italic_s start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ) terms, we get the final NNLO soft-plus-virtual correction of 0.0123 fb. Thus, we see that the virtual contribution to the NNLO soft-plus-virtual corrections is small, and the soft-gluon corrections alone are dominant. We also see that the leading-logarithmic contribution (i.e. from 𝒟 3 subscript 𝒟 3 {\cal D}_{3} caligraphic_D start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ) alone is quite close to the complete soft-plus-virtual contribution at NNLO.
For the N3 LO soft-gluon contributions, again with a 600 GeV charged-Higgs mass at 13.6 TeV energy, we find a contribution of 0.00993 fb from the 𝒟 5 subscript 𝒟 5 {\cal D}_{5} caligraphic_D start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT terms alone, 0.0176 fb from 𝒟 5 + 𝒟 4 subscript 𝒟 5 subscript 𝒟 4 {\cal D}_{5}+{\cal D}_{4} caligraphic_D start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT + caligraphic_D start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT , 0.0155 fb from 𝒟 5 + 𝒟 4 + 𝒟 3 subscript 𝒟 5 subscript 𝒟 4 subscript 𝒟 3 {\cal D}_{5}+{\cal D}_{4}+{\cal D}_{3} caligraphic_D start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT + caligraphic_D start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT + caligraphic_D start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT , 0.00898 fb from 𝒟 5 + 𝒟 4 + 𝒟 3 + 𝒟 2 subscript 𝒟 5 subscript 𝒟 4 subscript 𝒟 3 subscript 𝒟 2 {\cal D}_{5}+{\cal D}_{4}+{\cal D}_{3}+{\cal D}_{2} caligraphic_D start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT + caligraphic_D start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT + caligraphic_D start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT + caligraphic_D start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , 0.00200 from
𝒟 5 + 𝒟 4 + 𝒟 3 + 𝒟 2 + 𝒟 1 subscript 𝒟 5 subscript 𝒟 4 subscript 𝒟 3 subscript 𝒟 2 subscript 𝒟 1 {\cal D}_{5}+{\cal D}_{4}+{\cal D}_{3}+{\cal D}_{2}+{\cal D}_{1} caligraphic_D start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT + caligraphic_D start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT + caligraphic_D start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT + caligraphic_D start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + caligraphic_D start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , and finally 0.00133 fb from 𝒟 5 + 𝒟 4 + 𝒟 3 + 𝒟 2 + 𝒟 1 + 𝒟 0 subscript 𝒟 5 subscript 𝒟 4 subscript 𝒟 3 subscript 𝒟 2 subscript 𝒟 1 subscript 𝒟 0 {\cal D}_{5}+{\cal D}_{4}+{\cal D}_{3}+{\cal D}_{2}+{\cal D}_{1}+{\cal D}_{0} caligraphic_D start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT + caligraphic_D start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT + caligraphic_D start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT + caligraphic_D start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + caligraphic_D start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + caligraphic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT which is the total N3 LO soft-gluon correction. Thus, at N3 LO, the leading-logarithmic contribution (i.e. from 𝒟 5 subscript 𝒟 5 {\cal D}_{5} caligraphic_D start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT ) is very different from the complete soft-gluon correction, in contrast to what we found above for the NNLO corrections. We also see that the N3 LO soft-gluon corrections are an order of magnitude smaller than those at NNLO.
As mentioned in the previous section, for small values of tan β 𝛽 \tan\beta roman_tan italic_β , the additional contributions from b b ¯ → H + H − → 𝑏 ¯ 𝑏 superscript 𝐻 superscript 𝐻 b{\bar{b}}\to H^{+}H^{-} italic_b over¯ start_ARG italic_b end_ARG → italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT are rather small and less than the aNNLO corrections of the Drell-Yan type q q ¯ → H + H − → 𝑞 ¯ 𝑞 superscript 𝐻 superscript 𝐻 q{\bar{q}}\to H^{+}H^{-} italic_q over¯ start_ARG italic_q end_ARG → italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT italic_H start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT diagrams. One can calculate these b b ¯ 𝑏 ¯ 𝑏 b{\bar{b}} italic_b over¯ start_ARG italic_b end_ARG contributions at NLO for any choice of the Higgs-potential parameters, and then add them to our aN3 LO result for the Drell-Yan type q q ¯ 𝑞 ¯ 𝑞 q{\bar{q}} italic_q over¯ start_ARG italic_q end_ARG contributions. For example, for tan β = 3 𝛽 3 \tan\beta=3 roman_tan italic_β = 3 , m H + = 600 subscript 𝑚 superscript 𝐻 600 m_{H^{+}}=600 italic_m start_POSTSUBSCRIPT italic_H start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT end_POSTSUBSCRIPT = 600 GeV, m H = 400 subscript 𝑚 𝐻 400 m_{H}=400 italic_m start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT = 400 GeV, cos ( β − α ) = 0 𝛽 𝛼 0 \cos(\beta-\alpha)=0 roman_cos ( italic_β - italic_α ) = 0 , and m 12 = 0 subscript 𝑚 12 0 m_{12}=0 italic_m start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT = 0 GeV, the missing b b ¯ 𝑏 ¯ 𝑏 b{\bar{b}} italic_b over¯ start_ARG italic_b end_ARG contribution at 13 TeV energy is 0.0035 fb. Looking at the first line of Table 1, we see that this is much smaller than the aNNLO corrections of 0.011 fb for q q ¯ 𝑞 ¯ 𝑞 q{\bar{q}} italic_q over¯ start_ARG italic_q end_ARG as well as the total uncertainty of the overall q q ¯ 𝑞 ¯ 𝑞 q{\bar{q}} italic_q over¯ start_ARG italic_q end_ARG cross section. Adding the b b ¯ 𝑏 ¯ 𝑏 b{\bar{b}} italic_b over¯ start_ARG italic_b end_ARG contribution to the aN3 LO q q ¯ 𝑞 ¯ 𝑞 q{\bar{q}} italic_q over¯ start_ARG italic_q end_ARG cross section gives a total of 0.232 fb. Clearly, in this case it is not necessary to calculate the aNNLO and aN3 LO corrections for the b b ¯ 𝑏 ¯ 𝑏 b{\bar{b}} italic_b over¯ start_ARG italic_b end_ARG contribution since they would be much smaller than the precision of the numbers displayed.