• ### The Current Ability to Test Theories of Gravity with Black Hole Shadows(1804.05812)

April 16, 2018 gr-qc, astro-ph.GA, astro-ph.HE
Our Galactic Center, Sagittarius A* (Sgr A*), is believed to harbour a supermassive black hole (BH), as suggested by observations tracking individual orbiting stars. Upcoming sub-millimetre very-long-baseline-interferometry (VLBI) images of Sgr A* carried out by the Event-Horizon-Telescope Collaboration (EHTC) are expected to provide critical evidence for the existence of this supermassive BH. We assess our present ability to use EHTC images to determine if they correspond to a Kerr BH as predicted by Einstein's theory of general relativity (GR) or to a BH in alternative theories of gravity. To this end, we perform general-relativistic magnetohydrodynamical (GRMHD) simulations and use general-relativistic radiative transfer (GRRT) calculations to generate synthetic shadow images of a magnetised accretion flow onto a Kerr BH. In addition, and for the first time, we perform GRMHD simulations and GRRT calculations for a dilaton BH, which we take as a representative solution of an alternative theory of gravity. Adopting the VLBI configuration from the 2017 EHTC campaign, we find that it could be extremely difficult to distinguish between BHs from different theories of gravity, thus highlighting that great caution is needed when interpreting BH images as tests of GR.
• ### Antenna Design and Implementation for the Future Space Ultra-Long Wavelength Radio Telescope(1802.07640)

Feb. 22, 2018 astro-ph.IM
In radio astronomy, the Ultra-Long Wavelengths (ULW) regime of longer than 10 m (frequencies below 30 MHz), remains the last virtually unexplored window of the celestial electromagnetic spectrum. The strength of the science case for extending radio astronomy into the ULW window is growing. However, the opaqueness of the Earth's ionosphere makes ULW observations by ground-based facilities practically impossible. Furthermore, the ULW spectrum is full of anthropogenic radio frequency interference (RFI). The only radical solution for both problems is in placing an ULW astronomy facility in space. We present a concept of a key element of a space-borne ULW array facility, an antenna that addresses radio astronomical specifications. A tripole-type antenna and amplifier are analysed as a solution for ULW implementation. A receiver system with a low power dissipation is discussed as well. The active antenna is optimized to operate at the noise level defined by the celestial emission in the frequency band 1 - 30 MHz. Field experiments with a prototype tripole antenna enabled estimates of the system noise temperature. They indicated that the proposed concept meets the requirements of a space-borne ULW array facility.
• ### RAPTOR I: Time-dependent radiative transfer in arbitrary spacetimes(1801.10452)

Jan. 31, 2018 astro-ph.HE
Observational efforts to image the immediate environment of a black hole at the scale of the event horizon benefit from the development of efficient imaging codes that are capable of producing synthetic data, which may be compared with observational data. We aim to present RAPTOR, a new public code that produces accurate images, animations, and spectra of relativistic plasmas in strong gravity by numerically integrating the equations of motion of light rays and performing time-dependent radiative transfer calculations along the rays. The code is compatible with any analytical or numerical spacetime. It is hardware-agnostic and may be compiled and run both on GPUs and CPUs. We describe the algorithms used in RAPTOR and test the code's performance. We have performed a detailed comparison of RAPTOR output with that of other radiative-transfer codes and demonstrate convergence of the results. We then applied RAPTOR to study accretion models of supermassive black holes, performing time-dependent radiative transfer through general relativistic magneto-hydrodynamical (GRMHD) simulations and investigating the expected observational differences between the so-called fast-light and slow-light paradigms. Using RAPTOR to produce synthetic images and light curves of a GRMHD model of an accreting black hole, we find that the relative difference between fast-light and slow-light light curves is less than 5%. Using two distinct radiative-transfer codes to process the same data, we find integrated flux densities with a relative difference less than 0.01%. For two-dimensional GRMHD models, such as those examined in this paper, the fast-light approximation suffices as long as errors of a few percent are acceptable. The convergence of the results of two different codes demonstrates that they are, at a minimum, consistent.
• ### Imaging black holes: past, present and future(1801.03298)

