M87 Jet Research Articles

The Structure and Propagation of the Misaligned Jet M87

Due to its proximity, M87 is a prime target for next-generation high-resolution VLBI at short millimeter wavelengths, by which the jet launching region and the black hole shadow are expected to be resolved and imaged sometime soon. Along with this situation, high-quality VLBI imaging and monitoring at lower frequencies play an important role in complementing the high-frequency data. Here, we present our recent and ongoing observational studies of the M87 jet on pc-to-subpc scales based on ultra-deep VLBI imaging programs at 86 GHz and 15 GHz. The high-dynamic-range images have allowed us to obtain some remarkably improved views on this jet. We also introduce the KVN and VERA Array (KaVA), a new regularly-operating VLBI network in East Asia, which is quite suitable for studying the structure and propagation of relativistic jets. Some early results from our pilot study for M87—including the detection of superluminal motions near the jet base—implying an efficient magnetic-to-kinetic conversion at these scales, are reported.

Conical Stream of the Two-Sided Jets in NGC 4261 over the Range of 103–109 Schwarzschild Radii

We report the jet width profile of of the nearby ( ∼ 30 Mpc ) AGN NGC 4261 for both the approaching jet and the counter jet at radial distances ranging from ∼ 10 3 – 10 9 Schwarzschild radius ( R S ) from the central engine. Our Very Large Array (VLA) and Very Long Baseline Array (VLBA) observations reveal that the jets maintain a conical structure on both sides over the range 10 3 – 10 9 R S without any structural transition (i.e., parabolic to conical) like in the approaching jet in M87. Thus, NGC 4261 will provide a unique opportunity to examine the conical jet hypothesis in blazars, while it may require some additional consideration on the acceleration and collimation process in AGN jets.

STRUCTURAL TRANSITION IN THE NGC 6251 JET: AN INTERPLAY WITH THE SUPERMASSIVE BLACK HOLE AND ITS HOST GALAXY

ABSTRACT The structure of the NGC 6251 jet on the milliarcsecond scale is investigated using images taken with the European VLBI Network and the Very Long Baseline Array. We detect a structural transition of the jet from a parabolic to a conical shape at a distance of (1–2) × 105 times the Schwarzschild radius from the central engine, which is close to the sphere of gravitational influence of the supermassive black hole (SMBH). We also examine the jet pressure profiles with the synchrotron minimum energy assumption to discuss the physical origin of the structural transition. The NGC 6251 jet, together with the M87 jet, suggests a fundamental process of structural transition in the jets of active galactic nuclei (AGNs). Collimated AGN jets are characterized by their external galactic medium, showing that AGN jets interplay with the SMBH and its host galaxy.

Reconfinement of highly magnetized jets: implications for HST-1 in M87

Stationary features are occasionally observed in AGN jets. A notable example is the HST-1 knot in M87. Such features are commonly interpreted as re-confinement shocks in hydrodynamic jets or focusing nozzles in Poynting jets. In this paper we compute the structure and Lorentz factor of a highly magnetized jet confined by external pressure having a profile that flattens abruptly at some radius. We find the development of strong oscillations upon transition from the steeper to the flatter pressure profile medium. Analytic formula is derived for the location of the nodes of these oscillations. We apply the model to the M87 jet and show that if the jet remains magnetically dominated up to sub-kpc scales, then focusing is expected. The location of the HST-1 knot can be reconciled with recent measurements of the pressure profile around the Bondi radius if the jet luminosity satisfies $L_j\simeq10^{43}$ erg/s. However, we find that magnetic domination at the collimation break implies a Lorentz factor in excess of $10^2$, atypical to FRI sources. A much lower value of the asymptotic Lorentz factor would require substantial loading close to the black hole. In that case HST-1 may be associated with a collimation nozzle of a hydrodynamic flow.

MULTI-WAVELENGTH POLARIMETRY AND SPECTRAL STUDY OF THE M87 JET DURING 2002–2008*

