arxiv: v1 [hep-ph] 28 Aug 2013

Transcription

1 Prepared for submission to JHEP Light-Cone Distribution Amplitudes for Heavy-Quark Hadrons arxiv: v1 [hep-ph] 8 Aug 13 Guido Bell, a Thorsten Feldmann, b Yu-Ming Wang, c and Matthew W Y Yip d a Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford, UK b Theoretische Elementarteilchenphysik, Universität Siegen, Siegen, Germany c Institut für Theoretische Physik E, RWTH Aachen, Aachen, Germany; Physik Department T31, Technische Universität München, Garching, Germany d Institute for Particle Physics Phenomenology, Durham University, Durham, UK guido.bell@physics.ox.ac.uk, thorsten.feldmann@uni-siegen.de, yuming.wang@tum.de, m.w.yip@durham.ac.uk Abstract: We construct parametrizations of light-cone distribution amplitudes LCDAs for B-mesons and Λ b -baryons that obey various theoretical constraints, and which are simple to use in factorization theorems relevant for phenomenological applications in heavyflavour physics. In particular, we find the eigenfunctions of the Lange-Neubert renormalization kernel, which allow for a systematic implementation of renormalization-group evolution effects for both B-meson and Λ b -baryon decays. We also present a new strategy to construct LCDA models from momentum-space projectors, which can be used to implement Wandzura-Wilczek like relations, and which allow for a comparison with theoretical approaches that go beyond the collinear limit for the light-quark momenta in energetic heavy-hadron decays. Keywords: Heavy Quarks, Light-Cone Distribution Amplitudes, Renormalization 1 October 8, 18. Oxford: OUTP-13-15P; Siegen: SI-HEP-13-5, QFET-13-5; Aachen: TTK-13-19, SFB/CPP-13-6; Munich: TUM-HEP-899/13; Durham: IPPP/13/55, DCPT/13/11.

2 Contents 1 Motivation B Mesons 3.1 Light-Cone Distribution Amplitudes Light-Cone Projector for -Particle Fock State 3.1. Wandzura-Wilczek Approximation.1.3 Renormalization-Group Evolution 5.1. Application in Factorization Theorems 7. Construction from Momentum Space Particle Wave Functions 9.. Corrections to Wandzura-Wilczek Relation Higher Fock States 13 3 Λ b Baryons Light-Cone Distribution Amplitudes Light-Cone Projectors for 3-Particle Fock State The Chiral-Odd Projector M The Chiral-Even Projector M Wandzura-Wilczek Approximation Construction from Momentum Space 3..1 The Chiral-Odd Projector M 3.. The Chiral-Even Projector M Renormalization-Group Evolution Analytic Solution Numerical Examples and Asymptotic Form 7 Summary 31 A B-Meson Projectors in arbitrary Frame 3 A.1 -Particle Projector 3 A. 3-Particle Projector 33 B The ERBL Term for the dual LCDA of the Λ b Baryon 33 C Some Relations with Bessel Functions 35 1

3 1 Motivation Light-cone distribution amplitudes LCDAs play a key role in the factorization of shortand long-distance dynamics entering exclusive transition amplitudes in Quantum Chromodynamics QCD. Defined as hadronic matrix elements of composite QCD operators with fields separated along the light-cone, they represent genuine non-perturbative quantities, whose renormalization-scale dependence follows from the anomalous dimension of the defining operators which can be computed in perturbation theory. While early analyses focused on light hadrons pions, kaons, nucleons, etc. [1 ], with the running b-physics program at B-Factories or high-luminosity hadron colliders, the interest in exclusive decay channels of B-mesons and Λ b -baryons revealed the importance of LCDAs for hadrons containing a heavy quark. For instance, the former provide the dominant hadronic input for the theoretical description of radiative leptonic B-decays [5 1], they enter the QCD-factorization approach for non-leptonic and semi-leptonic B-decays [11 16], and they naturally appear in the context of soft-collinear effective theory SCET [17, 18]. They can also be used in correlation functions that form the basis of light-cone sum rules for heavy-to-light transition form factors [19, ], and they appear as limiting cases in the k T -factorization approach for heavy-meson [1, ] or heavy-baryon [3] decays. Theoretical properties of LCDAs have been classified for heavy mesons see e.g. [1,, 5] and heavy baryons see e.g. [6 8]. The renormalization-group RG evolution of the LCDAs has been derived in the heavy-quark limit, where the b-quark fields are treated in heavy-quark effective theory HQET, for both, mesons [9 35] and baryons [6]. Available model parametrizations are typically inspired by sum-rule analyses, for instance [3] for B-mesons or [6] for Λ b -baryons. Certain inverse moments of the B-meson LCDA can also be constrained from experimental data, most notably from the radiative leptonic decay B γlν see [9, 1] for recent analyses. In this article, our aim is to improve the theoretical modelling of LCDAs for B-mesons and Λ b -baryons in two aspects: First, we introduce a new representation for LCDAs in terms of a convolution of Bessel functions and a spectral function, which behaves in a dual way compared to the original LCDAs. This representation diagonalizes the Lange- Neubert LN evolution kernel, and therefore greatly simplifies the analytic solution of the RG equation. Our formalism also provides a systematic method to include the Efremov- Radyushkin Brodsky-Lepage ERBL kernel, which arises from gluon exchange among the light spectator quarks in the Λ b -baryon. Second, we develop a general procedure for an efficient modelling of LCDAs, starting from momentum-space projectors, which satisfy the equations of motion for the partonic Fock-state components. In this approach, the LCDAs follow as simple integrals over a set of fewer on-shell wave functions, which allows for a comparison between LCDAs used in the collinear factorization approach and transversemomentum-dependent wave-functions employed in the k T -factorization approach. For B- mesons, our formalism reproduces the so-called Wandzura-Wilczek WW relations established in [1], and for Λ b -baryons, we derive analogous relations between the various 3-particle LCDAs.

4 Our paper is organized as follows. In the first part, we focus on B-meson LCDAs for which we first recollect the main definitions and results, before we introduce our new representation in terms of dual spectral functions and give a detailed discussion on the RG properties. As a sample application of our formalism, we briefly reconsider the RGimproved factorization formula for the radiative leptonic B γlν decay. We proceed with the definition of on-shell momentum-space wave functions, and show how to construct the LCDAs and the corresponding momentum-space projectors to be used in applications of collinear QCD factorization. We also discuss how to incorporate corrections to the WW relations from 3-particle Fock states, and discuss two simple wave function models in detail. In the second part we present an analogous analysis for baryon LCDAs. Here, we show that the momentum-space construction together with some simplifying assumptions leads to an enormous reduction of independent hadronic functions. The RG equations for the baryon LCDAs are more complicated than for the meson LCDAs, because the LN kernel and the ERBL kernel depend differently on the two light-quark momenta. Nevertheless, we show that the RG equations can be solved systematically in the dual space in an expansion in Gegenbauer polynomials. This also establishes the general properties of leading-twist baryon LCDAs in the asymptotic limit. Our results are shortly summarized in Section. In the appendix we generalize our results for the B-meson momentum-space projectors to an arbitrary frame, where the heavy quark has a transverse velocity component. We also provide some details about the ERBL evolution of the baryon LCDAs, and collect some general integral relations for Bessel functions. B Mesons.1 Light-Cone Distribution Amplitudes We first recapitulate the properties of B-meson LCDAs in the heavy-quark limit, see e.g. [] and references therein..1.1 Light-Cone Projector for -Particle Fock State The -particle LCDAs of the B-meson in HQET which differ from the QCD definition, see [36, 37] can be obtained from the coordinate-space matrix elements q β z [z, ] h α v Bv = i f [ { B m B 1 + /v Φ + B t, z + Φ B t, z Φ + B t, z t /z } γ 5 ] αβ,.1 in the limit of light-like separation z. Here t = v z, and [z, ] denotes a gauge-link represented by a straight Wilson line along z µ. The heavy b-quark field in HQET is denoted as h v for a B-meson that moves with velocity v µ. The decay constant f B in the heavyquark limit includes the non-trivial scale dependence from the non-vanishing anomalous dimension of the local heavy-to-light current in HQET. 3

