Hadron energy spectrum in polarized top-quark decays considering the effects of hadron and bottom quark masses

We present the analytical expressions for the next-to-leading order corrections to the partial decay width t(↑)→bW+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t(\uparrow ) \rightarrow bW^+$$\end{document}, followed by b→HbX\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b\rightarrow H_bX$$\end{document}, for nonzero b-quark mass (mb≠0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_b\ne 0$$\end{document}) in the fixed-flavor-number scheme (FFNs). To make the predictions for the energy distribution of outgoing hadrons Hb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_b$$\end{document}, as a function of the normalized Hb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_b$$\end{document}-energy fraction xH\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x_H$$\end{document}, we apply the general-mass variable-flavor-number scheme (GM-VFNs) in a specific helicity coordinate system where the polarization of top quark is evaluated relative to the b-quark momentum. We also study the effects of gluon fragmentation and finite hadron mass on the hadron energy spectrum so that hadron masses are responsible for the low-xH\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x_H$$\end{document} threshold. In order to describe both the b-quark and the gluon hadronizations in top decays we apply realistic and nonperturbative fragmentation functions extracted through a global fit to the e+e-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^+e^-$$\end{document} annihilation data from CERN LEP1 and SLAC SLC by relying on their universality and scaling violations.


Introduction
Ever since the top-quark discovery in pp collisions in 1995 by the CDF [1] and D0 experiments [2] at the Fermilab Tevatron, it has been in or near the center of attention in high-energy physics. Its characteristics such as its mass m t , total decay width t , branching fractions, and elements V tq of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix, have not yet been determined precisely. Its property which is likely most central in many aspects of top physics is its mass. Since, for example, the top interacts with the Higgs boson through the potential y t htt where y t is proportional to the top mass m t , therefore the top mass has a fundamental role in the issue of stability of the Higgs potential (see Eq. (5)). Using the full sample of pp collision data collected a e-mails: moosavi@ipm.ir;mmoosavi@yazd.ac.ir by the D0 experiment in the Tevatron Run II [3] at a centerof-mass energy of 1.96 TeV and for an integrated luminosity of up to 8.7 fb −1 the mass of top quark is measured as m t = 174.98 ± 0.76 GeV. Due to its remarkably large mass, the top couples strongly to the agents of electroweak symmetry breaking and this makes it both an object of interest itself, and a tool to investigate that mechanism in detail.
Among other things, the CERN LHC is a genuine top factory, in particular in Run II, producing about 90 million top-quark pairs per year of running at design energy 14 TeV. The existing and upcoming data will allow us to study the top quark and its behavior in LHC collisions in great detail, if also the theoretical descriptions and simulations are of proportionate quality.
The top decay width itself is very difficult to determine in hadron colliders, though a recent experimental inference of the width was performed by D0 [4] and found t = 2.00 ± 0.47 GeV in the context of single top t-channel production, and CDF [5] also reported 1.10 < t < 4.05 GeV at the 68 % confidence level. The top decay characteristics play an important role in studying the top quark at colliders. Since the width-to-mass ratio t /m t of the top quark is small enough then, for many purposes, the notion of top quark as a stable particle makes sense, so that its production and decay processes can be factorized through the narrow width approximation [6]. In fact, if it were not for the confinement of color, the top could be considered as a free particle. This property allows it to behave like a real particle and one can safely describe its decay in perturbative theory.
The top decay width is largely due to decays to a Wboson and a bottom quark (with the mass m b ), as it is represented in the element |V tb | ≈ 1 of the CKM matrix [7,8]. Since m t m W + m b , its width is sufficiently large to pre-empt top-quark hadronization, then this rapid decay (τ t ≈ 0.5 × 10 −24 s [9]) enables transmission of top-quark spin information to final states. Part of the top attractiveness is its power to self-analyze its spin, through its purely left-handed standard model (SM) weak decay. The interplay between the top spin and its mass is of crucial importance in studying the SM. The top quark polarization can be studied by the angular correlations between the top-quark spin and its decay product momenta, and these spin-momentum correlations will allow for detailed studies of the top decay mechanism. In [10,11], we showed that these correlations depend on the choice of the possible helicity coordinate systems.
Since produced b-quarks hadronize before they decay, then each b-jet contains a bottom-flavored hadron which, most of the times, is a B-meson. They are identified by a displaced decay vertex associated which charged lepton tracks. These nonperturbative transitions are described by realistic, nonperturbative fragmentation functions (FFs) that are usually obtained through a global fit to e + e − data. At LHC, the decay process t → BW + + X is of prime importance, and it is an urgent task to predict its partial decay width as reliably as possible, specifically the distribution in the scaledenergy of B-mesons (x B ) in the top quark rest frame is of particular interest. In [10,11], using the zero-mass variableflavor-number scheme (ZM-VFNs) in which the mass of bquarks are set to zero at the parton level, we studied the NLO angular distribution of the scaled-energy of B-hadrons through polarized top decays. For that, we calculated the polar angular correlation in the rest frame decay of a polarized top quark into a stable W + -boson and the B-hadron, i.e., t (↑) → W + + b(→ B + X ). We analyzed this correlation in a helicity coordinate system where the event plane, including the top quark and its decay products, is defined in the (x, z)-plane with the Z-axes along the b-quark momentum. Here, the top polarization vector was evaluated with respect to the b-momentum direction.
Here, using the same frame, we revisit B-hadron production from polarized top decays by working at NLO in the general-mass variable-flavor-number scheme (GM-VFNs), where b-quark masses are preserved from the beginning. This makes the calculations more complicated. Being manifestly based on Collin's QCD factorization theorem [12] convenient for massive quarks, this factorization scheme allows us to resum the large logarithms in m b , to retain the finite-m b effects and to preserve the universality of the FFs, whose scaling violations remain to be subject to DGLAP evolution [13][14][15][16][17]. In this way, it combines the virtues of the fixed-flavor-number scheme (FFNs) and the ZM-VFN scheme and also avoids their flaws. In fact, it is an elaborate tool for global analyses of experimental data on the inclusive production of heavy flavored hadrons, allowing one to transfer nonperturbative information on the hadronization of partons from one type of experiment to another and from one energy scale μ F to another, without the restriction μ F m b , which is essential for the ZM-VFNs. Our analysis is supposed to enhance our previous result [10,11] in the ZM-VFN scheme by retaining all nonlogarithmic m bterms.
Moreover, we also include finite-m B effects, which modify the relations between hadronic and partonic scaling variables and reduce the available phase space. In comparison with the finite-m b effects, these effects turn out to be moderate. Their study is necessary to fully exploit the enormous statistics of the LHC data to be taken in the long run for a high-precision determination of the top properties.
This paper is organized as follows. In Sect. 2, we provide a compact review of some of the key aspects of top-quark physics, largely from a SM point of view. In Sect. 3 our formalism is explained by studying the real and virtual corrections. In Sect. 5, we explain the effects of finite-m B corrections in the evaluation of energy spectrum of B-mesons. Our numerical analysis is presented in Sect. 6 and in Sect. 7, we summarize our conclusions.

