The Galactic Center black hole Sgr A* shows significant variability and flares in the submillimeter, infrared, and X-ray wavelengths. Owing to its exquisite resolution in the IR bands, the GRAVITY experiment for the first time spatially resolved the locations of three flares and showed that a bright region moves in ellipse-like trajectories close to but offset from the black hole over the course of each event. We present a model for plasmoids that form during reconnection events and orbit in the coronal region around a black hole to explain these observations. We utilize general-relativistic radiative transfer calculations that include effects from finite light travel time, plasmoid motion, particle acceleration, and synchrotron cooling and obtain a rich structure in the flare lightcurves. This model can naturally account for the observed motion of the bright regions observed by the GRAVITY experiment and the offset between the center of the centroid motion and the position of the black hole. It also explains why some flares may be double-peaked while others have only a single peak and uncovers a correlation between the structure in the lightcurve and the location of the flare. Finally, we make predictions for future observations of flares from the inner accretion flow of Sgr A* that will provide a test of this model.

We introduce a new Markov Chain Monte Carlo (MCMC) algorithm with parallel tempering for fitting theoretical models of horizon-scale images of black holes to the interferometric data from the Event Horizon Telescope (EHT). The algorithm implements forms of the noise distribution in the data that are accurate for all signal-to-noise ratios. In addition to being trivially parallelizable, the algorithm is optimized for high performance, achieving 1 million MCMC chain steps in under 20 seconds on a single processor. We use synthetic data for the 2017 EHT coverage of M87 that are generated based on analytic as well as General Relativistic Magnetohydrodynamic (GRMHD) model images to explore several potential sources of biases in fitting models to sparse interferometric data. We demonstrate that a very small number of data points that lie near salient features of the interferometric data exert disproportionate influence on the inferred model parameters. We also show that the preferred orientations of the EHT baselines introduce significant biases in the inference of the orientation of the model images. Finally, we discuss strategies that help identify the presence and severity of such biases in realistic applications.

We show that the occultation of Sagittarius A* by stars can be detected with space-based or space-ground very-long-baseline-interferometers (SVLBIs), with an expected event rate that is high due to relativistic precession. We compute the tell-tale signal of an occultation event, and describe methods to flag non-occultation events that can masquerade as the signal.

Interferometric visibilities, reflecting the complex correlations between signals recorded at antennas in an interferometric array, carry information about the angular structure of a distant source. While unknown antenna gains in both amplitude and phase can prevent direct interpretation of these measurements, certain combinations of visibilities called closure phases and closure amplitudes are independent of antenna gains and provide a convenient set of robust observables. However, these closure quantities have subtle noise properties and are generally both linearly and statistically dependent. These complications have obstructed the proper use of closure quantities in interferometric analysis, and they have obscured the relationship between analysis with closure quantities and other analysis techniques such as self calibration. We review the statistics of closure quantities, noting common pitfalls that arise when approaching low signal-to-noise due to the nonlinear propagation of statistical errors. We then develop a strategy for isolating and fitting to the independent degrees of freedom captured by the closure quantities through explicit construction of linearly independent sets of quantities along with their noise covariance in the Gaussian limit, valid for moderate signal-to-noise, and we demonstrate that model fits have biased posteriors when this covariance is ignored. Finally, we introduce a unified procedure for fitting to both closure information and partially calibrated visibilities, and we demonstrate both analytically and numerically the direct equivalence of inference based on closure quantities to that based on self calibration of complex visibilities with unconstrained antenna gains.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

Closure phases are critical in astronomical interferometry. However, their uncertainties are difficult to compute numerically. We provide a method to efficiently compute interferometric closure phase distributions in terms of an approximate distribution that is valid in the low signal-to-noise ratio regime. This is done by first showing that the true phase distribution is well approximated by the von Mises distribution, then performing a convolution of three von Mises distributions. The resulting approximation is superior than the normal distribution for all signal-to-noise ratios and, being fully analytic, allow for fast computations in statistical algorithms.