Jan. 11, 2018 astro-ph.IM, astro-ph.HE
This paper briefly reviews past, current, and future efforts to image black holes in the radio regime. Black holes seem like mystical objects, but they are an integral part of current astrophysics and are at the center of attempts to unify quantum physics and general relativity. Yet, nobody has ever seen a black hole. What do they look like? Initially, this question seemed more of an academic nature. However, this has changed over the past two decades. Observations and theoretical considerations suggest that the supermassive black hole, Sgr A*, in the center of our Milky Way is surrounded by a compact, foggy emission region radiating at and above 230 GHz. It has been predicted that the event horizon of Sgr A* should cast its shadow onto that emission region, which could be detectable with a global VLBI array of radio telescopes. In contrast to earlier pictures of black holes, that dark feature is not supposed to be due to a hole in the accretion flow, but would represent a true negative image of the event horizon. Currently, the global Event Horizon Telescope consortium is attempting to make such an image. In the future those images could be improved by adding more telescopes to the array, in particular at high sites in Africa. Ultimately, a space array at THz frequencies, the Event Horizon Imager, could produce much more detailed images of black holes. In combination with numerical simulations and precise measurements of the orbits of stars - ideally also of pulsars - these images will allow us to study black holes with unprecedented precision.
• ### General relativistic magnetohydrodynamical $\kappa$-jet models for Sgr A*(1712.02266)

Dec. 6, 2017 astro-ph.HE
The observed spectral energy distribution of an accreting supermassive black hole typically forms a power-law spectrum in the Near Infrared (NIR) and optical wavelengths, that may be interpreted as a signature of accelerated electrons along the jet. However, the details of acceleration remain uncertain. In this paper, we study the radiative properties of jets produced in axisymmetric GRMHD simulations of hot accretion flows onto underluminous supermassive black holes both numerically and semi-analytically, with the aim of investigating the differences between models with and without accelerated electrons inside the jet. We assume that electrons are accelerated in the jet regions of our GRMHD simulation. To model them, we modify the electrons' distribution function in the jet regions from a purely relativistic thermal distribution to a combination of a relativistic thermal distribution and the $\kappa$-distribution function. Inside the disk, we assume a thermal distribution for the electrons. We calculate jet spectra and synchrotron maps by using the ray tracing code {\tt RAPTOR}, and compare the synthetic observations to observations of Sgr~A*. Finally, we compare numerical models of jets to semi-analytical ones. We find that in the $\kappa$-jet models, the radio-emitting region size, radio flux, and spectral index in NIR/optical bands increase for decreasing values of the $\kappa$ parameter, which corresponds to a larger amount of accelerated electrons. The model with $\kappa = 3.5$, $\eta_{\rm acc}=5-10\%$ (the percentage of electrons that are accelerated), and observing angle $i = 30^{\rm o}$ fits the observed Sgr~A* emission in the flaring state from the radio to the NIR/optical regimes, while $\kappa = 3.5$, $\eta_{\rm acc}< 1\%$, and observing angle $i = 30^{\rm o}$ fit the upper limits in quiescence.
• ### Large magnetic field variations towards the Galactic Centre magnetar, PSR J1745-2900(1711.10323)

Nov. 28, 2017 astro-ph.HE
Polarised radio emission from PSR J1745-2900 has already been used to investigate the strength of the magnetic field in the Galactic Centre, close to Sagittarius A*. Here we report how persistent radio emission from this magnetar, for over four years since its discovery, has revealed large changes in the observed Faraday rotation measure, by up to 3500 rad m$^{-2}$ (a five per cent fractional change). From simultaneous analysis of the dispersion measure, we determine that these fluctuations are dominated by variations in the projected magnetic field, rather than the integrated free electron density, along the changing line of sight to the rapidly moving magnetar. From a structure function analysis of rotation measure variations, and a recent epoch of rapid change of rotation measure, we determine a minimum scale of magnetic fluctuations of size ~ 2 au at the Galactic Centre distance, inferring PSR J1745-2900 is just ~ 0.1 pc behind an additional scattering screen.
• ### Quantifying Intrinsic Variability of Sagittarius A* using Closure Phase Measurements of the Event Horizon Telescope(1708.01056)