Top quark in the standard model
At first, we briefly review the various interactions of the top quark field t (x μ ) in the SM Lagrangian; a topic needed for the calculation of top decay widths.
The charged weak interaction of the top quark is lefthanded and flavor-changing, so expressed as where f (x) stands for the fields of down, strange and bottom quarks and the weak coupling factor g w is related to the Fermi coupling constant as g 2 w = 4 √ 2m 2 W G F , while its neutral weak interaction is flavor-conserving and parity violating where θ W is the Weinberg mixing angle, so that sin 2 θ W = 0.23124 [18]. Its interaction with gluons is a vector-like coupling, involving an SU(3) generator (T a ) in the fundamental representation where g s is the strong coupling constant, a = 1, 2, . . . , 8, is the QCD color index so T r(T a T a )/3 = C F . The top interaction with photons is also simply vector-like as that is, proportional to the top-quark electric charge. Finally, the interaction of the top quark with the Higgs field h(x) is of the Yukawa type, with a coupling constant y t = √ 2m t /v, where v is the Higgs vacuum expectation value. The Yukawa coupling is almost y t ≈ 1 in the SM.
In many extensions of the SM such as minimal supersymmetric standard model (MSSM), the Higgs sector of the SM is enlarged by considering an extra doublet of complex Higgs field [19][20][21][22][23]. In [24,25], we studied the top decay in the general two Higgs doublet model (2HDM).
Moreover, beyond the interactions above, effective interactions such as for flavor-changing neutral currents occur due to loop corrections. However, they are generally very small in comparison with those above. All these interactions could be modified in structure and strength by virtual effects due to new interactions associated with the physics beyond the SM. This is of interest to investigate that if the top quark, evidently, has a large coupling to the electroweak symmetry breaking sector. Therefore, it is so important to test these structures in detail, and indeed this is the thrust behind the field of top physics.
One proposed way to study the properties of top quarks is to consider the scaled-energy distribution of outgoing hadrons. In the next section we shall study this approach in detail, using the GM-VFN scheme where the mass of b-quark is preserved from the beginning.

