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.

Millimetre-wavelength very long baseline interferometric (mm-VLBI) observations of M87 by the Event Horizon Telescope (EHT) should provide a unique opportunity to observe and characterize the origins of jet variability already seen at longer wavelengths. Synchrotron spot models have been used to model variability near black holes; this work extends these by allowing spots to shear and deform in the jet velocity field. Depending on the position of the spot, shearing forces can significantly alter the structure of the spot, producing distinct signals in reconstructed images and light curves. The maximum intensity of the shearing spot can vary by as much as a factor of five depending on the spot azimuthal launch position, but the intensity decay time depends most significantly on the spot radial launch position. Spots launched by a black hole driven jet exhibit distinct arc structures in reconstructed images, and exhibit brighter and shorter-lived enhancements of the light curve. Spots launched by a wind-driven jet have exhibit much simpler structures in the image, and longer-lived light curve enhancements than spots launched by a black hole driven jet.

We have now entered the new era of high-resolution imaging astronomy with the beginning of the Event Horizon Telescope (EHT). The EHT can resolve the dynamics of matter in the immediate vicinity around black holes at and below the horizon scale. One of the candidate black holes, Sagittarius A* flares 1\-4 times a day depending on the wavelength. A possible interpretation of these flares could be hotspots generated through magnetic reconnection events in the accretion flow. In this paper, we construct a semi-analytical model for hotspots that include the effects of shearing as a spot moves along the accretion flow. We then explore the ability of the EHT to recover these hotspots. Even including significant systematic uncertainties, such as thermal noise, diffractive scattering, and background emission due to an accretion disk, we were able to recover the hotspots and spacetime structure to sub-percent precision. Moreover, by observing multiple flaring events we show how the EHT could be used to tomographically map spacetime. This provides new avenues for testing relativistic fluid dynamics and general relativity near the event horizon of supermassive black holes.

High-resolution imaging of supermassive black holes is now possible, with new applications to testing general relativity and horizon-scale accretion and relativistic jet formation processes. Over the coming decade, the EHT will propose to add new strategically placed VLBI elements operating at 1.3mm and 0.87mm wavelength. In parallel, development of next-generation backend instrumentation, coupled with high throughput correlation architectures, will boost sensitivity, allowing the new stations to be of modest collecting area while still improving imaging fidelity and angular resolution. The goal of these efforts is to move from imaging static horizon scale structure to dynamic reconstructions that capture the processes of accretion and jet launching in near real time.

We set out to bridge the gap between regular black-hole spacetimes and observations of a black-hole shadow by the Event Horizon Telescope. We explore modifications of spinning and non-spinning black-hole spacetimes inspired by asymptotically safe quantum gravity which features a scale dependence of the Newton coupling. As a consequence, the predictions of our model, such as the shadow shape and size, depend on one free parameter determining the curvature scale at which deviations from General Relativity set in. In more general new-physics settings, it can also depart substantially from the Planck scale. In this case, the free parameter is constrained by observations, since the corresponding curvature scale is significantly below the Planck-scale. The leading new-physics effect can be recast as a scale-dependent black-hole mass, resulting in distinct observational signatures of our model. As a concrete example, we show that two mass-measurements, extracted from the size of the shadow and from Keplerian orbital motion of stars, allow to distinguish the classical from the modified, regular black-hole spacetime, yielding a bound on the free parameter. For spinning black holes, we further find that the singularity-resolving new physics puts a characteristic dent in the shadow. Finally, we argue, based on the underlying physical mechanism, that the effects we derive could be generic consequences of a large class of quantum-gravity theories.