ABSTRACT We present a multi-wavelength polarimetric and spectral study of the M87 jet obtained at sub-arcsecond resolution between 2002 and 2008. The observations include multi-band archival VLA polarimetry data sets along with Hubble Space Telescope (HST) imaging polarimetry. These observations have better angular resolution than previous work by factors of 2–3 and in addition, allow us to explore the time domain. These observations envelop the huge flare in HST-1 located 0.″86 from the nucleus. The increased resolution enables us to view more structure in each knot, showing several resolved sub-components. We also see apparent helical structure in the polarization vectors in several knots, with polarization vectors turning either clockwise or counterclockwise near the flux maxima in various places as well as showing filamentary undulations. Some of these characteristics are correlated with flux and polarization maxima while others are not. We also examine the total flux and fractional polarization and look for changes in both radio and optical since the observations of Perlman et al. (1999) and test them against various models based on shocks and instabilities in the jet. Our results are broadly consistent with previous spine-sheath models and recollimation shock models; however, they require additional combinations of features to explain the observed complexity, e.g., shearing of magnetic field lines near the jet surface and compression of the toroidal component near shocks. In particular, in many regions we find apparently helical features both in total flux and polarization. We discuss the physical interpretation of these features.

Kinematics of the jet in M 87 on scales of 100–1000 Schwarzschild radii

Very long baseline interferometry (VLBI) imaging of radio emission from extragalactic jets provides a unique probe of physical mechanisms governing the launching, acceleration, and collimation of relativistic outflows. The two-dimensional structure and kinematics of the jet in M\,87 (NGC\,4486) have been studied by applying the Wavelet-based Image Segmentation and Evaluation (WISE) method to 11 images obtained from multi-epoch Very Long Baseline Array (VLBA) observations made in January-August 2007 at 43 GHz ($\lambda = 7$ mm). The WISE analysis recovers a detailed two-dimensional velocity field in the jet in M\,87 at sub-parsec scales. The observed evolution of the flow velocity with distance from the jet base can be explained in the framework of MHD jet acceleration and Poynting flux conversion. A linear acceleration regime is observed up to $z_{obs} \sim 2$\,mas. The acceleration is reduced at larger scales, which is consistent with saturation of Poynting flux conversion. Stacked cross correlation analysis of the images reveals a pronounced stratification of the flow. The flow consists of a slow, mildly relativistic layer (moving at $\beta \sim 0.5\,c$), associated either with instability pattern speed or an outer wind, and a fast, accelerating stream line (with $\beta \sim 0.92$, corresponding to a bulk Lorentz factor $\gamma \sim 2.5$). A systematic difference of the apparent speeds in the northern and southern limbs of the jet is detected, providing evidence for jet rotation. The angular velocity of the magnetic field line associated with this rotation suggests that the jet in M87 is launched in the inner part of the disk, at a distance $r_0 \sim 5\, R_\mathrm{s}$ from the central engine. The combined results of the analysis imply that MHD acceleration and conversion of Poynting flux to kinetic energy play the dominant roles in collimation and acceleration of the flow in M\,87.

Observations of the Structure and Dynamics of the Inner M87 Jet

M87 is the best source in which to study a jet at high resolution in gravitational units because it has a very high mass black hole and is nearby. The angular size of the black hole is second only to Sgr A*, which does not have a strong jet. The jet structure is edge brightened with a wide opening angle base and a weak counterjet. We have roughly annual observations for 17 years plus intensive monitoring at three week intervals for a year and five day intervals for 2.5 months made with the Very Long Baseline Array (VLBA) at 43 GHz. The inner jet shows very complex dynamics, with apparent motions both along and across the jet. Speeds from zero to over 2c are seen, with acceleration observed over the first 3 milli-arcseconds. The counterjet decreases in brightness much more rapidly than the main jet, as is expected from relativistic beaming in an accelerating jet oriented near the line-of-sight. Details of the structure and dynamics are discussed. The roughly annual observations show side-to-side motion of the whole jet with a characteristic time scale of about 9 years.

Resolving the Base of the Relativistic Jet in M87 at 6Rsch Resolution with Global mm-VLBI

