arxiv: v1 [hep-ph] 4 Aug 2014

Transcription

1 Prepared for sumission to JHEP HRI-RECAPP Higgs oson production through annihilation at threshold in N LO QCD arxiv: v1 [hep-ph] 4 Aug 014 Taushif Ahmed, a Narayan Rana a and V. Ravindran a Regional Centre for Accelerator-ased Particle Physics, Harish-Chandra Research Institute, Allahaad, India The Institute of Mathematical Sciences, Chennai, India taushif@hri.res.in, narayan@hri.res.in, ravindra@imsc.res.in Astract: We present threshold enhanced N LO QCD corrections to inclusive Higgs production through ottom anti-ottom annihilation at hadron colliders using threshold resummed cross section. The resummed cross section is otained using factorization properties and Sudakov resummation of the inclusive cross section. We use the recent results on threshold N LO corrections in QCD for Drell-Yan production and three loop QCD corrections to Higgs form factor with ottom anti-ottom quark to achieve this task. This is the first step towards the evaluation of complete N LO result. We have numerically demonstrated the importance of such corrections at the LHC. Keywords: QCD, Higgs, Threshold corrections

2 Contents 1 Introduction 1 Threshold resummation Results 1 Introduction The discovery of Higgs oson y ATLAS [1] and CMS [] collaorations of the Large Hadron Collider LHC at CERN has not only shed the light on the dynamics ehind the electroweak symmetry reaking ut also put the Standard Model SM of particle physics on a firmer ground. In the SM, the elementary particles such as quarks, leptons and gauge osons, Z, W ± acquire their masses through spontaneous symmetry reaking SSB. The Higgs mechanism provides the framework for SSB. The SM predicts the existence of a Higgs oson whose mass is a parameter of the model. The recent discovery of Higgs oson provides a valuale information on this, namely on its mass which is aout 15.5 GeV. The searches for the Higgs oson have een going on for several decades in various experiments. Earlier experiments such as LEP [] and Tevatron [4] played an important role in the discovery y the LHC collaorations through narrowing down its possile mass range. LEP excluded Higgs oson of mass elow GeV and their precision electroweak measurements [5] hinted the mass less than 15 GeV at 5% confidence level CL, while Tevatron excluded Higgs oson of mass in the range GeV at 5% CL. Higgs osons are produced dominantly at the LHC via gluon gluon fusion through top quark loop, while the su-dominant ones are vector oson fusion, associated production of Higgs oson with vector osons, with top anti-top pairs and also in ottom anti-ottom annihilation. The inclusive productions of Higgs oson in gluon gluon [6], vector oson fusion processes [7] and associated production with vector osons [8] are known to next to next to leading order NNLO accuracy in QCD. Higgs production in ottom anti-ottom annihilation is also known to NNLO accuracy in the variale flavour scheme VFS [, 10], while it is known to NLO in the fixed flavour scheme FFS [11]. In the minimal super symmetric standard model MSSM, the coupling of ottom quarks to Higgs ecomes large in the large tan β region, where tan β is the ratio of vacuum expectation values of up and down type Higgs fields. This can enhance contriutions from ottom anti-ottom annihilation suprocesses. While the theoretical predictions of NNLO [6] and next to next to leading log NNLL [1] QCD corrections and of two loop electroweak effects [1] played an important role in the Higgs discovery, the theoretical uncertainties resulting from factorization and renormalization scales are not fully under control. Hence, the efforts to go eyond NNLO are 1