The Galactic Center supermassive black hole Sagittarius A* (Sgr A*) is one of the most promising targets to study the dynamics of black hole accretion and outflow via direct imaging with very long baseline interferometry (VLBI). At 3.5 mm (86 GHz), the emission from Sgr A* is resolvable with the Global Millimeter VLBI Array (GMVA). We present the first observations of Sgr A* with the phased Atacama Large Millimeter/submillimeter Array (ALMA) joining the GMVA. Our observations achieve an angular resolution of ~87{\mu}as, improving upon previous experiments by a factor of two. We reconstruct a first image of the unscattered source structure of Sgr A* at 3.5 mm, mitigating effects of interstellar scattering. The unscattered source has a major axis size of 120 $\pm$ 34{\mu}as (12 $\pm$ 3.4 Schwarzschild radii), and a symmetrical morphology (axial ratio of 1.2$^{+0.3}_{-0.2}$), which is further supported by closure phases consistent with zero within 3{\sigma}. We show that multiple disk-dominated models of Sgr A* match our observational constraints, while the two jet-dominated models considered are constrained to small viewing angles. Our long-baseline detections to ALMA also provide new constraints on the scattering of Sgr A*, and we show that refractive scattering effects are likely to be weak for images of Sgr A* at 1.3 mm with the Event Horizon Telescope. Our results provide the most stringent constraints to date for the intrinsic morphology and refractive scattering of Sgr A*, demonstrating the exceptional contribution of ALMA to millimeter VLBI.

Gravitational wave and electromagnetic observations can provide new insights into the nature of matter at supra-nuclear densities inside neutron stars. Improvements in electromagnetic and gravitational wave sensing instruments continue to enhance the accuracy with which they can measure the masses, radii, and tidal deformability of neutron stars. These better measurements place tighter constraints on the equation of state of cold matter above nuclear density. In this article, we discuss a complementary approach to get insights into the structure of neutron stars by providing a model prediction for non-linear fundamental eigenmodes (f-modes) and their decay over time, which are thought to be induced by time-dependent tides in neutron star binaries. Building on pioneering studies that relate the properties of f-modes to the structure of neutron stars, we systematically study this link in the non-perturbative regime using models that utilize numerical relativity. Using a suite of fully relativistic numerical relativity simulations of oscillating TOV stars, we establish blueprints for the numerical accuracy needed to accurately compute the frequency and damping times of f-mode oscillations, which we expect to be a good guide for the requirements in the binary case. We show that the resulting f-mode frequencies match established results from linear perturbation theory, but the damping times within numerical errors depart from linear predictions. This work lays the foundation for upcoming studies aimed at a comparison of theoretical models of f-mode signatures in gravitational waves, and their uncertainties with actual gravitational wave data, searching for neutron star binaries on highly eccentric orbits, and probing neutron star structure at high densities.

With the recent advent of multi-messenger gravitational-wave astronomy and in anticipation of more sensitive, next-generation gravitational-wave detectors, we investigate the dynamics, gravitational-wave emission, and nucleosynthetic yields of numerous eccentric binary neutron-star mergers having different equations of state. For each equation of state we vary the orbital properties around the threshold of immediate merger, as well as the binary mass ratio. In addition to a study of the gravitational-wave emission including $f$-mode oscillations before and after merger, we couple the dynamical ejecta output from the simulations to the nuclear-reaction network code \texttt{SkyNet} to compute nucleosynthetic yields and compare to the corresponding results in the case of a quasi-circular merger. We find that the amount and velocity of dynamically ejected material is always much larger than in the quasi-circular case, reaching maximal values of $M_{\rm ej, max} \sim 0.1 \, M_{\odot}$ and $v_{\rm max}/c \sim 0.75$. At the same time, the properties of this material are rather insensitive to the details of the orbit, such as pericenter distance or post-encounter apoastron distance. Furthermore, while the composition of the ejected matter depends on the orbital parameters and on the equation of state, the relative nucleosynthetic yields do not, thus indicating that kilonova signatures could provide information on the orbital properties of dynamically captured neutron-star binaries.

