Moments Of Spectral Functions Research Articles

Reconciling the contour-improved and fixed-order approaches for τ hadronic spectral moments. Part II. Renormalon norm and application in αs determinations

In a previous article, we have shown that the discrepancy between the fixed-order (FOPT) and contour-improved (CIPT) perturbative expansions for τ hadronic spectral function moments, which had affected the precision of αs determinations for many years, may be reconciled by employing a renormalon-free (RF) scheme for the gluon condensate (GC) matrix element. In addition, the perturbative convergence of spectral function moments with a sizeable GC correction can be improved. The RF GC scheme depends on an IR factorization scale R and the normalization Ng of the GC renormalon. In the present work, we use three different methods to determine Ng, yielding a result with an uncertainty of 40%. Following two recent state-of-the-art strong coupling determination analyses at mathcal{O} ( {alpha}_s^5 ), we show that using the renormalon-free GC scheme successfully reconciles the results for αs( {m}_{tau}^2 ) based on CIPT and FOPT. The uncertainties due to variations of R and the uncertainty of Ng only lead to a small or moderate increase of the final uncertainty of αs( {m}_{tau}^2 ), and affect mainly the CIPT expansion method. The FOPT and CIPT results obtained in the RF GC scheme may be consistently averaged. The RF GC scheme thus constitutes a powerful new ingredient for future analyses of τ hadronic spectral function moments.

Reconciling the contour-improved and fixed-order approaches for τ hadronic spectral moments. Part I. Renormalon-free gluon condensate scheme

We propose a simple and easy-to-implement scheme for a renormalon-free gluon condensate (GC) matrix element, which is analogous to implementations of short-distance heavy-quark mass renormalization schemes existing in the literature already for a long time. Because the scheme is based on a perturbative subtraction at the level of the matrix element, with a freely adaptable infrared factorization scale, it can be implemented with little effort for any observable where the GC is relevant. The scheme depends on the renormalon norm of the GC which has to be supplemented independently. We apply the scheme to the fixed-order (FOPT) and contour-improved (CIPT) perturbative expansions of τ hadronic spectral function moments. These expansions exhibit a long-standing discrepancy for moments used in high-precision determinations of the strong coupling in the commonly used GC scheme that is not renormalon-free. We show that the scheme is capable of resolving the FOPT-CIPT discrepancy problem. At the same time, the perturbative behaviour of the moments that previously showed bad convergence properties and for which the non-perturbative corrections from the GC are sizeable, is substantially improved. The new GC scheme may provide a powerful theoretical tool for future phenomenological applications.

On the discrepancy between the FOPT and CIPT approaches for hadronic tau spectral function moments

The discrepancy between the FOPT and CIPT approaches for hadronic \tauτ spectral function moments constitutes the major theoretical uncertainty for strong coupling determinations from tau decay data. We show the discrepancy can be analytically understood since the Borel representations – which have been assumed to be identical for both approaches previously – differ in the presence of IR renormalons. This implies that the OPE condensate corrections are different for both approaches and that the discrepancy may eventually be reconciled. In the talk we explain the difference and some mathematical aspects of of the FOPT and CIPT Borel representations and show numerical results.

Borel representation of τ hadronic spectral function moments in contour-improved perturbation theory

We show that the Borel representations of $\ensuremath{\tau}$ hadronic spectral function moments based on contour-improved perturbation theory (CIPT) in general differ from those obtained within fixed-order perturbation theory (FOPT) in the presence of IR renormalons in the underlying Adler function. The Borel sums obtained from both types of Borel representations in general differ as well, and the apparently conflicting behavior of the FOPT and CIPT spectral function moment series at intermediate orders, which has been subject to many studies in the past literature, can be understood quantitatively using concrete Borel function models. The difference between the CIPT and FOPT Borel sums, which we call the ``asymptotic separation,'' can be computed analytically for any Borel function model and is proportional to inverse exponential terms in the strong coupling. Even though moments can be designed where the asymptotic separation is strongly suppressed, it is as a matter of principle unavoidable. If the Borel function of the Euclidean Adler function has a sizeable gluon condensate renormalon cut, the asymptotic separation can explain the observed disparity of the CIPT and FOPT spectral function moments at the 5-loop level. The existence of the asymptotic separation implies that the power corrections in the operator product expansion for the spectral function moments in the CIPT expansion approach do not have the commonly assumed analytic standard form.

Reconciling the FOPT and CIPT Predictions for τ Hadronic Spectral Function Moments

