University of Illinois Urbana-Champaign

Time-symmetric integration of partial differential equations with applications to black hole physics

O'Boyle, Michael F.

Content Files
OBOYLE-DISSERTATION-2022.pdf

Permalink

https://hdl.handle.net/2142/116185

Description

Title
Time-symmetric integration of partial differential equations with applications to black hole physics
Author(s)
O'Boyle, Michael F.
Issue Date
2022-07-07
Director of Research (if dissertation) or Advisor (if thesis)
  • Gammie, Charles F
  • Markakis, Charalampos M
Doctoral Committee Chair(s)
Ricker, Paul M
Committee Member(s)
Witek, Helvi
Department of Study
Physics
Discipline
Physics
Degree Granting Institution
University of Illinois at Urbana-Champaign
Degree Name
Ph.D.
Degree Level
Dissertation
Keyword(s)
  • Geometric integrators
  • numerical analysis
  • black hole perturbation theory
  • hyperbolic partial differential equations
  • hyperboloidal compactification
  • gravitational waves
  • neutron star equation of state
Abstract
The detection of gravitational waves from merging binary black holes has heralded a new era of astronomy: gravitational radiation has joined electromagnetic radiation and elementary particles as a channel through which we may observe the universe around us. This has created the need for accurate, robust models and calculations describing the predictions that general relativity and other theories make for a whole new class of exotic, complicated phenomena. One such phenomenon is the ``extreme mass ratio inspiral,'' wherein a star or stellar-mass black hole plunges into a supermassive black hole found in a galactic nucleus and emits gravitational waves. These events present the possibility to learn about the objects found under such conditions, which would reveal information about galaxy formation and stellar evolution and probe the validity of general relativity in these regions' strong gravitational fields. Such phenomena are incredibly complicated, with only the most trivial problems possessing closed-form solutions. The aim of this work is to present a series of tools to facilitate the numerical study of these phenomena. First, we will emphasize the need for numerical methods which accurately track the physics of the problem: i.e. respect the qualitative features of the evolution equations--symmetry under time-reversal and symplecticity--and accurately track conserved quantities of interest in gravitational wave astronomy: energy and angular momentum. Next, we will discuss the problem of representing a distributional source term, i.e. a point mass, on a discrete grid, and consider the effect this has on discrete operations like numerical differentiation and integration. After providing a brief overview of the relevant aspects of black hole theory, we will discuss how to extract gravitational waves at infinite distances from the black hole using numerical methods. Finally, we will demonstrate the methods' utility in working with problems of black hole theory, showing that conservative, geometric evolution is possible. Second, we will discuss a problem related to the astrophysics of neutron stars. A current goal of gravitational wave astronomy is the constraint of the neutron star thermodynamic equation of state, which will inform the fundamental theory of nuclear interactions. It is necessary to parametrize the equation of state to compare to observation. We present an extension of the popular piecewise polytrope formalism that allows for a continuous sound speed in the nuclear matter. We found that in addition to capturing the thermodynamics of many microscopic models, it more accurately predicts the astrophysically observable properties of neutron stars -- mass, radius, and tidal deformability -- when compared to the original model.
Graduation Semester
2022-08
Type of Resource
Thesis
Handle URL
https://hdl.handle.net/2142/116185
Copyright and License Information
Copyright 2022 Michael F O'Boyle

Owning Collections

Graduate Dissertations and Theses at Illinois PRIMARY
Graduate Theses and Dissertations at Illinois