As gravitational wave instrumentation becomes more sensitive, it is interesting to speculate about subtle effects that could be analyzed using upcoming generations of detectors. One such effect that has great potential for revealing the properties of very dense matter is fluid oscillations of neutron stars. These have been found in numerical simulations of the hypermassive remnants of double neutron star mergers and of highly eccentric neutron star orbits. Here we focus on the latter and sketch out some ideas for the production, gravitational-wave detection, and analysis of neutron star oscillations. These events will be rare (perhaps up to several tens per year could be detected using third-generation detectors such as the Einstein Telescope or the Cosmic Explorer), but they would have unique diagnostic power for the analysis of cold, catalyzed, dense matter. Furthermore, these systems are unusual in that analysis of the tidally excited f-modes of the stars could yield simultaneous measurements of their masses, moments of inertia, and tidal Love numbers, using the frequency, damping time, and amplitude of the modes. They would thus present a nearly unique opportunity to test observationally the I-Love-Q relation. The analysis of such events will require significant further work in nuclear physics and general relativistic nonlinear mode coupling, and thus we discuss further directions that will need to be pursued. For example, we note that for nearly-grazing encounters, numerical simulations show that the energy delivered to the f-modes may be up to two orders of magnitude greater than predicted in the linear theory.

Magnetic fields are believed to drive accretion and relativistic jets in black hole accretion systems, but the magnetic-field structure that controls these phenomena remains uncertain. We perform general relativistic (GR) polarized radiative transfer of time-dependent three-dimensional GR magnetohydrodynamical (MHD) simulations to model thermal synchrotron emission from the Galactic Center source Sagittarius A$^\ast$ (Sgr A*). We compare our results to new polarimetry measurements by the Event Horizon Telescope (EHT) and show how polarization in the visibility (Fourier) domain distinguishes and constrains accretion flow models with different magnetic field structures. These include models with small-scale fields in disks driven by the magnetorotational instability (MRI) as well as models with large-scale ordered fields in magnetically-arrested disks (MAD). We also consider different electron temperature and jet mass-loading prescriptions that control the brightness of the disk, funnel-wall jet, and Blandford-Znajek-driven funnel jet. Our comparisons between the simulations and observations favor models with ordered magnetic fields near the black hole event horizon in Sgr A*, although both disk- and jet-dominated emission can satisfactorily explain most of the current EHT data. We show that stronger model constraints should be possible with upcoming circular polarization and higher frequency ($349~{\rm GHz}$) measurements.

Near a black hole, differential rotation of a magnetized accretion disk is thought to produce an instability that amplifies weak magnetic fields, driving accretion and outflow. These magnetic fields would naturally give rise to the observed synchrotron emission in galaxy cores and to the formation of relativistic jets, but no observations to date have been able to resolve the expected horizon-scale magnetic-field structure. We report interferometric observations at 1.3-millimeter wavelength that spatially resolve the linearly polarized emission from the Galactic Center supermassive black hole, Sagittarius A*. We have found evidence for partially ordered fields near the event horizon, on scales of ~6 Schwarzschild radii, and we have detected and localized the intra-hour variability associated with these fields.

We report results from simulations in general relativity of magnetized disks accreting onto merging black hole binaries, starting from relaxed disk initial data. The simulations feature an effective, rapid radiative cooling scheme as a limiting case of future treatments with radiative transfer. Here we evolve the systems after binary-disk decoupling through inspiral and merger, and analyze the dependence on the binary mass ratio with $q\equiv m_{\rm bh}/M_{\rm BH}=1,1/2,$ and $1/4$. We find that the luminosity associated with local cooling is larger than the luminosity associated with matter kinetic outflows, while the electromagnetic (Poynting) luminosity associated with bulk transport of magnetic field energy is the smallest. The cooling luminosity around merger is only marginally smaller than that of a single, non-spinning black hole. Incipient jets are launched independently of the mass ratio, while the same initial disk accreting on a single non-spinning black hole does not lead to a jet, as expected. For all mass ratios we see a transient behavior in the collimated, magnetized outflows lasting $2-5 ( M/10^8M_\odot ) \rm days$ after merger: the outflows become increasingly magnetically dominated and accelerated to higher velocities, boosting the Poynting luminosity. These sudden changes can alter the electromagnetic emission across the jet and potentially help distinguish mergers of black holes in AGNs from single accreting black holes based on jet morphology alone.

