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 study the longitudinal stability of beam-plasma systems in the presence of a density inhomogeneity in the background plasma. Previous works have focused on the non-relativistic regime where hydrodynamical models are used to evolve pre-existing Langmuir waves within inhomogeneous background plasmas. Here, for the first time we study the problem with kinetic equations in a fully-relativistic way. We do not assume the existence of Langmuir waves, and we focus on the rate and the mechanism by which waves are excited in such systems from an initial perturbation. We derive the structure of the unstable modes and compute an analytical approximation for their growth rates. Our computation is limited to dilute and cold beams, and shows an excellent agreement with particle-in-cell simulations performed using the SHARP code. We show that, due to such an inhomogeneity, the virulent beam-plasma instabilities in the intergalactic medium are not suppressed but their counterparts in the solar wind can be suppressed as evidenced by propagating type-III solar radio bursts.

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.

Pair creation on the cosmic infrared background and subsequent inverse-Compton scattering on the CMB potentially reprocesses the TeV emission of blazars into faint GeV halos with structures sensitive to intergalactic magnetic fields (IGMF). Previous work has shown that these halos are then highly-anisotropic and extended. If the coherence length of the IGMF is greater than the inverse-Compton cooling length of the pairs, then the orientation of the gamma-ray halo will be correlated with the direction of the magnetic field which is unknown and expected to change for each source. In order to constructively add each source we then use angular power spectra which are insensitive to the jet orientation. By looking at known GeV blazars detected by Fermi, we exclude the existence of an IGMF with coherence lengths $>100$ Mpc at greater than $3.9\sigma$ with current-day strengths in the range $10^{-16}$ to $10^{-15}$ G, and at 2$\sigma$ from $10^{-17}$ to $10^{-14}$ G. This provides a direct measurement of the non-existence of gamma-ray halos, providing an important check on previous results.

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.

At low redshift, the majority of supermassive black hole (SMBH) mass estimates are obtained from modeling stellar kinematics or ionized gas dynamics in the vicinity of the galaxy nucleus. For large early type galaxies, stellar kinematics models predict higher masses than gas-dynamical models. In the case of M87, this discrepancy is larger than 2 $\sigma$. Critical to gas-dynamical modeling is the assumed underlying dynamical state of the gas: that it lies on circular Keplerian orbits, potentially with some additional turbulent pressure support. This is inconsistent with models of the gas flow about low-accretion-rate SMBHs and at odds with observations of the Galactic Center. We present a simple model for non-Keplerian gas disks and explore their implications for SMBH mass measurements. We show that a larger central black hole with gas experiencing small amounts of sub-Keplerian motion can produce velocity curves similar to models that just contain circular Keplerian motions and a lower black hole mass. However, these non-Keplerian models are distinguishable from low-mass Keplerian models primarily through measurements of the velocity dispersion, wherein non-Keplerian models produce higher and narrower peak dispersions. Away from the galaxy center, but still within the circumnuclear gas disk, non-Keplerian models also become distinguishable from Keplerian models via a shift in the velocity curve. The velocity model presented in this paper is capable of resolving the discrepancy between the ionized gas dynamics and stellar kinematics mass estimates, and is applicable to gas-dynamical mass estimates of SMBHs in general.

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.

An intergalactic magnetic field stronger than $3\times10^{-13}$~G would explain the lack of a bright, extended degree-scale, GeV-energy inverse Compton component in the gamma-ray spectra of TeV-blazars. A robustly predicted consequence of the presence of such a field is the existence of degree-scale GeV-energy gamma-ray halos -- gamma-ray bow ties -- about TeV-bright active galactic nuclei, corresponding to more than half of all radio galaxies. However, the emitting regions of these halos are confined to and aligned with the direction of the relativistic jets associated with gamma-ray sources. Based on the orientation of radio jets, we align and stack corresponding degree-scale gamma-ray images of isolated Fanaroff-Riley class I and II objects and exclude the existence of these halos at overwhelming confidence, limiting the intergalactic field strength to $<10^{-15}$~G for large-scale fields and progressively larger in the diffusive regime when the correlation length of the field becomes small in comparison to 1 Mpc. When combined with prior limits on the strength of the intergalactic magnetic field, this excludes a purely magnetic explanation for the absence of halos. Thus, it requires the existence of novel physical processes that preempt the creation of halos, e.g., the presence of beam-plasma instabilities in the intergalactic medium or a drastic cutoff of the very high energy spectrum of these sources.

Millimeter polarimetry of Sgr A* probes the linearly polarized emission region on a scale of $\sim 10$ Schwarzschild radii ($R_S$) as well as the dense, magnetized accretion flow on scales out to the Bondi radius ($\sim 10^5 R_S$) through Faraday rotation. We present here multi-epoch ALMA Band 6 (230 GHz) polarimetry of Sgr A*. The results confirm a mean rotation measure, ${\rm RM} \approx -5 \times 10^5 {\rm\ rad\ m^{-2}}$, consistent with measurements over the past 20 years and support an interpretation of the RM as originating from a radiatively inefficient accretion flow (RIAF) with $\dot{M} \approx 10^{-8} { \rm\ M_{\odot}\ y^{-1} }$. Variability is observed for the first time in the RM on time scales that range from hours to months. The long-term variations may be the result of changes in the line of sight properties in a turbulent accretion flow. Short-term variations in the apparent RM are not necessarily the result of Faraday rotation and may be the result of complex emission and propagatation effects close to the black hole, some of which have been predicted in numerical modeling. We also confirm the detection of circular polarization at a mean value of $-1.1 \pm 0.2 \%$. It is variable in amplitude on time scales from hours to months but the handedness remains unchanged from that observed in past centimeter- and millimeter-wavelength detections. These results provide critical constraints for the analysis and interpretation of Event Horizon Telescope data of Sgr A*, M87, and similar sources.

Upcoming Event Horizon Telescope (EHT) observations will provide a unique opportunity to reveal the innermost region of the radiative inefficient accretion flow (RIAF) around the Galactic black hole, Sgr A*. Depending on the flow dynamics and accumulated magnetic flux, the innermost region of an RIAF could have a quasi-spherical or disk-like geometry. Here we present a phenomenological model to investigate the characteristics of the black hole shadow images with different flow geometries, together with the effect of black hole spin and flow dynamics. The resulting image consists in general of two major components: a crescent, which may surround the funnel region of the black hole or the black hole itself, and a photon ring, which may be partially luminous and overlapped with the crescent component. Compared to a quasi-spherical flow case, a disk-like flow in the vicinity of a black hole exhibits the following image features: (i) due to less material near the funnel region, the crescent structure has a smaller size, and (ii) due to the combination of emission from the flow beside and behind the black hole, the crescent structure has a more irregular shape, and a less smooth brightness distribution. How these features can result in different observables for EHT observations is discussed.

Radio images of the Galactic Center supermassive black hole, Sagittarius A* (Sgr A*), are dominated by interstellar scattering. Previous studies of Sgr A* have adopted an anisotropic Gaussian model for both the intrinsic source and the scattering, and they have extrapolated the scattering using a purely $\lambda^2$ scaling to estimate intrinsic properties. However, physically motivated source and scattering models break all three of these assumptions. They also predict that refractive scattering effects will be significant, which have been ignored in standard model fitting procedures. We analyze radio observations of Sgr A* using a physically motivated scattering model, and we develop a prescription to incorporate refractive scattering uncertainties when model fitting. We show that an anisotropic Gaussian scattering kernel is an excellent approximation for Sgr A* at wavelengths longer than 1cm, with an angular size of $(1.380 \pm 0.013) \lambda_{\rm cm}^2\,{\rm mas}$ along the major axis, $(0.703 \pm 0.013) \lambda_{\rm cm}^2\,{\rm mas}$ along the minor axis, and a position angle of $81.9^\circ \pm 0.2^\circ$. We estimate that the turbulent dissipation scale is at least $600\,{\rm km}$, with tentative support for $r_{\rm in} = 800 \pm 200\,{\rm km}$, suggesting that the ion Larmor radius defines the dissipation scale. We find that the power-law index for density fluctuations in the scattering material is $\beta < 3.47$, shallower than expected for a Kolmogorov spectrum ($\beta=11/3$), and we estimate $\beta = 3.38^{+0.08}_{-0.04}$ in the case of $r_{\rm in} = 800\,{\rm km}$. We find that the intrinsic structure of Sgr A* is nearly isotropic over wavelengths from 1.3mm to 1.3cm, with a size that is roughly proportional to wavelength. We discuss implications for models of Sgr A*, for theories of interstellar turbulence, and for imaging Sgr A* with the Event Horizon Telescope.

We report results from very long baseline interferometric (VLBI) observations of the supermassive black hole in the Galactic center, Sgr A*, at 1.3 mm (230 GHz). The observations were performed in 2013 March using six VLBI stations in Hawaii, California, Arizona, and Chile. Compared to earlier observations, the addition of the APEX telescope in Chile almost doubles the longest baseline length in the array, provides additional {\it uv} coverage in the N-S direction, and leads to a spatial resolution of $\sim$30 $\mu$as ($\sim$3 Schwarzschild radii) for Sgr A*. The source is detected even at the longest baselines with visibility amplitudes of $\sim$4-13% of the total flux density. We argue that such flux densities cannot result from interstellar refractive scattering alone, but indicate the presence of compact intrinsic source structure on scales of $\sim$3 Schwarzschild radii. The measured nonzero closure phases rule out point-symmetric emission. We discuss our results in the context of simple geometric models that capture the basic characteristics and brightness distributions of disk- and jet-dominated models and show that both can reproduce the observed data. Common to these models are the brightness asymmetry, the orientation, and characteristic sizes, which are comparable to the expected size of the black hole shadow. Future 1.3 mm VLBI observations with an expanded array and better sensitivity will allow a more detailed imaging of the horizon-scale structure and bear the potential for a deep insight into the physical processes at the black hole boundary.

