Subleading power corrections to the pion-photon transition form factor in QCD

We reconsider QCD factorization for the leading power contribution to the γ*γ → π0 form factor Fγ*γ→π0(Q2) at one loop using the evanescent operator approach, and demonstrate the equivalence of the resulting factorization formulae derived with distinct prescriptions of γ5 in dimensional regularization. Applying the light-cone QCD sum rules (LCSRs) with photon distribution amplitudes (DAs) we further compute the subleading power contribution to the pion-photon form factor induced by the “hadronic” component of the real photon at the next-to-leading-order in Oαs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{O}\left({\alpha}_s\right) $$\end{document}, with both naive dimensional regularization and ’t Hooft-Veltman schemes of γ5. Confronting our theoretical predictions of Fγ*γ→π0 (Q2) with the experimental measurements from the BaBar and the Belle Collaborations implies that a reasonable agreement can be achieved without introducing an “exotic” end-point behaviour for the twist-2 pion DA.


JHEP12(2017)037
as well as lattice QCD simulations (see for instance [48][49][50]). In particular, an "exotic" twist-two pion DA with the non-vanishing end-point behaviour was proposed [51,52] to accommodate the anomalous BaBar data at high Q 2 , but was soon critically examined in [22] concluding that a reasonable description of the BaBar data in [51,52] is achieved rather due to the introduction of a sizable nonperturbative soft correction from the TMD pion wavefunction.
The outline of this paper is as follows: in section 2 we recalculate the one-loop hard function entering the factorization formula for the pion-photon form factor at leading power in 1/Q 2 with both the NDR and HV schemes of γ 5 , and demonstrate the renormalizationscheme independence of the factorization formulae for physical quantities explicitly. It will be shown that our expression of the NLO hard-scattering kernel in the NDR scheme reproduces the classical result obtained by Braaten [6] and the renormalization-prescription dependence of the short-distance coefficient at O(α s ) will be cancelled precisely by the scheme dependent twist-2 pion DA at one loop. We then establish QCD factorization for the vacuum-to-photon correlation function defined with a pseudoscalar interpolating current for the pion state and an electromagnetic current carrying a space-like momentum q µ (q 2 = −Q 2 ) at one loop in section 3. It will be further proved that the resulting hard matching coefficients obtained in the NDR and HV schemes are related by the finite renormalization constant term, which is introduced in the HV scheme in order to fulfill the Adler-Bardeen theorem for the non-renormalization of the axial anomaly [11,53,54]. The next-to-leading-logarithmic (NLL) resummation improved LCSR for the hadronic photon correction to F γ * γ→π 0 (Q 2 ) will be also presented with the aid of the parton-hadron duality ansatz. Taking advantage of the newly derived subleading power correction and the twistfour effect from both the two-particle and three-particle pion DAs at tree level [22,26], we will provide updated theoretical predictions for the pion-photon form factor in section 4 with distinct nonperturbative models for the twist-2 pion DA. A summary of our observations and the concluding remarks are presented in section 5. We collect the two-loop evolution functions for the leading twist DAs of the pion and the photon in appendix A and display the spectral representations of the convolution integrals for the construction of the NLL LCSR of the hadronic photon contribution in appendix B.

Factorization of the leading power contribution
We will first compute the leading power contribution to the pion-photon form factor at one loop π(p)|j em µ |γ(p ) = g 2 em µναβ q α p β ν (p )F γ * γ→π 0 (Q 2 ) , (2.1) with both the NDR and HV schemes for the γ 5 matrix in D dimensions, applying the technique of soft-collinear effective theory (SCET) and the evanescent operator approach.
The main purpose of this section is to provide an alternative way to establish the QCD factorization theorem for the leading twist contribution to the π 0 γ * γ form factor at the operator level. Here, q = p − p , p refers to the four-momentum of the pion, the on-shell JHEP12(2017)037 x p x p ν µ q p ′ x p x p µ ν q p ′ (a) (b) Figure 1. Diagrammatical representation of the tree-level contribution to the partonic amplitude γ γ * → qq induced by two electromagnetic currents.
photon carries the four-momentum p and j em µ = q g em Q qq γ µ q , 0123 = −1 . (2.2) We further introduce a light-cone vectorn µ parallel to the photon momentum p , define another light-cone vector n µ along the direction of the momentum p in the massless pion limit, and employ the following power counting scheme at large momentum transfer 2.1 QCD factorization of F γ * γ→π 0 (Q 2 ) at tree level QCD factorization for the leading-twist contribution to the γ * γ → π 0 form factor at tree level can be established by inspecting the four-point QCD matrix element Π µ = q(x p)q(x p)|j em µ |γ(p ) (2.4) at leading order (LO) in α s , where x indicates the momentum fraction carried by the collinear quark of the pion andx ≡ 1 − x. Computing the two diagrams displayed in figure 1 yields Here, O A, µν (0) and O B, µν (0) denote the tree-level partonic matrix elements of the collinear operators O A, µν and O B, µν in SCET O j, µν (x, x ) ≡ q(x p)q(x p)|O j, µν (x )|0 =ξ(x p) Γ j, µν ξ(x p) δ(x − x ) + O(α s ) , (2.6) and the convolution integration is represented by an asterisk. The definition of the SCET operator O j, µν in the momentum space is given by dλn · A c (λn) (2.8) and Γ A, µν = γ µ,⊥ n γ ν,⊥ , Γ B, µν = γ ν,⊥ n γ µ,⊥ . (2.9) To facilitate the determination of the hard function entering the leading power factorization formula of F γ * γ→π 0 (Q 2 ), we employ the SCET operator basis {O 1, µν , O 2, µν , O E, µν } with where O E, µν is an evanescent operator vanishing in four dimensions and It needs to be pointed out that the -tensor is always defined in four dimensions, since its indices will be never contracted with a D-dimensional index. The effective operator O 1, µν cannot couple with a collinear pion state due to the parity conservation. Taking advantages of the operator identities we observe that the two tree-level diagrams in figure 1 give rise to the identical contribution to the pion-photon transition form factor and such observation can be further generalized to all orders in QCD applying the charge-conjugation transformation. We now employ the operator matching equation with the evanescent operator 13) and expand all quantities to the tree level, yielding Utilizing the definition of the leading twist pion DA on the light cone it is straightforward to derive the tree-level factorization formula of the π 0 γ * γ form factor Figure 2. Diagrammatical representation of the one-loop contribution to the partonic amplitude γ γ * → qq induced by two electromagnetic currents. The corresponding symmetric diagrams obtained by exchanging the two photon states are not shown.