Aug. 3, 2017 astro-ph.HE
General relativistic magnetohydrodynamic (GRMHD) simulations of accretion disks and jets associated with supermassive black holes show variability on a wide range of timescales. On timescales comparable to or longer than the gravitational timescale $t_G=GM/c^3$, variation may be dominated by orbital dynamics of the inhomogeneous accretion flow. Turbulent evolution within the accretion disk is expected on timescales comparable to the orbital period, typically an order of magnitude larger than $t_G$. For Sgr A*, $t_G$ is much shorter than the typical duration of a VLBI experiment, enabling us to study this variability within a single observation. Closure phases, the sum of interferometric visibility phases on a triangle of baselines, are particularly useful for studying this variability. In addition to a changing source structure, variations in observed closure phase can also be due to interstellar scattering, thermal noise, and the changing geometry of projected baselines over time due to Earth rotation. We present a metric that is able to distinguish the latter two from intrinsic or scattering variability. This metric is validated using synthetic observations of GRMHD simulations of Sgr A*. When applied to existing multi-epoch EHT data of Sgr A*, this metric shows that the data are most consistent with source models containing intrinsic variability from source dynamics, interstellar scattering, or a combination of those. The effects of black hole inclination, orientation, spin, and morphology (disk or jet) on the expected closure phase variability are also discussed.
• ### Simultaneous Monitoring of X-ray and Radio Variability in Sagittarius A*(1707.01937)

July 6, 2017 astro-ph.HE
Monitoring of Sagittarius A* from X-ray to radio wavelengths has revealed structured variability --- including X-ray flares --- but it is challenging to establish correlations between them. Most studies have focused on variability in the X-ray and infrared, where variations are often simultaneous, and because long time series at sub-millimeter and radio wavelengths are limited. Previous work on sub-mm and radio variability hints at a lag between X-ray flares and their candidate sub-millimeter or radio counterparts, with the long wavelength data lagging the X-ray. However, there is only one published time lag between an X-ray flare and a possible radio counterpart. Here we report 9 contemporaneous X-ray and radio observations of Sgr A*. We detect significant radio variability peaking $\gtrsim$176 minutes after the brightest X-ray flare ever detected from Sgr A*. We also report other potentially associated X-ray and radio variability, with the radio peaks appearing $\lesssim$80 minutes after these weaker X-ray flares. Taken at face value, these results suggest that stronger X-ray flares lead to longer time lags in the radio. However, we also test the possibility that the variability at X-ray and radio wavelengths is not temporally correlated. We cross-correlate data from mismatched X-ray and radio epochs and obtain comparable correlations to the matched data. Hence, we find no overall statistical evidence that X-ray flares and radio variability are correlated, underscoring a need for more simultaneous, long duration X-ray--radio monitoring of Sgr A*.
• ### The Black Hole Accretion Code(1611.09720)

June 18, 2017 gr-qc, astro-ph.HE
We present the black hole accretion code (BHAC), a new multidimensional general-relativistic magnetohydrodynamics module for the MPI-AMRVAC framework. BHAC has been designed to solve the equations of ideal general-relativistic magnetohydrodynamics in arbitrary spacetimes and exploits adaptive mesh refinement techniques with an efficient block-based approach. Several spacetimes have already been implemented and tested. We demonstrate the validity of BHAC by means of various one-, two-, and three-dimensional test problems, as well as through a close comparison with the HARM3D code in the case of a torus accreting onto a black hole. The convergence of a turbulent accretion scenario is investigated with several diagnostics and we find accretion rates and horizon-penetrating fluxes to be convergent to within a few percent when the problem is run in three dimensions. Our analysis also involves the study of the corresponding thermal synchrotron emission, which is performed by means of a new general-relativistic radiative transfer code, BHOSS. The resulting synthetic intensity maps of accretion onto black holes are found to be convergent with increasing resolution and are anticipated to play a crucial role in the interpretation of horizon-scale images resulting from upcoming radio observations of the source at the Galactic Center.
• ### Overview of lunar detection of ultra-high energy particles and new plans for the SKA(1704.05336)