Formalism
We consider the decay process of polarized on-shell top quark at NLO, as where X stands for the unobserved final-state particles. We wish to study the angular distribution of the scaled energy of B-hadrons by considering the contribution of bottom and gluon fragmentations into the heavy meson B, so that the gluon contributes to the real radiation at NLO. To obtain this energy spectrum, we need to have the parton-level differential width of the process (6). The LO contribution results from t (↑) → bW + . We define the partonic scaled-energy fraction where p i stands for the gluon or bottom quark momenta. In the top-quark rest frame where p t = (m t , 0), one has x i = 2E i /m t where E i refers to the energy of outgoing partons; gluon or bottom at NLO. By preserving the bottom quark mass, one has As in [26], throughout this paper, we shall make use of the normalized energy fraction of the partons: The allowed values of x b and x g shall be discussed in Sect. 7. We analyse the decay t (↑) → bW + in the rest frame of the top quark where the 3-momentum of the b-quark points to the direction of the positive Z-axis. For a polarized top quark, the general angular distribution of differential decay width d /dx is given by This form clarifies the correlations between the top decay products and the spin of the top quark.
In (8), P is the magnitude of the top-quark polarization with 0 ≤ P ≤ 1, so that P = 0 corresponds to an unpolarized top quark and P = 1 is for the 100 % polarization. Here, θ P is defined as the polar angle between the top-quark polarization vector P and the Z-axis (b-quark momentum direction).
In (8), d˜ unpol /dx i refers to the unpolarized differential widths which studied in [27], both in ZM-and GM-VFN schemes. In the following, we discuss the evaluation of the quantities d˜ pol /dx i in the GM-VFN scheme.

Born term result
Considering the charged weak interaction Lagrangian (1), the dynamics of the current and at the Born level and O(α s ) one-loop contributions the intermediate state is |X b = |b . Also, this hadronic tensor depends on the top spin s t . It is straightforward to compute the Born term contribution to the decay (6). In the top rest frame, the fourmomentum of the bottom quark is set to p b = (E b ; 0, 0, p b ) and the polarization four-vector of the top quark is set as s t = P(0; sin θ P cos φ P , sin θ P sin φ P , cos θ P ). Considering the general distribution (8), the Born term helicity structure of the partial rates reads where the LO polarized (˜ (0) B ) and unpolarized (˜ (0) A ) total decay widths read Here, we used the following kinematic variables, in the notations of Ref. [26]: In the limit of vanishing bottom quark mass, the tree-level decay widths converted to our results in [10,11] and [27], respectively.

