High-precision computation of uniform asymptotic expansions for special functions

In this dissertation, we investigate new methods to obtain uniform asymptotic expansions for the numerical evaluation of special functions to high-precision. We shall first present the theoretical and computational fundamental aspects required for the development and ultimately implementation of such methods. Applying some of these methods, we obtain efficient new convergent and uniform expansions for numerically evaluating the confluent hypergeometric functions and the Lerch transcendent at high-precision. In addition, we also investigate a new scheme of computation for the generalized exponential integral, obtaining on the fastest and most robust implementations in double-precision floating-point arithmetic. In this work, we aim to combine new developments in asymptotic analysis with fast and effective open-source implementations. These implementations are comparable and often faster than current open-source and commercial stateof-the-art software for the evaluation of special functions.Esta tesis presenta nuevos métodos para obtener expansiones uniformes asintóticas, para la evaluación numérica de funciones especiales en alta precisión. En primer lugar, se introducen fundamentos teóricos y de carácter computacional necesarios para el desarrollado y posterior implementación de tales métodos. Aplicando varios de dichos métodos, se obtienen nuevas expansiones uniformes convergentes para la evaluación numérica de las funciones hipergeométricas confluentes y de la función transcendental de Lerch. Por otro lado, se estudian nuevos esquemas de computo para evaluar la integral exponencial generalizada, desarrollando una de las implementaciones más eficientes y robustas en aritmética de punto flotante de doble precisión. En este trabajo, se combinan nuevos desarrollos en análisis asintótico con implementaciones rigurosas, distribuidas en código abierto. Las implementaciones resultantes son comparables, y en ocasiones superiores, a las soluciones comerciales y de código abierto actuales, que representan el estado de la técnica en el campo de la evaluación de funciones especiales

High-precision computation of uniform asymptotic expansions for special functions