3 going on intensively. Some of the ingredients to otain N LO QCD corrections are already availale. For example, quark and gluon form factors [14 17], the mass factorization kernels [18] and the renormalization constant [1] for the effective operator descriing the coupling of Higgs oson with the SM fields in the infinite top quark mass limit up to three loop level in dimensional regularization are known for some time. In addition, NNLO soft contriutions are known [0] to all orders in ɛ for oth DY and Higgs productions using dimensional regularization with space time dimension eing d = 4 + ɛ. They were used to otain the partial N LO threshold effects [1, ] to Drell-Yan production of di-leptons and inclusive productions of Higgs oson through gluon gluon fusion and in ottom antiottom annihilation. Threshold contriution to the inclusive production cross section is expanded in terms of δ1 z and D i z where ln1 z D i z = z with the scaling parameter z = m H /ŝ for Higgs and z = m l + l /ŝ for DY. Here m H, m l + l and ŝ are mass of the Higgs oson, invariant mass of the di-leptons and center of mass energy of the partonic reaction responsile for production mechanism respectively. The missing δ1 z terms for the complete N LO threshold contriutions to the Higgs production through gluon gluon fusion are now availale due to the seminal work y Anastasiou et al [] where the relevant soft contriutions were otained from the real radiations at N LO level. In [4], we exploited the universal structure of the soft radiations to otain the corresponding soft gluon contriutions to DY production, which led to the evaluation of missing δ1 z part of the N LO threshold corrections, later confirmed in [5]. For the Higgs production through annihilation, till date, only partial N LO threshold corrections are known []. It was not possile to determine the δ1 z at N LO due to the lack of information on three loop finite part of ottom anti-ottom higgs form factor in QCD and the soft gluon radiation at N LO level. The recent results on Higgs form factor with ottom anti-ottom y Gehrmann and Kara [6] and on the universal soft distriution otained for the Drell-Yan production [4] can now e used to otain δ1 z part of the threshold N LO contriution. For the soft gluon radiations in the annihilation, the results from [4] can e used as they do not depend on the flavour of the incoming quark states. We have set ottom quark mass to e zero throughout except in the Yukawa coupling. In the next section, we present the details of the threshold resummation and in section we present our results for threshold N LO QCD contriutions to Higgs production through annihilation at hadron colliders and their numerical impact. The numerical impact of threshold enhanced N LO contriutions is demonstrated for the LHC energy s = 14 TeV y studying the staility of the perturation theory under factorization and renormalization scales. Finally we give a rief summary of our findings. Threshold resummation The interaction of ottom quarks and Higgs oson is given y the action SI = λ d 4 x φxψ xψ x.1 +

4 where ψ x denotes the ottom quark field and φx the scalar field. λ is the Yukawa coupling given y m /ν, with the ottom quark mass m and the vacuum expectation value ν 46 GeV. In MSSM, for the pseudoscalar Higgs oson, we need to replace λφxψ xψ x y λ φxψ xγ 5 ψ x in the aove equation. The MSSM couplings are given y m sin α ν cos β, φ = h, λ = m cos α ν cos β, φ = H, m tan β ν, φ = A respectively. The angle α measures the mixing of weak and mass eigenstates of neutral Higgs osons. We use VFS scheme throughout, hence except in the Yukawa coupling, m is taken to e zero like other light quarks in the theory. The inclusive Higgs production through ottom anti-ottom annihilation can e computed using σ s, q = σ 0 µ R ac dx 1 dx f a x 1, µ F f c x, µ F acŝ, q, µ F, µ R,. where f a x 1, µ F and f cx, µ F are parton distriution functions with momentum fractions x 1 and x respectively. µ F is the factorization scale and ŝ = x 1 x s where s ŝ is the square of hadronic partonic center of mass energy. The orn cross section is given y σ 0 = πλ µ R 1m H.. The orn normalized partonic suprocesses after mass factorization are denoted y ac where initial state partons are a and c. Here, q = m H, µ F results from mass factorization and the renormalization scale µ R is due to UV renormalization. ac can e decomposed into two parts denoted y sv and,r ac. acz, q, µ R, µ F = sv z, q, µ R, µ F +,R ac z, q, µ R, µ F.4 where sv contains only distriutions such as δ1 z and D i, often called threshold contriutions and the second term contains regular terms in z. Following [], the threshold resumed cross section at the partonic level in d = 4 + ɛ dimensions is given y sv z, q, µ R, µ F = C exp Ψ z, q, µ R, µ F, ɛ.5 ɛ=0 where the scaling variale z = q /ŝ and Ψ z, q, µ R, µ F, ɛ is a finite distriution. The aove resummed expression follows from factorization properties of the inclusive cross section and Sudakov resummation of soft gluons in the QCD amplitudes. The symol C implies convolution with the following expansion Ce fz = δ1 z + 1 1! fz + 1 fz fz ! Here means Mellin convolution and fz is a distriution of the kind δ1 z and D i. We drop all the regular terms in z in the evaluation of threshold contriution sv.