We explore how inhomogeneity in the background plasma number density alters the growth of electrostatic unstable wavemodes of beam plasma systems. This is particularly interesting for blazar-driven beam-plasma instabilities, which may be suppressed by inhomogeneities in the intergalactic medium as was recently claimed in the literature. Using high resolution Particle-In-Cell simulations with the SHARP code, we show that the growth of the instability is local, i.e., regions with almost homogeneous background density will support the growth of the Langmuir waves, even when they are separated by strongly inhomogeneous regions, resulting in an overall slower growth of the instability. We also show that if the background density is continuously varying, the growth rate of the instability is lower; though in all cases, the system remains within the linear regime longer and the instability is not extinguished. In all cases, the beam loses approximately the same fraction of its initial kinetic energy in comparison to the uniform case at non-linear saturation. Thus, inhomogeneities in the intergalactic medium are unlikely to suppress the growth of blazar-driven beam-plasma instabilities.

In principle, the angular anisotropy in the extragalactic gamma-ray background (EGRB) places severe constraints upon putative populations of unresolved gamma-ray point sources. Existing estimates of the EGRB anisotropy have been constructed by excising known point sources, e.g., taken from the First or 2 Year Fermi-LAT Source Catalog (1FGL or 2FGL, respectively) and statistically analyzing the residual gamma-ray sky maps. We perform an independent check of the EGRB anisotropy limits by comparing the values obtained from the 1FGL-masked sky maps to the signal implied by sources that lie below the 1FGL detection threshold in the more sensitive 2FGL and 1FHL (First Fermi-LAT catalog of >10 GeV sources). As such, our analysis provides an internal consistency check of implications for source counts and spectral index distributions of gamma-ray bright active galactic nuclei obtained from Fermi-LAT data. Based on this, we find evidence for substantially larger anisotropies than those previously reported at energies above 5 GeV, where BL Lac objects are likely to provide the bulk of their contribution to the EGRB. This uncertainty in the EGRB anisotropy cautions against using it as an independent constraint for the high-redshift gamma-ray universe. Moreover, this would suggest that contrary to previous claims, smooth extensions of the resolved point-source population may be able to simultaneously explain both the isotropic and anisotropic components of the EGRB.

We demonstrate the unprecedented capabilities of the Event Horizon Telescope (EHT) to image the innermost dark matter profile in the vicinity of the supermassive black hole at the center of the M87 radio galaxy. We present the first model of the synchrotron emission induced by dark matter annihilations from a spiky profile in the close vicinity of a supermassive black hole, accounting for strong gravitational lensing effects. Our results show that the EHT should readily resolve dark matter spikes if present. Moreover, the photon ring surrounding the silhouette of the black hole is clearly visible in the spike emission, which introduces observable small-scale structure into the signal. We find that the dark matter-induced emission provides an adequate fit to the existing EHT data, implying that in addition to the jet, a dark matter spike may account for a sizable portion of the millimeter emission from the innermost (subparsec) region of M87. Regardless, our results show that the EHT can probe very weakly annihilating dark matter. Current EHT observations already constrain very small cross sections, typically down to a few 10^{-31} cm^{3} s^{-1} for a 10 GeV candidate, close to characteristic values for p-wave-suppressed annihilation. Future EHT observations will further improve constraints on the DM scenario.

Many astrophysical plasmas are prone to beam-plasma instabilities. For relativistic and dilute beams, the {\it spectral} support of the beam-plasma instabilities is narrow, i.e., the linearly unstable modes that grow with rates comparable to the maximum growth rate occupy a narrow range of wave numbers. This places stringent requirements on the box-sizes when simulating the evolution of the instabilities. We identify the implied lower limits on the box size imposed by the longitudinal beam plasma instability, i.e., typically the most stringent condition required to correctly capture the linear evolution of the instabilities in multidimensional simulations. We find that sizes many orders of magnitude larger than the resonant wavelength are typically required. Using one-dimensional particle-in-cell simulations, we show that the failure to sufficiently resolve the spectral support of the longitudinal instability yields slower growth and lower levels of saturation, potentially leading to erroneous physical conclusion.

Numerical heating in particle-in-cell (PIC) codes currently precludes the accurate simulation of cold, relativistic plasma over long periods, severely limiting their applications in astrophysical environments. We present a spatially higher-order accurate relativistic PIC algorithm in one spatial dimension, which conserves charge and momentum exactly. We utilize the smoothness implied by the usage of higher-order interpolation functions to achieve a spatially higher-order accurate algorithm (up to fifth order). We validate our algorithm against several test problems -- thermal stability of stationary plasma, stability of linear plasma waves, and two-stream instability in the relativistic and non-relativistic regimes. Comparing our simulations to exact solutions of the dispersion relations, we demonstrate that SHARP can quantitatively reproduce important kinetic features of the linear regime. Our simulations have a superior ability to control energy non-conservation and avoid numerical heating in comparison to common second-order schemes. We provide a natural definition for convergence of a general PIC algorithm: the complement of physical modes captured by the simulation, i.e., those that lie above the Poisson noise, must grow commensurately with the resolution. This implies that it is necessary to simultaneously increase the number of particles per cell and decrease the cell size. We demonstrate that traditional ways for testing for convergence fail, leading to plateauing of the energy error. This new PIC code enables us to faithfully study the long-term evolution of plasma problems that require absolute control of the energy and momentum conservation.

We propose a new technique for radio interferometry to obtain super-resolution full polarization images in all four Stokes parameters using sparse modeling. The proposed technique reconstructs the image in each Stokes parameter from the corresponding full-complex Stokes visibilities by utilizing two regularization functions: the $\ell _1$-norm and total variation (TV) of the brightness distribution. As an application of this technique, we present simulated linear polarization observations of two physically motivated models of M87 with the Event Horizon Telescope (EHT). We confirm that $\ell _1$+TV regularization can achieve an optimal resolution of $\sim 25-30$\% of the diffraction limit $\lambda/D_{\rm max}$, which is the nominal spatial resolution of a radio interferometer for both the total intensity (i.e. Stokes $I$) and linear polarizations (i.e. Stokes $Q$ and $U$). This optimal resolution is better than that obtained from the widely used Cotton-Schwab CLEAN algorithm or from using $\ell _1$ or TV regularizations alone. Furthermore, we find that $\ell _1$+TV regularization can achieve much better image fidelity in linear polarization than other techniques over a wide range of spatial scales, not only in the super-resolution regime, but also on scales larger than the diffraction limit. Our results clearly demonstrate that sparse reconstruction is a useful choice for high-fidelity full-polarimetric interferometric imaging.

We propose a new imaging technique for radio and optical/infrared interferometry. The proposed technique reconstructs the image from the visibility amplitude and closure phase, which are standard data products of short-millimeter very long baseline interferometers such as the Event Horizon Telescope (EHT) and optical/infrared interferometers, by utilizing two regularization functions: the $\ell_1$-norm and total variation (TV) of the brightness distribution. In the proposed method, optimal regularization parameters, which represent the sparseness and effective spatial resolution of the image, are derived from data themselves using cross validation (CV). As an application of this technique, we present simulated observations of M87 with the EHT based on four physically motivated models. We confirm that $\ell_1$+TV regularization can achieve an optimal resolution of $\sim 20-30$% of the diffraction limit $\lambda/D_{\rm max}$, which is the nominal spatial resolution of a radio interferometer. With the proposed technique, the EHT can robustly and reasonably achieve super-resolution sufficient to clearly resolve the black hole shadow. These results make it promising for the EHT to provide an unprecedented view of the event-horizon-scale structure in the vicinity of the super-massive black hole in M87 and also the Galactic center Sgr A*.

Many-degree-scale gamma-ray halos are expected to surround extragalactic high-energy gamma ray sources. These arise from the inverse Compton emission of an intergalactic population of relativistic electron/positron pairs generated by the annihilation of >100 GeV gamma rays on the extragalactic background light. These are typically anisotropic due to the jetted structure from which they originate or the presence of intergalactic magnetic fields. Here we propose a novel method for detecting these inverse-Compton gamma-ray halos based upon this anisotropic structure. Specifically, we show that by stacking suitably defined angular power spectra instead of images it is possible to robustly detect gamma-ray halos with existing Fermi Large Area Telescope (LAT) observations for a broad class of intergalactic magnetic fields. Importantly, these are largely insensitive to systematic uncertainties within the LAT instrumental response or associated with contaminating astronomical sources.

Microlensing events provide a unique capacity to study the stellar remnant population of the Galaxy. Optical microlensing suffers from a near complete degeneracy between the mass, the velocity and the distance. However, a subpopulation of lensed stars, Mira variable stars, are also radio bright, exhibiting strong SiO masers. These are sufficiently bright and compact to permit direct imaging using existing very long baseline interferometers such as the Very Long Baseline Array (VLBA). We show that these events are relatively common, occurring at a rate of $\approx 2~{\rm yr^{\ -1}}$ of which $0.1~{\rm yr^{-1}}$ are associated with Galactic black holes. Features in the associated images, e.g., the Einstein ring, are sufficiently well resolved to fully reconstruct the lens properties, enabling the measurement of mass, distance, and tangential velocity of the lensing object to a precision better than 15%. Future radio microlensing surveys conducted with upcoming radio telescopes combined with modest improvements in the VLBA could increase the rate of Galactic black hole events to roughly 10~${\rm yr}^{-1}$, sufficient to double the number of known stellar mass black holes in a couple years, and permitting the construction of distribution functions of stellar mass black hole properties.

