Bayesian Model Comparison and Analysis of the Galactic Disk Population of Gamma-Ray Millisecond Pulsars

Pulsed emission from almost one hundred millisecond pulsars (MSPs) has been detected in $\gamma$-rays by the Fermi Large-Area Telescope. The global properties of this population remain relatively unconstrained despite many attempts to model their spatial and luminosity distributions. We perform here a self-consistent Bayesian analysis of both the spatial distribution and luminosity function simultaneously. Distance uncertainties, arising from errors in the parallax measurement or Galactic electron-density model, are marginalized over. We provide a public Python package for calculating distance uncertainties to pulsars derived using the dispersion measure by accounting for the uncertainties in Galactic electron-density model YMW16. Finally, we use multiple parameterizations for the MSP population and perform Bayesian model comparison, finding that a broken power law luminosity function with Lorimer spatial profile are preferred over multiple other parameterizations used in the past. The best-fit spatial distribution and number of $\gamma$-ray MSPs is consistent with results for the radio population of MSPs.Comment: 13 pages, 8 figures, 3 tables + Appendix. Public code and source list available from http://github.com/tedwards2412/MSPDis

Representation of a complex Green function on a real basis: I. General Theory

When the Hamiltonian of a system is represented by a finite matrix, constructed from a discrete basis, the matrix representation of the resolvent covers only one branch. We show how all branches can be specified by the phase of a complex unit of time. This permits the Hamiltonian matrix to be constructed on a real basis; the only duty of the basis is to span the dynamical region of space, without regard for the particular asymptotic boundary conditions that pertain to the problem of interest.Comment: about 40 pages with 5 eps-figure

Acceleration of generalized hypergeometric functions through precise remainder asymptotics

We express the asymptotics of the remainders of the partial sums {s_n} of the generalized hypergeometric function q+1_F_q through an inverse power series z^n n^l \sum_k c_k/n^k, where the exponent l and the asymptotic coefficients {c_k} may be recursively computed to any desired order from the hypergeometric parameters and argument. From this we derive a new series acceleration technique that can be applied to any such function, even with complex parameters and at the branch point z=1. For moderate parameters (up to approximately ten) a C implementation at fixed precision is very effective at computing these functions; for larger parameters an implementation in higher than machine precision would be needed. Even for larger parameters, however, our C implementation is able to correctly determine whether or not it has converged; and when it converges, its estimate of its error is accurate.Comment: 36 pages, 6 figures, LaTeX2e. Fixed sign error in Eq. (2.28), added several references, added comparison to other methods, and added discussion of recursion stabilit

Mathematical Properties of a New Levin-Type Sequence Transformation Introduced by \v{C}\'{\i}\v{z}ek, Zamastil, and Sk\'{a}la. I. Algebraic Theory

\v{C}\'{\i}\v{z}ek, Zamastil, and Sk\'{a}la [J. Math. Phys. \textbf{44}, 962 - 968 (2003)] introduced in connection with the summation of the divergent perturbation expansion of the hydrogen atom in an external magnetic field a new sequence transformation which uses as input data not only the elements of a sequence $\{s_n \}_{n=0}^{\infty}$ of partial sums, but also explicit estimates $\{\omega_n \}_{n=0}^{\infty}$ for the truncation errors. The explicit incorporation of the information contained in the truncation error estimates makes this and related transformations potentially much more powerful than for instance Pad\'{e} approximants. Special cases of the new transformation are sequence transformations introduced by Levin [Int. J. Comput. Math. B \textbf{3}, 371 - 388 (1973)] and Weniger [Comput. Phys. Rep. \textbf{10}, 189 - 371 (1989), Sections 7 -9; Numer. Algor. \textbf{3}, 477 - 486 (1992)] and also a variant of Richardson extrapolation [Phil. Trans. Roy. Soc. London A \textbf{226}, 299 - 349 (1927)]. The algebraic theory of these transformations - explicit expressions, recurrence formulas, explicit expressions in the case of special remainder estimates, and asymptotic order estimates satisfied by rational approximants to power series - is formulated in terms of hitherto unknown mathematical properties of the new transformation introduced by \v{C}\'{\i}\v{z}ek, Zamastil, and Sk\'{a}la. This leads to a considerable formal simplification and unification.Comment: 41 + ii pages, LaTeX2e, 0 figures. Submitted to Journal of Mathematical Physic