5 The terms in curly brackets can be expanded around z =, using the power-counting induced by the convolution with generic hard-scattering kernels in factorization theorems see the discussion in [1]. Introducing light-like vectors n µ ± with n n + = and v µ = nµ +nµ +, and taking n z z n + z, we obtain 1 /v 1 Φ + B t, z + Φ B t, z Φ + B t, z /z t φ + B τ /n + + φ B τ /n + φ B τ φ + B τ τ /z + Oz, n z,. where t τ = n +z can be interpreted as the Fourier-conjugated variable to the momentum component = n k associated with the light anti-quark field. Defining the LCDAs in momentum space from the Fourier transform, dτ φ ± B π eiτ φ± B τ,.3 the light-cone expansion in. corresponds to the momentum-space projector [1] M B v, = i f [ { B m B 1 + /v φ + B /n + + φ B /n dη ] φ B η φ+ B η γ µ }γ 5,. to be used in factorization theorems with hard-scattering kernels where the k and n +k- dependence can be neglected see below..1. Wandzura-Wilczek Approximation In the approximation where 3-particle contributions to the B-meson wave function are neglected, the equations of motion for the massless light-quark field yield WW relations between the -particle LCDAs [1], k µ d φ B τ + 1 dτ τ φ B τ φ + B τ = dη φ B η φ+ B η = φ B φ+ B = dφ B d..5 Including the contributions from 3-particle LCDAs as defined in.5 below, this is generalized to φ B dη φ B η φ+ B η dξ = dη η ξ ξ Ψ Aη, ξ Ψ V η, ξ..6 1 The generalization to frames where v µ n µ + + nµ / can be found in Appendix A.1.

6 .1.3 Renormalization-Group Evolution The RG evolution equation for the LCDA φ + B, µ reads with a slight change of notation compared to [9] dφ + B, µ [ = Γ cusp α s ln µ d ln µ + γ +α s ] φ + B, µ dη Γ +, η, α s φ + B η, µ,.7 where the leading terms in the various contributions to the anomalous dimensions in units of α s C F /π are [ ] Γ 1 cusp =, γ 1 θη θ η + =, Γ1 +, η = Γ1 cusp +,.8 ηη η + with the usual definition of the plus distribution [9]. The RG equation can be solved in closed form [3]. Starting from the Fourier transform 3 with respect to the variable ln /µ, ϕ + B θ, µ = d iθ φ + B µ, µ φ+ B, µ = dθ π one has an explicit solution of the RG equation, ϕ + B θ, µ = ev γ Eg µ µ where the RG functions are [3] V := V µ, µ = g := gµ, µ = iθ ϕ + B θ, µ,.9 µ iθ Γ1 iθ Γ1 + iθ g Γ1 + iθ Γ1 iθ + g ϕ+ B θ + ig, µ,.1 α sµ α sµ α sµ α sµ dα βα dα Γ cuspα βα Γ cusp α α α sµ dα βα + γ +α,..11 Transforming back to momentum space, the RG solution can be written as a convolution involving hypergeometric functions, φ + B, µ = ev γ Eg Γ g Γg min, η max, η F 1 dη η φ+ B η, µ 1 g, g,, max, η µ min, η max, η g,.1 which is valid for < g < 1 larger values of g do not appear in phenomenological applications and will not be considered further. The evolution kernel has been calculated at one-loop, and in [3] it has been conjectured that the structure of the evolution kernel remains a general feature at higher orders in the perturbative analysis. 3 This can also be viewed as Mellin moments N 1 + B for N = iθ. The corresponding formalism for the LCDA φ + B τ, µ in coordinate space has been worked out in [38]. 5

7 As one of the central new ideas of this paper, we suggest an alternative representation of the RG solution, which is obtained from the ansatz ϕ + Γ1 iθ d µ iθ B θ, µ := Γ1 + iθ ρ + B, µ..13 The particular parametrization for ϕ + B θ implies a simple RG behaviour for the spectral function ρ + B, ρ + B, µ = e V µ ˆ g ρ + B, µ = e V as a solution to the RG equation, Here we have defined dρ + B, µ d ln µ = g/ µµ ˆ ρ + B, µ,.1 [ Γ cusp α s ln µˆ ] + γ +α s ρ + B, µ..15 ˆ = e γ E, V µ, µ = 1 V µ, µ V µ, µ,.16 to write the solution in a manifestly symmetric form with respect to µ µ. The relation between the original LCDA φ + B and the spectral function ρ+ B then follows as 5 φ + B, µ = = = e V dθ Γ1 iθ π Γ1 + iθ d J 1 d J 1 d iθ ρ + B, µ ρ + B, µ µ ˆ g ρ + B, µ,.17 where J 1 is the Bessel function of the first kind. The expansion of the LCDA φ + B, µ in terms of Bessel functions and a spectral weight ρ + B, µ with simple RG properties can be viewed as the analogue of the Gegenbauer expansion for light mesons, where the Gegenbauer polynomials diagonalize the corresponding one-loop RG kernel. Therefore, instead of defining a model for the input function φ + B, µ, one can equivalently define a model for the dual spectral function ρ + B, µ at a given hadronic input scale, and determine φ B, µ at a different scale from a relatively simple convolution integral. Alternatively, one can express ρ + B, µ in terms of φ + B, µ as ρ + B, µ = d J 1 φ + B, µ,.18 such that the solution of the RG equation for φ + B, µ can also be obtained from a double convolution, φ + B, µ = ev d J 1 η dη η η µ J 1 ˆ g φ + B η, µ In the following, we implicitly assume that ρ + B, µ 1/ ɛ with ɛ > in the limit, such that the integrals in.17 converge for g ɛ < 1. 6