The annihilation of TeV photons from extragalactic TeV sources and the extragalactic background light produces ultrarelativistic $e^{\pm}$ beams, which are subject to powerful plasma instabilities that sap their kinetic energy. Here we study the linear phase of the plasma instabilities that these pair beams drive. To this end, we calculate the linear growth rate of the beam plasma and oblique instability in the electrostatic approximation in both the reactive and kinetic regimes, assuming a Maxwell-J{\"u}ttner distribution for the pair beam. We reproduce the well-known reactive and kinetic growth rates for both the beam plasma and oblique mode. We demonstrate for the oblique instability that there is a broad spectrum of unstable modes that grow at the maximum rate for a wide range of beam temperatures and wave vector orientations relative to the beam. We also delineate the conditions for applicability for the reactive and kinetic regimes and find that the beam plasma mode transitions to the reactive regime at a lower Lorentz factor than the oblique mode due to a combination of their different scalings and the anisotropy of the velocity dispersions. Applying these results to the ultrarelativistic $e^{\pm}$ beams from TeV blazars, we confirm that these beams are unstable to both, the kinetic oblique mode and the reactive beam-plasma mode. These results are important in understanding how powerful plasma instabilities may sap the energy of the ultrarelativistic $e^{\pm}$ beams as they propagate through intergalactic space.

The Event Horizon Telescope is a global very-long baseline interferometer capable of probing potential deviations from the Kerr metric, which is believed to provide the unique description of astrophysical black holes. Here we report an updated constraint on the quadrupolar deviation of Sagittarius A* within the context of a radiatively inefficient accretion flow model in a quasi-Kerr background. We also simulate near-future constraints obtainable by the forthcoming eight-station array and show that in this model already a one-day observation can measure the spin magnitude to within $0.005$, the inclination to within $0.09^\circ$, the position angle to within $0.04^\circ$, and the quadrupolar deviation to within $0.005$ at $3\sigma$ confidence. Thus, we are entering an era of high-precision strong gravity measurements.

Extended inverse Compton halos are generally anticipated around extragalactic sources of gamma rays with energies above 100 GeV. These result from inverse Compton scattered cosmic microwave background photons by a population of high-energy electron/positron pairs produced by the annihilation of the high-energy gamma rays on the infrared background. Despite the observed attenuation of the high-energy gamma rays, the halo emission has yet to be directly detected. Here, we demonstrate that in most cases these halos are expected to be highly anisotropic, distributing the up-scattered gamma rays along axes defined either by the radio jets of the sources or oriented perpendicular to a global magnetic field. We present a pedagogical derivation of the angular structure in the inverse Compton halo and provide an analytic formalism that facilitates the generation of mock images. We discuss exploiting this fact for the purpose of detecting gamma-ray halos in a set of companion papers.

An initial three-station version of the Event Horizon Telescope, a millimeter-wavelength very-long baseline interferometer, has observed Sagittarius A* (Sgr A*) repeatedly from 2007 to 2013, resulting in the measurement of a variety of interferometric quantities. Of particular importance, there is now a large set of closure phases, measured over a number of independent observing epochs. We analyze these observations within the context of a realization of semi-analytic radiatively inefficient disk models, implicated by the low luminosity of Sgr A*. We find a broad consistency among the various observing epochs and between different interferometric data types, with the latter providing significant support for this class of models of Sgr A*. The new data significantly tighten existing constraints on the spin magnitude and its orientation within this model context, finding a spin magnitude of $a=0.10^{+0.30+0.56}_{-0.10-0.10}$, an inclination with respect to the line of sight of $\theta={60^\circ}^{+5^\circ+10^\circ}_{-8^\circ-13^\circ}$, and a position angle of $\xi={156^\circ}^{+10^\circ+14^\circ}_{-17^\circ-27^\circ}$ east of north. These are in good agreement with previous analyses. Notably, the previous $180^\circ$ degeneracy in the position angle has now been conclusively broken by the inclusion of the closure phase measurements. A reflection degeneracy in the inclination remains, permitting two localizations of the spin vector orientation, one of which is in agreement with the orbital angular momentum of the infrared gas cloud G2 and the clockwise disk of young stars. This possibly supports a relationship between Sgr A*'s accretion flow and these larger-scale features.

The Galactic Center black hole Sagittarius A* (Sgr A*) is a prime observing target for the Event Horizon Telescope (EHT), which can resolve the 1.3 mm emission from this source on angular scales comparable to that of the general relativistic shadow. Previous EHT observations have used visibility amplitudes to infer the morphology of the millimeter-wavelength emission. Potentially much richer source information is contained in the phases. We report on 1.3 mm phase information on Sgr A* obtained with the EHT on a total of 13 observing nights over 4 years. Closure phases, the sum of visibility phases along a closed triangle of interferometer baselines, are used because they are robust against phase corruptions introduced by instrumentation and the rapidly variable atmosphere. The median closure phase on a triangle including telescopes in California, Hawaii, and Arizona is nonzero. This result conclusively demonstrates that the millimeter emission is asymmetric on scales of a few Schwarzschild radii and can be used to break 180-degree rotational ambiguities inherent from amplitude data alone. The stability of the sign of the closure phase over most observing nights indicates persistent asymmetry in the image of Sgr A* that is not obscured by refraction due to interstellar electrons along the line of sight.

In general relativity, the angular radius of the shadow of a black hole is primarily determined by its mass-to-distance ratio and depends only weakly on its spin and inclination. If general relativity is violated, however, the shadow size may also depend strongly on parametric deviations from the Kerr metric. Based on a reconstructed image of Sagittarius A* (Sgr A*) from a simulated one-day observing run of a seven-station Event Horizon Telescope (EHT) array, we employ a Markov chain Monte Carlo algorithm to demonstrate that such an observation can measure the angular radius of the shadow of Sgr A* with an uncertainty of ~1.5 uas (6%). We show that existing mass and distance measurements can be improved significantly when combined with upcoming EHT measurements of the shadow size and that tight constraints on potential deviations from the Kerr metric can be obtained.

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 here on key science topics for the Next Generation Very Large Array in the areas of time domain, fundamental physics, and cosmology. Key science cases considered are pulsars in orbit around the Galactic Center massive black hole, Sagittarius A*, electromagnetic counterparts to gravitational waves, and astrometric cosmology. These areas all have the potential for ground-breaking and transformative discovery. Numerous other topics were discussed during the preparation of this report and some of those discussions are summarized here, as well. There is no doubt that further investigation of the science case will reveal rich and compelling opportunities.

TeV-blazars potentially heat the intergalactic medium (IGM) as their gamma rays interact with photons of the extragalactic background light to produce electron-positron pairs, which lose their kinetic energy to the surrounding medium through plasma instabilities. This results in a heating mechanism that is only weakly sensitive to the local density, and therefore approximately spatially uniform, naturally producing an inverted temperature-density relation in underdense regions. In this paper we go beyond the approximation of uniform heating and quantify the heating rate fluctuations due to the clustered distribution of blazars and how this impacts on the thermal history of the IGM. We analytically compute a filtering function that relates the heating rate fluctuations to the underlying dark matter density field. We implement it in the cosmological code GADGET-3 and perform large scale simulations to determine the impact of inhomogeneous heating. We show that, because of blazar clustering, blazar heating is inhomogeneous for z>= 2. At high redshift, the temperature-density relation shows an important scatter and presents a low temperature envelope of unheated regions, in particular at low densities and within voids. However, the median temperature of the IGM is close to that in the uniform case, albeit slightly lower at low redshift. We find that blazar heating is more complex than initially assumed and that the temperature-density relation is not unique. Our analytic model for the heating rate fluctuations couples well with large scale simulations and provides a cost-effective alternative to subgrid models.

We report on 230 GHz (1.3 mm) VLBI observations of M87 with the Event Horizon Telescope using antennas on Mauna Kea in Hawaii, Mt. Graham in Arizona and Cedar Flat in California. For the first time, we have acquired 230 GHz VLBI interferometric phase information on M87 through measurement of closure phase on the triangle of long baselines. Most of the measured closure phases are consistent with 0$^{\circ}$ as expected by physically-motivated models for 230 GHz structure such as jet models and accretion disk models. The brightness temperature of the event-horizon-scale structure is $\sim 1 \times 10^{10}$ K derived from the compact flux density of $\sim 1$ Jy and the angular size of $\sim 40 $ $\rm \mu$as $\sim$ 5.5 $R_{{\rm s}}$, which is broadly consistent with the peak brightness of the radio cores at 1-86 GHz located within $\sim 10^2$ $R_{{\rm s}}$. Our observations occurred in the middle of an enhancement in very-high-energy (VHE) $\rm \gamma$-ray flux, presumably originating in the vicinity of the central black hole. Our measurements, combined with results of multi-wavelength observations, favor a scenario in which the VHE region has an extended size of $\sim$20-60 $R_{{\rm s}}$.

It has now become clear that the radio jet in the giant elliptical galaxy M87 must turn on very close to the black hole. This implies the efficient acceleration of leptons within the jet at scales much smaller than feasible by the typical dissipative events usually invoked to explain jet synchrotron emission. Here we show that the stagnation surface, the separatrix between material that falls back into the black hole and material that is accelerated outward forming the jet, is a natural site of pair formation and particle acceleration. This occurs via an inverse-Compton pair catastrophe driven by unscreened electric fields within the charge-starved region about the stagnation surface and substantially amplified by a post-gap cascade. For typical estimates of the jet properties in M87, we find excellent quantitive agreement between the predicted relativistic lepton densities and those required by recent high-frequency radio observations of M87. This mechanism fails to adequately fill a putative jet from Sagittarius A* with relativistic leptons, which may explain the lack of an obvious radio jet in the Galactic center. Finally, this process implies a relationship between the kinetic jet power and the gamma-ray luminosity of blazars, produced during the post-gap cascade.

Chandra X-ray observations of the nearby brightest cluster galaxy M87 resolve the hot gas structure across the Bondi accretion radius of the central supermassive black hole, a measurement possible in only a handful of systems but complicated by the bright nucleus and jet emission. By stacking only short frame-time observations to limit pileup, and after subtracting the nuclear PSF, we analysed the X-ray gas properties within the Bondi radius at 0.12-0.22 kpc (1.5-2.8 arcsec), depending on the black hole mass. Within 2 kpc radius, we detect two significant temperature components, which are consistent with constant values of 2 keV and 0.9 keV down to 0.15 kpc radius. No evidence was found for the expected temperature increase within ~0.25 kpc due to the influence of the SMBH. Within the Bondi radius, the density profile is consistent with $\rho\propto r^{-1}$. The lack of a temperature increase inside the Bondi radius suggests that the hot gas structure is not dictated by the SMBH's potential and, together with the shallow density profile, shows that the classical Bondi rate may not reflect the accretion rate onto the SMBH. If this density profile extends in towards the SMBH, the mass accretion rate onto the SMBH could be at least two orders of magnitude less than the Bondi rate, which agrees with Faraday rotation measurements for M87. We discuss the evidence for outflow from the hot gas and the cold gas disk and for cold feedback, where gas cooling rapidly from the hot atmosphere could feed the cirumnuclear disk and fuel the SMBH. At 0.2 kpc radius, the cooler X-ray temperature component represents ~20% of the total X-ray gas mass and, by losing angular momentum to the hot gas component, could provide a fuel source of cold clouds within the Bondi radius.

The 6 billion solar mass supermassive black hole at the center of the giant elliptical galaxy M87 powers a relativistic jet. Observations at millimeter wavelengths with the Event Horizon Telescope have localized the emission from the base of this jet to angular scales comparable to the putative black hole horizon. The jet might be powered directly by an accretion disk or by electromagnetic extraction of the rotational energy of the black hole. However, even the latter mechanism requires a confining thick accretion disk to maintain the required magnetic flux near the black hole. Therefore, regardless of the jet mechanism, the observed jet power in M87 implies a certain minimum mass accretion rate. If the central compact object in M87 were not a black hole but had a surface, this accretion would result in considerable thermal near-infrared and optical emission from the surface. Current flux limits on the nucleus of M87 strongly constrain any such surface emission. This rules out the presence of a surface and thereby provides indirect evidence for an event horizon.

It has been clear for some time now that super-critical surface magnetic fields, exceeding 4 x 10^13 G, exist on a subset of neutron stars. These magnetars may harbor interior fields many orders of magnitude larger, potentially reaching equipartition values. However, the impact of these strong fields on stellar structure has been largely ignored, potentially complicating attempts to infer the high density nuclear equation of state. Here we assess the effect of these strong magnetic fields on the mass-radius relationship of neutron stars. We employ an effective field theory model for the nuclear equation of state that includes the impact of hyperons, anomalous magnetic moments, and the physics of the crust. We consider two magnetic field geometries, bounding the likely magnitude of the impact of magnetic fields: a statistically isotropic, tangled field and a force-free configuration. In both cases even equipartition fields have at most a 30% impact on the maximum mass. However, the direction of the effect of the magnetic field depends on the geometry employed - force-free fields leading to reductions in the maximum neutron star mass and radius while tangled fields increase both - challenging the common intuition in the literature on the impact of magnetic fields.

Very-high energy gamma-rays from extragalactic sources pair-produce off of the extragalactic background light, yielding an electron-positron pair beam. This pair beam is unstable to various plasma instabilities, especially the "oblique" instability, which can be the dominant cooling mechanism for the beam. However, recently, it has been claimed that nonlinear Landau damping renders it physically irrelevant by reducing the effective damping rate to a low level. Here, we show with numerical calculations that the effective damping rate is $8\times 10^{-4}$ of the growth rate of the linear instability, which is sufficient for the "oblique" instability to be the dominant cooling mechanism of these pair beams. In particular, we show that previous estimates of this rate ignored the exponential cutoff in the scattering amplitude at large wavenumber and assumed that the damping of scattered waves entirely depends on collisions, ignoring collisionless processes. We find that the total wave energy eventually grows to approximate equipartition with the beam by increasingly depositing energy into long wavelength modes. As we have not included the effect of nonlinear wave-wave interactions on these long wavelength modes, this scenario represents the "worst-case" scenario for the oblique instability. As it continues to drain energy from the beam at a faster rate than other processes, we conclude that the "oblique" instability is sufficiently strong to make it the physically dominant cooling mechanism for high-energy pair beams in the intergalactic medium.

Observations of the black hole in the center of the Milky Way with the Event Horizon Telescope at 1.3 mm have revealed a size of the emitting region that is smaller than the size of the black-hole shadow. This can be reconciled with the spectral properties of the source, if the accretion flow is seen at a relatively high inclination (50-60 degrees). Such an inclination makes the angular momentum of the flow, and perhaps of the black hole, nearly aligned with the angular momenta of the orbits of stars that lie within 3 arcsec from the black hole. We discuss the implications of such an alignment for the properties of the black hole and of its accretion flow. We argue that future Event-Horizon-Telescope observations will not only refine the inclination of Sgr A* but also measure precisely its orientation on the plane of the sky.

The image of the emission surrounding the black hole in the center of the Milky Way is predicted to exhibit the imprint of general relativistic (GR) effects, including the existence of a shadow feature and a photon ring of diameter ~50 microarcseconds. Structure on these scales can be resolved by millimeter-wavelength very long baseline interferometry (VLBI). However, strong-field GR features of interest will be blurred at lambda >= 1.3 mm due to scattering by interstellar electrons. The scattering properties are well understood over most of the relevant range of baseline lengths, suggesting that the scattering may be (mostly) invertible. We simulate observations of a model image of Sgr A* and demonstrate that the effects of scattering can indeed be mitigated by correcting the visibilities before reconstructing the image. This technique is also applicable to Sgr A* at longer wavelengths.

We demonstrate that polarimetric interferometry can be used to extract precise spatial information about compact polarized flares of Sgr A*. We show that, for a faint dynamical component, a single interferometric baseline suffices to determine both its polarization and projected displacement from the quiescent intensity centroid. A second baseline enables two-dimensional reconstruction of the displacement, and additional baselines can self-calibrate using the flare, enhancing synthesis imaging of the quiescent emission. We apply this technique to simulated 1.3-mm wavelength observations of a "hot spot" embedded in a radiatively inefficient accretion disk around Sgr A*. Our results indicate that, even with current sensitivities, polarimetric interferometry with the Event Horizon Telescope can achieve ~5 microarcsecond relative astrometry of compact flaring structures near Sgr A* on timescales of minutes.

The Event Horizon Telescope (EHT) is a project to assemble a Very Long Baseline Interferometry (VLBI) network of mm wavelength dishes that can resolve strong field General Relativistic signatures near a supermassive black hole. As planned, the EHT will include enough dishes to enable imaging of the predicted black hole "shadow", a feature caused by severe light bending at the black hole boundary. The center of M87, a giant elliptical galaxy, presents one of the most interesting EHT targets as it exhibits a relativistic jet, offering the additional possibility of studying jet genesis on Schwarzschild radius scales. Fully relativistic models of the M87 jet that fit all existing observational constraints now allow horizon-scale images to be generated. We perform realistic VLBI simulations of M87 model images to examine detectability of the black shadow with the EHT, focusing on a sequence of model images with a changing jet mass load radius. When the jet is launched close to the black hole, the shadow is clearly visible both at 230 and 345 GHz. The EHT array with a resolution of 20-30$\mu$as resolution ($\sim$2-4 Schwarzschild radii) is able to image this feature independent of any theoretical models and we show that imaging methods used to process data from optical interferometers are applicable and effective for EHT data sets. We demonstrate that the EHT is also capable of tracing real-time structural changes on a few Schwarzschild radii scales, such as those implicated by VHE flaring activity of M87. While inclusion of ALMA in the EHT is critical for shadow imaging, generally the array is robust against loss of a station.

The advent of the Event Horizon Telescope (EHT), a millimeter-wave very-long baseline interferometric array, has enabled spatially-resolved studies of the sub-horizon-scale structure for a handful of supermassive black holes. Among these, the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), presents the largest angular cross section. Thus far, these studies have focused upon measurements of the black hole spin and the validation of low-luminosity accretion models. However, a critical input into the analysis of EHT data is the structure of the black hole spacetime, and thus these observations provide the novel opportunity to test the applicability of the Kerr metric to astrophysical black holes. Here we present the first simulated images of a radiatively inefficient accretion flow (RIAF) around Sgr A* employing a quasi-Kerr metric that contains an independent quadrupole moment in addition to the mass and spin that fully characterize a black hole in general relativity. We show that these images can be significantly different from the images of a RIAF around a Kerr black hole with the same spin and demonstrate the feasibility of testing the no-hair theorem by constraining the quadrupolar deviation from the Kerr metric with existing EHT data. Equally important, we find that the disk inclination and spin orientation angles are robust to the inclusion of additional parameters, providing confidence in previous estimations assuming the Kerr metric based upon EHT observations. However, at present the limits upon potential modifications of the Kerr metric remain weak.