We report on simulations in general relativity of magnetized disks onto black hole binaries. We vary the binary mass ratio from 1:1 to 1:10 and evolve the systems when they orbit near the binary-disk decoupling radius. We compare (surface) density profiles, accretion rates (relative to a single, non-spinning black hole), variability, effective $\alpha$-stress levels and luminosities as functions of the mass ratio. We treat the disks in two limiting regimes: rapid radiative cooling and no radiative cooling. The magnetic field lines clearly reveal jets emerging from both black hole horizons and merging into one common jet at large distances. The magnetic fields give rise to much stronger shock heating than the pure hydrodynamic flows, completely alter the disk structure, and boost accretion rates and luminosities. Accretion streams near the horizons are among the densest structures; in fact, the 1:10 no-cooling evolution results in a refilling of the cavity. The typical effective temperature in the bulk of the disk is $\sim 10^5 (M/10^8 M_\odot)^{-1/4} (L/L_{\rm edd})^{1/4} {\rm K}$ yielding characteristic thermal frequencies $\sim 10^{15} (M/10^8 M_\odot)^{-1/4} (L/L_{\rm edd})^{1/4}(1+z)^{-1}{\rm Hz} $. These systems are thus promising targets for many extragalactic optical surveys, such as LSST, WFIRST, and PanSTARRS.

Between the launch of the \textit{GGS Wind} spacecraft in 1994 November and the end of 2010, the Konus-\textit{Wind} experiment detected 296 short-duration gamma-ray bursts (including 23 bursts which can be classified as short bursts with extended emission). During this period, the IPN consisted of up to eleven spacecraft, and using triangulation, the localizations of 271 bursts were obtained. We present the most comprehensive IPN localization data on these events. The short burst detection rate, $\sim$18 per year, exceeds that of many individual experiments.

TwoPunctures is perhaps the most widely-adopted code for generating binary black hole "puncture" initial data and interpolating these (spectral) data onto evolution grids. In typical usage, the bulk of this code's run time is spent in its spectral interpolation routine. We announce a new publicly-available spectral interpolation routine that improves the performance of the original interpolation routine by a factor of ~100, yielding results consistent with the original spectral interpolation routine to roundoff precision. This note serves as a guide for installing this routine both in the original standalone TwoPunctures code and the Einstein Toolkit supported version of this code.

Neutron star binaries offer a rich phenomenology in terms of gravitational waves and merger remnants. However, most general relativistic studies have been performed for nearly circular binaries, with the exception of head-on collisions. We present the first numerical relativity investigation of mergers of eccentric equal-mass neutron-star binaries that probes the regime between head-on and circular. In addition to gravitational waves generated by the orbital motion, we find that the signal also contains a strong component due to stellar oscillations (f-modes) induced by tidal forces, extending a classical result for Newtonian binaries. The merger can lead to rather massive disks on the order of 10% of the total initial mass.

We present results from the first fully general relativistic, magnetohydrodynamic (GRMHD) simulations of an equal-mass black hole binary (BHBH) in a magnetized, circumbinary accretion disk. We simulate both the pre and post-decoupling phases of a BHBH-disk system and both "cooling" and "no-cooling" gas flows. Prior to decoupling, the competition between the binary tidal torques and the effective viscous torques due to MHD turbulence depletes the disk interior to the binary orbit. However, it also induces a two-stream accretion flow and mildly relativistic polar outflows from the BHs. Following decoupling, but before gas fills the low-density "hollow" surrounding the remnant, the accretion rate is reduced, while there is a prompt electromagnetic (EM) luminosity enhancement following merger due to shock heating and accretion onto the spinning BH remnant. This investigation, though preliminary, previews more detailed GRMHD simulations we plan to perform in anticipation of future, simultaneous detections of gravitational and EM radiation from a merging BHBH-disk system.