QCD factorization of
We proceed to compute the NLO QCD correction to the four-point partonic amplitude Π (1) µ at leading power in 1/Q 2 for the determination of the hard function T 2 at O(α s ). It needs to be stressed that the QCD matrix element Π µ defined by two electromagnetic currents is independent of the prescription of γ 5 in the D-dimensional space and the renormalization scheme dependence of the perturbative matching coefficient T 2 comes solely from the radiative correction to the twist-2 pion DA φ π (x, µ), whose definition depends on the precise treatment of the Dirac matrix γ 5 in dimensional regularization, due to the infrared subtraction. Evaluating the hard contribution to the one-loop diagrams displayed in figure 2 with the method of regions [55] leads to 2,hard (x ) + . . . , (2.21) where the γ 5 -prescription independent amplitude A 2,hard reads Expanding the matching equation (2.13) to the one-loop order yields Now we are in a position to derive the master formula for the one-loop perturbative matching coefficient T (1) i by implementing both the ultraviolet (UV) renormalization and the infrared (IR) subtraction. Following the strategy presented in [18] the UV renormalized matrix element of the SCET operator O i, µν at O(α s ) is given by where the bare matrix element M (1) ij,bare depends on the IR regularization scheme R. Applying the dimensional regularization for both the UV and IR divergences, the bare matrix element M (1) ij,bare vanishes due to scaleless integrals entering the relevant one-loop computation. Inserting (2.24) into (2.23) and comparing the coefficient of O 2, µν (0) give rise to We then need to evaluate the above-mentioned renormalization constants Z indicating the absence of the operator mixing between O 1 µν and O 2 µν under renormalization. The renormalization constant Z 11 governs the renormalization-scale evolution of the leading twist (longitudinally polarized) ρ-meson DA. In addition, the IR subtraction term T 22 will remove the collinear contribution to the QCD amplitude Π µ at one loop so that the matching coefficient T We are now ready to discuss the renormalization constant Z E2 of the evanescent operator O E, µν for the derivation of the final result of the matching coefficient T (1) 2 . Applying the renormalization prescription that the IR finite matrix element of the evanescent operator O E, µν vanishes with dimensional regularization applied only to the UV divergences and with the IR singularities regularized by any parameter other than the dimensions of spacetime [19,20] and making use of the identity (2.24) yield (2.27) The one-loop matching coefficient of the physical operator O 2, µν can be readily obtained by substituting (2.27) into (2.25) . (2.28) The one-loop contribution to the matrix element of the evanescent operator O E, µν depends on the renormalization prescription of γ 5 in the D-dimensional space. We will employ both the NDR and HV schemes of γ 5 below for the illustration of the prescription independence of the factorization formula of F γ * γ→π 0 (Q 2 ) at O(α s ) and at leading power in 1/Q 2 . This amounts to showing that the renormalization scheme dependence of the short-distance coefficient function cancels against that of the twist-2 pion DA precisely. Evaluating the SCET diagrams displayed in figure 3 with the Wilson-line Feynman rules, we find that only a single diagram 3(a) can generate a nonvanishing contribution to M (1)off E2 using the NDR scheme of γ 5 , which turns out to be proportional to the spin-dependent term of the Brodsky-Lepage evolution kernel [1,2]. The resulting expression of the infrared subtraction term T (0) is given by However, computing the one-loop matrix element of the evanescent operator O E, µν with the HV scheme of γ 5 leads to   in the NDR scheme is identical to that presented in [6] using the trace formalism.
We now aim at demonstrating the renormalization prescription independence of the one-loop factorization formula for the pion-photon form factor 33) where the superscript "∆" indicates the γ 5 -scheme in dimensional regularization. Taking advantage of the relation of the twist-2 pion DA between the NDR and HV schemes where the finite renormalization kernel Z −1 HV is given by [56] It is then straightforward to show that which cancels against the renormalization scheme dependence of the NLO hard function T (1), ∆ 2 (x, µ) as displayed in (2.31). We emphasize again that the γ 5 -prescription independence of the leading power factorization formula for F γ * γ→π 0 (Q 2 ) stems from the fact that the QCD matrix element q(x p)q(x p)|j em µ |γ(p ) itself is free of the γ 5 ambiguity in dimensional regularization and both the NDR and HV prescriptions can be employed to construct QCD factorization theorems for hard processes provided that the corresponding matching coefficients are computed in an appropriate way without overlooking the potential evanescent operators.
The renormalization scale independence of the factorization formula (2.33) can be readily verified by making use of the evolution equation of the pion DA φ π (x, µ)