5 In d = 4 + ɛ dimensions, the distriution Ψ receives contriutions from UV renormalization constant Z for Yukawa coupling λ, the Higgs form factor ˆF from ottom anti-ottom, soft gluon distriution Φ from the real radiations in partonic suprocesses and the mass factorization kernels Γ that remove the collinear singularities from the initial state ottom quark states. That is [ Ψ z, q, µ R, µ F, ɛ = ln Z â s, µ R, µ, ɛ] + ln ˆF â s, Q, µ, ɛ δ1 z.7 + Φ â s, q, µ, z, ɛ C ln Γ â s, µ, µ F, z, ɛ..8 While the individual contriutions are divergent, the sum is a finite distriution. The scale µ is introduced to define the dimensionless coupling constant â s = ĝs/16π in dimensional regularization and Q = q. The renormalized strong coupling constant a s µ R is related to are â s through strong coupling constant renormalization Zµ R, µ ɛzµ â s = µ R Sɛ 1 a s µ R, S ɛ = exp γ E ln 4π ɛ.. R where γ E is Euler-Mascheroni constant. In 4 + ɛ dimensions, Zµ R can e expressed in terms of the coefficients β i of the β function of the strong coupling RG equation. Solving the RGE, we otain, up to three loop level Zµ R = 1 + a s µ R β 0 ɛ + a sµ R The coefficients β 0, β 1 and β are 4β0 ɛ + β 1 ɛ β 0 = 11 C A 4 T F n f, + a sµ R β 1 = 4 C A 4T F n f C F 0 T F n f C A, β = C A C AT F n f C AT F n f 8β0 ɛ + 14β 0β 1 ɛ + β ɛ C F T F n f 05 C F C A T F n f + C F T F n f.11 where the SUN QCD color factors are given y C A = N, C F = N 1 N, T F = 1 and n f is the numer of active flavours. The UV renormalization constant Z â s, µ R, µ, ɛ for Yukawa coupling λ satisfies d µ R dµ R ln Z â s, µ R, µ, ɛ = a i sµ Rγi 1. i=1.1 The aove renormalization group equation RGE can e solved in 4 + ɛ dimensions to otain up to Oa s level: µ Z â s, µ R, µ ɛ ], ɛ = 1 + â R [1 µ s Sɛ µ γ0 + â ɛ [ R s ɛ µ Sɛ 1 ɛ γ0 β0 γ0 4

6 + 1 ɛ γ 1 ] µ + â ɛ [ R s S 1 4 µ ɛ ɛ + 1 ɛ γ 0 γ 1 β 1 γ 0 8 β 0 γ 1 γ0 4 β0 γ ɛ ] γ β 0 γ 0.1 where the anomalous dimensions γ i can e otained from the quark mass anomalous dimensions [7] γ 0 = C F, γ1 = C F C F C A 10 C F T F n f, γ = 1 C F 1 4 C F C A ζ ζ C F T F n f 108 C F CA + C F C A T F n f C F TF n f..14 The are form factor ˆF â s, Q, µ, ɛ satisfies the following differential equation [8 1] Q d dq ln ˆF = 1 [ K â s, µ R µ, ɛ + G â s, Q µ R, µ R µ, ɛ ].15 where K contains all the poles in ɛ and G contains the terms finite in ɛ. Renormalization group invariance of ˆF â s, Q, µ, ɛ gives d µ R dµ R K = µ R d dµ R G = where A q i s are the cusp anomalous dimensions [18, ] given y A q 1 = 4C F, A q 67 = 8C F C A 18 ζ + 8C F n f 5, A q 45 = 16C F CA 4 67ζ + 11ζ + 11ζ C F C A n f ζ 7ζ a i sµ RA q i..16 i=1 + 16C F n f C F n f 1 7 Expanding µ R independent part of the solution of RG equation for G as 4 + ζ..17 G a s Q, 1, ɛ = a i sq G iɛ,.18 i=1 one finds that G i can e decomposed in terms of collinear B q i and soft f q i anomalous dimensions as follows [ 5] G iɛ = B q i γ i + f q i + C i + k=1 ɛ k g,k i..1 5

7 The collinear anomalous dimensions B q i [18] are given y B q 1 = C F, CF B q = 1 4ζ + 48ζ + C A C F ζ 4ζ + n f T F C F B q 1 = 16C A C F 8 ζ 81 7 ζ + 7 ζ 5 ζ ζ ζ + 15 ζ C A C F n f 0 ζ C F 5 ζ ζ ζ + 8 ζ ζ 15ζ C F n f 0 ζ 5 1 ζ ζ + 16 the soft anomalous dimensions f q i [] are + 16C A C F 47 4 ζ, 60 ζ + ζ ζ 54 ζ ζ C F n f 5 7 ζ + 1 ζ , f q 1 = 0, f q = C AC F ζ 8ζ C F n f T F ζ 4, 7 f q 5 = C A C F 5 ζ ζ ζ ζ 116 ζ + 1ζ C A C F n f 6 5 ζ ζ ζ C F n f 5 ζ + 4ζ + 04 ζ C F n f 40 7 ζ ζ and the constants C i are given y C1 = 0, C = β 0 g,1 1, C = β 1 g,1 1 β 0 g,1 + β 0 g, 1..1 Since G 1 ɛ and G ɛ are known to all orders in ɛ and G ɛ is known to Oɛ [6], the coefficients g,k i for i = 1,, can e readily otained from G i ɛ. The relevant one loop terms are found to e g,1 1 = C F + ζ, g, 1 = C F 7 ζ, g, 1 = C F ζ ζ,. the relevant two loop terms [10, ] are g,1 = C F n f 616 g, ζ 8 ζ + C F C A 8 + ζ 88 5 ζ 60ζ, + CF 7 = C F n f 1 ζ 55 7 ζ ζ C A C F 4 7 ζ 51ζ CF ζ ζ + 15 ζ 65 4 ζ + 8 ζ ζ ζ 6 5 ζ 8ζ ζ 44ζ + 116ζ + 1ζ 5 4 6

8 and finally the relevant three loop term [6] is g,1 = CAC F ζ + 78 ζ + 76 ζ ζ ζ 116 ζ ζ + 18 ζ C A CF ζ ζ 584 ζ ζ ζ + 6ζ ζ 140 ζ C A C F n f ζ + ζ ζ ζ ζ 7ζ C F 77 n f ζ 15 ζ ζ ζ ζ 68 ζ C F n f 40 4 ζ 8 81 ζ ζ 75 + CF ζ ζ + 64ζ ζ 75ζ + 48ζ 14ζ + 17ζ The mass factorization kernel Γ z, µ F, ɛ removes the collinear singularities resulting from massless partons. It satisfies the following RG equation µ F d dµ F Γ z, µ F, ɛ = 1 where P c z, µ F are Altarelli-Parisi splitting functions. expand them as P c z, µ F = P c z, µ F Γc z, µ F, ɛ,.4 c i=1 In perturative QCD, we can a i sµ F P i 1 c z..5 The off diagonal splitting functions are regular as z 1. On the other hand diagonal ones contain δ1 z and D 0 as well as regular terms i.e. P i z = B i+1δ1 z + A i+1d 0 + P i reg, z..6 We find that the regular part of the splitting function, P i reg,, does not contriute to threshold corrections. The fact that sv is finite as ɛ 0 implies that soft distriution function Φ â s, q, µ, z, ɛ also satisfies Sudakov type differential equations [] namely q d dq Φ = 1 [ K â s, µ R µ, z, ɛ + G â s, q µ, µ ] R R µ, z, ɛ. K and G take the forms similar to those of K and G of the form factors. This guarantees Ψ is finite order y order in perturation theory as ɛ 0. The solution to the aove equation is found to e q Φ = â i 1 z i ɛ s S i iɛ µ ɛ ˆφ,i ɛ.7 1 z i=1 7

9 where ˆφ,i ɛ = 1 iɛ [ K,i ɛ + G,i ɛ ].8 and the constants K can e otained using µ dk R dµ R = δ1 zµ dk R dµ.. R This implies that K,i ɛ can e written in terms of A q i and β i. Defining G iɛ through i=1 â i s q i ɛ z S i µ ɛ G,i ɛ = where qz = q 1 z and using the fact that sv as G 1 ɛ = f q 1 + k=1 ɛ k G,k 1, G ɛ = f q β 0G,1 1 + k=1 ɛ k G,k, a i sqzg iɛ.0 i=1 G ɛ = f q β 1G,1 1 β 0 G,1 + β 0 G, 1 + is finite as ɛ 0, we can express G iɛ k=1 ɛ k G,k..1 The constants G,j i arise from the soft part of the partonic reactions. Since the soft part does not depend on the hard process, G iɛ can e directly otained from G q i ɛ that contriutes to Drell-Yan production [] G iɛ = G q i ɛ.. We list the relevant ones that contriute up to N LO level. The terms required for one and two loops [] are G q,1 1 = C F ζ, G q, 1 = C F 7 ζ, G q, 1 = C F 16 ζ, G q,1 = C F n f ζ + ζ + C A C F ζ + 4ζ 176 G q, 11 = C A C F 40 ζ 0 ζ ζ ζ ζ + 4ζ C F n f 1 0 ζ 16 7 ζ 10 7 ζ ζ, 8

10 and for three loops [4] G q,1 15 = C F C A 6 ζ ζ ζ Results ζ ζ ζ ζ ζ 5648 ζ 7 + C A n f 5 ζ 108 ζ ζ ζ ζ C F n f 15 ζ 88 ζ ζ ζ ζ n f ζ ζ ζ ζ 5 Expanding eqn..5 in powers of a s µ R using eqn.6 and performing the convolutions we find sv z = i=0 a i sµ R sv,i z, µ R.1 where sv,i can e expressed in terms of the distriutions δ1 z and D i : sv,i = sv,i i 1 µ R δ δ1 z + sv,i µ R Dj D j.. The results up to N LO are given elow sv,1 = δ1 z G,1 1 + g,1 1 + ζ A 1 D 0 f1 + 4D 1 A 1, sv, = δ1 z G,1 + G,1 1 + g,1 + 4g,1 1 G,1 1 + g,1 1 + β0 G, 1 + g, 1 8ζ f1a 1 + ζ A f1,1 + 6G 1 A 1 + 6g,1 1 A 1 + ζ β 0 6γ0 + 6B1 + f ζ A 1 + D 0 4β 0 G, ζ A1 + ζ f1a 1 f 4G,1 1 f1 4g,1 1 f 1 + D 1 4β 0 f1 4ζ A1 + 4A + 4 f1,1 + 8G 1 A 1 + 8g,1 1 A 1 + D 4β 0 A 1 1f1A 1 + 8D A 1, sv, = δ1 z G,1 + G,1 1 G,1 + 4 G,1 1 + g,1 + g,1 G,1 1 + g,1 1 G,1 + 4g,1 1 G,1 1 + g,1 1 g,1 + 4g,1,1 1 G g, β 1 G, 1 + g, β 0 G, + G,1 1 G, g, + g, 1 G,1 1 + g,1 1 G, 1 + g,1 1 g, β 0 G, 1 + g, 1 ζ 5 f1 + β 0 A 1 8ζ fa 1 + f1a + 1 f1,1 + G 1 f 1A 1 + g,1 1 f 1A 1 8ζ β 0 f1,1 + G 1 A ζ A 1 + ζ A 4f1f + G,1 A 1 + 6G,1 1 A 4G,1 1 f1,1 + 6G 1 A 1 + g,1 A 1 + 6g,1 1 A j=0

11 4g,1 1 f1 + 1g,1 1 G,1 1 A 1 + 6g,1 1 A 1 ζ β 1 γ 0 B1 f1 + 6ζ β 0 γ1 + B + f + G, 1 A 1 G,1 1 γ0 + G,1 1 B1 1 G,1 1 f1 + g, 1 A 1 g,1 1 γ 0 + g,1 1 B 1 + g,1 1 f 1 1ζ β0g, ζ ζ f1 A 1 + ζ ζ β 0 A ζ A 1A 8 A f G,1 1 A 1 + g,1 1 A 1 + ζβ 0 18A 1B1 18A 1γ0 + f1a 1 ζβ 0A ζ A 1 + D 0 4β 1 G,1 1 4β 0 G,1 + G, 1 f1 + G,1 1 + g, 1 f 1 + g,1 1 G,1 1 8β0G, 1 + 1ζ 5 A 1 + ζ A 1A A + f1 1 + G,1 1 A 1 + g,1 1 A ζ β 0 f1a 1 + ζ fa 1 + f1a + f1,1 + G 1 f1a 1 + g,1 1 f 1A 1 + ζ β 0 6f1γ 0 6f1B 1 + f1 + G,1 1 A 1 80ζ ζ A1 + ζ f1 A ζ β 0 A 1 f + G,1 f1 + G,1 1 f + G,1 1 f 1 + g,1 f 1 + g,1 1 f + 4g,1 1 G,1 1 f1 + g,1 1 f 1 + D 1 4β 1 f1 + 8β 0 f + G, 1 A 1 + G,1 1 f1 + g, 1 A 1 + g,1 1 f β0G, ζ f1 A 1 6ζ β 0 A 1 8ζ A 1A + 7 A f1 1 + G,1 1 A 1 + g,1 1 A 1 4ζ β 0 A 1γ0 A 1B1 + f1a 1 46ζ A A + f1f + G,1 A 1 + G,1 1 A + G,1 1 f1,1 + G 1 A 1 + g,1 A 1 + g,1 1 A + g,1 1 f1 + 4g,1 1 G,1 1 A 1 + g,1 1 A 1 + D 4β 1 A 1 8β 0 A + f1,1 + 4G 1 A 1 + g,1 1 A 1 8β0f ζ A1 + 60ζ f1 A 1 + 6ζ β 0 A 1 1 fa 1 + f1a + 1 f1,1 + G 1 f 1A 1 + g, f 1A 1 + D β 0f1A β 0A 1 A 40ζ A A 1A + f1 1 + G,1 1 A 1 + g,1 1 A 1 + D 4 40 β 0 A 1 0f 1 + 8D 5 A1.. A 1 In the aove equation A i = Aq i, B i = Bq i, f i = f q i and we have set µ R = µ F = q. The finite sv,i Q depend on the anomalous dimensions A i, B i, f i and γi, the β function coefficients β i and ɛ expansion coefficients of G ɛ, g,i j s and of the corresponding G ɛ, G,i j s. The results for the sv,1 and sv, agree with those from the exact NLO and NNLO results [10]. At N LO level, only sv, Di s were known [] as the terms g,, g,1 of the form factor and G q,, G q,1 needed for sv, δ were not availale. The recent results for g, and g,1 from [6], G q, from [0] and the N LO G q,1 from [4] can e used to otain the missing δ1 z part namely sv, δ. This completes the evaluation of full threshold N LO contriutions in QCD for Higgs production through ottom anti- 10

12 ottom annihilation at hadron colliders. Below we present our results up to N LO level after sustituting A i, B i, f i, γ i and β i terms as well as the constants g,i j from the form factors and G,i j from the soft distriutions function in eqn... sv, = δ1 z 164 C A C F 15 ζ ζ 1064 ζ ζ 7ζ 400 ζ ζ 84ζ C A C F ζ 1040 ζ 7144 ζ ζ ζ 8ζ ζ ζ 1188ζ + 848ζ ζ ζ ζ D 0 15 ζ ζ ζ ζ ζ + C A C F n f ζ + 08 ζ ζ 15 ζ C F ζ 64ζ ζ C F n f 15 ζ C F n f 7 ζ 81 ζ ζ ζ ζ C A C F 15 ζ 5 ζ ζ ζ ζ 84ζ C A C F 1408 ζ 147ζ ζ ζ 480 ζ C F n f ζ C A C F n f ζ 7 ζ + 88 ζ C F 6144ζ ζ 104ζ + 188ζ C F n f 7 ζ ζ 578 ζ D 1 C A C F 7 5 ζ 10 ζ 704ζ C A C F 5 ζ 0864 ζ 5760ζ C A C F n f + C F 1408 ζ + 104ζ 60ζ D + C F n f ζ + 144ζ ζ + C A C F n f + D 4 + C A C F n f C F n f D 7 ζ 84ζ C F n f ζ + 180ζ C A C F ζ ζ C A C F ζ C F 1040ζ + C F n f C A C F C F 07ζ 51 + C F n f 7 C F n f C A C F D 5 51C F C A C F 51ζ C F n f 7..4 The numerical impact of our results can e studied using the exact LO, NLO, NNLO 11

13 0.7 LO NLO NNLO N LO SV 1.7 LO NLO NNLO N LO SV 0.6 σ [p] 0.5 σ [p] µ = M H /4 F µ = M R H µ /m R H µ /m F H Figure 1. Total cross section for Higgs production in annihilation at various orders in a s as a function of µ R /m H left panel and of µ F /m H right panel at the LHC with s = 14 TeV.,i ac, i = 0, 1, and the threshold N LO result sv,. We have used s = 14 TeV for the LHC, the Z oson mass M Z = GeV and Higgs oson mass m H = 15.5 GeV throughout. The strong coupling constant α s µ R is evolved using the 4-loop RG equations with α N LO s m Z = and for parton density sets we use MSTW 008NNLO [6]. The Yukawa coupling is evolved using 4 loop RG with λm = m m /ν and m m = 4. GeV. The renormalization scale dependence is studied y varying µ R etween 0.1 m H and 10 m H keeping µ F = m H /4 fixed. For the factorization scale, we have fixed µ R = m H and varied µ F etween 0.1 m H and 10 m H. We find that the perturation theory ehaves etter if we include more and more higher order terms see Fig.1. To summarize, we have systematically developed a framework to compute threshold contriutions in QCD to the production of Higgs oson in ottom anti-ottom annihilation suprocesses at the hadron colliders. Factorization of UV, soft and collinear singularities and exponentiation of their sum allow us to otain threshold corrections order y order in perturation theory. Using the recently otained N LO soft distriution function for Drell-Yan production and the three loop Higgs form factor with ottom anti-ottom quarks, we have otained threshold N LO corrections to Higgs production through ottom antiottom annihilation. We have also studied the staility of our result under renormalization and factorization scales. Acknowledgments We sincerely thank T. Gehrmann for constant encouragement and fruitful discussions. We also thank M. K. Mandal for useful discussions. The work of TA and NR has een partially supported y funding from RECAPP, Department of Atomic Energy, Govt. of India. 1

14 References [1] G. Aad et al. [ATLAS Collaoration], Phys. Lett. B [] S. Chatrchyan et al. [CMS Collaoration], Phys. Lett. B [] R. Barate et al. [ LEP Working Group for Higgs oson searches and ALEPH and DELPHI and L and OPAL Collaorations ], Phys. Lett. B [4] CDF and D0 Collaorations, Phys. Rev. Lett [5] ALEPH,CDF,D0, DELPHI,L,OPAL,SLD Collaorations, CERN PH-EP ,010. [6] H. M. Georgi, S. L. Glashow, M. E. Machacek, and D. V. Nanopoulos, Phys. Rev. Lett ; A. Djouadi, M. Spira, and P. M. Zerwas, Phys. Lett. B ; S. Dawson, Nucl. Phys. B ; M. Spira, A. Djouadi, D. Graudenz, and P. M. Zerwas, Nucl. Phys. B ; S. Catani, D. de Florian and M. Grazzini, JHEP ; R. V. Harlander and W. B. Kilgore, Phys. Rev. D , Phys. Rev. Lett ; C. Anastasiou and K. Melnikov, Nucl. Phys. B ; V. Ravindran, J. Smith, and W. L. van Neerven, Nucl. Phys. B [7] P. Bolzoni, F. Maltoni, S. Moch, M. Zaro, Phys. Rev. Lett [8] T. Han and S. Willenrock, Phys. Lett. B ; O. Brein, A. Djouadi and R. Harlander, Phys. Lett. B [] D. A. Dicus and S. Willenrock, Phys. Rev. D ; D. Dicus, T. Stelzer, Z. Sullivan and S. Willenrock, Phys. Rev. D ; F. Maltoni, Z. Sullivan and S. Willenrock, Phys. Rev. D ; F. I. Olness and W. -K. Tung, Nucl. Phys. B ; J. F. Gunion, H. E. Haer, F. E. Paige, W. -K. Tung and S. S. D. Willenrock, Nucl. Phys. B [10] R. V. Harlander and W. B. Kilgore, Phys. Rev. D [11] L. Reina and S. Dawson, Phys. Rev. Lett ; W. Beenakker, S. Dittmaier, M. Kramer, B. Plumper, M. Spira and P. M. Zerwas, Phys. Rev. Lett ; S. Dawson, L. H. Orr, L. Reina and D. Wackeroth, Phys. Rev. D ; W. Beenakker, S. Dittmaier, M. Kramer, B. Plumper, M. Spira and P. M. Zerwas, Nucl. Phys. B ; R. Raitio and W. W. Wada, Phys. Rev. D ; Z. Kunszt, Nucl. Phys. B [1] S. Catani, D. de Florian, M. Grazzini, and P. Nason, JHEP [1] U. Aglietti, R. Bonciani, G. Degrassi and A. Vicini, Phys. Lett. B ; S. Actis, G. Passarino, C. Sturm and S. Uccirati, Phys. Lett. B G. Altarelli, R. K. Ellis and G. Martinelli, Nucl. Phys. B 14, [Erratum-iid. B 146, ] ; T. Matsuura and W. L. van Neerven, Z. Phys. C 8, ; T. Matsuura, S. C. van der Marck and W. L. van Neerven, Nucl. Phys. B 1, ; R. Hamerg, W. L. van Neerven and T. Matsuura, Nucl. Phys. B 5, 4 11 [Erratum-iid. B 644, 40 00]. [14] S. Moch, J. A. M. Vermaseren and A. Vogt, JHEP ; Phys. Lett. B [15] P. A. Baikov, K. G. Chetyrkin, A. V. Smirnov, V. A. Smirnov and M. Steinhauser, Phys. Rev. Lett [16] T. Gehrmann, E. W. N. Glover, T. Huer, N. Ikizlerli and C. Studerus, JHEP

15 [17] T. Gehrmann, T. Huer and D. Maitre, Phys. Lett. B [18] S. Moch, J. A. M. Vermaseren and A. Vogt, Nucl. Phys. B 688, ; Nucl. Phys. B 61, [1] K. G. Chetyrkin, B. A. Kniehl and M. Steinhauser, Nucl. Phys. B 510, [0] D. de Florian and J. Mazzitelli, JHEP [1] S. Moch and A. Vogt, Phys. Lett. B ; E. Laenen and L. Magnea, Phys. Lett. B ; A. Idili, X. -d. Ji, J. -P. Ma and F. Yuan, Phys. Rev. D [] V. Ravindran, Nucl. Phys. B ; Nucl. Phys. B [] C. Anastasiou, C. Duhr, F. Dulat, E. Furlan, T. Gehrmann, F. Herzog and B. Mistlerger, arxiv: [hep-ph]. [4] T. Ahmed, M. Mahakhud, N. Rana and V. Ravindran, arxiv: [hep-ph]. [5] Y. Li, A. von Manteuffel, R. M. Schainger and H. X. Zhu, arxiv: [hep-ph]; S. Catani, L. Cieri, D. de Florian, G. Ferrera and M. Grazzini, arxiv: [hep-ph]. [6] T. Gehrmann and D. Kara, arxiv: [hep-ph]. [7] T. van Ritergen, J. A. M. Vermaseren and S. A. Larin, Phys. Lett. B ; M. Czakon, Nucl. Phys. B [8] V. V. Sudakov, Sov. Phys. JETP [Zh. Eksp. Teor. Fiz ]. [] A. H. Mueller, Phys. Rev. D [0] J. C. Collins, Phys. Rev. D [1] A. Sen, Phys. Rev. D [] S. Catani and L. Trentadue, Nucl. Phys. B 7 18 ; Nucl. Phys. B ; A. Vogt, Phys. Lett. B ; [] V. Ravindran, J. Smith and W. L. van Neerven, Nucl. Phys. B ; [4] T. Becher and M. Neuert, Phys. Rev. Lett ; [5] E. Gardi and L. Magnea, JHEP [6] A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Eur. Phys. J. C