8 It is also interesting to note, how the spectral function is connected with the function φ + B τ appearing in the defining light-cone matrix elements. From.3 and.18 we obtain [ dτ ρ + B, µ = 1 exp i ] φ + π B τ, µ.. τ The fact that the product τ appears as the inverse in the exponential as compared to τ in the conventional Fourier transform in.3 justifies the notion dual function. It also illustrates, why the function ρ + B cannot be reconstructed through its positive moments related to a local operator product expansion around τ see also the discussion in [3, 3, 39]. The -particle operator that defines the LCDA φ + B, µ does not mix with the contributions from 3-particle operators under RG evolution [3, 35]. In contrast, the evolution of the LCDA φ B, µ contains a part that is independent of the 3-particle contributions and can be reconstructed from the WW relation [33], and a part that describes the explicit mixing with the combination Ψ A Ψ V of 3-particle LCDAs that enters the corrections to the WW relations in.6. In the WW approximation, the RG equation for the LCDA, µ can be solved in closed form [33]. In terms of our new strategy for the RG solution, φ B we find that the transformation with inverse transformation φ B, µ = ρ B, µ = d J ρ B, µ.1 d J φ B, µ. diagonalizes the corresponding one-loop renormalization kernel. The WW relation.5 implies a particularly simple relation between the spectral functions, and the solution for the LCDA φ B, µ takes the form φ B, µ = ev ρ B = ρ + B,.3 d µ g J ρ ˆ B, µ.. We will give explicit examples for φ ± B, and ρ± B from different models further below..1. Application in Factorization Theorems The most relevant parameter in phenomenological applications of the factorization approach is the first inverse moment of the LCDA φ + B. Interestingly, equation.19 implies as a consequence of the completeness relation C.1 for Fourier-Bessel transforms that this moment is identical to the first inverse moment of the corresponding spectral 7