JHEP12(2017)037
where the evolution kernel V (x, y) can be expanded perturbatively in QCD with the "plus" function defined as It is appropriate to point out that the one-loop evolution kernel V 0 (x, y) is independent of the γ 5 prescription in the complex D-dimensional space, however, the two-loop evolution kernel V 1 (x, y) does depend on the renormalization scheme. Applying the renormalizationgroup (RG) equation (2.37) then leads to We further turn to sum the parametrically large logarithms of Q 2 /µ 2 in the shortdistance function at next-to-leading-logarithmic (NLL) accuracy employing the standard RG approach in the momentum space. Technically, the desired NLL resummation can be readily achieved by setting the factorization scale of order µ ∼ Q 2 and evolving the leading twist pion DA up to that scale at two loops. The NLO evolution kernel V 1 (x, y) was first obtained with the diagrammatic approach in the light-cone gauge [57,58], then in the Feynman gauge [59,60], and was further reconstructed [61] based upon the knowledge of the conformal anomalies and the available forward DGLAP splitting functions at O(α 2 s ). The two-loop evolution potential V 1 (x, y) can be organized as 42) where N f is the number of the active quark flavours. The explicit expressions of the kernels V N , V G and V F are given by [62]

JHEP12(2017)037
where we have introduced the functions F (x, y) and H(x, y) as follows In order to perform the NLL QCD resummation, it turns out to be convenient to adopt the Gegenbauer expansion of the pion DA where the LO coefficient a 0 (µ) = 1 is renormalization invariant due to the normalization condition. The exact solution to the two-loop RG equation (2.37) can be constructed from the forward anomalous dimensions and the special conformal anomaly matrix in the Gegenbauer moment space [63,64], and we obtain (see also [22]) where both n and k are non-negative even integers and the explicit expressions of E NLO V,n and d k V,n are collected in appendix A. Inserting (2.48) into (2.33) and employing the technique developed in [65], we obtain where the hard coefficient C n (Q 2 , µ) in the NDR scheme of γ 5 is given by with the harmonic number defined as H n = Σ k=n k=1 1/k.

The subleading-power correction from photon DAs
In this section we aim at evaluating the power suppressed contribution to the pion-photon form factor due to the hadronic photon effect at O(α s ) with the LCSR approach. To this end, we construct the following vacuum-to-photon correlation function

JHEP12(2017)037
defined with an electromagnetic current (2.2) carrying a four-momentum q µ and a pion interpolating current j π whose explicit structure is as follows Here we have introduced the convention ⊥ µναβ ≡ g ρ ⊥ µ ρναβ . Following the standard strategy, the primary task for the sum-rule construction consists in the demonstration of QCD factorization for the considered correlation function (3.1) at space-like interpolating momentum p = p + q. Due to the appearance of γ 5 in the QCD matrix element (3.1), the nonminimal renormalization of the axial current (dependent on the γ 5 prescription) will be performed to restore the Ward identities of chiral gauge theories order by order in perturbation theory (see the discussion of subsection 3.2 for the details). We will employ both the NDR and HV schemes of the Dirac matrix γ 5 to establish the QCD factorization formula of the transition amplitude (3.1) at O(α s ) and then derive the NLL resummation improved LCSR for the hadronic photon correction to the π 0 γ * γ form factor. Furthermore, the power counting rule for the external momenta will be adopted to determine the perturbative matching coefficient entering the factorization formula of G µ (p , q) to the one-loop order.

The hadronic photon effect at tree level
QCD factorization for the correlation function (3.1) at tree level can be established by investigating the following four-point QCD amplitude at LO in α s . Evaluating the two diagrams in figure 4 leads to where r = −p 2 /Q 2 , η u = 1 and η d = −1. The partonic matrix element of the (anti)collinear SCET operator O A,µ at tree level is given by The explicit definition of the (anti)-collinear operator O j,µ in the momentum space is where we have suppressed the flavour indices of O j,µ for brevity, Γ A, µ = γ 5 n γ µ,⊥ and the corresponding Wilson line is defined as To achieve the hard-collinear factorization for the correlation function (3.1), we introduce the SCET operator basis where O E,µ is an evanescent operator. Applying the operator matching equation including the evanescent operator and expanding all quantities to the tree level, we can readily find that Making use of the leading twist photon DA defined in [35] 0|χ(0) Wc(0, y) σ αβ χ(y)|γ(p ) the tree-level factorization formula of the form factor G(p 2 , Q 2 ) can be written as where the magnetic susceptibility of the quark condensate χ(µ) encodes the dynamical information of the QCD vacuum [66]. Applying the standard definition for the pion decay constant we can write down the hadronic dispersion relation of G(p 2 , Q 2 ) The final expression of the LCSR for the hadronic photon correction to the pion-photon form factor can then be derived by implementing the continuum subtraction and the Borel transformation with the aid of the parton-hadron duality with u 0 = Q 2 /(s 0 +Q 2 ). Employing the power counting scheme for the sum rule parameters we can obtain the scaling behaviour of the hadronic photon effect at large Q 2