April 18, 2017 astro-ph.IM, astro-ph.HE
The lunar technique is a method for maximising the collection area for ultra-high-energy (UHE) cosmic ray and neutrino searches. The method uses either ground-based radio telescopes or lunar orbiters to search for Askaryan emission from particles cascading near the lunar surface. While experiments using the technique have made important advances in the detection of nanosecond-scale pulses, only at the very highest energies has the lunar technique achieved competitive limits. This is expected to change with the advent of the Square Kilometre Array (SKA), the low-frequency component of which (SKA-low) is predicted to be able to detect an unprecedented number of UHE cosmic rays. In this contribution, the status of lunar particle detection is reviewed, with particular attention paid to outstanding theoretical questions, and the technical challenges of using a giant radio array to search for nanosecond pulses. The activities of SKA's High Energy Cosmic Particles Focus Group are described, as is a roadmap by which this group plans to incorporate this detection mode into SKA-low observations. Estimates for the sensitivity of SKA-low phases 1 and 2 to UHE particles are given, along with the achievable science goals with each stage. Prospects for near-future observations with other instruments are also described.
• ### Faraday rotation in GRMHD simulations of the jet launching zone of M87(1703.02390)

March 7, 2017 astro-ph.HE
Non-VLBI measurements of Faraday rotation at mm wavelengths have been used to constrain mass accretion rates ($\mdot$) onto supermassive black holes in the centre of the Milky Way and in the centre of M87. We constructed general relativistic magnetohydrodynamics models for these sources that qualitatively well describe their spectra and radio/mm images invoking a coupled jet-disk system. Using general relativistic polarized radiative transfer, we now also model the observed mm rotation measure (RM) of M87. The models are tied to the observed radio flux, however, electron temperature and accretion rate are degenerate parameters and are allowed to vary. For the inferred low viewing angles of the M87 jet, the RM is low even as the black hole $\mdot$ increases by a factor of $\simeq100$. In jet-dominated models, the observed linear polarization is produced in the forward-jet, while the dense accretion disk depolarizes the bulk of the near-horizon scale emission which originates in the counter-jet. In the jet-dominated models, with increasing $\mdot$ and increasing Faraday optical depth one is progressively sensitive only to polarized emission in the forward-jet, keeping the measured RM relatively constant. The jet-dominated model reproduces a low net-polarization of $\simeq1$ per cent and RMs in agreement with observed values due to Faraday depolarization, however, with $\mdot$ much larger than the previously inferred limit of $9\times10^{-4}\,\mdotu$. All jet-dominated models produce much higher RMs for inclination angles $i\gtrsim30^\circ$, where the line-of-sight passes through the accretion flow, thereby providing independent constraints on the viewing geometry of the M87 jet.
• ### Ultimate precision in cosmic-ray radio detection --- the SKA(1608.08869)

Feb. 6, 2017 hep-ex, astro-ph.IM, astro-ph.HE
As of 2023, the low-frequency part of the Square Kilometre Array will go online in Australia. It will constitute the largest and most powerful low-frequency radio-astronomical observatory to date, and will facilitate a rich science programme in astronomy and astrophysics. With modest engineering changes, it will also be able to measure cosmic rays via the radio emission from extensive air showers. The extreme antenna density and the homogeneous coverage provided by more than 60,000 antennas within an area of one km$^2$ will push radio detection of cosmic rays in the energy range around 10$^{17}$ eV to ultimate precision, with superior capabilities in the reconstruction of arrival direction, energy, and an expected depth-of-shower-maximum resolution of 6~g/cm${^2}$.
• ### Scale-invariant radio jets and varying black hole spin(1610.08652)

Oct. 27, 2016 astro-ph.HE
• ### Asymmetric structure in Sgr A* at 3mm from closure phase measurements with VLBA, GBT and LMT(1608.06515)