An international consortium is presently constructing a beamformer for the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile that will be available as a facility instrument. The beamformer will aggregate the entire collecting area of the array into a single, very large aperture. The extraordinary sensitivity of phased ALMA, combined with the extremely fine angular resolution available on baselines to the Northern Hemisphere, will enable transformational new very long baseline interferometry (VLBI) observations in Bands 6 and 7 (1.3 and 0.8 mm) and provide substantial improvements to existing VLBI arrays in Bands 1 and 3 (7 and 3 mm). The ALMA beamformer will have impact on a variety of scientific topics, including accretion and outflow processes around black holes in active galactic nuclei (AGN), tests of general relativity near black holes, jet launch and collimation from AGN and microquasars, pulsar and magnetar emission processes, the chemical history of the universe and the evolution of fundamental constants across cosmic time, maser science, and astrometry.

The extragalactic gamma-ray sky at TeV energies is dominated by blazars, a subclass of accreting super-massive black holes with powerful relativistic outflows directed at us. Only constituting a small fraction of the total power output of black holes, blazars were thought to have a minor impact on the universe at best. As we argue here, the opposite is true and the gamma-ray emission from TeV blazars can be thermalized via beam-plasma instabilities on cosmological scales with order unity efficiency, resulting in a potentially dramatic heating of the low-density intergalactic medium. Here, we review this novel heating mechanism and explore the consequences for the formation of structure in the universe. In particular, we show how it produces an inverted temperature-density relation of the intergalactic medium that is in agreement with observations of the Lyman-alpha forest. This suggests that blazar heating can potentially explain the paucity of dwarf galaxies in galactic halos and voids, and the bimodality of galaxy clusters. This also transforms our understanding of the evolution of blazars, their contribution to the extra-galactic gamma-ray background, and how their individual spectra can be used in constraining intergalactic magnetic fields.

Fermi has been instrumental in constraining the luminosity function and redshift evolution of gamma-ray bright blazars. This includes limits upon the spectrum and anisotropy of the extragalactic gamma-ray background (EGRB), redshift distribution of nearby Fermi active galactic nuclei (AGN), and the construction of a log(N)-log(S) relation. Based upon these, it has been argued that the evolution of the gamma-ray bright blazar population must be much less dramatic than that of other AGN. However, critical to such claims is the assumption that inverse Compton cascades reprocess emission above a TeV into the Fermi energy range, substantially enhancing the strength of the observed limits. Here we demonstrate that in the absence of such a process, due, e.g., to the presence of virulent plasma beam instabilities that preempt the cascade, a population of TeV-bright blazars that evolve similarly to quasars is consistent with the population of hard gamma-ray blazars observed by Fermi. Specifically, we show that a simple model for the properties and luminosity function is simultaneously able to reproduce their log(N)-log(S) relation, local redshift distribution, and contribution to the EGRB and its anisotropy without any free parameters. Insofar the naturalness of a picture in which the hard gamma-ray blazar population exhibits the strong redshift evolution observed in other tracers of the cosmological history of accretion onto halos is desirable, this lends support for the absence of the inverse Compton cascades and the existence of the beam plasma instabilities.