The hadronic photon effect at one loop
To construct the NLL LCSR for the hadronic photon effect, we first need to establish the one-loop factorization formula for the correlation function (3.1) at the leading power in 1/Q 2 . Following the strategy for demonstrating QCD factorization of the leading power contribution presented in section 2, the perturbative matching coefficient entering the factorization formula of the form factor G(p 2 , Q 2 ) can be determined by evaluating the oneloop diagrams for the QCD amplitude Π µ (3.4) in figure 5. We will compute the hard JHEP12(2017)037 contributions from these diagrams one-by-one with both the NDR and HV schemes applying the strategy of regions. The one-loop QCD correction to the electromagnetic vertex diagram displayed in figure 5(a) is free of the γ 5 ambiguity and a straightforward calculation yields where the term proportional to O E,µ (x, x ) is not shown explicitly. Due to the appearance of γ 5 in the pion interpolating current, the one-loop QCD correction to the pion vertex diagram depends on the γ 5 prescription employed in the reduction of the Dirac algebra.
Computing the hard contribution from the one-loop diagram 5(b) in both the NDR and HV schemes gives where Π (0a) µ represents the tree-level contribution to the diagram 4(a) and can be obtained from (3.5) by keeping only the first term in the square bracket. The self-energy correction to the intermediate hard propagator displayed in figure 5(c) is independent of the γ 5 prescription in the D-dimensional space and we obtain We finally turn to compute the hard contribution from the one-loop box diagram shown in figure 5(d), which depends on the actual prescription of γ 5 adopted in the calculation of the corresponding QCD amplitude. Evaluating the contribution from the box diagram with both the NDR and HV schemes, we find that the corresponding hard coefficients only contribute at O( ), vanishing in four-dimensional space. Explicitly,

JHEP12(2017)037
Collecting different pieces together, we derive the NLO QCD correction to the four-point amplitude Π µ as follows (3.24) where the renormalization prescription dependent hard amplitude A 1,hard is given by

25)
Applying the strategy to implement the IR subtraction for the four-point QCD amplitude Π µ discussed in section 2, the master formula for the one-loop hard coefficient of the physical SCET operator O 1,µ can be written as where the bare matrix element M represents the QCD mixing of the evanescent operator O E,µ into O 1,µ at one loop. It is obvious that the infrared subtraction term T (0) suffers from the γ 5 ambiguity in dimensional regularization. The corresponding SCET diagrams at one loop are in analogy to that displayed in figure 3, but with the vertex "⊗" indicating an insertion of O E,µ . Computing these effective diagrams with dimensional regularization applied to the UV divergences and with the IR singularities regularized by the fictitious gluon mass, we find that M (1)off E1 vanishes at one loop with the NDR scheme of γ 5 and it receives a nonvanishing contribution of O( ) with the HV scheme of γ 5 from the effective diagram with a collinear gluon exchange between two external quarks. We are then led to conclude that (3.28) for both the NDR and HV schemes of the γ 5 matrix, with A 1,hard presented in (3.25) and (3.26). We mention in passing that the γ 5 scheme dependence of the short-distance function T (1) 1 will not be cancelled by the one-loop QCD correction to the twist-2 photon DA defined by the light-cone matrix element of the tensor current, which is clearly free of the γ 5 ambiguity in dimensional regularization, and the γ 5 ambiguity of A 1,hard can be traced back to the renormalization prescription dependence of the QCD amplitude (3.4) itself.

JHEP12(2017)037
To preserve the one-loop character of the axial anomaly, an additional finite counterterm must be introduced [11] when performing the UV renormalization of the pseudoscalar current in the HV scheme. Making use of (3.11), (3.25), (3.26) and (3.29), it is then straightforward to verify that at one loop, which provides a nontrivial check to justify the obtained one-loop hard amplitude T 1 . The NLO factorization formula for the vacuum-to-photon correlation function can be further derived as follows at hand, we can also obtain the one-loop short-distance function entering the factorization formula of the H → J/ψ γ form factor at leading power in 1/m 2 H by taking the r → ∞ limit of T 1 and by performing the analytical continuation in the variable p 2 , which reproduces the expression displayed in (3.17) of [67] (see also [68,69]) computed from an alternative approach precisely.
We are now in a position to demonstrate the factorization-scale independence of (3.32) by employing the RG equation of the leading twist photon DA with the perturbative expansion of the evolution kernel The one-loop renormalization kernel V 0 (x, y) is given by [70] V 0 (x, y) Taking into account the factorization scale dependence of T 1 (x, µ), we can further deduce (3.36) The residual µ-dependence of G(p 2 , Q 2 ) originates from the UV renormalization of the QCD pseudoscalar current defining the correlation function (3.1). Taking advantage of the evolution equation of the QCD renormalization constant for the pseudoscalar current [71,72]