Aug. 23, 2016 astro-ph.GA
We present the results of a closure phase analysis of 3 mm very long baseline interferometry (VLBI) measurements performed on Sagittarius A* (Sgr A*). We have analyzed observations made in May 2015 using the Very Long Baseline Array, the Robert C. Byrd Green Bank Telescope and the Large Millimeter Telescope Alfonso Serrano and obtained non-zero closure phase measurements on several station triangles - indicative of a non-point-symmetric source structure. The data are fitted with an asymmetric source structure model in Sgr A*, represented by a simple two-component model, which favours a fainter component due East of the main source. This result is discussed in light of a scattering screen with substructure or an intrinsically asymmetric source.
• ### The lunar Askaryan technique with the Square Kilometre Array(1608.02408)

Aug. 8, 2016 astro-ph.IM, astro-ph.HE
The lunar Askaryan technique is a method to study the highest-energy cosmic rays, and their predicted counterparts, the ultra-high-energy neutrinos. By observing the Moon with a radio telescope, and searching for the characteristic nanosecond-scale Askaryan pulses emitted when a high-energy particle interacts in the outer layers of the Moon, the visible lunar surface can be used as a detection area. Several previous experiments, at Parkes, Goldstone, Kalyazin, Westerbork, the ATCA, Lovell, LOFAR, and the VLA, have developed the necessary techniques to search for these pulses, but existing instruments have lacked the necessary sensitivity to detect the known flux of cosmic rays from such a distance. This will change with the advent of the SKA. The Square Kilometre Array (SKA) will be the world's most powerful radio telescope. To be built in southern Africa, Australia and New Zealand during the next decade, it will have an unsurpassed sensitivity over the key 100 MHz to few-GHZ band. We introduce a planned experiment to use the SKA to observe the highest-energy cosmic rays and, potentially, neutrinos. The estimated event rate will be presented, along with the predicted energy and directional resolution. Prospects for directional studies with phase 1 of the SKA will be discussed, as will the major technical challenges to be overcome to make full use of this powerful instrument. Finally, we show how phase 2 of the SKA could provide a vast increase in the number of detected cosmic rays at the highest energies, and thus to provide new insight into their spectrum and origin.
• We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principle calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.
• ### Lorentz factor distribution of blazars from the optical Fundamental plane of black hole activity(1606.06147)

June 20, 2016 astro-ph.HE
Blazar radiation is dominated by a relativistic jet which can be modeled at first approximation using just two intrinsic parameters - the Lorentz factor $\Gamma$ and the viewing angle $\theta$. Blazar jet observations are often beamed due to relativistic effects, complicating the understanding of these intrinsic properties. The most common way to estimate blazar Lorentz factors needs the estimation of apparent jet speeds and Doppler beaming factors. We present a new and independent method of constructing the blazar Lorentz factor distribution, using the optical fundamental plane of black hole activity. The optical fundamental plane is a plane stretched out by both the supermassive black holes and the X-ray binaries, in the 3D space provided by their [OIII] line luminosity, radio luminosity and black hole mass. We use the intrinsic radio luminosity obtained from the optical fundamental plane to constrain the boosting parameters of the VLBA Imaging and Polarimetry Survey (VIPS) blazar sample. We find a blazar bulk Lorentz factor distribution in the form of a power law as $N(\Gamma) \propto \Gamma^{-2.1 \pm 0.4}$ for the $\Gamma$ range of 1 to 40. We also discuss the viewing angle distribution of the blazars and the dependence of our results on the input parameters.
• ### Linearly Polarized Millimeter and Submillimeter Continuum Emission of Sgr A* Constrained by ALMA(1605.05544)