Approximately 10% of active galactic nuclei exhibit relativistic jets, which are powered by accretion of matter onto super massive black holes. While the measured width profiles of such jets on large scales agree with theories of magnetic collimation, predicted structure on accretion disk scales at the jet launch point has not been detected. We report radio interferometry observations at 1.3mm wavelength of the elliptical galaxy M87 that spatially resolve the base of the jet in this source. The derived size of 5.5 +/- 0.4 Schwarzschild radii is significantly smaller than the innermost edge of a retrograde accretion disk, suggesting that the M87 jet is powered by an accretion disk in a prograde orbit around a spinning black hole.

The Universe is opaque to extragalactic very high-energy gamma rays (VHEGRs, E>100 GeV) because they annihilate and pair produce on the extragalactic background light. The resulting ultra-relativistic pairs are assumed to lose energy through inverse Compton scattering of CMB photons. In Broderick et al. (2011, Paper I of this three paper series), we argued that instead powerful plasma instabilities in the ultra-relativistic pair beam dissipate the kinetic energy of the TeV-generated pairs locally, heating the intergalactic medium (IGM). Here, we explore the effect of this heating upon the thermal history of the IGM. We collate the observed extragalactic VHEGR sources to determine a local VHEGR heating rate and correct for the pointed nature of VHEGR observations using Fermi observations of high and intermediate peaked BL Lacs. Because the local extragalactic VHEGR flux is dominated by TeV blazars, we tie the TeV blazar luminosity density to the quasar luminosity density, and produce a VHEGR heating rate as a function of redshift. This heating is relatively homogeneous for z<~4 with increasing spatial variation at higher redshift (order unity at z~6). This new heating process dominates photoheating at low redshift and the inclusion of TeV blazar heating qualitatively and quantitatively changes the structure and history of the IGM. TeV blazars produce a uniform volumetric heating rate that is sufficient to increase the temperature of the mean density IGM by nearly an order of magnitude, and at low densities by substantially more, naturally producing an inverted equation of state inferred by observations of the Ly-alpha forest, a feature that is difficult to reconcile with standard reionization models. Finally, we close with a discussion on the possibility of detecting this hot low-density IGM, but find that such measurements are currently not feasible. (abridged)

A subset of blazars emit TeV gamma rays which annihilate and pair produce on the extragalactic background light. We have argued in Broderick et al. (2011, Paper I) that plasma beam instabilities can dissipate the pairs' energy locally. This heats the intergalactic medium and dramatically increases its entropy after redshift z~2, with important implications for structure formation: (1) This suggests a scenario for the origin of the cool core (CC)/non-cool core (NCC) bimodality in galaxy clusters and groups. Early forming galaxy groups are unaffected because they can efficiently radiate the additional entropy, developing a CC. However, late forming groups do not have sufficient time to cool before the entropy is gravitationally reprocessed through successive mergers - counteracting cooling and raising the core entropy further. Hence blazar heating works different than feedback by active galactic nuclei, which balances radiative cooling but is unable to transform CC into NCC clusters due to the weak coupling to the cluster gas. (2) We predict a suppression of the Sunyaev-Zel'dovich power spectrum on angular scales smaller than 5' due to the globally reduced central pressure of groups and clusters forming after z~1. (3) Our redshift dependent entropy floor increases the characteristic halo mass below which dwarf galaxies cannot form by a factor of ~10 (50) at mean density (in voids) over that found in models that include photoionization alone. This prevents the formation of late forming dwarf galaxies (z<2) with masses ranging from 10^{10} to 10^{11} M_sun for redshifts z~2 to 0, respectively. This may help resolve the "missing satellite problem" in the Milky Way and the "void phenomenon" of the low observed abundances of dwarf satellites compared to cold dark matter simulations and may bring the observed early star formation histories into agreement with galaxy formation models. (abridged)

Inverse-Compton cascades initiated by energetic gamma rays (E>100 GeV) enhance the GeV emission from bright, extragalactic TeV sources. The absence of this emission from bright TeV blazars has been used to constrain the intergalactic magnetic field (IGMF), and the stringent limits placed upon the unresolved extragalactic gamma-ray background (EGRB) by Fermi has been used to argue against a large number of such objects at high redshifts. However, these are predicated upon the assumption that inverse-Compton scattering is the primary energy-loss mechanism for the ultra-relativistic pairs produced by the annihilation of the energetic gamma rays on extragalactic background light photons. Here we show that for sufficiently bright TeV sources (isotropic-equivalent luminosities >10^{42} erg/s) plasma beam instabilities, specifically the "oblique" instability, present a plausible mechanism by which the energy of these pairs can be dissipated locally, heating the intergalactic medium. Since these instabilities typically grow on timescales short in comparison to the inverse-Compton cooling rate, they necessarily suppress the inverse-Compton cascades. As a consequence, this places a severe constraint upon efforts to limit the IGMF from the lack of a discernible GeV bump in TeV sources. Similarly, it considerably weakens the Fermi limits upon the evolution of blazar populations. Specifically, we construct a TeV-blazar luminosity function from those objects presently observed and find that it is very well described by the quasar luminosity function at z~0.1, shifted to lower luminosities and number densities, suggesting that both classes of sources are regulated by similar processes. Extending this relationship to higher redshifts, we show that the magnitude and shape of the EGRB above ~10 GeV is naturally reproduced with this particular example of a rapidly evolving TeV-blazar luminosity function.

It has been realised only recently that TeV emission from blazars can significantly heat the intergalactic medium (IGM) by pair-producing high-energy electrons and positrons, which in turn excite vigorous plasma instabilities, leading to a local dissipation of the pairs' kinetic energy. In this work, we use cosmological hydrodynamical simulations to model the impact of this blazar heating on the Lyman-alpha forest at redshifts z~2-3. We find that blazar heating produces an inverted temperature-density relation in the IGM and naturally resolves many of the problems present in previous simulations of the forest that included photoheating alone. In particular, our simulations with blazar heating simultaneously reproduce the observed effective optical depth and temperature as a function of redshift, the observed probability distribution functions of the transmitted flux, and the observed flux power spectra, over the full redshift range 2<z<3 analysed here. Additionally, by deblending the Lyman-alpha forest into a sum of thermally broadened individual lines, we find superb agreement with the observed lower cutoff of the line-width distribution and abundances of neutral hydrogen column densities. Using the most recent constraints on the cosmic ultraviolet (UV) background, this excellent agreement with observations does not require rescaling the amplitude of the UV background; a procedure that was routinely used in the past to match the observed level of transmitted flux. We also show that our blazar-heated model matches the data better than standard simulations even when such a rescaling is allowed. This concordance between Lyman-alpha data and simulations, which are based on the most recent cosmological parameters, suggests that the inclusion of blazar heating alleviates previous tensions on constraints for sigma_8 derived from Lyman-alpha measurements and other cosmological data. [abridged]

Millimeter wave Very Long Baseline Interferometry (mm-VLBI) provides access to the emission region surrounding Sagittarius A*, the supermassive black hole at the center of the Milky Way, on sub-horizon scales. Recently, a closure phase of 0+-40 degrees was reported on a triangle of Earth-sized baselines (SMT-CARMA-JCMT) representing a new constraint upon the structure and orientation of the emission region, independent from those provided by the previously measured 1.3mm-VLBI visibility amplitudes alone. Here, we compare this to the closure phases associated with a class of physically motivated, radiatively inefficient accretion flow models, and present predictions for future mm-VLBI experiments with the developing Event Horizon Telescope (EHT). We find that the accretion flow models are capable of producing a wide variety of closure phases on the SMT-CARMA-JCMT triangle, and thus not all models are consistent with the recent observations. However, those models that reproduce the 1.3mm-VLBI visibility amplitudes overwhelmingly have SMT-CARMA-JCMT closure phases between +-30 degrees, and are therefore broadly consistent with all current mm-VLBI observations. Improving station sensitivity by factors of a few, achievable by increases in bandwidth and phasing together multiple antennas at individual sites, should result in physically relevant additional constraints upon the model parameters and eliminate the current 180 degree ambiguity on the source orientation. When additional stations are included, closure phases of order 45--90 degrees are typical. In all cases the EHT will be able to measure these with sufficient precision to produce dramatic improvements in the constraints upon the spin of Sgr A*.

With the advent of the Event Horizon Telescope (EHT), a millimeter/sub-millimeter very-long baseline interferometer (VLBI), it has become possible to image a handful of black holes with sub-horizon resolutions. However, these images do not translate into microarcsecond absolute positions due to the lack of absolute phase information when an external phase reference is not used. Due to the short atmospheric coherence time at these wavelengths, nodding between the source and phase reference is impractical. However, here we suggest an alternative scheme which makes use of the fact that many of the VLBI stations within the EHT are arrays in their own right. With this we show that it should be possible to absolutely position the supermassive black holes at the centers of the Milky Way (Sgr A*) and M87 relative to nearby objects with precisions of roughly 1 microarcsecond. This is sufficient to detect the perturbations to Sgr A*'s position resulting from interactions with the stars and stellar-mass black holes in the Galactic cusp on year timescales, and severely constrain the astrophysically relevant parameter space for an orbiting intermediate mass black hole, implicated in some mechanisms for producing the young massive stars in the Galactic center. For M87, it allows the registering of millimeter images, in which the black hole may be identified by its silhouette against nearby emission, and existing larger scale radio images, eliminating present ambiguities in the nature of the radio core and inclination, opening angle, and source of the radio jet.