JHEP12(2017)037
and distinguishing the renormalization scale of the QCD current from the factorization scale due to the IR subtraction (see [73] for more details), we can then find that the expression for the form factor G(p 2 , Q 2 ) (3.32) is indeed factorization-scale invariant at O(α s ). We proceed to perform the NLL resummation for the parametrically large logarithms in the short-distance function T 1 , which can be achieved alternatively by fixing the factorization scale as µ ∼ Q 2 and by evolving the twist-2 photon DA from the hadronic scale to that scale. To this end, we need the two-loop coefficient of the evolution kernel V (x, y) [61, 74, 75] where the explicit expressions of the kernel functions are [74] Applying the Gegenbauer expansion of the twist-2 photon DA [35] φ γ (x, µ) = 6 xx

JHEP12(2017)037
and implementing the conformal consistency relation discussed in [63], the two-loop evolution of the Gegenbauer moment b n (µ 0 ) can be constructed as follows with even k, n ≥ 0. The detailed expressions the evolution functions E NLO T,n and d k T,n can be found in appendix A. Combining everything together we arrive at the NLL resummation improved factorization formula where the perturbative matching coefficient C n (Q 2 , µ) is defined by 1 (x, µ) We will not present the analytical result of C n (Q 2 , µ) by evaluating the appeared convolution integral explicitly, since the continuum subtraction needs to be performed for the dispersion representation of (3.47) in order to construct the desired LCSRs for the hadronic photon correction to the pion-photon form factor. Employing the spectral representations of the convolution integrals displayed in appendix B, it is straightforward to derive the dispersion form of the NLL factorization formula where we have exploited the symmetric property of the photon DA φ γ (x, µ) = φ γ (x, µ) due to the charge-parity conservation. The resulting QCD spectral densities ρ (i) (s, Q 2 ) (i = 0, 1) can be written as

JHEP12(2017)037
where P indicates the principle-value prescription. The NLL LCSRs for the subleading power contribution to the π 0 γ * γ form factor can be further derived as Collecting different contributions together, we now present the final expression for the pion-photon form factor including the twist-4 correction computed in [22,26] where the explicit expressions of F LP γ * γ→π 0 and F NLP γ * γ→π 0 are displayed in (2.50) and (3.51), respectively. The obtained factorization formula of F tw−4 γ * γ→π 0 from both the two-particle and the three-particle pion DAs at tree level reads where the definition of the twist-4 pion DA F π can be found in (38) of [22] and keeping only the leading conformal spin (i.e., "S"-wave) contribution we obtain The nonperturbative parameter δ 2 π is defined by the local QCD matrix element 0|g sqGµν γ ν q|π(p) = i f π δ 2 π (µ) p µ , (3.55) with the renormalization-scale evolution at one loop Several comments on the general structure of the π 0 γ * γ form factor (3.52) are in order.
• It is apparent that the twist-four correction to the π 0 γ * γ form factor is suppressed by a factor of δ 2 π /Q 2 compared with that of the leading twist contribution. Such subleading power contribution turns out to be numerically significant at Q 2 ≤ 5 GeV 2 due to the large prefactor "80/3" entering the asymptotic expression of F π (x, µ), however, it is still far from sufficient to generate the scaling violation at Q 2 ∼ 40 GeV 2 indicated by the BaBar measurement [21]. Furthermore, it is of high interest to compute the NLO correction to the twist-four contribution in order to develop a better understanding of factorization properties of the high twist effects, where the infrared subtraction for constructing the factorization formula is complicated by the mixing of different twist-four pion DAs under the QCD renormalization.

JHEP12(2017)037
• The twist-six correction to the pion-photon form factor computed from the dispersion approach [22] is partially absorbed into the hadronic photon effect F NLP γ * γ→π 0 (Q 2 ) displayed in (3.51). The precise correspondence of distinct contributions in two frameworks cannot be established without identifying the operator-level definitions of the "soft" corrections in [22], which originate from the nonperturbative modification of the QCD spectral density appeared in the dispersion form of the π 0 γ * γ * form factor.
• The subleading power corrections from the yet higher twist pion/photon DAs, which are not taken into account in this work, are conjectured to be suppressed by only one power of Λ 2 /Q 2 due to the absent correspondence between the twist counting and the large-momentum expansion [22]. A direct calculation of the two-particle and three-particle corrections to the pion-photon form factor from the twist-three and twist-four photon DAs based upon the LCSR approach is in demand to verify this interesting hypothesis.
• We do not include the NNLO QCD correction to the leading power contribution in the large β 0 approximation [8] on account of the absence of a complete NNLO contribution, which also necessitates the three-loop evolution equation of the twisttwo pion DA [76] to obtain the factorization formula at the next-to-next-to-leadinglogarithmic accuracy. A recent discussion of the NNLO radiative corrections in the framework of the dispersion approach can be found in [30].