Virtual corrections and counterterms
The QCD one-loop vertex corrections arise from the emission and absorption of the virtual gluons, so an interaction Lagrangian as in (3) is needed to calculate the virtual radiative corrections. Here, we adopt the on-shell mass renormalization scheme and use dimensional regularization to regulate the ultraviolet (UV) and soft singularities which appear in one-loop corrections. For example, the UV-singularities appear when the integration region of the internal momentum of the virtual gluon goes to infinity. The singularities are regularized by dimensional regularization in D = 4 − 2 space-time dimensions to become single poles in , so that 0 < ≤ 1. In the massless case, all singularities are subtracted at factorization scale μ F and absorbed into the bare FFs in accordance with the modified minimal subtraction (M S) scheme; see [10,11]. In the massive case, all singularities are automatically canceled after summing all radiative corrections up. Considering the notations (11), the contribution of virtual corrections into the doubly differential decay width (8) is obtained: where M 0 stands for the Born term amplitude and the renormalized amplitude M 1−loop refers to the virtual-gluon corrections, presented in [10,11]. The virtual contributions include the counterterm and the one-loop vertex corrections. The counterterm of the vertex contains the wave-function renormalization constants of the top (δ Z t ) and the bottom quark (δ Z b ). These constants can be found in [10,11]. The wave-function renormalization and the one-loop vertex correction contain the UV and infrared (IR) singularities so that all UV-divergences are canceled after summing all virtual corrections up and, from now on, we label the remaining IR-singularities by . Therefore, the virtual decay width is given by where the unpolarized differential decay rate reads and the polarized one is expressed by wherẽ In the equations above, is the dilogarithmic function (or Spence function) and the term F includes the IR-singularity ( ) as Here, μ F stands for the factorization scale which will be removed after summing all corrections up in the GM-VFN scheme.

Real-gluon corrections
The O(α s ) real graph contributions result from the real-gluon emissions from the bottom and top quarks, individually. In the rest frame of a top quark decaying into a b-quark, a W + boson and a gluon, the outgoing particles define an event plane.
Relative to this plane one can, then, define the spin direction of the polarized top quark. As in [10,11], here we apply a specific helicity coordinate system where the momenta of the b-quark and the W + boson are defined by p b ||ẑ, ( p W ) x ≥ 0, and the polarization vector of top quark is evaluated relative to theẑ-axis. In the following, we explain a brief technical detail of our calculation for the NLO radiative corrections to the tree-level decay rate of t (↑) → bW . In [10,11], where we set the mass of b-quark to zero from the first, the IR-singularities arose from the soft-and collinear gluon emissions. Since here we preserve the mass of b-quark then all IR-singularities arise from the soft real-gluon emission and the collinear divergences would be absent. As before, to regularize the IR-singularities we work in D dimensions where the real differential rate is given by To compute the differential rate d˜ real /dx b , we fix the bquark momentum and integrate over the gluon energy which In the GM-VFN scheme the real and virtual differential widths include the pole ∝ 1/ , which shall disappear in the total NLO result. Due to the radiation of a soft gluon (E g → 0) in top decay, on integration over the phase space for the real-gluon radiation, terms of the form (1 − x b ) −1−2 arise, which are divergent when x b → 1. Therefore, for a massive scheme (m b = 0) where β ≤ x b ≤ 1, we shall make use of the following expression [26]: with the plus prescription defined as