Recently it has been clarified by Hoang and Regner that the longstanding discrepancy between the CIPT and FOPT expansion approaches in αs determinations from the τ hadronic spectral function moments has been caused by an inconsistency of CIPT with the standard OPE approach. This inconsistency arises in the presence of IR renormalons in the underlying Adler function and is numerically dominated by the dimension-4 gluon condensate renormalon. In this talk we report on an approach to reconcile the CIPT based on a perturbative definition of a renormalon-free and scale-invariant gluon condensate scheme, called RF GC scheme. The scheme implies perturbative subtractions which eliminate the CIPT inconsistency for all practical applications of the τ hadronic spectral function moments. The scheme depends on the gluon condensate renormalon norm Ng as an independent input and on an IR subtraction scale R. We discuss three different approaches to determine Ng which yield consistent results and we apply the RF GC scheme in two full-fledged phenomenological αs determinations based on the truncated OPE and the duality violation model approach. In the RF GC scheme the long-standing CIPT-FOPT discrepancy problem is gone and the CIPT and FOPT αs determinations can be consistently combined.

Renormalons in integrated spectral function moments and αs extractions

Precise extractions of $\alpha_s$ from $\tau\to {\rm (hadrons)}+\nu_\tau$ and from $e^+e^-\to {\rm (hadrons)}$ below the charm threshold rely on finite energy sum rules (FESRs) where the experimental side is given by integrated spectral function moments. Here we study the renormalons that appear in the Borel transform of polynomial moments in the large-$\beta_0$ limit and in full QCD. In large-$\beta_0$, we establish a direct connection between the renormalons and the perturbative behaviour of moments often employed in the literature. The leading IR singularity is particularly prominent and is behind the fate of moments whose perturbative series are unstable, while those with good perturbative behaviour benefit from partial cancellations of renormalon singularities. The conclusions can be extended to QCD through a convenient scheme transformation to the $C$-scheme together with the use of a modified Borel transform which make the results particularly simple; the leading IR singularity becomes a simple pole, as in large-$\beta_0$. Finally, for the moments that display good perturbative behaviour, we discuss an optimized truncation based on renormalisation scheme (or scale) variation. Our results allow for a deeper understanding of the perturbative behaviour of integrated spectral function moments and provide theoretical support for low-$Q^2$ $\alpha_s$ determinations.

Tau hadronic spectral function moments: perturbative expansion and αs extractions

In the extraction of αs from hadronic τ decays different moments of the spectral functions have been used. Furthermore, the two mainstream renormalization group improvement (RGI) frameworks, namely Fixed Order Perturbation Theory (FOPT) and Contour Improved Perturbation Theory (CIPT), lead to conflicting values of αs. In order to improve the strategy used in αs determinations, we have performed a systematic study of the perturbative behaviour of these spectral moments in the context of FOPT and CIPT. Higher order coefficients of the perturbative series, yet unknown, were modelled using available knowledge of the renormalon content of the QCD Adler function. We conclude that within these RGI frameworks some of the moments often employed in αs extractions should be avoided due to their poor perturbative behaviour. Finally, under reasonable assumptions about higher orders, we conclude that FOPT is the preferred method to perform the renormalization group improvement of the perturbative series.

Expansions ofτhadronic spectral function moments in a nonpower QCD perturbation theory with tamed large order behavior

The moments of the hadronic spectral functions are of interest for the extraction of the strong coupling $\alpha_s$ and other QCD parameters from the hadronic decays of the $\tau$ lepton. Motivated by the recent analyses of a large class of moments in the standard fixed-order and contour-improved perturbation theories, we consider the perturbative behavior of these moments in the framework of a QCD nonpower perturbation theory, defined by the technique of series acceleration by conformal mappings, which simultaneously implements renormalization-group summation and has a tame large-order behavior. Two recently proposed models of the Adler function are employed to generate the higher order coefficients of the perturbation series and to predict the exact values of the moments, required for testing the properties of the perturbative expansions. We show that the contour-improved nonpower perturbation theories and the renormalization-group-summed nonpower perturbation theories have very good convergence properties for a large class of moments of the so-called "reference model", including moments that are poorly described by the standard expansions. The results provide additional support for the plausibility of the description of the Adler function in terms of a small number of dominant renormalons.

STRONG COUPLING FROM THE TAU HADRONIC WIDTH BY NON-POWER QCD PERTURBATION THEORY

Starting from the divergent character of the perturbative expansions in QCD and using the technique of series acceleration by the conformal mappings of the Borel plane, I define a novel, non-power perturbative expansion for the Adler function, which simultaneously implements renormalisation-group summation and has a tamed large-order behaviour. The new expansion functions, which replace the standard powers of the coupling, are singular at the origin of the coupling plane and have divergent perturbative expansions, resembling the expanded function itself. Confronting the new perturbative expansions with the standard ones on specific models investigated recently in the literature, I show that they approximate in an impressive way the exact Adler function and the spectral function moments. Applied to the τ hadronic width, the contour-improved and the renormalisation-group summed non-power expansions in the [Formula: see text] scheme lead to the prediction [Formula: see text], which translates to [Formula: see text].