Numerical analysis
We are now ready to explore the phenomenological consequences of the hadronic photon correction to the pion-photon form factor applying the master formula (3.52). In doing so, we will first need to specify the non-perturbative models for the twist-2 pion DA, the magnetic susceptibility χ(µ), the Gegenbauer moments of the photon DA, and to determine the "internal" sum rule parameters entering the expression (3.51).

Theory input parameters
The fundamental ingredients entering the NLL factorization formula of the leading power contribution are the Gegenbauer moments of the twist-2 pion DA. Tremendous efforts have been devoted to the determinations of the lowest moment a 2 (µ) from the direct calculations with the QCD sum rules pioneered by Chernyak and Zhitnitsky (CZ) [77] and with the lattice simulations, and from the indirect calculations by matching the LCSR predictions with the experimental data. To quantify the systematic uncertainty from the Gegenbauer moments, we will consider the following four models for the leading twist pion DA

JHEP12(2017)037
The obtained Gegenbauer coefficients in the Bakulev-Mikhailov-Stefanis (BMS) model [30,78] are computed from the QCD sum rules with non-local condensates absorbing the high-order terms in the operator-product-expansion (OPE) partially (see, however, [79]). The first and second nontrivial Gegenbauer moments of the KMOW model [80] are determined by comparing the LCSR predictions for the pion electromagnetic form factor, including the NLO correction to the twist-2 effect and the subleading terms up to twist-6, with the intermediate-Q 2 data from the JLab experiment. The holographic model of the twist-2 pion DA [81] φ Hol is motivated by the correspondence between the string theory in the five-dimensional antide Sitter space and conformal field theories in the physical space-time (see also [82] for a similar end-point behaviour of the pion DA) and implementing the Gegenbauer expansion of φ Hol π (x, µ 0 ) leads to the expression of a n displayed in (4.1). For the phenomenological analysis of the π 0 γ * γ form factor, we will truncate the expansion of the "holographic" model at n = 12, which was demonstrated to be a good approximation in [22]. It needs to be pointed our that the values of the second Gegenbauer coefficient in the first three models of (4.1) are in line with the recent lattice determinations [83] within the theory uncertainties and the CZ model is introduced for the illustration purpose to understand the model dependence of the predictions for the pion-photon form factor.
The normalization parameter for the twist-four pion DAs will be taken as δ 2 π (1 GeV) = (0.2 ± 0.04) GeV 2 computed from the QCD sum rules [84] (see also [85]). We further adopt the value of the quark condensate density qq (1 GeV) = − 256 +14 −16 MeV 3 determined in [80]. A key nonperturbative quantity appearing in the twist-2 photon DA is the magnetic susceptibility of the quark condensate χ(µ) describing a response of the QCD vacuum in the presence of an external photon field. Different QCD-based approaches have been proposed to evaluate χ(µ) (see, e.g., [35,86,87]) with the aid of the resonance information from the experimental data and the interval χ(1 GeV) = (3.15 ± 0.3) GeV −2 [35] will be employed in the numerical calculations. In contrast, our understanding of the higher Gegenbauer moments of the leading twist photon DA is rather limited, even for the leading non-asymptotic correction due to b 2 (µ 0 ). The available information of the second Gegenbauer coefficient mainly comes from the QCD sum rules constructed from the correlation function with a light-ray tensor operator and a local vector current, which unfortunately give rise to the theory predictions sensitive to the choice of the input parameters. The crude estimate b 2 (1 GeV) = 0.07 ± 0.07 from [85] will be used in our numerical analysis and an independent determination from the lattice QCD calculation will be very welcome in the future. A natural choice of the factorization scale in the leading power factorization formula (2.50) is µ 2 = x Q 2 with 1/4 ≤ x ≤ 3/4 corresponding to the characteristic virtuality of the intermediate quark displayed in figure 1(a), and it will be frozen at µ = 1 GeV for x Q 2 < 1 GeV 2 at low Q 2 in order not to run into the nonperturbative QCD regime (see [8] for the discussion about the BLM proposal). Along similar lines, the factoriza-JHEP12(2017)037 Figure 6. Left: distinct contributions to the π 0 γ * γ form factor from the twist-two pion DA ("LP") with the BMS model at NLL, from the hadronic photon effect ("NLP") at NLL and from the twistfour pion DAs ("tw-4") at LO. The solid curve is obtained by adding up the above-mentioned three pieces together with the central inputs. Right: dependence of perturbative QCD corrections to the leading power contribution and to the hadronic photon effect with the BMS model, at LL and NLL accuracy, on the momentum transfer accessible at the BaBar and Belle experiments. tion scale entering the NLL LCSRs for the hadronic photon effect (3.51) will be taken as µ 2 = x M 2 + x Q 2 as widely employed in the sum rule calculations [22].
Finally, the determination of the Borel mass M 2 and the threshold parameter s 0 can be achieved by applying the standard strategies described in [88,89], and we can readily obtain in agreement with the intervals adopted in [90].