An efficient algorithm for accelerating the convergence of oscillatory series, useful for computing the polylogarithm and Hurwitz zeta functions

This paper sketches a technique for improving the rate of convergence of a general oscillatory sequence, and then applies this series acceleration algorithm to the polylogarithm and the Hurwitz zeta function. As such, it may be taken as an extension of the techniques given by Borwein's "An efficient algorithm for computing the Riemann zeta function", to more general series. The algorithm provides a rapid means of evaluating Li_s(z) for general values of complex s and the region of complex z values given by |z^2/(z-1)|<4. Alternatively, the Hurwitz zeta can be very rapidly evaluated by means of an Euler-Maclaurin series. The polylogarithm and the Hurwitz zeta are related, in that two evaluations of the one can be used to obtain a value of the other; thus, either algorithm can be used to evaluate either function. The Euler-Maclaurin series is a clear performance winner for the Hurwitz zeta, while the Borwein algorithm is superior for evaluating the polylogarithm in the kidney-shaped region. Both algorithms are superior to the simple Taylor's series or direct summation. The primary, concrete result of this paper is an algorithm allows the exploration of the Hurwitz zeta in the critical strip, where fast algorithms are otherwise unavailable. A discussion of the monodromy group of the polylogarithm is included.Comment: 37 pages, 6 graphs, 14 full-color phase plots. v3: Added discussion of a fast Hurwitz algorithm; expanded development of the monodromy v4:Correction and clarifiction of monodrom

Lamm, Valluri, Jentschura and Weniger comment on "A Convergent Series for the QED Effective Action" by Cho and Pak [Phys. Rev. Lett. vol. 86, pp. 1947-1950 (2001)]

Complete results were obtained by us in [Can. J. Phys. 71, 389 (1993)] for convergent series representations of both the real and the imaginary part of the QED effective action; these derivations were based on correct intermediate steps. In this comment, we argue that the physical significance of the "logarithmic correction term" found by Cho and Pak in [Phys. Rev. Lett. 86, 1947 (2001)] in comparison to the usual expression for the QED effective action remains to be demonstrated. Further information on related subjects can be found in Appendix A of hep-ph/0308223 and in hep-th/0210240.Comment: 1 page, RevTeX; only "meta-data" update

Global analysis of the pMSSM in light of the Fermi GeV excess: prospects for the LHC Run-II and astroparticle experiments

We present a new global fit of the 19-dimensional phenomenological Minimal Supersymmetric Standard Model (pMSSM-19) that comply with all the latest experimental results from dark matter indirect, direct and accelerator dark matter searches. We show that the model provides a satisfactory explanation of the excess of gamma-rays from the Galactic centre observed by the Fermi~Large Area Telescope, assuming that it is produced by the annihilation of neutralinos in the Milky Way halo. We identify two regions that pass all the constraints: the first corresponds to neutralinos with a mass ~80-100 GeV annihilating into WW with a branching ratio of 95% ; the second to heavier neutralinos, with mass ~180-200 GeV annihilating into t tbar with a branching ratio of 87%. We show that neutralinos compatible with the Galactic centre GeV excess will soon be within the reach of LHC run-II -- notably through searches for charginos and neutralinos, squarks and light smuons -- and of Xenon1T, thanks to its unprecedented sensitivity to spin-dependent cross-section off neutrons.Comment: Minor changes following referee reports. Main conclusions unchanged. Matches version published in JCA