Millimeter very-long baseline interferometry (mm-VLBI) provides the novel capacity to probe the emission region of a handful of supermassive black holes on sub-horizon scales. For Sagittarius A* (Sgr A*), the supermassive black hole at the center of the Milky Way, this provides access to the region in the immediate vicinity of the horizon. Broderick et al. (2009) have already shown that by leveraging spectral and polarization information as well as accretion theory, it is possible to extract accretion-model parameters (including black hole spin) from mm-VLBI experiments containing only a handful of telescopes. Here we repeat this analysis with the most recent mm-VLBI data, considering a class of aligned, radiatively inefficient accretion flow (RIAF) models. We find that the combined data set rules out symmetric models for Sgr A*'s flux distribution at the 3.9-sigma level, strongly favoring length-to-width ratios of roughly 2.4:1. More importantly, we find that physically motivated accretion flow models provide a significantly better fit to the mm-VLBI observations than phenomenological models, at the 2.9-sigma level. This implies that not only is mm-VLBI presently capable of distinguishing between potential physical models for Sgr A*'s emission, but further that it is sensitive to the strong gravitational lensing associated with the propagation of photons near the black hole. Based upon this analysis we find that the most probable magnitude, viewing angle, and position angle for the black hole spin are a=0.0(+0.64+0.86), theta=68(+5+9)(-20-28) degrees, and xi=-52(+17+33)(-15-24) east of north, where the errors quoted are the 1-sigma and 2-sigma uncertainties.

For the first time it has become possible to compare global 3D general relativistic magnetohydrodynamic (GRMHD) jet formation simulations directly to very-long baseline interferometric multi-frequency polarization observations of the pc-scale structure of active galactic nucleus (AGN) jets. Unlike the jet emission, which requires post hoc modeling of the non-thermal electrons, the Faraday rotation measures (RMs) depend primarily upon simulated quantities and thus provide a robust way in which to confront simulations with observations. We compute RM distributions of 3D global GRMHD jet formation simulations, with which we explore the dependence upon model and observational parameters, emphasizing the signatures of structures generic to the theory of MHD jets. With typical parameters, we find that it is possible to reproduce the observed magnitudes and many of the structures found in AGN jet RMs, including the presence of transverse RM gradients. In our simulations the RMs are generated within a smooth extension of the jet itself, containing ordered toroidally dominated magnetic fields. This results in a particular bilateral morphology that is unlikely to arise due to Faraday rotation in distant foreground clouds. However, critical to efforts to probe the Faraday screen will be resolving the transverse jet structure. Therefore, the RMs of radio cores may not be reliable indicators of the properties of the rotating medium. Finally, we are able to constrain the particle content of the jet, finding that at pc-scales AGN jets are electromagnetically dominated, with roughly 2% of the comoving energy in nonthermal leptons and much less in baryons.

For a number of starless cores, self-absorbed molecular line and column density observations have implied the presence of large-amplitude oscillations. We examine the consequences of these oscillations on the evolution of the cores and the interpretation of their observations. We find that the pulsation energy helps support the cores and that the dissipation of this energy can lead toward instability and star formation. In this picture, the core lifetimes are limited by the pulsation decay timescales, dominated by non-linear mode-mode coupling, and on the order of ~few x 10^5--10^6 yr. Notably, this is similar to what is required to explain the relatively low rate of conversion of cores into stars. For cores with large-amplitude oscillations, dust continuum observations may appear asymmetric or irregular. As a consequence, some of the cores that would be classified as supercritical may be dynamically stable when oscillations are taken into account. Thus, our investigation motivates a simple hydrodynamic picture, capable of reproducing many of the features of the progenitors of stars without the inclusion of additional physical processes, such as large-scale magnetic fields.

Faraday rotation measurements have provided an invaluable technique with which to measure the properties of astrophysical magnetized plasmas. Unfortunately, typical observations provide information only about the density-weighted average of the magnetic field component parallel to the line of sight. As a result, the magnetic field geometry along the line of sight, and in many cases even the location of the rotating material, is poorly constrained. Frequently, interpretations of Faraday rotation observations are dependent upon underlying models of the magnetic field being probed (e.g., uniform, turbulent, equipartition). However, we show that at sufficiently low frequencies, specifically below roughly 13 (RM/rad m^-2)^(1/4) (B/G)^(1/2) MHz, the character of Faraday rotation changes, entering what we term the ``super-adiabatic regime'' in which the rotation measure is proportional to the integrated absolute value of the line-of-sight component of the field. As a consequence, comparing rotation measures at high frequencies with those in this new regime provides direct information about the geometry of the magnetic field along the line of sight. Furthermore, the frequency defining the transition to this new regime, nu_SA, depends directly upon the local electron density and magnetic field strength where the magnetic field is perpendicular to the line of sight, allowing the unambiguous distinction between Faraday rotation within and in front of the emission region. Typical values of nu_SA range from 10 kHz to 10 GHz, depending upon the details of the Faraday rotating environment. In particular, for resolved AGN, including the black holes at the center of the Milky Way (Sgr A*) and M81, nu_SA ranges from roughly 10 MHz to 10 GHz, and thus can be probed via existing and up-coming ground-based radio observatories.

Sagittarius A* is the source of near infrared, X-ray, radio, and (sub)millimeter emission associated with the supermassive black hole at the Galactic Center. In the submillimeter regime, Sgr A* exhibits time-variable linear polarization on timescales corresponding to <10 Schwarzschild radii of the presumed 4 million solar mass black hole. In previous work, we demonstrated the potential for total-intensity (sub)millimeter-wavelength VLBI to detect time-variable -- and periodic -- source structure changes in the Sgr A* black hole system using nonimaging analyses. Here we extend this work to include full polarimetric VLBI observations. We simulate full-polarization (sub)millimeter VLBI data of Sgr A* using a hot-spot model that is embedded within an accretion disk, with emphasis on nonimaging polarimetric data products that are robust against calibration errors. Although the source-integrated linear polarization fraction in the models is typically only a few percent, the linear polarization fraction on small angular scales can be much higher, enabling the detection of changes in the polarimetric structure of Sgr A* on a wide variety of baselines. The shortest baselines track the source-integrated linear polarization fraction, while longer baselines are sensitive to polarization substructures that are beam-diluted by connected-element interferometry. The detection of periodic variability in source polarization should not be significantly affected even if instrumental polarization terms cannot be calibrated out. As more antennas are included in the (sub)mm-VLBI array, observations with full polarization will provide important new diagnostics to help disentangle intrinsic source polarization from Faraday rotation effects in the accretion and outflow region close to the black hole event horizon.

Polarization has proved an invaluable tool for probing magnetic fields in relativistic jets. Maps of the intrinsic polarization vectors have provided the best evidence to date for uniform, toroidally dominated magnetic fields within jets. More recently, maps of the rotation measure (RM) in jets have for the first time probed the field geometry of the cool, moderately relativistic surrounding material. In most cases, clear signatures of toroidal magnetic field are detected, corresponding to gradients in RM profiles transverse to the jet. However, in many objects these profiles also display marked asymmetries which are difficult to explain in simple helical jet models. Furthermore, in some cases the RM profiles are strongly frequency and/or time dependent. Here we show that these features may be naturally accounted for by including relativistic helical motion in the jet model. In particular, we are able to reproduce bent RM profiles observed in a variety of jets, frequency dependent RM profile morphologies and even the time dependence of the RM profiles of knots in 3C 273. Finally, we predict that some sources may show reversals in their RM profiles at sufficiently high frequencies, depending upon the the ratio of the components of jet sheath velocity transverse and parallel to the jet. Thus, multi-frequency RM maps promise a novel way in which to probe the velocity structure of relativistic outflows.

A long standing goal in astrophysics is to directly observe the immediate environment of a black hole with angular resolution comparable to the event horizon. Realizing this goal would open a new window on the study of General Relativity in the strong field regime, accretion and outflow processes at the edge of a black hole, the existence of an event horizon, and fundamental black hole physics (e.g., spin). Steady long-term progress on improving the capability of Very Long Baseline Interferometry (VLBI) at short wavelengths has now made it extremely likely that this goal will be achieved within the next decade. The most compelling evidence for this is the recent observation by 1.3mm VLBI of Schwarzschild radius scale structure in SgrA*, the compact source of radio, submm, NIR and xrays at the center of the Milky Way. SgrA* is thought to mark the position of a ~4 million solar mass black hole, and because of its proximity and estimated mass presents the largest apparent event horizon size of any black hole candidate in the Universe. Over the next decade, existing and planned mm/submm facilities will be combined into a high sensitivity, high angular resolution "Event Horizon Telescope" that will bring us as close to the edge of black hole as we will come for decades. This white paper describes the science case for mm/submm VLBI observations of both SgrA* and M87 (a radio loud AGN of a much more luminous class that SgrA*). We emphasize that while there is development and procurement involved, the technical path forward is clear, and the recent successful observations have removed much of the risk that would normally be associated with such an ambitious project.

Black hole event horizons, causally separating the external universe from compact regions of spacetime, are one of the most exotic predictions of General Relativity (GR). Until recently, their compact size has prevented efforts to study them directly. Here we show that recent millimeter and infrared observations of Sagittarius A* (Sgr A*), the supermassive black hole at the center of the Milky Way, all but requires the existence of a horizon. Specifically, we show that these observations limit the luminosity of any putative visible compact emitting region to below 0.4% of Sgr A*'s accretion luminosity. Equivalently, this requires the efficiency of converting the gravitational binding energy liberated during accretion into radiation and kinetic outflows to be greater than 99.6%, considerably larger than those implicated in Sgr A*, and therefore inconsistent with the existence of such a visible region. Finally, since we are able to frame this argument entirely in terms of observable quantities, our results apply to all geometric theories of gravity that admit stationary solutions, including the commonly discussed f(R) class of theories.

The silhouette cast by the horizon of the supermassive black hole in M87 can now be resolved with the emerging millimeter very-long baseline interferometry (VLBI) capability. Despite being ~2000 times farther away than SgrA* (the supermassive black hole at the center of the Milky-Way and the primary target for horizon-scale imaging), M87's much larger black hole mass results in a horizon angular scale roughly half that of SgrA*'s, providing another practical target for direct imaging. However, unlike SgrA*, M87 exhibits a powerful radio jet, providing an opportunity to study jet formation physics on horizon scales. We employ a simple, qualitatively correct force-free jet model to explore the expected high-resolution images of M87 at wavelengths of 1.3mm and 0.87mm (230GHz and 345GHz), for a variety of jet parameters. We show that future VLBI data will be able to constrain the size of the jet footprint, the jet collimation rate, and the black hole spin. Polarization will further probe the structure of the jet's magnetic field and its effect on the emitting gas. Horizon-scale imaging of M87 and SgrA* will enable for the first time the empirical exploration of the relationship between the mass and spin of a black hole and the characteristics of the gas inflow/outflow around it.