Parton-level results for angular distribution of partial decay rates in FFN scheme
Considering the tree-level, the real and virtual contributions, we present our analytic expression for the angular distribution of the partial decay rate in the FFN scheme. According to the Lee-Nauenberg theorem, after summing all corrections up the singularities cancel each other and the final result is free from IR-singularities. Therefore, the complete NLO results read where d˜ unpol N L O /dx b is given in [26,27], and d˜ pol N L O /dx b in the M S scheme is presented, for the first time, as where Since the B-mesons can also be produced through the fragmentation of emitted real gluons, then to obtain the most accurate result for the energy spectrum of mesons one needs the doubly differential distribution d 2˜ /(dx g d cos θ P ) where the scaled-variable x g is defined in (7). As we will show in Fig. 1, the contribution of the gluon fragmentation into the B-meson is negative and leads to a significant reduction in size in the threshold region, so that this contribution would be important at a low energy of the observed meson.
To get the d 2˜ /(dx g d cos θ P ), one has to integrate over the b-quark energy by fixing the gluon momentum in the phase space so that the b-quark energy ranges as m t Sφ  The d˜ unpol N L O /dx g in the FFN scheme is given in [27], and by considering the following notations: and 1+3b−4S x g , the polarized contribution reads

General-mass variable-flavor-number scheme
In [10,11], for obtaining the parton-level results for angular distribution of partial decay rates we used the ZM-VFN scheme, where m b = 0 was put right from the beginning and all collinear singularities were absorbed into the bare FFs according to the M S scheme. This approach renormalizes the FFs and produces finite terms of the form (α s /π ) ln(μ 2 F /m 2 t ) in the partial decay rates dˆ /dx a , which are rendered perturbatively small by choosing μ F = O(m t ). In this scheme, the b-quark mass m b just sets the initial scale μ ini F = O(m b ) of the DGLAP evolution equations, where ansaetze for the z-dependences of the FFs D a (z, μ ini F ) are injected by some proposed models [28]. The DGLAP evolution from μ ini F to a higher scale μ F then effectively resums the problematic logarithms (α s /π ) ln(m 2 t /m 2 b ) of the FFN scheme, however, all information on the m b -dependence of dˆ /dx a is wasted.
The GM-VFN scheme provides an ideal theoretical framework to study the effects of heavy quark masses, so it combines the virtues of the ZM-VFN and FFN schemes and, at the same time, avoids their flaws. In the GM-VFN scheme, the perturbative fragmentation functions enter the formalism via subtraction terms for the hard-scattering decay rates, so that the actual FFs are truly nonperturbative and may be assumed to have some smooth forms which can be specified through global data fits. In opposition with the FFN scheme, the GM-VFN scheme also accommodates FFs for light quarks and gluons, as in the ZM-VFN scheme. In our present work, the GM-VFNs is applied to resum the large logarithms in m b and to retain the entire nonlogarithmic m bdependence at the same time. This is reached by introducing convenient subtraction terms in the O(α s ) FFN expressions for d˜ /dx i , so that the O(α s ) ZM-VFN results are exactly recovered in the limit m b /m t → 0. These subtraction terms are universal and so are the FFs in the FFN scheme, as is guaranteed by Collin's hard-scattering factorization theorem [12].
As explained above, the GM-VFN results for the angular decay distributions are obtained by matching the FFN results (22,26) to the ZM-VFN ones [10,11] by subtraction, as where the subtraction terms are obtained: Taking the limit m b → 0 in (22) and (26), we recover the results presented in [10,11] up to the terms and As we have already shown in [27], for the unpolarized top decay in the SM, i.e. t → bW + , and also for the top decay in the theories beyond the SM including the two Higgs doublet where t → bH + [24,25], Eq. (29) coincides with the perturbative FF of the transition b → b. This is in agreement with the Collin factorization theorem which guarantees that the subtraction terms are universal. Thus the results presented in (29) and (30) ensure the correctness of our calculations shown in (22) and (26).