M87 is one of the nearest radio galaxies with a central Super-Massive Black Hole (SMBH) and a prominent relativistic jet. Due to its close distance to the observer and the large SMBH mass, the source is one of the best laboratories to obtain strong observational constraints on the theoretical models for the formation and evolution of the AGN jets. In this article, we present preliminary results from our ongoing observational study about the innermost jet of M87 at an ultra-high resolution of ∼50 μ as achieved by the Global Millimeter-Very Long Baseline Interferometry Array (GMVA). The data obtained between 2004 and 2015 clearly show limb-brightened jets at extreme resolution and sensitivity. Our preliminary analysis reveals that the innermost jet expands in an edge-brightened cone structure (parabolic shape) but with the jet expansion profile slightly different from the outer regions of the jet. Brightness temperatures of the Very Long Baseline Interferometry (VLBI) core obtained from cm- to mm-wavelengths show a systematic evolution, which can be interpreted as the evolution as a function of distance from the BH. We also adopt an alternative imaging algorithm, Bi-Spectrum Maximum Entropy Method (BSMEM), to test reliable imaging at higher angular resolution than provided by the standard CLEAN method (i.e., super-resolution). A demonstration with a VLBA 7 mm example data set shows consistent results with a near-in-time 3 mm VLBI image. Application of the method to the 2009 GMVA data yields an image with remarkable fine-scale structures that have been never imaged before. We present a brief interpretation of the complexity in the structure.

Understanding the light curves of the HST-1 knot in M87 with internal relativistic shock waves along its jet

Knots or blobs observed in astrophysical jets are commonly interpreted as shock waves moving along them. Long time observations of the HST-1 knot inside the jet of the galaxy M87 have produced detailed multi-wavelength light curves. In this article, we model these light curves using the semi-analytical approach developed by Mendoza et al. (2009). This model was developed to account for the light curves of working surfaces moving along relativistic jets. These working surfaces are generated by periodic oscillations of the injected flow velocity and mass ejection rates at the base of the jet. Using genetic algorithms to fit the parameters of the model, we are able to explain the outbursts observed in the light curves of the HST-1 knot with an accuracy greater than a 2 sigma statistical confidence level.

RESOLVING THE ROTATION MEASURE OF THE M87 JET ON KILOPARSEC SCALES

We investigate the distribution of Faraday rotation measure (RM) in the M87 jet at arc-second scales by using archival polarimetric VLA data at 8, 15, 22 and 43 GHz. We resolve the structure of the RM in several knots along the jet for the first time. We derive the power spectrum in the arcsecond scale jet and find indications that the RM cannot be associated with a turbulent magnetic field with 3D Kolmogorov spectrum. Our analysis indicates that the RM probed on jet scales has a significant contribution of a Faraday screen associated with the vicinity of the jet, in contrast with that on kiloparsec scales, typically assumed to be disconnected from the jet. Comparison with previous RM analyses suggests that the magnetic fields giving rise to the RMs observed in jet scales have different properties and are well less turbulent than these observed in the lobes.

Detection of multiple velocity components in partially overlapping emitting regions

Velocity measurements made from multiple-epoch astronomical images of evolving objects with optically thin continuum emission (e.g. as relativistic jets or expanding supernova shells) may be confused as a result of the overlap of semi-transparent features moving at different speeds. Multi-scale wavelet decomposition can be effectively applied to identify and track such overlapping features, provided that their respective structural responses can be separated over the spatial scales used for the decomposition. We developed a new method that combines the stacked cross-correlation with the wavelet-based image segmentation and evaluation (WISE) technique of decomposition of two-dimensional structures, to separate and track dominant spatial responses of overlapping evolving features. The method is tested on a set of simulated images of a stratified relativistic jet, demonstrating the robust detection of both the faster spine and the slower sheath speeds. The method is applied to mutliple-epoch images from the MOJAVE survey, revealing two different superluminal streams inside the jet in 3C 273 and the acceleration of the flow in 3C 120. The method can be applied to densely monitored objects with composite structural evolution such as the parsec-scale jet in M87 or heavily resolved expanding supernova shells.