May 18, 2016 astro-ph.GA, astro-ph.HE
Our aim is to characterize the polarized continuum emission properties including intensity, polarization position angle, and polarization percentage of Sgr A* at $\sim$100 (3.0 mm), $\sim$230 (1.3 mm), $\sim$345 (0.87 mm), $\sim$500 (0.6 mm), and $\sim$700 GHz (0.43 mm). We report continuum emission properties of Sgr A* at the above frequency bands, based on the Atacama Large Millimeter Array (ALMA) observations. We measured flux densities of Sgr A* from ALMA single pointing and mosaic observations. We performed sinusoidal fittings to the observed (XX-YY)/I intensity ratios, to derive the polarization position angles and polarization percentages. We successfully detect polarized continuum emission from all observed frequency bands. We observed lower Stokes I intensity at $\sim$700 GHz than that at $\sim$500 GHz, which suggests that emission at $\gtrsim$500 GHz is from optically thin part of a synchrotron emission spectrum. Both the Stokes I intensity and the polarization position angle at our highest observing frequency of $\sim$700 GHz, may be varying with time. However, we do not yet detect variation in the polarization percentage at $>$500 GHz. The polarization percentage at $\sim$700 GHz is likely lower than that at $\sim$500 GHz. By comparing the $\sim$500 GHz and $\sim$700 GHz observations with the observations at lower frequency bands, we suggest that the intrinsic polarization position angle of Sgr A* is varying with time. This paper also reports the measurable polarization properties from the observed calibration quasars. The future simultaneous multi-frequency polarization observations are required for clarifying the time and frequency variation of polarization position angle and polarization percentage.
• ### The 492 GHz emission of Sgr A* constrained by ALMA(1604.00599)

May 6, 2016 astro-ph.GA, astro-ph.HE
We report linearly polarized continuum emission properties of Sgr A* at $\sim$492 GHz, based on the Atacama Large Millimeter Array (ALMA) observations. We used the observations of the likely unpolarized continuum emission of Titan, and the observations of C\textsc{i} line emission, to gauge the degree of spurious polarization. The Stokes I flux of 3.6$\pm$0.72 Jy during our run is consistent with extrapolations from the previous, lower frequency observations. We found that the continuum emission of Sgr A* at $\sim$492 GHz shows large amplitude differences between the XX and the YY correlations. The observed intensity ratio between the XX and YY correlations as a function of parallactic angle may be explained by a constant polarization position angle of $\sim$158$^{\circ}$$\pm3^{\circ}. The fitted polarization percentage of Sgr A* during our observational period is 14\%\pm1.2\%. The calibrator quasar J1744-3116 we observed at the same night can be fitted to Stokes I = 252 mJy, with 7.9\%\pm0.9\% polarization in position angle P.A. = 4.1^{\circ}$$\pm$4.2$^{\circ}$. The observed polarization percentage and polarization position angle in the present work appear consistent with those expected from longer wavelength observations in the period of 1999-2005. In particular, the polarization position angle at 492 GHz, expected from the previously fitted 167$^{\circ}$$\pm$7$^{\circ}$ intrinsic polarization position angle and (-5.6$\pm$0.7)$\times$10$^{5}$ rotation measure, is 155$^{+9}_{-8}$, which is consistent with our new measurement of polarization position angle within 1$\sigma$. The polarization percentage and the polarization position angle may be varying over the period of our ALMA 12m Array observations, which demands further investigation with future polarization observations.
• ### General relativistic magnetohydrodynamical simulations of the jet in M87(1510.07243)

April 29, 2016 astro-ph.HE
(abridged) The connection between black hole, accretion disk, and radio jet can be best constrained by fitting models to observations of nearby low luminosity galactic nuclei, in particular the well studied sources Sgr~A* and M87. There has been considerable progress in modeling the central engine of active galactic nuclei by an accreting supermassive black hole coupled to a relativistic plasma jet. However, can a single model be applied to a range of black hole masses and accretion rates? Here we want to compare the latest three-dimensional numerical model, originally developed for Sgr A* in the center of the Milky Way, to radio observations of the much more powerful and more massive black hole in M87. We postprocess three-dimensional GRMHD models of a jet-producing radiatively inefficient accretion flow around a spinning black hole using relativistic radiative transfer and ray-tracing to produce model spectra and images. As a key new ingredient to these models, we allow the proton-electron coupling in these simulations depend on the magnetic properties of the plasma. We find that the radio emission in M87 is well described by a combination of a two-temperature accretion flow and a hot single-temperature jet. The model fits the basic observed characteristics of the M87 radio core. The best fit model has a mass-accretion rate of Mdot approx 9x10^{-3} MSUN/YR and a total jet power of P_j \sim 10^{43} erg/s. Emission at 1.3mm is produced by the counter jet close to the event horizon. Its characteristic crescent shape surrounding the black hole shadow could be resolved by future millimeter-wave VLBI experiments. The model was successfully derived from one for the supermassive black hole in center of the Milky Way by appropriately scaling mass and accretion rate. This suggests the possibility that this model could also apply to a larger range of low-luminosity black holes.
• ### What are the megahertz peaked-spectrum sources?(1604.00171)