Recent millimeter-VLBI observations of Sagittarius A* (Sgr A*) have, for the first time, directly probed distances comparable to the horizon scale of a black hole. This provides unprecedented access to the environment immediately around the horizon of an accreting black hole. We leverage both existing spectral and polarization measurements and our present understanding of accretion theory to produce a suite of generic radiatively inefficient accretion flow (RIAF) models of Sgr A*, which we then fit to these recent millimeter-VLBI observations. We find that if the accretion flow onto Sgr A* is well described by a RIAF model, the orientation and magnitude of the black hole's spin is constrained to a two-dimensional surface in the spin, inclination, position angle parameter space. For each of these we find the likeliest values and their 1-sigma & 2-sigma errors to be a=0(+0.4+0.7), inclination=50(+10+30)(-10-10) degrees, and position angle=-20(+31+107)(-16-29) degrees, when the resulting probability distribution is marginalized over the others. The most probable combination is a=0(+0.2+0.4), inclination=90(-40-50) degrees and position angle=-14(+7+11)(-7-11) degrees, though the uncertainties on these are very strongly correlated, and high probability configurations exist for a variety of inclination angles above 30 degrees and spins below 0.99. Nevertheless, this demonstrates the ability millimeter-VLBI observations, even with only a few stations, to significantly constrain the properties of Sgr A*.

The recent detection of Sagittarius A* at lambda = 1.3 mm on a baseline from Hawaii to Arizona demonstrates that millimeter wavelength very long baseline interferometry (VLBI) can now spatially resolve emission from the innermost accretion flow of the Galactic center region. Here, we investigate the ability of future millimeter VLBI arrays to constrain the spin and inclination of the putative black hole and the orientation of the accretion disk major axis within the context of radiatively inefficient accretion flow (RIAF) models. We examine the range of baseline visibility and closure amplitudes predicted by RIAF models to identify critical telescopes for determining the spin, inclination, and disk orientation of the Sgr A* black hole and accretion disk system. We find that baseline lengths near 3 gigalambda have the greatest power to distinguish amongst RIAF model parameters, and that it will be important to include new telescopes that will form north-south baselines with a range of lengths. If a RIAF model describes the emission from Sgr A*, it is likely that the orientation of the accretion disk can be determined with the addition of a Chilean telescope to the array. Some likely disk orientations predict detectable fluxes on baselines between the continental United States and even a single 10-12 m dish in Chile. The extra information provided from closure amplitudes by a four-antenna array enhances the ability of VLBI to discriminate amongst model parameters.

The super massive black hole candidate, Sagittarius A*, exhibits variability from radio to X-ray wavelengths on time scales that correspond to < 10 Schwarzschild radii. We survey the potential of millimeter-wavelength VLBI to detect and constrain time variable structures that could give rise to such variations, focusing on a model in which an orbiting hot spot is embedded in an accretion disk. Non-imaging algorithms are developed that use interferometric closure quantities to test for periodicity, and applied to an ensemble of hot-spot models that sample a range of parameter space. We find that structural periodicity in a wide range of cases can be detected on most potential VLBI arrays using modern VLBI instrumentation. Future enhancements of mm/sub-mm VLBI arrays including phased array processors to aggregate VLBI station collecting area, increased bandwidth recording, and addition of new VLBI sites all significantly aid periodicity detection. The methods described herein can be applied to other models of Sagittarius A*, including jet outflows and Magneto-Hydrodynamic accretion simulations.

Multiwavelength monitoring observations of Sagittarius A* exhibit variability on timescales of minutes to hours, indicating emission regions localized near the event horizon. (Sub)Millimeter-wavelength VLBI is uniquely suited to probe the environment of the assumed black hole on these scales. We consider a range of orbiting hot-spot and accretion-disk models and find that periodicity in Sgr A* flares is detectable using closure quantities. Our methods are applicable to any model producing source structure changes near the black hole, including jets and magnetohydrodynamic disk instabilities, and suggest that (sub)millimeter VLBI will play a prominent role in investigating Sgr A* near the event horizon.

In a previous paper we demonstrated that non-radial hydrodynamic oscillations of a thermally-supported (Bonnor-Ebert) sphere embedded in a low-density, high-temperature medium persist for many periods. The predicted column density variations and molecular spectral line profiles are similar to those observed in the Bok globule B68 suggesting that the motions in some starless cores may be oscillating perturbations on a thermally supported equilibrium structure. Such oscillations can produce molecular line maps which mimic rotation, collapse or expansion, and thus could make determining the dynamical state from such observations alone difficult. However, while B68 is embedded in a very hot, low-density medium, many starless cores are not, having interior/exterior density contrasts closer to unity. In this paper we investigate the oscillation damping rate as a function of the exterior density. For concreteness we use the same interior model employed in Broderick et al. (2007), with varying models for the exterior gas. We also develop a simple analytical formalism, based upon the linear perturbation analysis of the oscillations, which predicts the contribution to the damping rates due to the excitation of sound waves in the external medium. We find that the damping rate of oscillations on globules in dense molecular environments is always many periods, corresponding to hundreds of thousands of years, and persisting over the inferred lifetimes of the globules.

We present measurements with the VLBA of the variability in the centroid position of Sgr A* relative to a background quasar at 7-mm wavelength. We find an average centroid wander of 71 +/- 45 micro-arcsec for time scales between 50 and 100 min and 113 +/- 50 micro-arcsec for timescales between 100 and 200 min, with no secular trend. These are sufficient to begin constraining the viability of the hot-spot model for the radio variability of Sgr A*. It is possible to rule out hot spots with orbital radii above 15GM_SgrA*/c^2 that contribute more than 30% of the total 7-mm flux. However, closer or less luminous hot spots remain unconstrained. Since the fractional variability of Sgr A* during our observations was ~20% on time scales of hours, the hot-spot model for Sgr A*'s radio variability remains consistent with these limits. Improved monitoring of Sgr A*'s centroid position has the potential to place significant constraints upon the existence and morphology of inhomogeneities in a supermassive black hole accretion flow.

In the absence of an active dynamo, purely poloidal magnetic field configurations are unstable to large-scale dynamical perturbations, and decay via reconnection on an Alfvenic timescale. Nevertheless, a number of classes of dynamo-free stars do exhibit significant, long-lived, surface magnetic fields. Numerical simulations suggest that the large-scale poloidal field in these systems is stabilized by a toroidal component of the field in the stellar interior. Using the principle of conservation of total helicity, we develop a variational principle for computing the structure of the magnetic field inside a conducting sphere surrounded by an insulating vacuum. We show that, for a fixed total helicity, the minimum energy state corresponds to a force-free configuration. We find a simple class of axisymmetric solutions, parametrized by angular and radial quantum numbers. However, these solutions have a discontinuity in the toroidal magnetic field at the stellar surface which will exert a toroidal stress on the surface of the star. We then describe two other classes of solutions, the standard spheromak solutions and ones with fixed surface magnetic fields, the latter being relevant for neutron stars with rigid crusts. We discuss the implications of our results for the structure of neutron star magnetic fields, the decay of fields, and the origin of variability and outbursts in magnetars.

In a previous paper, we modeled the oscillations of a thermally-supported (Bonnor-Ebert) sphere as non-radial, linear perturbations following a standard analysis developed for stellar pulsations. The predicted column density variations and molecular spectral line profiles are similar to those observed in the Bok globule B68 suggesting that the motions in some starless cores may be oscillating perturbations on a thermally supported equilibrium structure. However, the linear analysis is unable to address several questions, among them the stability, and lifetime of the perturbations. In this paper we simulate the oscillations using a three-dimensional numerical hydrodynamic code. We find that the oscillations are damped predominantly by non-linear mode-coupling, and the damping time scale is typically many oscillation periods, corresponding to a few million years, and persisting over the inferred lifetime of gobules.

The gravastar model, which postulates a strongly correlated thin shell of anisotropic matter surrounding a region of anti-de Sitter space, has been proposed as an alternative to black holes. We discuss constraints that present-day observations of well-known black hole candidates place on this model. We focus upon two black hole candidates known to have extraordinarily low luminosities: the supermassive black hole in the Galactic Center, Sagittarius A*, and the stellar-mass black hole, XTE J1118+480. We find that the length scale for modifications of the type discussed in Chapline et al. (2003) must be sub-Planckian.

Submilliarcsecond astrometry and imaging of the black hole Sgr A* at the Galactic Center may become possible in the near future at infrared and submillimetre wavelengths. This resolution is sufficient to observe the silhouette the supermassive black hole in the Galactic center casts upon background emission. However, more exciting is the prospect of observing ``hot spots'' in the accretion flow. Here we discuss how such measurements may be used to test not only the consistency of General Relativity, but also the validity of the Kerr metric in particular.

We consider a model in which Sgr A*, the 3.5x10^6 M_sun supermassive black hole candidate at the Galactic Center, is a compact object with a surface. Given the very low quiescent luminosity of Sgr A* in the near infrared, the existence of a hard surface, even in the limit in which the radius approaches the horizon, places severe constraints upon the steady mass accretion rate in the source, requiring dM/dt < 10^-12 M_sun/yr. This limit is well below the minimum accretion rate needed to power the observed submillimeter luminosity of Sgr A*. We thus argue that Sgr A* does not have a surface, i.e., it must have an event horizon. The argument could be made more restrictive by an order of magnitude with microarcsecond resolution imaging, e.g., with submillimeter VLBI.

If the split, asymmetric molecular spectral line profiles that are seen in many starless cores are interpreted as indicative of global collapse or expansion of the core then one possible implication is that most starless cores have short lifetimes on the order of the collapse or sound crossing time scale. An alternative interpretation of the line profiles as indicative of perturbations on an underlying equilibrium structure leads to the opposite implication, that many cores have long lifetimes. While evidence suggests that some cores are collapsing on a free-fall time scale, we show that observations of some other starless cores can be reproduced by a model of non-radial oscillations about the equilibrium configuration of a pressure-bounded, thermally-supported sphere (Bonnor-Ebert sphere). We model the oscillations as linear perturbations following a standard analysis developed for stellar pulsations and compare the column densities and molecular spectral line profiles predicted from a particular model to observations of the Bok globule B68.

