QED self-energy contribution to highly-excited atomic states

We present numerical values for the self-energy shifts predicted by QED (Quantum Electrodynamics) for hydrogenlike ions (nuclear charge $60 \le Z \le 110$) with an electron in an $n=3$, 4 or 5 level with high angular momentum ($5/2\le j \le 9/2$). Applications include predictions of precision transition energies and studies of the outer-shell structure of atoms and ions.Comment: 20 pages, 5 figure

O(1) Computation of Legendre polynomials and Gauss-Legendre nodes and weights for parallel computing

A self-contained set of algorithms is proposed for the fast evaluation of Legendre polynomials of arbitrary degree and argument is an element of [-1, 1]. More specifically the time required to evaluate any Legendre polynomial, regardless of argument and degree, is bounded by a constant; i.e., the complexity is O(1). The proposed algorithm also immediately yields an O(1) algorithm for computing an arbitrary Gauss-Legendre quadrature node. Such a capability is crucial for efficiently performing certain parallel computations with high order Legendre polynomials, such as computing an integral in parallel by means of Gauss-Legendre quadrature and the parallel evaluation of Legendre series. In order to achieve the O(1) complexity, novel efficient asymptotic expansions are derived and used alongside known results. A C++ implementation is available from the authors that includes the evaluation routines of the Legendre polynomials and Gauss-Legendre quadrature rules

A fast, simple, and stable Chebyshev-Legendre transform using an asymptotic formula

A fast, simple, and numerically stable transform for converting between Legendre and Chebyshev coefficients of a degree $N$ polynomial in $O(N(\log N)^{2}/ \log \log N)$ operations is derived. The basis of the algorithm is to rewrite a well-known asymptotic formula for Legendre polynomials of large degree as a weighted linear combination of Chebyshev polynomials, which can then be evaluated by using the discrete cosine transform. Numerical results are provided to demonstrate the efficiency and numerical stability. Since the algorithm evaluates a Legendre expansion at an $N+1$ Chebyshev grid as an intermediate step, it also provides a fast transform between Legendre coefficients and values on a Chebyshev grid

Design of quadrature rules for Müntz and Müntz-logarithmic polynomials using monomial transformation

A method for constructing the exact quadratures for Müntz and Müntz-logarithmic polynomials is presented. The algorithm does permit to anticipate the precision (machine precision) of the numerical integration of Müntz-logarithmic polynomials in terms of the number of Gauss-Legendre (GL) quadrature samples and monomial transformation order. To investigate in depth the properties of classical GL quadrature, we present new optimal asymptotic estimates for the remainder. In boundary element integrals this quadrature rule can be applied to evaluate singular functions with end-point singularity, singular kernel as well as smooth functions. The method is numerically stable, efficient, easy to be implemented. The rule has been fully tested and several numerical examples are included. The proposed quadrature method is more efficient in run-time evaluation than the existing methods for Müntz polynomial