Megahertz peaked-spectrum (MPS) sources have spectra that peak at frequencies below 1 GHz in the observer's frame and are believed to be radio-loud active galactic nuclei (AGN). We recently presented a new method to search for high-redshift AGN by identifying unusually compact MPS sources. In this paper, we present European VLBI Network (EVN) observations of 11 MPS sources which we use to determine their sizes and investigate the nature of the sources with ~10 mas resolution. Of the 11 sources, we detect nine with the EVN. Combining the EVN observations with spectral and redshift information, we show that the detected sources are all AGN with linear sizes smaller than 1.1 kpc and are likely young. This shows that low-frequency colour-colour diagrams are an easy and efficient way of selecting small AGN and explains our high detection fraction (82%) in comparison to comparable surveys. Finally we argue that the detected sources are all likely compact symmetric objects and that none of the sources are blazars.
• ### The Proper Motion of the Galactic Center Pulsar Relative to Sagittarius A*(1411.0399)

March 16, 2016 astro-ph.HE
We measure the proper motion of the pulsar PSR J1745-2900 relative to the Galactic Center massive black hole, Sgr A*, using the Very Long Baseline Array (VLBA). The pulsar has a transverse velocity of 236 +/- 11 km s^-1 at position angle 22 +/- 2 deg East of North at a projected separation of 0.097 pc from Sgr A*. Given the unknown radial velocity, this transverse velocity measurement does not conclusively prove that the pulsar is bound to Sgr A*; however, the probability of chance alignment is very small. We do show that the velocity and position is consistent with a bound orbit originating in the clockwise disk of massive stars orbiting Sgr A* and a natal velocity kick of <~ 500 km s^-1. An origin among the isotropic stellar cluster is possible but less probable. If the pulsar remains radio-bright, multi-year astrometry of PSR J1745-2900 can detect its acceleration and determine the full three-dimensional orbit. We also demonstrate that PSR J1745-2900 exhibits the same angular broadening as Sgr A* over a wavelength range of 3.6 cm to 0.7 cm, further confirming that the two sources share the same interstellar scattering properties. Finally, we place the first limits on the presence of a wavelength-dependent shift in the position of Sgr A*, i.e., the core shift, one of the expected properties of optically-thick jet emission. Our results for PSR J1745-2900 support the hypothesis that Galactic Center pulsars will originate from the stellar disk and deepens the mystery regarding the small number of detected Galactic Center pulsars.
• ### Persistent Asymmetric Structure of Sagittarius A* on Event Horizon Scales(1602.05527)

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.
• Neutrinos in the cosmic ray flux with energies near 1 EeV and above are detectable with the Surface Detector array of the Pierre Auger Observatory. We report here on searches through Auger data from 1 January 2004 until 20 June 2013. No neutrino candidates were found, yielding a limit to the diffuse flux of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound predictions. Neutrino identification is attempted using the broad time-structure of the signals expected in the SD stations, and is efficiently done for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for "Earth-skimming" neutrino interactions in the case of tau neutrinos. In this paper the searches for downward-going neutrinos in the zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an $E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times 10^{19}$ eV is $E_\nu^2 dN_\nu/dE_\nu < 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~ s^{-1}~ sr^{-1}}$.