9 function ρ + B. This result can be generalized to moments with one or two additional powers of ln, = e V d lnn d µ φ + n=,1, B, µ = ˆ ln n µ µ d ˆ ln n µ ρ + B, µ ˆ g ρ + B, µ..5 The logarithmic moments in the dual space, d ˆ L n µ ln n ρ + B µ, µ,.6 obey a RG equation, dl n µ d ln µ = Γ cuspα s L n+1 µ γ + α s L n µ n L n 1 µ,.7 which is simpler than the corresponding expressions for the logarithmic moments of φ + B, see [33]. Its solution can be explicitly written as L n µ = e V m= g m m! n j= n! n j! j! lnn j µ L m+j µ..8 µ However, to make use of this relation in practice, one would have to truncate the infinite sum, i.e. expand the result for g 1. Evidently, the solution for the spectral function.1 which contains the same information as.8 is much simpler and very economic. Moreover, in terms of the spectral function ρ + B, µ, the RG-improved factorization formulas in exclusive B-decays take a particularly simple form. For example, using the results from [7 9], the form factor relevant for the radiative leptonic B γlν decay in the heavy-quark limit, can be written as F E γ = HE γ, µ d je γ ˆ, µ ρ + B, µ..9 Here HE γ, µ contains the hard-matching coefficients from QCD onto SCET and HQET, and satisfies the RG equation [ dhe γ, µ = Γ cusp α s ln E ] γ d ln µ µ γ hα s HE γ, µ,.3 and js, µ is the hard-collinear function which has a perturbative expansion js, µ = 1 + α sc F ln s π µ 1 π and obeys the simple RG equation djs ], µ [Γ d ln µ = cusp α s ln s µ + γ hcα s js, µ,.3 8

10 with γ hc = Oαs. Notice that the hard-collinear function in dual space depends on ˆ via the combination lne γ ˆ /µ, where ˆ is defined in.16. The resulting form factor is RG-invariant, df E γ /d ln µ =, when γ h + γ hc + γ + =, which has been checked at one-loop accuracy. The RG-improved factorization formula for the form factor can thus be written as [ ] gµ,µh F E γ = e V hµ,µ h µh HE γ, µ h d E γ [ e V hcµ,µ hc µ gµ,µhc ] hc E γ ˆ je γ ˆ, µ hc [ e V µ,µ µ ] gµ,µ ρ + ˆ B, µ,.33 where the evolution factor for the hard coefficient, V h µ, µ h, is defined as V µ, µ in.11 with γ + replaced by γ h, and similarly for V hc µ, µ hc with γ + replaced by γ hc. Here, µ hc has to be chosen as a hard-collinear scale of order E γ ˆ, and µ h is the hard scale of order m b E γ. 6. Construction from Momentum Space As the second main subject of this paper, we are now going to present a general framework to construct momentum-space projectors from on-shell wave functions with definite relations to the LCDAs as defined above...1 -Particle Wave Functions An alternative method to construct momentum-space projectors that obey the WW relations starts from a generic Dirac matrix with the correct behaviour under space-time transformations as defined by the hadronic bound state that can be constructed from on-shell momenta k µ with k = for the light anti-quark, and m b v µ with v = 1 for the heavy quark, M B v, k = i f B m B [1 + /v /k γ 5 ] ψ B v k,.3 and that also fulfills the free Dirac equations for both constituents, /v 1 M B = M B /k =..35 The interpretation as a wave function for a -particle Fock state requires to consider a particular gauge, in this case light-cone gauge with n Ax =. We define a Lorentzinvariant integration measure dk for an on-shell massless particle, for which we can choose 6 For simplicity, we do not disentangle the matching and running for the B-meson decay constant in HQET. 9

11 a representation that reflects the light-cone kinematics of a hard-scattering process with the azimuthal angle in the transverse plane integrated out, dk := d k d = d3 k π v k, with kµ = nµ + + kµ + k n µ..36 To make contact with the general definition of LCDAs, we consider the convolution with a hard-scattering kernel that is at most linear in k, and obtain [ T dk tr + k µ T µ ] M B v, k = i f B m B i f B m B [ tr [ tr d T d T µ d k ψ Bx 1 + /v d k ψ Bx 1 + /v with x = v k = + k /. Comparison with. implies together with φ + B = d k ψ Bx = φ B = d k k ψ B x = dη φ B η φ+ B η = dx ψ B x, /n + + k ] /n γ 5 k ] γµ γ 5,.37 dx x ψ B x,.38 dx x ψ B x = φ B,.39 which is in line with the WW relations. Interestingly, in this approximation, the LCDAs φ + B and φ B can be obtained by simple k -integrals of a light-cone wave function ψ B x in particular, this procedure could be used to match calculations in the k T - factorization approach [1 3] to collinear QCD factorization at tree level. Notice that in the WW approximation the wave function ψ B x can be reconstructed from φ + B or as follows, φ B ψ B x = φ+ B x x 1 x d φ + B x dx = d φ B x dx WW.. This also implies that the RG evolution for ψ B x, µ can be easily obtained from that of ρ + B x, µ as in.17, ψ B x, µ = 1 d 1 x x J ρ + B, µ..1 Using that ρ + B, µ = = d J 1 dx x J x dx ψ B x, µ ψ B x, µ,. 1

12 we obtain an explicit solution for the RG evolution of the wave function ψ B x, µ, ψ B x, µ = e V d dx x µ g x ˆ x x J J ψ B x, µ,.3 which is similar to the one for the LCDA φ + B, µ in.19. We stress that the scale dependence of ψ B x, µ, according to our definition, follows directly from the renormalization of the LCDAs. In the k T -factorization approach, the RG properties of the transversemomentum-dependent wave functions will in general be different, depending on the exact definition of the gauge-link operator beyond the collinear limit we refer to [] for more detailed discussions. The ansatz for the momentum-space projector is thus consistent with the RG behaviour within the WW approximation. As an example, exponentially decreasing LCDAs which are often used in phenomenological applications can be obtained from the model ψ B x e x/ 3 φ + B e /, φ B e /.. They correspond to a spectral function ρ + B = ρ B e /..5 Note that the spectral function ρ + B shows dual behaviour compared to the function φ + B, i.e. it is exponentially suppressed at small values of and vanishes linearly with 1/ for, whereas φ + B decreases linearly at small and vanishes exponentially for. For comparison, a free parton picture with v k = M B m b = Λ cf. [5] would correspond to ψ B x δx Λ Λ φ + B Λ θ Λ, φ B Λ Λ θ Λ,.6 with a spectral function ρ + B = ρ B 1 Λ J Λ..7 Numerical examples for the two models with a sample RG evolution are plotted in Fig. 1. As expected, the RG evolution tends to wash out the differences between the shapes of the input functions at higher scales. In particular, the LCDA of the free parton model has become a smooth function after RG evolution the oscillatory behaviour of ρ + B at small values of is a relic from the singular behaviour of ρ + B at = Λ. 11

13 Figure 1. Left: exponential model from. showing the LCDAs φ + B top, φ B center and the spectral function ρ + B bottom, using µ = 1 GeV and =.3 GeV as input solid line and evolving to another scale corresponding to g =.3 dashed line and neglecting the overall factor e V γeg for simplicity. Right: the same for the free parton model.6 with Λ =.3 GeV. 1

14 .. Corrections to Wandzura-Wilczek Relation Corrections to the WW approximation can be incorporated by abandoning the constraint /k = from the light quark s Dirac equation. We thus write a more general ansatz, M B M B v, k = i f B m B [1 + /v {/k ψ B1 x + x ψ B x} γ 5 ]..8 Considering again the convolution with an appropriate hard-scattering kernel, one has [ T dk tr + k µ T µ ] M B v, k = i f B m B i f B m B i f B m B [ tr [ tr [ tr d T d T d T µ dx ψ B1 x 1 + /v dx ψ B x 1 + /v dx ψ B1 x 1 + /v The LCDAs for this ansatz follow as φ + B = dx { ψ B1 x + x ψ B x}, together with /n+ + x /n γ5 ] x /n+ + x /n γ5 ] ] x γ µ γ5..9 φ B = dx {x ψ B1 x + x ψ B x},.5 φ B dη φ B η φ+ B η = dx x ψ B x..51 The wave function ψ B x can thus be reconstructed from integrals of the 3-particle LCDAs Ψ A Ψ V using.6, and vice versa...3 Higher Fock States Our formalism can also be applied to construct on-shell wave functions for higher Fock states. As an example, we will consider the quark-antiquark-gluon contribution in the WW approximation. We recall that the partonic Fock-state interpretation refers to light-cone gauge n Ax =, which we will assume in the following. For the 3-particle contributions, it is for practical purposes often sufficient to consider convolutions with a hard-scattering kernel that does not depend on any transverse momenta. To keep the discussion simple, we will therefore focus on the strict collinear limit, i.e. we will assume a corresponding hard-scattering kernel without the linear terms in partonic transverse momenta. The conventional definition of 3-particle LCDAs starts from the position-space matrix element [5] see also [1, ] z ν q β z [z, uz] gg µν uz [uz, ] h α v Bv = f [ { B m B 1 + /v v µ /z tγ µ ΨA t, u Ψ V t, u } ] αβ iσ µν z ν ΨV t, u z µ XA t, u + zµ t /zỹat, u γ 5,.5 13

15 with t = v z. In order to relate this representation to the wave-function approach, we need the momentum-space projector corresponding to the non-local matrix element q β z g A µ uz h α v Bv.53 in light-cone gauge. In the collinear approximation, we further have z µ = τn µ, such that t = τ in a frame 7 where v n = 1. We then obtain M 3 B µ v,, ξ = i [ { f B m B 1 + /v n µ ξ /n n µ γµ Ψ A, ξ Ψ V, ξ ] +n µ /n n µ + γµ /n Ψ V, ξ n µ X A, ξ + n µ /n Y A, ξ }γ 5,.5 where we introduced the Fourier-transformed LCDAs dτ du Ψ A, ξ = π π τ ei+ξuτ ΨA τ, u, etc..55 Proceeding in analogy to the -particle construction, we may formulate an equivalent representation of the momentum-space projector starting from on-shell momenta for the Fockstate components k µ with k = for the light anti-quark, l µ with l = for the gluon, and m b v µ with v = 1 for the heavy quark. The most general ansatz that fulfills the equations of motion for all constituents, is given by /v 1 M 3 B µ = M 3 B µ /k = M 3 B µ l µ =,.56 M 3 B µ v, k, l = i f { B m B 1 + /v yk µ Φ 1 x, y + xl µ Φ x, y + x yγ µ v µ /l Φ 3 x, y + k µ /lφ x, y + l µ /lφ 5 x, y + x γ µ /l /lγ µ } Φ 6 x, y /kγ In general, the wave functions Φ i depend on three invariants, x = v k, y = v l and z = k l, but in the above decomposition we have neglected the invariant mass of the antiquark-gluon subsystem, z = k + l, for simplicity, which is based on the assumption that the wave functions only depend on the total invariant mass of the partonic configuration, m b v + k + l m b + m bx + y. Writing k µ = nµ + + kµ + k n µ, lµ = ξ nµ + + lµ + l ξ n µ,.58 which implies x = + k / and y = ξ + l /ξ, we may neglect any odd powers of k and l in the transverse-momentum integrals with a hard-scattering kernel. In light-cone 7 In Appendix A., we show the result for a general Lorentz frame. 1

16 coordinates, the momentum-space projector.57 may then be rewritten in terms of six independent structures, which we choose as n µ, nµ +, γµ, nµ /n, n µ +/n, γ µ /n..59 The light-cone gauge condition, M 3 B µ n µ =, further eliminates the nµ + and nµ +/n structures we use these constraints to determine Φ 1 and Φ. The remaining four Lorentz structures are of the form.5, and we can read off { Ψ A, ξ Ψ V, ξ = ξ dx dy ξ xx Φ x, y + y ξ x x Φ 3 x, y ξ + ξ xξ + ξ y Φ 5 x, y + xξ + ξ y } x Φ 6 x, y, Ψ V, ξ = ξ { dx dy xxy ξ + ξ y Φ 3 x, y ξ } + ξ3y ξ + x3ξ y Φ 5 x, y + xy ξ Φ 6 x, y, { Ψ A, ξ + X A, ξ = ξ dx dy xξx y Φ x, y y ξy ξx Φ 5 x, y ξ + x ξx + y ξ } Φ 3 x, y Φ 6 x, y, Ψ A, ξ + Y A, ξ = ξ { dx dy xx ξx y Φ x, y ξ + ξx y Φ5 x, y x ξx + y ξ Φ 3 x, y Φ 6 x, y }..6 In the collinear approximation, the four 3-particle LCDAs Ψ A, Ψ V, X A and Y A can thus be expressed through four independent on-shell wave functions. 3 Λ b Baryons 3.1 Light-Cone Distribution Amplitudes Light-cone distribution amplitudes for Λ b -baryons in HQET have been classified in [6] for related work, see also [8]. They contain the hadronic information entering factorization theorems for exclusive Λ b transitions in the heavy-quark limit see e.g. [7, 3]. In this work, we will focus on the LCDAs related to the leading 3-particle operators [6] gaugelinks are understood implicitly, but not shown for simplicity, ɛ abc u a τ 1 n Cγ 5 /n d b τ n ɛ abc u a τ 1 n Cγ 5 /n + d b τ n h c v Λ b v, s = f Λ b φ τ 1, τ u Λb v, s, h c v Λ b v, s = f Λ b φ τ 1, τ u Λb v, s,

17 for the chiral-odd part i.e. with an odd number of Dirac matrices in the light-diquark current, and ɛ abc u a τ 1 n Cγ 5 d b τ n h c v Λ b v, s = f 1 Λ φs b 3 τ 1, τ u Λb v, s, ɛ abc u a τ 1 n Cγ 5 iσ µν n µ + nν d b τ n h c v Λ b v, s = f 1 Λ φσ b 3 τ 1, τ u Λb v, s, 3. for the chiral-even part. Here u Λb v, s denotes the on-shell Dirac spinor for the Λ b -baryon, and the prefactors f 1, Λ b are defined by the normalization of the matrix elements of the corresponding local operators in HQET at a given renormalization scale Light-Cone Projectors for 3-Particle Fock State The above definitions can be cast into a manifestly Lorentz-invariant form by defining the most general non-local matrix elements in coordinate space as [7] ɛ abc u a αz 1 d b β z h c v Λ b v, s 1 { [ ] [ ] } f 1 Λ b M 1 v, z 1, z γ 5 C T + f Λ βα b M v, z 1, z γ 5 C T u Λb v, s, 3.3 βα with a part that contains an odd number of Dirac matrices t i = v z i, M v, z 1, z = /v Φ t 1, t, z 1, z, z 1 z + Φ X t 1, t, z 1, z, z 1 z t 1 t /z /v/z 1 /z 1 /v/z + i Φ t 1, t, z1, z, z 1 z /z 1 + t 1 Φ ii t 1, t, z 1, z, z 1 z t /z, 3. and a part that contains an even number of Dirac matrices, M 1 v, z 1, z = Φ 3 t 1, t, z1, z, z 1 z + Φ Y t 1, t, z1, z, z 1 z /z /z 1 /z 1 /z t 1 t + i Φ 3 t 1, t, z1, z, z 1 z /v/z 1 + t 1 Φ ii 3 t 1, t, z 1, z, z 1 z t /z /v. 3.5 Considering isospin invariance for the light-quark fields exchanging z 1 z and taking care of the charge-conjugation properties of the Dirac matrices, one obtains the following relations between the individual functions, Φ t 1, t, z 1, z, z 1 z = Φ t, t 1, z, z 1, z 1 z, Φ i t 1, t, z1, z, ii z 1 z = Φ t, t 1, z, z1, z 1 z, Φ X t 1, t, z1, z, z 1 z = Φ X t, t 1, z, z1, z 1 z,

18 and Φ Φ 3 t 1, t, z1, z, z 1 z = 3 t, t 1, z, z1, z 1 z, 3 t 1, t, z1, z, ii z 1 z = Φ 3 t, t 1, z, z1, z 1 z, Φ Y t 1, t, z1, z, z 1 z = Φ Y t, t 1, z, z1, z 1 z. 3.7 Φ i 3.1. The Chiral-Odd Projector M In the same way as we argued for B-mesons, we may again expand the arguments z 1 and z around the light-cone, such that n z i zi n + z i, to obtain the projector in coordinate-space as M v, z 1, z /n + φ τ 1, τ + /n i φ τ 1, τ + φ τ ii 1, τ + φ τ 1, τ i φ + τ ii 1, τ φ /z 1 + τ 1, τ /z τ 1 τ + φ X τ 1, τ /z 1 τ 1 /z τ /n /n + /n +/n + Ozi, n z i. 3.8 Here again we denote with τ i = n +z i the Fourier-conjugate variables to the momentum components i = n k i of the associated light-quark states in the heavy baryon, such that φ 1, dτ1 π ei 1τ 1 Comparison with the definition in 3.1 yields the relation dτ π ei τ φ τ 1, τ etc. 3.9 φ i τ ii 1, τ + φ τ 1, τ = φ τ 1, τ φ τ 1, τ, 3.1 φ i ii while the asymmetric combination of and φ, as well as φ X do not contribute in the collinear limit zi. After Fourier transformation, the general momentum-space representation for 3.8 including the first-order terms off the light-cone becomes M 1, = /n + φ 1, + /n φ 1, 1 γ µ 1 γ µ 1 γ µ 1 γ µ 1 1 dη 1 φ i η 1, φ X η 1, dη 1 φ i η 1, + φ X η 1, dη φ ii 1, η φ X 1, η dη φ ii 1, η + φ X 1, η /n+ /n /n /n + /n /n + /n+ /n k 1µ k 1µ k µ k µ Compared to the mesonic analogue, we observe that a larger number of independent terms that are sensitive to the transverse momenta of the light quarks in the hard-scattering kernel appear. 17

19 3.1.3 The Chiral-Even Projector M 1 Similarly, for the chiral-even projector we obtain the expansion around the collinear limit as M 1 v, z 1, z + φ 3 τ i 1, τ + φ 3 τ 1, τ /n +/n ii + φ 3 τ 1, τ /n /n + φ i 3 τ 1, τ /v/z 1 + φ Y τ 1, τ ii φ 3 τ 1, τ /z /v τ + τ 1 /z /n + /n /z 1 τ τ 1 where now from the comparison with 3. one has the relations φ s 3τ 1, τ = φ σ 3 τ 1, τ = φ 3 τ i 1, τ + φ 3 τ 1, τ + φ ii 3 τ 1, τ i φ 3 τ 1, τ + Oz i, n z i, 3.1 φ ii 3 τ 1, τ, It is sometimes more convenient to define symmetric and antisymmetric combinations of the functions φ s 3 and φσ 3 which will be denoted as [6] φ + 3 τ 1, τ = φs 3 τ 1, τ + φ σ 3 τ 1, τ φ + 3 τ 1, τ = φs 3 τ 1, τ φ σ 3 τ 1, τ = = φ φ ii 3 τ 1, τ + φ 3 τ i 1, τ + φ 3 τ 1, τ 3 τ 1, τ,. 3.1 The expansion of the corresponding momentum-space projector then takes the general form M 1 1, = /n /n φ , + /n +/n 8 dη 1 φ i 3 η 1, /v γ µ dη 1 φ Y η 1, /n γ µ φ + 3 1, k 1µ 1 k 1µ 1 dη φ ii 3 1, η γ µ /v dη φ Y 1, η γ µ /n k µ k µ Again, in the general case, it involves four independent structures related to transverse momenta of the light quarks Wandzura-Wilczek Approximation In the WW approximation, the matrices free quark fields, M 1, z 1, z fulfill the equations of motion for γ µ i µ M 1, v, z 1, z = i µ 1 M 1, v, z 1, z γ µ =, 3.16 which translates into differential equations for the LCDAs in the collinear limit. These can be obtained by expanding the above equations around the light-cone, and solving for the 18

20 derivatives with respect to the arguments z i, z 1 z off the light cone. Alternatively, one can start from the expanded form of M 1, and consider the projected equations /n + /n γ µ i µ M 1, v, z 1, z z 1, = = i µ M 1 1, /n /n + v, z 1, z γ µ z 1, = = In both cases, this yields the following WW relations for the LCDAs in M, φ i τ 1, τ + φ X τ 1, τ + τ 1 τ 1 φ τ 1, τ =, φ ii τ 1, τ + φ X τ 1, τ + τ τ φ τ 1, τ = For the Fourier-transformed LCDAs this implies φ i 1, + φ X 1, 1 1 φ 1, =, φ ii 1, + φ X 1, φ 1, = Equivalently, by considering linear combinations of the above, the following relations hold φ i 1, φ ii 1, = 1 1 φ 1, φ X 1, + φ 1, φ 1, = 1 1 φ 1, + φ 1,, φ 1,. 3. This reveals that, once the functions φ and φ which are the relevant LCDAs in the collinear limit are given, the function φ X and the asymmetric combination of φ i,ii can be calculated from the WW approximation. In a similar way, for the LCDAs in M 1 we obtain the relations or, in momentum space, φ i 3 τ 1, τ + τ φ 1 3 τ τ i 1, τ + φ 3 τ 1, τ =, 1 φ ii 3 τ 1, τ + τ φ 3 τ τ ii 1, τ + φ 3 τ 1, τ =, 3.1 φ i 3 1, 1 φ 3 1, + 1 φ i 3 1, =, 1 φ ii 3 1, φ 3 1, + φ ii 3 1, =. 3. Notice that in this case, the function φ Y does not appear in the WW relations, and therefore remains independent, whereas the functions φ + 3 and φ + 3 that are relevant in the collinear limit are related by 1 1 φ + 3 1, = φ + 3 1, = φ 3 1,

21 3. Construction from Momentum Space We are now going to apply the same formalism that we have developed for B-mesons to construct momentum-space projectors for Λ b -baryons from 3-particle wave functions. To keep the discussion simple, we will ignore corrections to the WW relation in the rest of the paper. The most general form of the momentum-space projectors can then be written as M 1 v, k 1, k = ψ s x 1, x, K /k /k 1, M v, k 1, k = ψ v x 1, x, K /k /v /k 1, 3. with x i = v k i and K = k 1 + k and two independent wave functions ψ s and ψ v. The equations of motion, /k M 1, v, k 1, k = M 1, v, k 1, k /k 1 =, are again trivially fulfilled for on-shell quarks, ki =. On the other hand, the invariant mass of the diquark system in principle can be arbitrary, K. For simplicity, we will ignore the potential K dependence, which would correspond to the case where the wave function only depends on the total invariant mass of the three quarks in the Λ b -baryon, m b v + k 1 + k m b + m bx 1 + x The Chiral-Odd Projector M To compare with the general definition of LCDAs, we again consider the convolution with a hard-scattering kernel that is at most linear in k i. For the chiral-odd projector M, we obtain dk 1 = d 1 d [ T dk tr 1, + k µ i T µ i 1, ] M v, k 1, k { dx 1 dx 1 [ /n + tr T 1, 1 + x 1 1 x /n ] [ tr Tµ 1 1, 1 x 1 1 /n +/n [ tr } Tµ 1, γµ + 1 x 1 1 x /n /n + γ µ ] 1 x /n /n + + x 1 1 x /n +/n ] ψ v x 1, x, 3.5 where the momentum integrations for the light quarks are defined as in the mesonic case. Comparison with the momentum-space projector shown in 3.11 above, yields φ 1, = φ 1, = 1 dx 1 1 dx 1 dx 1 ψ v x 1, x, dx x 1 1 x ψ v x 1, x, 3.6

22 together with φ i 1, = 1 φ ii 1, = 1 1 dx 1 1 dx 1 dx x x 1 1 ψ v x 1, x, dx x 1 x ψ v x 1, x, 3.7 and φ X 1, = 1 dx 1 dx x 1 1 x ψ v x 1, x It can easily be checked that the LCDAs constructed in this way satisfy the WW relations as derived above. Notice that our simplified ansatz relates all LCDAs to x i moments of only two fundamental wave functions ψ v and ψ s below. The functional form of ψ v can be reconstructed, for instance, from ψ v x 1, x = d φ x 1, x dx 1 dx x 1 x = d φ x 1, x dx 1 dx WW. 3.9 In a more general ansatz, these relations would be modified by the non-trivial K -dependence of the wave functions. In the simplest case, we could again model the wave functions by assuming an exponential dependence of ψ v on x 1 + x with a single hadronic parameter measuring the average energy of the light quarks, exp x 1+x ψ v x 1, x. 3.3 This ansatz yields the following exponential model for the various LCDAs defined above, and and with φ 1, 1 6 e 1+ /, φ 1, 1 e 1+ /, 3.31 φ i 1, 1 + e 1+ /, φ ii 1, + 1 e 1+ /, 3.3 φ X 1, 1 e 1+ / φ i 1, φ ii 1, 1 3 e 1+ /

23 For comparison, an alternative model based on a free parton picture with the constraint x 1 + x = Λ = M Λb m b would correspond to a wave function ψ v x 1, x 15 Λ 5 δx 1 + x Λ From this, our construction immediately yield the corresponding expressions for the LCDAs in terms of θ-functions, φ 1, 15 1 Λ 1 Λ 5 θ Λ 1, φ 1, 5 Λ Λ 5 θ Λ 1, etc To illustrate these model results, we compare the LCDA φ 1, following from the exponential ansatz in 3.31, the free-parton approximation 3.36 and the model from Eq. 38 in [6]. For that purpose, we disentangle the dependence on the total lightcone momentum and the momentum fractions u of the light quarks by considering the projections 1 3 e / with = Λ 5 =. GeV h := du φ u, ū = 3 e /ɛ 6ɛ [6] with ɛ =. GeV 5 3 Λ 8 Λ 5 θ Λ 3.36 with Λ = 1 GeV 3.37 and g u := d φ u, ū = uū uū ɛ + 3a 5u ū 1 ɛ 1 5uū Λ 3.31 with =. GeV ɛ =. GeV [6] with ɛ 1 =.65 GeV a = 1/ with Λ = 1 GeV 3.38 The parameter in the first case has been related to the value of Λ in the third case, such that the 1 moment of h is identical in both cases. Notice that the two models for the LCDA φ which are based on a wave function ψ v that only depends on the sum of the light-quark energies, lead to ψ v = ψ v x 1 + x φ u, ū uū In contrast, the model in [6] takes into account a non-trivial shape from the next-toleading term in a Gegenbauer expansion. That model also prefers a smaller value for the parameter ɛ and a corresponding larger value for the inverse moment 1 than in the other two models. The numerical comparison between the three models is shown in Fig..

24 Figure. The functions h left and g u right for different models of the LCDA φ 1,. Solid lines correspond to the model 3.31, dashed lines to 3.36 and light-dotted lines to [6] in the right figure the dashed line is on top of the solid line. 3.. The Chiral-Even Projector M 1 For the chiral-even projector M 1, we again consider the convolution with a hard-scattering kernel, and obtain [ T dk 1 dk tr 1, + k µ i T µ i 1, ] M 1 v, k 1, k = d 1 d dx 1 1 [ tr T 1, { dx x 1 1 /n +/n [ tr Tµ 1 1, 1 x 1 1 /n + [ tr Tµ 1, γµ } + 1 x /n ] /n x 1 1 x /n γ µ ] 1 x /n + + x 1 1 x /n ] ψ s x 1, x. 3. Comparison with the momentum-space expression 3.15 yields φ + 3 1, = φ 3 1, + φ i 3 1, = dx 1 1 dx x 1 1 ψ s x 1, x, 3

25 φ + 3 1, = φ 3 1, + φ ii 3 1, = dx 1 1 and φ 3 1, = dx 1 dx 1 ψ s x 1, x, 1 dx 1 x ψ s x 1, x, 3.1 φ Y 1, = 1 dx 1 dx 1 x 1 x ψ s x 1, x Again, the wave function ψ s in our approximation can be reconstructed from ψ s x 1, x = d φ 3 x 1, x WW. 3.3 dx 1 dx x 1 x With the exponential model for the wave function, we now obtain exp x 1+x ψ s x 1, x, 3. which yields and φ + 3 1, 3 e 1+ /, φ + 3 1, 1 3 e 1+ /, 3.5 The free parton picture now yields 6 φ 3 1, 1 e 1+ /, φ Y 1, 1 e 1+ /. 3.6 etc. φ + 3 1, 15 Λ 1 Λ 5 θ Λ 1, φ + 3 1, 15 1 Λ 1 Λ 5 θ Λ 1, Renormalization-Group Evolution In the following, we will focus on the twist- Λ b LCDA φ 1, which enters the leading terms in factorization theorems for exclusive heavy-to-light decay amplitudes in the heavyquark limit, see e.g. []. Its one-loop RG equation has been extensively discussed in [6],

26 and reads 8 dφ 1,, µ d ln µ [ = Γ cusp α s ln 1 + α sc F π 1 ] µ + γ + α s φ 1,, µ 1 dη 1 Γ + 1, η 1, α s φ η 1,, µ dη Γ +, η, α s φ 1, η, µ dv V ERBL u, v φ v, v, µ, 3.8 where = 1 +, and u = 1 ū = 1 /, and v = 1 v. Here the first three lines correspond to the LN kernel for heavy baryons with the same anomalous dimensions as in.8, whereas the last term is the ERBL kernel [1, ], which arises from gluon exchange among the light quarks in the heavy baryon Analytic Solution We follow a similar strategy as for the B-meson LCDA φ +, and as a first step introduce the logarithmic Fourier transform which is in almost one-to-one correspondence to the Mellin moments discussed in [6], ϕ θ 1, θ, µ = d 1 1 Next, we introduce the ansatz ϕ θ 1, θ, µ := Γ1 iθ 1 Γ1 iθ Γ1 + iθ 1 Γ1 + iθ d φ 1,, µ d 1 1 in complete analogy to the mesonic case, such that d 1 d 1 φ 1,, µ = 1 1 J µ d iθ1 iθ. 3.9 µ µ iθ1 µ iθ ρ 1,, µ J 1 ρ 1,, µ The inverse transformation that expresses the dual spectral function ρ 1, in terms of the momentum-space LCDA is then given by ρ 1,, µ = d 1 1 d 1 1 J J 1 φ 1,, µ. 3.5 The one-loop RG equation 3.8 can be rewritten for the spectral function ρ 1, in a straightforward manner. In the absence of the ERBL kernel, the LN terms alone would take an analogous factorized form as in the case of the B-meson spectral function ρ + B, µ, dρ 1,, µ d ln µ [ = Γ cusp α s ln LN µ ˆ 1 ˆ + γ + α s ] ρ 1,, µ With a slight abuse of notation, we write a colour factor C F in the baryon case, although more precisely, the colour factor arises as 1 + 1/N C which only coincides with C F for N C = 3. It should be noted, however, that in the LN kernel for N C 3 one has to add up the contributions from N C 1 light spectators in a color singlet baryon, such that the net result in front of Γ cusp ln µ would be proportional to C F again. 5

27 In this approximation the RG equation would simply be solved by g ρ 1,, µ = e V µ ρ LN ˆ 1 ˆ 1,, µ, 3.5 with ˆ i = e γ E i, and the RG functions V and g from.11. The derivation of the ERBL term for the evolution of ρ 1,, µ, however, is more complicated the details can be found in Appendix B. Interestingly, the final result takes a simple form when written in terms of the reduced dual momentum and dual momentum fractions, r and u = 1 ū = Writing we obtain dˆρ r, u, µ d ln µ ρ 1, ˆρ r, u 3.56 [ = Γ cusp α s ln µ ] u ū + γ + α s ˆρ LN+ERBL r, u, µ + α sc F π 1 ˆ r dv V ERBL u, v ˆρ r, v, µ If we expand the spectral function in terms of Gegenbauer polynomials C 3/ n u 1, which are the eigenfunctions of the ERBL kernel, ˆρ r, u, µ := n=,,,... the coefficients f n r, µ satisfy the RG equation df n r, µ d ln µ = u ū f n r, µ C 3/ n u 1, 3.58 [ Γ cusp α s ln µˆ r + γ + α s + α ] sc F π γerbl n f n r, µ Γ nm α s f m r, µ, 3.59 with γn ERBL given in B., and a non-diagonal contribution from the substitution of variables in the LN kernel, given by Γ nm α s = Γ cusp α s n + 3 n + 1n + 1 du u ū lnu ū C 3/ n u 1 C 3/ m u = Γ cusp α s , 3.6 6

28 where the different lines refer to n =,,, etc and the columns to m =,,, etc. As one can see, the off-diagonal terms are typically smaller than the diagonal ones, and therefore, as a first approximation could be neglected. In that case, the particular form of the leadingtwist baryon LCDA in 3.39, which translates into ˆρ r, u, µ u ū f r, µ, 3.61 and and φ 1,, µ ψ v x 1 + x, µ 1 x 1 + x 5/ d r r d r r 1 + J 3 1 r 5/ J r x1 + x r f r, µ, 3.6 f r, µ, 3.63 would be stable under evolution. Diagonalizing the r.h.s. of the RG equation, truncated to a finite number of Gegenbauer coefficients, is now also a straightforward task, which will be illustrated below. As already discussed in [6] the numerical effect of the ERBL term is in any case expected to be sub-leading, and for practical applications it should be sufficient to treat it in an approximate way. We finally note that the connection between the function φ τ 1, τ, appearing in the light-cone matrix elements in coordinate space, and the spectral function ˆρ r, u in the baryonic case is given by ˆρ r, u, µ = dτ1 π dτ π ] ] 1 exp [ iū rτ 1 exp [ iu φ 1 rτ τ 1, τ, µ Numerical Examples and Asymptotic Form In the following, we study the coefficient functions f, r, µ and their RG behaviour, starting from different models for the LCDA φ defined at some input scale µ. For a given LCDA, making use of C.7 in the Appendix, we find using 1 = u, = ū f n r, µ = n + 3 n + 1n + 1 d r J n+3 r du C n 3/ u 1 φ u, ū, µ 3.65 Notice that the Gegenbauer expansion of the original LCDA φ u, ū directly translates to the Gegenbauer expansion of the spectral function ˆρ r, u. For the models discussed 7

29 above 3.37, this leads to model 1: f r, µ = 1 r e / r, f n r, µ = ; model : f r, µ = 1 r e ɛ / r, f n> r, µ =, f r, µ = a { r e 1/xr 7 x 5 r + 6 x r + x 3 r + 6 x r + 1 x r + 1 } + 1 x r 1 6 x r, x r r/ɛ 1 ; model 3: f r, µ = 15 r Λ Λ J Λ 5 r, f n r, µ = The qualitative behaviour of model 1 and model 3 is similar to what we have discussed for the corresponding functions in the B-meson case. Notice that the contribution of the coefficient function f r to the spectral function ˆρ r, u in model is concentrated at very low values of r, while for generic values of r its contribution is practically negligible. In order to study systematic deviations from the particular form of ˆρ r, u in 3.61, it would therefore be more convenient to define a modified version, for which we propose model : f r, µ 1 r e ɛ / r, f r, µ a 6 ɛ 1 r e ɛ 1/ r, 3.67 which has the same functional form as model at large r except for a different normalization such that the 1/ r moment of f r in model and model coincide. In the original momentum space, this corresponds to a LCDA model : φ 1,, µ = 1 ɛ e 1+ /ɛ 1 + a ɛ 6 e 1+ /ɛ Concerning the RG behaviour, we first note that the explicit solution for the functions f, r, µ in the absence of higher Gegenbauer coefficients reads g f r, µ = e V µ { c1 e.877g.378 c e.g}, ˆ r g f r, µ = e V µ {.17 c1 ˆ r e.877g + c e.g}, 3.69 where the two integration constants are related to the initial condition of the evolution via c 1 =.97 f r, µ f r, µ, c =.139 f r, µ +.97 f r, µ. 3.7 In the asymptotic limit, i.e. for large renormalization scales and large values of g, the first exponential in the curly brackets dominates, and the ratio of the two coefficients approaches 8

30 a constant, f /f.17. This is illustrated for model and model in Fig. 3. For both models, the asymptotic value for f /f is reached 9 for g > 3. When higher Gegenbauer moments are included, the asymptotic form is similarly determined by the largest eigenvalue of the RG equation 3.59 after subtracting the LN terms. We then find that the ratio f /f converges to about %, and that the admixtures of the higher Gegenbauer moments are less important, with f /f 9%, f 6 /f 5%, etc. The resulting asymptotic u -dependence of the spectral function, corresponding to the different levels of truncation in the Gegenbauer expansion, is illustrated in Fig.. The functional form that is approached asymptotically is well approximated by ˆρ r, u, µ u 1 u 1/3 f r, µ. Figure 3. Numerical examples for the evolution of the ratio of Gegenbauer coefficient functions f r, µ/f r, µ. Left: for model from [6]. Right: for model, with a = 1/3, ɛ =. GeV, ɛ 1 =.65 GeV, see text. Solid lines: input model at µ ; dashed lines: RG evolution corresponding to g =.3; dotted line: with g = 3. 9 In practice, this is limited by the fact that the evolution of the LCDAs within HQET has to be replaced by the standard QCD evolution above µ m b. 9

31 Figure. Asymptotic u -dependence of the spectral function ˆρ r, u, µ/f r, µ, for different levels of truncation in the Gegenbauer expansion: n = thin dotted, n = thin dashed, n = thick dotted, n = 6 thick dashed, n = 8 solid. The gray band shows the approximation ˆρ r, u, µ/f r, µ u 1 u 1/3. 3

32 Summary We have investigated light-cone distribution amplitudes LCDAs as defined in heavy-quark effective theory for B-mesons and Λ b -baryons. On the one hand, we have constructed easyto-use momentum-space representations for the leading Fock states, which reduce to an expansion in terms of conventional LCDAs when convoluted with a hard-scattering kernel in the collinear QCD factorization approach, but also allows for a comparison with models for transverse-momentum-dependent wave functions. In the simplest case, our construction automatically implements so-called Wandzura-Wilczek relations, which connect different LCDAs in the limit where higher Fock-state contributions to the equations of motion are neglected. We have also illustrated how corrections to the Wandzura-Wilczek approximation can be taken into account consistently within our approach. For the baryonic case in particular, our ansatz leads to a significant reduction of independent hadronic functions which appear at sub-leading order in the collinear expansion. The sub-leading functions are needed, for instance, in SCET sum rules analyses of exclusive Λ b transitions for cases where the standard QCD factorization approach would lead to endpoint-sensitive formally ill-defined convolution integrals see e.g. the discussion in [7]. Furthermore, we have found a new representation of LCDAs in terms of dual spectral functions, which are the eigenfunctions of the Lange-Neubert renormalization kernel. The connection between the LCDAs and their dual representations is via convolution integrals with Bessel functions. In the dual space, the solutions to the renormalization-group equations are extremely simple. In the mesonic case, they are local in the dual momentum variable of the light quark. We have demonstrated the simplifications that arise when reformulating the factorization theorem for radiative leptonic B γlν decays in terms of the new spectral function and an associated new hard-collinear function, which evolve both by a multiplicative factor. Again, the baryonic case is more complicated, because the Lange- Neubert kernel and the Efremov-Radyushkin Brodsky-Lepage kernel cannot be separated in the usual momentum space. In the dual space, however, a separation of the reduced dual momentum r = 1 / 1 + and the dual momentum fraction u = 1 / 1 + is possible, which opens the way for a systematic solution of the baryonic renormalizationgroup equation in terms of an expansion in Gegenbauer polynomials as known from the pion LCDAs. In summary, our results should be helpful in calculations of exclusive decay amplitudes for hadrons containing a heavy quark in the framework of QCD factorization, softcollinear effective theory or light-cone sum rules. Specifically for applications involving heavy baryons, we expect more transparent and efficient estimates of theoretical hadronic uncertainties, which are needed, for instance, to constrain physics beyond the Standard Model from rare decays like Λ b Λµ + µ which are currently studied at hadron colliders. 31