KEY SCIENCE OBSERVATIONS OF AGNs WITH THE KaVA ARRAY

KaVA (KVN and VERA Array) is a new combined VLBI array composed of KVN (Korean VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). Here, we report the following two issues. (1) We review the initial results of imaging observations of M87 at 23 GHz following Niinuma et al. (2014). The KaVA images reveal extended outflows including complex substructures such as knots and limb-brightening, in agreement with previous VLBI observations. KaVA achieves a high dynamic range of ~ 1000, more than three times better than that achieved by VERA alone. (2) Based on subsequent observations and discussions led by the KaVA AGN SubWorking Group, we set monitoring observations of Sgr A * and M87 as our Key Science Project (hereafter KSP) because of the closeness and largeness of their central super-massive black holes. The main science goals of the KSP are (i) testing the magnetically-driven-jet paradigm by mapping velocity elds of the M87 jet, and (ii) obtaining tight constraints on physical properties of the radio emitting region in Sgr A * . Towards KSP, we show the first preliminary images of M87 at 23 GHz and Sgr A * at 43 GHz with the bandwidth of 256 MHz.

HORIZON-SCALE LEPTON ACCELERATION IN JETS: EXPLAINING THE COMPACT RADIO EMISSION IN M87

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.

Variation of TeV emission from M87

In 2005, TeV flare from M87 was observed for the first time, which shows that the TeV flares are in a very short time scale, and detached. The viewing angle of the jet of M87 is very large, it is the first non-blazar like extragalactic TeV source, which might indicate a new production process for TeV emission. Possibly, the TeV source might be the hot plasma clouds from the accretion flow around black hole, these clouds might initially be in Keplerian orbits. In this work, the variation of the TeV emission is investigated in the framework of general relativity. In our model, both the effects of the Keplerian motion of the source and the different absorption optical depth of the TeV photons from the different locations of the source. The absorption of TeV photons is due to its collision with soft photons from the accretion disk. Our results show that variation of TeV emission is mainly dominated by their optical depth, the effect of the dynamical motion of the source could be neglected. When the TeV source is close to the central black hole, the optical depth significantly decreases with the increase of the black hole spin. If the black hole is non-rotating, the flux variation is obviously asymmetric. As the TeV source closes to black hole, their effects are enhanced.

Multi-wavelength polarimetry and variability study of the M87 jet

AbstractWe present a high resolution polarimetry and variability study of the M87 jet using VLA and HST data taken during 2002 to 2008. Both data-sets have an angular resolution as high as 0.06”, which is 2-3 times better than previous observations. New morphological details are revealed in both the optical and radio, which can help to reveal the energetic and magnetic field structure of the jet. By comparing the data with previously published HST and VLA observations, we show that the jet's morphology in total and polarized light is changing significantly on timescales of ~1 decade. We compare the evolution of the inner jet (particularly the nucleus and knot HST-1), when our observations overlap with the multi-wavelength monitoring campaigns conducted with HST and Chandra. We use these data to comment on particle acceleration and main emission processes.

Toy Model of Frame-Dragging Magnetosphere for the M87 Jet

We make a toy model for M87 jet to interpret its parabolic structure and acceleration in the apparent speeds, according to observations in milli-arcsecond to arcsecond scales upstream of HST-1. The outermost layer of jet is driven by the frame-dragging effect in the Kerr spacetime with a slowly to moderately spinning black hole. The corresponding magnetosphere has a foot-point R-0 in the vicinity of event-horizon, and rotating at a frequency Omega(F) equal to that of the frame-dragging omega(R-0).

A strong radio brightening at the jet base of M87 during the elevated very-high-energy γ-ray state in 2012