We perform a parameter study of non-spinning, equal and unequal mass black hole binaries on generic, eccentric orbits in numerical relativity. The linear momentum considered ranges from that of a circular orbit to ten times that value. We discuss the different manifestations of zoom-whirl behavior in the hyperbolic and the elliptic regime. The hyperbolic data set applies to dynamical capture scenarios (e.g. in globular clusters). Evolutions in the elliptic regime correspond to possible end states of supermassive black hole binaries. We spot zoom-whirl behavior for eccentricities as low as $e\sim0.5$, i.e. within the expected range of eccentricities in massive black hole binaries from galaxy mergers and binaries near galactic centers. The resulting gravitational waveforms reveal a rich structure, which will effectively break degeneracies in parameter space improving parameter estimation.

We have analyzed Cs-137 decay data, obtained from a small sample onboard the MESSENGER spacecraft en route to Mercury, with the aim of setting limits on a possible correlation between nuclear decay rates and solar activity. Such a correlation has been suggested recently on the basis of data from Mn-54 decay during the solar flare of 13 December 2006, and by indications of an annual and other periodic variations in the decay rates of Si-32, Cl-36, and Ra-226. Data from five measurements of the Cs-137 count rate over a period of approximately 5.4 years have been fit to a formula which accounts for the usual exponential decrease in count rate over time, along with the addition of a theoretical solar contribution varying with MESSENGER-Sun separation. The indication of solar influence is then characterized by a non-zero value of the calculated parameter \xi, and we find \xi=(2.8+/-8.1)x10^{-3} for Cs-137. A simulation of the increased data that can hypothetically be expected following Mercury orbit insertion on 18 March 2011 suggests that the anticipated improvement in the determination of \xi could reveal a non-zero value of \xi if present at a level consistent with other data.

Between 2000 November and 2006 May, one or more spacecraft of the interplanetary network (IPN) detected 226 cosmic gamma-ray bursts that were also detected by the FREGATE experiment aboard the HETE-II spacecraft. During this period, the IPN consisted of up to nine spacecraft, and using triangulation, the localizations of 157 bursts were obtained. We present the IPN localization data on these events.

Between 1996 July and 2002 April, one or more spacecraft of the interplanetary network detected 787 cosmic gamma-ray bursts that were also detected by the Gamma-Ray Burst Monitor and/or Wide-Field X-Ray Camera experiments aboard the BeppoSAX spacecraft. During this period, the network consisted of up to six spacecraft, and using triangulation, the localizations of 475 bursts were obtained. We present the localization data for these events.

We report near-infrared and optical follow-up observations of the afterglow of the Gamma-Ray Burst 000418 starting 2.5 days after the occurrence of the burst and extending over nearly seven weeks. GRB 000418 represents the second case for which the afterglow was initially identified by observations in the near-infrared. During the first 10 days its R-band afterglow was well characterized by a single power-law decay with a slope of 0.86. However, at later times the temporal evolution of the afterglow flattens with respect to a simple power-law decay. Attributing this to an underlying host galaxy we find its magnitude to be R=23.9 and an intrinsic afterglow decay slope of 1.22. The afterglow was very red with R-K=4 mag. The observations can be explained by an adiabatic, spherical fireball solution and a heavy reddening due to dust extinction in the host galaxy. This supports the picture that (long) bursts are associated with events in star-forming regions.

Soft gamma repeaters are high-energy transient sources associated with neutron stars in young supernova remnants. They emit sporadic, short (~ 0.1 s) bursts with soft energy spectra during periods of intense activity. The event of March 5, 1979 was the most intense and the only clearly periodic one to date. Here we report on an even more intense burst on August 27, 1998, from a different soft gamma repeater, which displayed a hard energy spectrum at its peak, and was followed by a ~300 s long tail with a soft energy spectrum and a dramatic 5.16 s period. Its peak and time integrated energy fluxes at Earth are the largest yet observed from any cosmic source. This event was probably initiated by a massive disruption of the neutron star crust, followed by an outflow of energetic particles rotating with the period of the star. Comparison of these two bursts supports the idea that magnetic energy plays an important role, and that such giant flares, while rare, are not unique, and may occur at any time in the neutron star's activity cycle.