Spectral moment sum rules for the retarded Green's function and self-energy of the inhomogeneous Bose-Hubbard model in equilibrium and nonequilibrium

We derive expressions for the zeroth and the first three spectral moment sum rules for the retarded Green's function and for the zeroth and the first spectral moment sum rules for the retarded self-energy of the inhomogeneous Bose-Hubbard model in nonequilibrium, when the local on-site repulsion and the chemical potential are time-dependent, and in the presence of an external time-dependent electromagnetic field. We also evaluate these moments for equilibrium, where all time dependence and external fields vanish. Unlike similar sum rules for the Fermi-Hubbard model, in the Bose-Hubbard model case, the sum rules often depend on expectation values that cannot be determined simply from parameters in the Hamiltonian like the interaction strength and chemical potential, but require knowledge of equal time many-body expectation values from some other source. We show how one can approximately evaluate these expectation values for the Mott-insulating phase in a systematic strong-coupling expansion in powers of the hopping divided by the interaction. We compare the exact moments to moments of spectral functions calculated from a variety of different approximations and use them to benchmark their accuracy.

Perturbative expansion of τ hadronic spectral function moments and α s extractions

Various moments of the hadronic spectral functions have been employed in the determination of the strong coupling α s from tau decays. In this work we study the behaviour of their perturbative series under different assumptions for the large-order behaviour of the Adler function, extending previous work on the tau hadronic width. We find that the moments can be divided into a small number of classes, whose characteristics depend only on generic features of the moment weight function and Adler function series. Some moments that are commonly employed in α s analyses from τ decays should be avoided because of their perturbative instability. This conclusion is corroborated by a simplified α s extraction from individual moments. Furthermore, under reasonable assumptions for the higher-order behaviour of the perturbative series, fixed-order perturbation theory (FOPT) provides the preferred framework for the renormalization group improvement of all moments that show good perturbative behaviour. Finally, we provide further evidence for the plausibility of the description of the Adler function in terms of a small number of leading renormalon singularities.

Duality violations in τ hadronic spectral moments

Evidence is presented for the necessity of including duality violations in a consistent description of spectral function moments employed in the precision determination of $\alpha_s$ from $\tau$ decay. A physically motivated ansatz for duality violations in the spectral functions enables us to perform fits to spectral moments employing both pinched and unpinched weights. We describe our analysis strategy and provide some preliminary findings. Final numerical results await completion of an ongoing re-determination of the ALEPH covariance matrices incorporating correlations due to the unfolding procedure which are absent from the currently posted versions. To what extent this issue affects existing analyses and our own work will require further study.

Recent progress in hadronic τ decays

The determination of alpha_s from hadronic tau decays is impeded by the fact that two choices for the renormalisation group resummation, namely fixed-order (FOPT) and contour-improved perturbation theory (CIPT), yield systematically differing results. On the basis of a model for higher-order terms in the perturbative series, which incorporates well-known structure from renormalons, it is found that FOPT smoothly approaches the Borel sum for the tau hadronic width, while CIPT is unable to account for the resummed series. An example for the behaviour of QCD spectral function moments, displaying a similar behaviour, is presented as well.

Moments of spectral functions: Monte Carlo evaluation and verification

The subject of the present study is the Monte Carlo path-integral evaluation of the moments of spectral functions. Such moments can be computed by formal differentiation of certain estimating functionals that are infinitely differentiable against time whenever the potential function is arbitrarily smooth. Here, I demonstrate that the numerical differentiation of the estimating functionals can be more successfully implemented by means of pseudospectral methods (e.g., exact differentiation of a Chebyshev polynomial interpolant), which utilize information from the entire interval . The algorithmic detail that leads to robust numerical approximations is the fact that the path-integral action and not the actual estimating functional are interpolated. Although the resulting approximation to the estimating functional is nonlinear, the derivatives can be computed from it in a fast and stable way by contour integration in the complex plane, with the help of the Cauchy integral formula (e.g., by Lyness' method). An interesting aspect of the present development is that Hamburger's conditions for a finite sequence of numbers to be a moment sequence provide the necessary and sufficient criteria for the computed data to be compatible with the existence of an inversion algorithm. Finally, the issue of appearance of the sign problem in the computation of moments, albeit in a milder form than for other quantities, is addressed.

The Exponential Accuracy of Fourier and Chebyshev Differencing Methods

It is shown that when differencing analytic functions using the pseudospectral Fourier or Chebyshev methods, the error committed decays to zero at an exponential rate.