AbstractThe nearby radio galaxy M87 offers a unique opportunity for exploring the connection between γ-ray production and jet formation at an unprecedented linear resolution. However, the origin and location of the γ-rays in this source is still elusive. Based on previous radio/TeV correlation events, the unresolved jet base (radio core) and the peculiar knot HST-1 at >120 pc from the nucleus are proposed as candidate site(s) of γ-ray production. Here we report our intensive, high-resolution radio monitoring observations of the M87 jet with the VLBI Exploration of Radio Astrometry (VERA) and the European VLBI Network (EVN) from February 2011 to October 2012. During this period, an elevated level of the M87 flux is reported at TeV with VERITAS. We detected a remarkable flux increase in the radio core with VERA at 22/43 GHz coincident with the VHE activity. Meanwhile, HST-1 remained quiescent in terms of its flux density and structure at radio. These results strongly suggest that the TeV γ-ray activity in 2012 originates in the jet base within 0.03 pc (projected) from the central supermassive black hole.

A STRONG RADIO BRIGHTENING AT THE JET BASE OF M 87 DURING THE ELEVATED VERY HIGH ENERGY GAMMA-RAY STATE IN 2012

We report our intensive radio monitoring observations of the jet in M87 with the VLBI Exploration of Radio Astrometry (VERA) and the European VLBI Network (EVN) from February 2011 to October 2012, together with contemporaneous high-energy gamma-ray light curves obtained by the Fermi-LAT. During this period, an elevated level of the M87 flux is reported at VHE gamma rays. We detected a remarkable increase of the radio flux density from the unresolved jet base (radio core) with VERA at 22 and 43GHz coincident with the VHE activity. Meanwhile, we confirmed with EVN at 5GHz that HST-1 (an alternative gamma-ray production candidate site) remained quiescent in terms of its flux density and structure. These results in the radio bands strongly suggest that the VHE gamma-ray activity in 2012 originates in the jet base within 0.03pc or 56 Schwarzschild radii from the central supermassive black hole. We further conducted VERA astrometry for the M87 core during the flaring period, and detected core shifts between 22 and 43GHz. We also discovered a clear frequency-dependent evolution of the radio core flare at 43, 22 and 5GHz; the radio flux density increased more rapidly at higher frequencies with a larger amplitude, and the light curves clearly showed a time-lag between the peaks at 22 and 43GHz. This indicates that a new radio-emitting component was created near the black hole in the period of the VHE event, and then propagated outward with progressively decreasing synchrotron opacity. By combining these results, we estimated an apparent speed of the newborn component, and derived a sub-luminal speed of less than ~0.2c. This value is significantly slower than the super-luminal (~1.1c) features that appeared from the core during the prominent VHE flaring event in 2008, suggesting that the stronger VHE activity can be associated with the production of the higher Lorentz factor jet.

IMAGING THE SUPERMASSIVE BLACK HOLE SHADOW AND JET BASE OF M87 WITH THE EVENT HORIZON TELESCOPE

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.

RELATIVISTIC ELECTRONS AND MAGNETIC FIELDS OF THE M87 JET ON THE ∼10 SCHWARZSCHILD RADII SCALE

We explore energy densities of magnetic field and relativistic electrons in the M87 jet. Since the radio core at the jet base is identical to the optically thick surface against synchrotron self absorption (SSA), the observing frequency is identical to the SSA turnover frequency. As a first step, we assume the radio core as a simple uniform sphere geometry. Using the observed angular size of the radio core measured by the Very Long Baseline Array at 43 GHz, we estimate the energy densities of magnetic field ($U_{B}$) and relativistic electrons ($U_{e}$) based on the standard SSA formula. Imposing the condition that the Poynting power and relativistic electron one should be smaller than the total power of the jet, we find that (i) the allowed range of the magnetic field strength ($B_{tot}$) is from 1 G to 15 G, and that (ii) $1 times 10^{-5} < U_{e}/U_{B} < 6 times 10^{2}$ holds. The uncertainty of $U_{e}/U_{B}$ comes from the strong dependence on the angular size of the radio core and the minimum Lorentz factor of non-thermal electrons ($gamma_{e,min}$) in the core. It is still open that the resultant energetics is consistent with either the magnetohydrodynamic jet or with kinetic power dominated jet even on ~10 Schwarzschild radii scale.