We investigate the ability of ground based gravitational wave observatories to detect gravitational wave lensing events caused by stellar mass lenses. We show that LIGO and Virgo possess the sensitivities required to detect lenses with masses as small as $\sim 30 M_\odot$ provided that the gravitational wave is observed with a signal-to-noise ratio of $\sim30$. Third generation observatories will allow detection of gravitational wave lenses with masses of $\sim 1 M_\odot$. Finally, we discuss the possibility of lensing by multiple stars, as is the case if the gravitational radiation is passing through galactic nucleus or a dense star cluster.

We investigate the effects of black hole mergers in star clusters on the black hole mass function. As black holes are not produced in pair-instability supernovae, it is suggested that there is a dearth of high mass stellar black holes. This dearth generates a gap in the upper end of the black hole mass function. Meanwhile, parameter fitting of X-ray binaries suggests the existence of a gap in the mass function under 5 solar masses. We show, through evolving a coagulation equation, that black hole mergers can appreciably fill the upper mass gap, and that the lower mass gap generates potentially observable features at larger mass scales. We also explore the importance of ejections in such systems and whether dynamical clusters can be formation sites of intermediate mass black hole seeds.

We discuss the consequences of violating the no-hair theorem on magnetic fields surrounding a black hole. This is achieved by parametrically deforming the Kerr spacetime and studying the effects of such deformations on asymptotically uniform magnetic fields around the black hole. We compute the deformed electromagnetic field for slow spins and small deformation parameter, and show that the correction is of order the deformation parameter and mimics the angular structure of a quadrupole.

Using the black hole merger rate inferred from LIGO, we calculate the abundance of tightly bound binary black holes in the Milky Way galaxy. Binaries with a small semimajor axis ($\lesssim 10 R_\odot$) originate at larger separations through conventional formation mechanisms and evolve as a result of gravitational wave emission. We find that LISA could detect them in the Milky Way. We also identify possible X-ray signatures of such binaries.

We show that an interferometer moving at a relativistic speed relative to a point source of light offers a sensitive probe of acceleration. Such an accelerometer contains no moving parts, and is thus more robust than conventional "mass-on-a-spring" accelerometers. In an interstellar mission to Alpha-Centauri, such an accelerometer could be used to measure the masses of planets around other stars as well as the mass distribution of the Milky Way Galaxy.

One avenue for testing the no-hair theorem is obtained through timing a pulsar orbiting close to a black hole and fitting for quadrupolar effects on the time-of-arrival of pulses. If deviations from the Kerr quadrupole are measured, then the no-hair theorem is invalidated. To this end, we derive an expression for the light travel time delay for a pulsar orbiting in a black-hole spacetime described by the Butterworth-Ipser metric, which has an arbitrary spin and quadrupole moment. We consider terms up to the quadrupole order in the black-hole metric and derive the time-delay expression in a closed analytic form. This allows for fast computations that are useful in fitting time-of-arrival observations of pulsars orbiting close to astrophysical black holes.

The presence of a circumnuclear stellar disk around Sgr A* and megamaser systems near other black holes indicates that dense neutral disks can be found in galactic nuclei. We show that depending on their inclination angle, optical depth, and spin temperature, these disks could be observed spectroscopically through 21 cm absorption. Related spectroscopic observations of Sgr A* can determine its HI disk parameters and the possible presence of gaps in the disk. Clumps of dense gas similar to the G2 could could also be detected in 21 cm absorption against Sgr A* radio emission.

Millimeter-wavelength VLBI observations of the supermassive black holes in Sgr A* and M87 by the Event Horizon Telescope could potentially trace the dynamics of ejected plasma blobs in real time. We demonstrate that the trajectory and tidal stretching of these blobs can be used to test general relativity and set new constraints on the mass and spin of these black holes.