Radiative transfer in curved spacetimes has become increasingly important to understanding high-energy astrophysical phenomena and testing general relativity in the strong field limit. The equations of radiative transfer are physically equivalent to the Boltzmann equation, where the latter has the virtue of being covariant. We show that by a judicious choice of the basis of the phase space, it is generally possible to make the momentum derivatives in the Boltzmann equation vanish along an arbitrary (including nongeodesic) path, thus reducing the problem of radiative transfer along a ray to a path integral in coordinate space.

Sub milli-arcsecond astrometry and imaging of the black hole Sgr A* at the Galactic centre may become possible in the near future at infrared and sub-millimetre wavelengths. Motivated by observations of short-term infrared and X-ray variability of Sgr A*, in a previous paper we computed the expected images and light curves, including polarization, associated with an compact emission region orbiting the central black hole. We extend this work, using a more realistic hot-spot model and including the effects of opacity in the underlying accretion flow. We find that at infrared wavelengths the qualitative features identified by our earlier work are present, namely it is possible to extract the black hole mass and spin from spot images and light curves of the observed flux and polarization. At radio wavelengths, disk opacity produces significant departures from the infrared behaviour, but there are still generic signatures of the black hole properties. Detailed comparison of these results with future data can be used to test general relativity and to improve existing models for the accretion flow in the immediate vicinity of the black hole.

The inferred black hole in the Galactic center spans the largest angle on the sky among all known black holes. Forthcoming observational programs plan to localize or potentially resolve the image of Sgr A* to an exquisite precision, comparable to the scale of the black hole horizon. Here we show that the location of the image centroid of Sgr A* should depend on observing frequency because of relativistic and radiative transfer effects. The same effects introduce a generic dependence of the source polarization on frequency. Future detection of the predicted centroid shift and the polarization dependence on frequency can be used to determine the unknown black hole spin and verify the validity of General Relativity.

Helium star--compact object binaries, and helium star--neutron star binaries in particular, are widely believed to be the progenitors of the observed double neutron star systems. In these, the second neutron star is presumed to be the compact remnant of the helium star supernova. Here, the observational implications of such a supernova are discussed, and in particular are explored as a candidate gamma-ray burst mechanism. In this scenario the supernova results in a transient period of rapid accretion onto the compact object, extracting via magnetic torques its rotational energy at highly super-Eddington luminosities in the form of a narrowly beamed, strongly electromagnetically dominated jet. Compton scattering of supernova photons advected within the ejecta, and photons originating at shocks driven into the ejecta by the jet, will cool the jet and can produce the observed prompt emission characteristics, including the peak--inferred isotropic energy relation, X-ray flash characteristics, subpulse light curves, energy dependent time lags and subpulse broadening, and late time spectral softening. The duration of the burst is limited by the rate of Compton cooling of the jet, eventually creating an optically thick, moderately relativistically expanding fireball which can produce the afterglow emission. If the black hole or neutron star stays bound to a compact remnant, late term light curve variability may be observed as in SN 2003dh.

Images from the vicinity of the black hole horizon at the Galactic centre (Sgr A*) could be obtained in the near future with a Very Large Baseline Array of sub-millimetre telescopes. The recently observed short-term infrared and X-ray variability of the emission from Sgr A* implies that the accretion flow near the black hole is clumpy or unsteady. We calculate the appearance of a compact emission region (bright spot) in a circular orbit around a spinning black hole as a function of orbital radius and orientation. We find that the mass and spin of the black hole can be extracted from their generic signatures on the spot image as well as on the lightcurves of its observed flux and polarisation. The strong-field distortion remains noticeable even when the spot image is smoothed over the expected ~20 microarcsecond resolution of future sub-millimetre observations.

We have developed a parallel, simple, and fast hydrodynamics code for multi-dimensional, self-gravitating, adiabatic flows. Our primary motivation is the study of the non-linear evolution of white dwarf oscillations excited via tidal resonances, typically over hundreds of stellar dynamical times. Consequently, we require long term stability, low diffusivity, and high algorithmic efficiency. An explicit, Eulerian, finite-difference scheme on a regular Cartesian grid fulfills these requirements. It provides uniform resolution throughout the flow, as well as simplifying the computation of the self-gravitational potential, which is done via spectral methods. In this paper, we describe the numerical scheme and present the results of some diagnostic problems. We also demonstrate the stability of a cold white dwarf in three dimensions over hundreds of dynamical times. Finally, we compare the results of the numerical scheme to the linear theory of adiabatic oscillations, finding numerical quality factors on the order of 6000, and excellent agreement with the oscillation frequency obtained by the linear analysis.

We consider the evolution of white dwarfs with compact object companions (specifically black holes with masses up to 10^6 solar masses, neutron stars, and other white dwarfs). We suppose that the orbits are initially quite elliptical and then shrink and circularise under the action of gravitational radiation. During this evolution, the white dwarfs will pass through resonances when harmonics of the orbital frequency match the stellar oscillation eigenfrequencies. As a star passes through these resonances, the associated modes will be excited and can be driven to amplitudes that are so large that there is a back reaction on the orbit which, in turn, limits the growth of the modes. A formalism is presented for describing this dynamical interaction for a non-rotating star in the linear approximation when the orbit can be treated as non-relativistic. A semi-analytical expression is found for computing the resonant energy transfer as a function of stellar and orbital parameters for the regime where back reaction may be neglected. This is used to calculate the results of passage through a sequence of resonances for several hypothetical systems. It is found that the amplitude of the l=m=2 f-mode can be driven into the non-linear regime for appropriate initial conditions. We also discuss where the no back reaction approximation is expected to fail, and the qualitative effects of back reaction.

Accretion onto compact objects plays a central role in high energy astrophysics. In these environments, both general relativistic and plasma effects may have significant impacts upon the spectral and polarimetric properties of the accretion flow. In paper I we presented a fully general relativistic magnetoionic theory, capable of tracing rays in the geometric optics approximation through a magnetised plasma in the vicinity of a compact object. In this paper we discuss how to perform polarised radiative transfer along these rays. In addition we apply the formalism to a barotropic thick disk model, appropriate for low luminosity active galactic nuclei. We find that it is possible to generate large fractional polarisations over the innermost portions of the accretion flow, even when the emission mechanism is unpolarised. This has implications for accreting systems ranging from pulsars and X-ray binaries to AGN.

Accretion onto compact objects plays a central role in high energy astrophysics. In these environments, both general relativistic and plasma effects may have significant impacts upon the propagation of photons. We present a full general relativistic magnetoionic theory, capable of tracing rays in the geometric optics approximation through a magnetised plasma in the vicinity of a compact object. We consider both the cold and warm, ion and pair plasmas. When plasma effects become large the two plasma eignemodes follow different ray trajectories resulting in a large observable polarisation. This has implications for accreting systems ranging from pulsars and X-ray binaries to AGN.

We have developed methods for tracing rays and performing radiative transfer through a magnetoactive plasma in a general relativistic environment. The two electromagnetic plasma modes propagate differently due to a combination of dispersive and gravitational effects. We have found that, when given an appropriate environment surrounding the central black hole, it is indeed possible to generate a significant degree of circular polarisation without an appreciable amount of linear polarisation due to these effects alone.

We present a variational formalism for describing the dynamical evolution of an oscillating star with a point-mass companion in the linear, non-relativistic regime. This includes both the excitation of normal modes and the back-reaction of the modes on the orbit. The general formalism for arbitrary fluid configurations is presented, and then specialized to a homentropic potential flow. Our formalism explicitly identifies and conserves both energy and angular momentum. We also consider corrections to the orbit up to 7/2 post-Newtonian order.

We investigate the effects of very strong magnetic fields upon the equation of state of dense bayonic matter in which hyperons are present. In the presence of a magnetic field, the equation of state above nuclear density is significantly affected both by Landau quantization and magnetic moment interactions, but only for field strengths $B>5\times10^{18}$ G. The former tends to soften the EOS and increase proton and lepton abundances, while the latter produces an overall stiffening of the EOS. Each results in a supression of hyperons relative to the field-free case. The structure of a neutron star is, however, primarily determined by the magnetic field stress. We utilize existing general relativistic magneto-hydrostatic calculations to demonstrate that maximum average fields within a stable neutron are limited to values $B\le 1-3 \times10^{18}$ G. This is not large enough to significantly influence particle compositions or the matter pressure, unless fluctuations dominate the average field strengths in the interior or configurations with significantly larger field gradients are considered.

Recent developments in the spectropolarimetric study of compact objects, specifically black holes (stellar and massive) and neutron stars are reviewed. The lectures are organized around five topics: disks, jets, outflows, neutron stars and black holes. They emphasize physical mechanisms and are intended to bridge the gap between the fundamentals of polarimetry and the phenomenology of observed cosmic sources of polarized radiation, as covered by the other lecturers. There has been considerable recent progress in spectropolarimetry from radio through optical frequencies and this is producing some unique diagnostics of the physical conditions around compact objects. It is argued that there is a great need to develop a correspondingly sensitive polarimetric capability at ultraviolet through gamma-ray energies.

We study the effects of very strong magnetic fields on the equation of state (EOS) in multicomponent, interacting matter by developing a covariant description for the inclusion of the anomalous magnetic moments of nucleons. For the description of neutron star matter, we employ a field-theoretical approach which permits the study of several models which differ in their behavior at high density. Effects of Landau quantization in ultra-strong magnetic fields ($B>10^{14}$ Gauss) lead to a reduction in the electron chemical potential and a substantial increase in the proton fraction. We find the generic result for $B>10^{18}$ Gauss that the softening of the EOS caused by Landau quantization is overwhelmed by stiffening due to the incorporation of the anomalous magnetic moments of the nucleons. In addition, the neutrons become completely spin polarized. The inclusion of ultra-strong magnetic fields leads to a dramatic increase in the proton fraction, with consequences for the direct Urca process and neutron star cooling. The magnetization of the matter never appears to become very large, as the value of $|H/B|$ never deviates from unity by more than a few percent. Our findings have implications for the structure of neutron stars in the presence of large frozen-in magnetic fields.