Predictions
for the π 0 γ * γ form factor Now we will turn to investigate the phenomenological significance of distinct terms contributing to the pion-photon form factor. Taking the BMS model for the twist-two pion DA as an example, it is apparent from figure 6 that the twist-four correction and the hadronic photon contribution generate the destructive and constructive interference with the leading power effect (a similar observation for the high twist corrections already made in [22]) and there appears to be a strong cancellation between these two mechanisms in the whole Q 2 ≤ 40 GeV 2 region. 1 However, both subleading power effects become rapidly suppressed with the growing of the momentum transfer squared in contrast to the numerically sizeable soft power correction estimated from the dispersion approach [22]. Such discrepancy may be ascribed to the very definition of the "soft" effect in the formalism of [26], roughly corresponding to the ρ-resonance contribution to the π 0 γ * γ form factor with the partonhadron duality approximation, which has no transparent counterpart in the framework of JHEP12(2017)037 perturbative QCD factorization. In addition, the NLL radiative corrections are observed to give rise to approximately O(15 %) (almost Q 2 -independent) shift to the LL predictions for both the leading power contribution and the hadronic photon effect.
To understand the phenomenological impact of the QCD resummation for the large logarithms appearing in the factorization formula for the leading power contribution and in the LCSRs for the hadronic photon correction, we further present in figure 7 our predictions for the π 0 γ * γ form factor, at LL, NLO and NLL accuracy, with the BMS model. The NLO QCD corrections are found to induce O (25 %) reduction of the tree-level results at 10 GeV 2 ≤ Q 2 ≤ 40 GeV 2 , however, the NLL resummation effect will enhance the NLO predictions by an amount of approximately O (10 %), in accordance with the pattern for the perturbative QCD corrections observed in [32,88]. Inspecting the model dependence of pion-photon form factor on the leading twist pion DA displayed in figure 7 implies that the theory predictions with both the holographic and KMOW models can reasonably balance the BaBar and Belle data at high Q 2 without resorting to the "exotic" end-point behaviour as advocated in [51,52]. In fact, we have checked that the predicted π 0 γ * γ form factor with the flat pion DA will overshoot both the BaBar and Belle data, in most Q 2 region of interest, at least in our framework. Given the fact that the end-point behaviour of the twisttwo pion DA in the holographic model differs from the standard postulation, motivated by the conformal expansion analysis, as employed for the KMOW model, we conclude that the local information of the pion DA cannot be extracted from the experimental measurements of the pion-photon form factor even in the leading power approximation. It needs further to be pointed out that our predictions with the holographic and KMOW models do not match the experimental data at 2 GeV 2 ≤ Q 2 ≤ 8 GeV 2 well, where the JHEP12(2017)037 power suppressed contributions from the yet higher-twist pion and photon DAs will become more pronounced and actually the large-momentum expansion applied for the construction of the factorization formula also becomes questionable. By contrast, the theory predictions from the dispersion approach [22,23] can result in a satisfactory description of the BaBar and Belle data in the whole Q 2 region by introducing the nonperturbative modification of the QCD spectral density function. Moreover, it becomes apparent that the computed pionphoton form factor with the BMS model and the asymptotic pion DA are less favorable by the experimental measurements at high Q 2 , albeit with the reasonable agreement achieved at low Q 2 . Also, confronting the theory predictions from the CZ model with the BaBar and Belle data indicates a large value of the second Gegenbauber moment a 2 (µ 0 ) is not favored, in agreement with the recent lattice QCD calculations [83,92]. It is worthwhile to mention that the pion-photon transition form factor depends on the strong interaction dynamics encoded by the higher-twist pion DAs and the photon DAs (including the uncalculated effects from higher Fock states), hence measuring F γ * γ→π 0 (Q 2 ) cannot be simply considered as a way of determining the higher Gegenbauer moments of the leading twist pion DA in the kinematic region accessible at the current experiments.
We present our final predictions for the π 0 γ * γ form factor from the expression (3.52) with three different models of the twist-two pion DA in figure 8, including the theory uncertainties due to the variations of the input parameters discussed before. To account for the yet higher-twist corrections discussed in the end of section 3, we have doubled the theory uncertainties of the twist-4 term F tw−4 γ * γ→π 0 (Q 2 ) for a crude estimate. In addition, we assume that the missing NNLO QCD correction can be (partially) taken into account by varying the factorization scale µ in the interval as specified above, and a direct calculation of the two-loop hard function T (2) 2 (x, µ) from the operator matching equation (2.13) will be in demand for an unambiguous analysis. We further assigned 20% uncertainty for the first six nontrivial Gegenbauer coefficients of the holographic model in the numerical estimation for the illustration purpose. It turns out that the dominant theory uncertainties originate from the shape parameters of the leading twist pion and photon DAs instead of the variations of the factorization scales. Precision determinations of the higher Gegenbauer coefficients for both two DAs along the lines of [83,92] will be essential to pin down the presently sizeable theory uncertainty in order to meet the challenge of the (potentially) more accurate experimental measurements at the BEPCII collider [93] and the SuperKEKB accelerator.
Concluding this section we make a few comments on the earlier QCD calculations of the pion-photon form factor in the literature.
• Applying the dispersion approach the subleading power "soft" contribution to the π 0 γ * γ form factor was demonstrated in [22,23] to be essential to accommodate the experimental data at low Q 2 . However, the absence of the operator-level definition for the numerically prominent soft correction makes it difficult to trace back the corresponding dynamical mechanism in the SCET framework. It was further argued in [22,23] that measuring the pion-photon form factor alone is not sufficient to determine the shape of the leading-twist pion DA unambiguously, in accordance with our observation made above.  [21] (orange circles) and Belle [47] (brown spades) are also displayed here.
• Implementing the threshold resummation for the enhanced logarithms (α s (µ) ln x) n in TMD factorization, the resulting theory predictions provide a reasonable description of the CLEO and BaBar data at Q 2 ≤ 40 GeV 2 on account of the Q 2 -dependent threshold function [94]. However, the importance of the small x region would invalidate QCD factorization for the π 0 γ * γ form factor, since the power counting scheme for the momentum fraction x ∼ O(1) cannot be justified consequently.
• Employing the LCSR approach developed in [26], a systematic analysis of the π 0 γ * γ form factor at 1 GeV 2 ≤ Q 2 ≤ 5 GeV 2 with an emphasis on reducing the theoretical uncertainties from a wide variety of sources was carried out in [30], and the obtained JHEP12(2017)037 theory predictions were found to comply with the available experimental data. Further development along this line can be achieved by working out a complete NNLO correction to the leading twist contribution and by performing an NLO computation for the twist-4 correction on the theoretical aspects.