We consider the time derivatives of the period $P$ of pulsars at the Galactic Center due to variations in their orbital Doppler shifts. We show that in conjunction with a measurement of a pulsar's proper motion and its projected separation from the supermassive black hole, Sgr A*,measuring two of the three derivatives $\dot{P}$, $\ddot{P}$, or $\dddot{P}$ sets a constraint that allows for the recovery of the complete six phase space coordinates of the pulsar's orbit, as well as the enclosed mass within the orbit. Thus, one can use multiple pulsars at different distances from Sgr A* to determine the radial mass distribution of stars and stellar remnants at the Galactic center. Furthermore, we consider the effect of passing stars on the pulsar's period derivatives and show how it can be exploited to measure the characteristic stellar mass in the Galactic Center.

The X-ray background during the epoch of reionization is currently poorly constrained. We demonstrate that it is possible to use first generation 21 cm experiments to calibrate it. Using the semi-numerical simulation, 21cmFAST, we calculate the dependence of the 21 cm power spectrum on the X-ray background flux. Comparing the signal to the sensitivity of the Murchison Widefield Array (MWA) we find that in the redshift interval z=8-14 the 21 cm signal is detectable based on the upper limit set by the present-day unresolved soft X-ray background. We show that there is no degeneracy between the X-ray production efficiency and the Lyman-Alpha production efficiency and that the degeneracy with the ionization fraction of the intergalactic medium can be broken.

We present broadband multi-wavelength observations of GRB 080310 at redshift z = 2.43. This burst was bright and long-lived, and unusual in having extensive optical and near IR follow-up during the prompt phase. Using these data we attempt to simultaneously model the gamma-ray, X-ray, optical and IR emission using a series of prompt pulses and an afterglow component. Initial attempts to extrapolate the high energy model directly to lower energies for each pulse reveal that a spectral break is required between the optical regime and 0.3 keV to avoid over predicting the optical flux. We demonstrate that afterglow emission alone is insufficient to describe all morphology seen in the optical and IR data. Allowing the prompt component to dominate the early-time optical and IR and permitting each pulse to have an independent low energy spectral indices we produce an alternative scenario which better describes the optical light curve. This, however, does not describe the spectral shape of GRB 080310 at early times. The fit statistics for the prompt and afterglow dominated models are nearly identical making it difficult to favour either. However one enduring result is that both models require a low energy spectral index consistent with self absorption for at least some of the pulses identified in the high energy emission model.

We present early-time optical through infrared photometry of the bright gamma-ray burst GRB 080607, starting only 6 s following the initial trigger in the rest frame. Complemented by our previously published spectroscopy, this high-quality photometric dataset allows us to solve for the extinction properties of the redshift 3.036 sightline, giving perhaps the most detailed information on the ultraviolet continuum absorption properties of any sightline outside our Local Group to date. The extinction properties are not adequately modeled by any ordinary extinction template (including the average Milky Way, Large Magellanic Cloud, and Small Magellanic Cloud curves), partially because the 2175-Angstrom feature (while present) is weaker by about a factor of two than when seen under similar circumstances locally. However, the spectral energy distribution is exquisitely fitted by the more general Fitzpatrick & Massa (1990) parameterization of Local-Group extinction, putting it in the same family as some peculiar Milky Way extinction curves. After correcting for this (considerable, A_V = 3.3 +/- 0.4 mag) extinction, GRB 080607 is revealed to have been among the most optically luminous events ever observed, comparable to the naked-eye burst GRB 080319B. Its early peak time (t_rest < 6 s) indicates a high initial Lorentz factor (Gamma > 600), while the extreme luminosity may be explained in part by a large circumburst density. Only because of its early high luminosity could the afterglow of GRB 080607 be studied in such detail in spite of the large attenuation and great distance, making this burst an excellent prototype for the understanding of other highly obscured extragalactic objects, and of the class of "dark" GRBs in particular.