In this dissertation, we investigate new methods to obtain uniform asymptotic expansions for the numerical evaluation of special functions to high-precision. We shall first present the theoretical and computational fundamental aspects required for the development and ultimately implementation of such methods. Applying some of these methods, we obtain efficient new convergent and uniform expansions for numerically evaluating the confluent hypergeometric functions and the Lerch transcendent at high-precision. In addition, we also investigate a new scheme of computation for the generalized exponential integral, obtaining on the fastest and most robust implementations in double-precision floating-point arithmetic. In this work, we aim to combine new developments in asymptotic analysis with fast and effective open-source implementations. These implementations are comparable and often faster than current open-source and commercial stateof-the-art software for the evaluation of special functions.Esta tesis presenta nuevos métodos para obtener expansiones uniformes asintóticas, para la evaluación numérica de funciones especiales en alta precisión. En primer lugar, se introducen fundamentos teóricos y de carácter computacional necesarios para el desarrollado y posterior implementación de tales métodos. Aplicando varios de dichos métodos, se obtienen nuevas expansiones uniformes convergentes para la evaluación numérica de las funciones hipergeométricas confluentes y de la función transcendental de Lerch. Por otro lado, se estudian nuevos esquemas de computo para evaluar la integral exponencial generalizada, desarrollando una de las implementaciones más eficientes y robustas en aritmética de punto flotante de doble precisión. En este trabajo, se combinan nuevos desarrollos en análisis asintótico con implementaciones rigurosas, distribuidas en código abierto. Las implementaciones resultantes son comparables, y en ocasiones superiores, a las soluciones comerciales y de código abierto actuales, que representan el estado de la técnica en el campo de la evaluación de funciones especiales.Postprint (published version

Numerical and analytical models of self-force effects in Kerr extreme-mass-ratio inspirals

Ground-based detectors now regularly observe the merger of stellar-mass compact objects and their gravitational waves. Building on this success, ESA, in partnership with NASA, will launch the space-based LISA observatory to detect milli-Hertz gravitational wave signals. Extreme-mass-ratio inspirals (EMRIs)---binaries composed of a stellar-mass compact object orbiting a massive black hole---are ideal gravitational wave sources for LISA. Because of their unique properties, EMRIs can provide new insights concerning the growth of massive black holes (and their host galaxies) and enable the most precise tests of general relativity. To achieve this science, LISA will rely on accurate EMRI models to search for and analyze gravitational wave signals. The most accurate EMRI models rely on a mechanism known as the gravitational self-force to calculate an EMRI inspiral and the resulting gravitational waveform. For EMRIs with rotating (Kerr) massive black holes, current gravitational self-force calculations are too computationally demanding to be incorporated into full EMRI models. For my dissertation, I built a developmental scalar self-force code to devise and implement new numerical and analytical techniques for calculating self-force effects in Kerr spacetime. I introduce spectral techniques for numerically evaluating Kerr geodesics and the sources of scalar perturbations. I discuss how these methods can be extended to gravitational self-force calculations. With this code, I produced the first calculations of the scalar self-force along resonant and non-resonant inclined, eccentric orbits in Kerr spacetime. With these new resonant calculations I provided one of the first tests of the integrability conjecture, which holds for these scalar self-force results. I also uncovered the existence of a physical effect in EMRI waveforms, now referred to as quasinormal bursts. Quasinormal bursts are periodic high-frequency oscillations in EMRI waveforms which may aid in the characterization of EMRI gravitational wave sources.Doctor of Philosoph

Probing Particle Physics with Gravitational Waves

The direct detection of gravitational waves offers an exciting new window onto our Universe. At the same time, multiple observational evidence and theoretical considerations motivate the presence of physics beyond the Standard Model. In this thesis, we explore new ways of probing particle physics in the era of gravitational-wave astronomy. We focus on the signatures of ultralight bosons on the gravitational waves emitted by binary systems, demonstrating how binary black holes are novel detectors of this class of dark matter. We also discuss probes of other types of new physics through their finite-size imprints on gravitational waveforms, and examine the extent to which current template-bank searches could be used to detect these signals. In the first two chapters of this thesis, we review several aspects of gravitational-wave physics and particle physics at the weak coupling frontier; we hope the reader would find these reviews helpful in delving further into the literature and in their research.Comment: PhD Thesis, University of Amsterdam, 2020; 298 page

Mathematical Economics

This book is devoted to the application of fractional calculus in economics to describe processes with memory and non-locality. Fractional calculus is a branch of mathematics that studies the properties of differential and integral operators that are characterized by real or complex orders. Fractional calculus methods are powerful tools for describing the processes and systems with memory and nonlocality. Recently, fractional integro-differential equations have been used to describe a wide class of economical processes with power law memory and spatial nonlocality. Generalizations of basic economic concepts and notions the economic processes with memory were proposed. New mathematical models with continuous time are proposed to describe economic dynamics with long memory. This book is a collection of articles reflecting the latest mathematical and conceptual developments in mathematical economics with memory and non-locality based on applications of fractional calculus

Dynamics of a class of vortex rings

The contour dynamics method is extended to vortex rings with vorticity varying linearly from the symmetry axis. An elliptic core model is also developed to explain some of the basic physics. Passage and collisions of two identical rings are studied focusing on core deformation, sound generation and stirring of fluid elements. With respect to core deformation, not only the strain rate but how rapidly it varies is important and accounts for greater susceptibility to vortex tearing than in two dimensions. For slow strain, as a passage interaction is completed and the strain relaxes, the cores return to their original shape while permanent deformations remain for rapidly varying strain. For collisions, if the strain changes slowly the core shapes migrate through a known family of two-dimensional steady vortex pairs up to the limiting member of the family. Thereafter energy conservation does not allow the cores to maintain a constant shape. For rapidly varying strain, core deformation is severe and a head-tail structure in good agreement with experiments is formed. With respect to sound generation, good agreement with the measured acoustic signal for colliding rings is obtained and a feature previously thought to be due to viscous effects is shown to be an effect of inviscid core deformation alone. For passage interactions, a component of high frequency is present. Evidence for the importance of this noise source in jet noise spectra is provided. Finally, processes of fluid engulfment and rejection for an unsteady vortex ring are studied using the stable and unstable manifolds. The unstable manifold shows excellent agreement with flow visualization experiments for leapfrogging rings suggesting that it may be a good tool for numerical flow visualization in other time periodic flows