Conclusion
Applying the standard OPE technique with the evanescent operator(s) we revisited the demonstration of QCD factorization for the pion-photon transition form factor at leading power in 1/Q 2 with both the NDR and HV schemes for γ 5 in the D-dimensional space. It has been shown explicitly at one loop that the renormalization scheme dependence of the short-distance matching coefficient and the twist-two pion DA are cancelled out precisely rendering the γ 5 -prescription independence of the factorization formula for the leading power contribution to F γ * γ→π 0 (Q 2 ). This can be readily understood from the fact that the QCD matrix element defined by two electromagnetic currents are free of the γ 5 ambiguity and the renormalization scheme dependence of the hard function arises from the infrared subtraction term completely. In the same vein, we established QCD factorization of the desired correlation function at one loop for the construction of the LCSRs for the hadronic photon contribution to the pion-photon form factor. By contrast, the corresponding QCD matrix element defined with an interpolating current for the pion and an electromagnetic current suffers from the γ 5 ambiguity, which can be resolved by applying the Ward identities of chiral gauge theories. The finite renormalization term introduced in the HV scheme to restore the appropriate Ward-Takahashi identities was found to provide the very transformation function to construct the hard matching coefficient in the NDR scheme. The NLL resummation of the parametrically large logarithms was also implemented by solving the relevant two-loop evolution equations in momentum space. Taking into account the leading power contribution and the hadronic photon effect at NLL as well as the twist-four correction from the pion DAs at tree level, we further explored the phenomenological consequence of the perturbative QCD corrections and the subleading power contributions. Interestingly, the observed strong cancellation between the two power suppressed mechanisms leads to the insignificant correction to the leading power contribution (almost) in the whole Q 2 region accessible at the current experiments. In addition, we paid a particular attention to the model dependence of the theory predictions for F γ * γ→π 0 (Q 2 ) on the twist-two pion DA. Both the holographic and KMOW models turned out to balance the BaBar and Belle data reasonably well at high Q 2 , despite the visible discrepancy at low Q 2 which could be compensated by the unaccounted subleading power corrections of both perturbative and nonperturbative origins. It was also demonstrated that the end-point behaviour of the pion DA cannot be extracted by matching the theory predictions for the pion-photon form factor with the experimental measurements.
Aiming at a better confrontation with the BaBar and Belle data, further improvements of our calculations can be made by first carrying out the perturbative correction to the twist-4 contribution from both the two-particle and three-particle pion DAs, which is also of conceptual interest in the framework of perturbative QCD factorization, and then JHEP12(2017)037 by evaluating the high twist contributions from the photon DAs with the LCSR approach. Phenomenological applications of the techniques discussed in this work can be also pursued in the context of the γ * γ → (η ( ) , η c ) transition form factors [33,95] for understanding the quark-gluon structure of eta mesons and heavy quarkonium states, the radiative leptonic B-meson decays for the determination of the inverse moment λ B , the radiative penguin decays of B-mesons for the precision test of the quark-flavour structure of the Standard Model, and the radiative heavy-hadron decays for constraining the magnetic susceptibility of the quark condensate [86]. To conclude, the anatomy of the subleading power contributions for the exclusive hadronic reactions is of high interest for understanding the general structures of the large momentum/mass expansion in QCD and for hunting new physics in the quark-flavour sector as indicated by the various flavour "anomalies" observed at the ongoing experiments.

Acknowledgments
We are grateful to Martin
(A. 16) The explicit expressions of the RG functions E NLO T,n and d k T,n can be obtained from that of E NLO V,n and d k V,n given above with the replacement rule γ

B Spectral representations
We present the dispersion representations of convolution integrals entering the NLL QCD factorization formula (3.32) in order to construct the sum rules for the hadronic photon correction to the pion-photon form factor. We have verified the spectral representations in what follows numerically by checking the corresponding dispersion integrals.

JHEP12(2017)037
Here, the parameter p 2 in the definition of r should be understood as s in the above convolution integrals and P represents the principle-value prescription.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.