Hadron mass effects and Hadron level results
Our main purpose is to obtain the scaled-energy (x B ) distribution of bottom-flavored hadrons (B) inclusively produced in polarized top decays at NLO. Here, the scaled-energy fraction of B-hadrons is defined as 1+b−ω)), as in (7). In [10,11], to obtain the partial decay width of the process (6) in the ZM-VFN scheme, we used the Collin's factorization theorem [12]. According to this theorem, the energy distribution of B-hadrons might be expressed as the convolution of the partonic hard-scattering decay rates d /dx a , with the nonperturbative FFs which describe the transition a → B, as Here, μ F and μ R are the factorization and the renormalization scales, respectively. The μ R is associated with the renormalization of the strong coupling constant and a choice often applied is μ R = μ F . In the massless (or ZM-VFN) scheme where one sets m b = 0, the b-quark, gluon, and B-hadron (with the mass m B ) have energies 0 As we demonstrated in [27], the relation (31) is convenient for the case m b = 0 = m B . To calculate the d /dx B when passing from the ZM-VFN scheme to the GM-VFN scheme by taking into account the finite-m B corrections, one should apply the following improved relation: where (d /dx a ) GM−VFN is given in (27), and with , g, B), . Now, the kinematically allowed scaling variables are for gluon transition.
where n f is the number of active quark flavors and (36) Considering n f = 5, we adopt the asymptotic scale parameter (5) M S = 231 MeV, adjusted such that α (5) s (m Z ) = 0.1184 for m Z = 91.1876 GeV [29]. To include the B-meson and the b-quark masses, we apply Eq. (32) in which for the transitions (b, g) → B, from [30] we employ the related realistic and nonperturbative FFs. In [30], a power model as at the initial scale μ F = 4.5 GeV of fragmentation, while the FFs of gluon and light quarks are set to zero at the starting scale and are evolved to higher scales via the DGLAP equations [13][14][15][16][17]. The fit parameters are obtained at NLO in the ZM-VFN scheme through a global fit to e + e − annihilation data from the ALEPH and OPAL collaborations at CERN LEP1 and by SLD at SLAC SLC and the results are N = 4684.1, α = 16.87 and β = 2.628.
In Fig. 1, our predictions for the scaled-energy (x B ) distribution of B-hadrons are shown by considering the corresponding quantity d (t (↑) → B X)/dx B in the GM-VFN scheme. For this study, we considered the size of the NLO corrections, by comparing the LO (dot-dashed line) and NLO (solid line) results, and the relative importance of the b → B (dotted line) and g → B (dashed line) fragmentation channels at NLO. To compare the size of the NLO corrections at the parton level, we evaluate the LO result using the same NLO FFs. As is seen, the g → B contribution into the NLO energy spectrum of the B-meson is negative and appreciable only in the low-x B region and for higher values of x B , the NLO result is practically exhausted by the b → B contribution, as expected in [26]. In fact, the contribution of the gluon is evaluated to see where it contributes to d /dx B and cannot be discriminated in the meson spectrum as an experimental quantity. In the scaled-energy of mesons, all contributions including the bottom quark, gluon and light quarks contribute.
From Fig. 1, it is also seen that the NLO corrections lead to a significant enhancement of the partial decay width in the peak region and above, by as much as 30 %, at the expense of a depletion in the lower-x B range. Moreover, the peak position is shifted toward higher values of x B .
In (31), the factorization (μ F ) and the renormalization (μ R ) scales are arbitrary and, in principle, one can use two different values for them. However, a choice often made consists of setting μ R = μ F and we shall adopt this convention for most of the results which we shall show. In Fig. 2, we show the dependence of the meson energy spectrum on the factorization scales by considering μ F = m t (dot-dashed line), μ F = m t /2 (solid line) and μ F = 2m t (dotted line). In [26], the dependence of the x b spectrum on the factorization scales μ 0F and μ F are studied in detail. Their results  show that the x b dependence on the initial scales μ 0F is small when one resums soft logarithms in the initial condition of the b-quark perturbative FF. According to their results, as a whole conclusion, one can states that resumming soft logarithms yields a reduction of the theoretical uncertainty, as the dependence on factorization scales is indeed an estimate of effects of higher order contributions which we have been neglecting. In Fig. 3, we study the improvement over our previous predictions [10,11], at NLO in the ZM-VFN scheme by neglecting g → B fragmentation and finite-m B corrections (dashed For comparison, we also plot the full ZM-VFN prediction, including g → B fragmentation, for m B = 0 (solid line). Note that our most reliable prediction for the energy spectrum of B-meson is made at NLO in the GM-VFN scheme by including finite-m B corrections. It is observed from Fig. 3 that the improvement in our new work is twofold. First, the finite-m B corrections are responsible for the appearance of the threshold at x B = ρ B = 0.083, and second, the finitem b corrections lead to a moderate reduction in size throughout the whole x B range allowed, specially in the peak by as much as 7 %. This improvement is important in order to fully exploit the enormous statistics of the LHC data to be taken in the long run for a high-precision determination of the top-quark properties. For a more quantitative description of Fig. 3, we consider in Fig. 4 the GM-VFN result (m b = 0) including finitem B corrections without g → B fragmentation (solid line) and the full GM-VFN result (dot-dashed line) using Eq.(32), both normalized to the full ZM-VFN result for m B = 0. It is observed that the omission of g → B fragmentation causes an excess by a factor of up to 2 close to threshold, while the finite-m b and m B corrections amount less.
In Fig. 5, to study the angular dependence of energy distributions we plot the ratio of the NLO unpolarized [27] and polarized results for the d (t → B + X )/dx B in the GM-VFN (m b = 0) scheme including finite-m B corrections.

Conclusions
Studying the fundamental properties of the top quark is an object of interest in theoretical and experimental particle physics. Among other things, the LHC is a superlative top factory which allows one to study the top characteristics in great detail, if also the theoretical descriptions and simulations are of commensurate quality. In particular, the LHC will allow for the study of the dominant decay mode t → BW + + X with unprecedented precision in the long run. As an application, these studies will enable us to deepen our conception of the nonperturbative aspects of B-hadron formation by hadronization. The key quantity for this purpose is the x B distribution d /dx B of t → B + X . Therefore, the distributions in the scaled B-hadron energy x B through the polarized (6) or unpolarized [27] top decays are of particular interest at the LHC. The top-quark decays rapidly so that has no enough time to hadronize and then passes on its full spin information to its decay products. This allows one to study the top spin state using the angular distributions of its decay products, so that in this work we studied the O(α s ) spin-dependent energy spectrum of hadrons produced from polarized topquark decays. For this purpose, we studied the observable d /dx B at NLO in the GM-VFN scheme [30]. This allowed us to investigate, for the first time, finite-m b corrections to the d /dx B . We also analyzed the size of finite-m B effects. Specifically, our analysis is supposed to enhance our previous results presented in [10,11] by retaining all nonlogarithmic m b terms of the result in the FFN scheme; see Fig. 3.
In this work, we analyzed the size of finite-m B effects and their theoretical uncertainty due to the freedom in the choice of the scaling variable. As we showed, the effects of infinite-m B corrections and the light-cone-momentum scaling variable lead to very similar results over most of the x B range. Appreciable differences only appear close to the kinematical bounds.
Here, we also studied the uncertainties raised from the choice of factorization scales (Fig. 2), but there are still other uncertainties such as the finite width effects, higher order corrections and etc, which should be taken into the calculations. Nevertheless, within the remaining theoretical uncertainties one can assert that the results of GM-and ZM-VFN schemes are equivalent.
Comparing future measurements of the polarized width d (t (↑) → BW + + X )/dx B at the LHC with our NLO predictions, one will also be able to test the universality and scaling violations of the B-hadron FFs.
Our formalism elaborated here is also applicable to the production of hadron species other than B-hadrons, such as pions, kaons and protons, etc., using the b, g → π/K /P FFs presented in our recent paper [31][32][33], relying on their universality and scaling violations.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .