Kaikki aineistot
Lisää
Context. Supermassive black holes can launch highly relativistic jets with velocities reaching Lorentz factors of as high as Gamma > 50. How the jets accelerate to such high velocities and where along the jet they reach terminal velocity are open questions that are tightly linked to their structure as well as their launching and dissipation mechanisms. Aims. Changes in the beaming factor along the jets could potentially reveal jet acceleration, deceleration, or bending. We aim to (1) quantify the relativistic effects in multiple radio frequencies and (2) study possible jet velocity-viewing angle variations at parsec scales. Methods. We used the state-of-the-art code Magnetron to model light curves from the University of Michigan Radio Observatory and the Metsahovi Radio Observatory's monitoring programs in five frequencies covering about 25 years of observations in the 4.8 to 37 GHz range for 61 sources. We supplement our data set with high-frequency radio observations in the 100-340 GHz range from ALMA, CARMA, and SMA. For each frequency we estimate the Doppler factor which we use to quantify possible changes in the relativistic effects along the jets. Results. The majority of our sources do not show any statistically significant difference in their Doppler factor across frequencies. This is consistent with constant velocity in a conical jet structure, as expected at parsec scales. However, our analysis reveals 17 sources where relativistic beaming changes as a function of frequency. In the majority of cases, the Doppler factor increases towards lower frequencies. Only 1253-053 shows the opposite behavior. By exploring their jet properties we find that the jet of 0420-014 is likely bent across the 4.8-340 GHz range. For 0212+735, the jet is likely parabolic, and still accelerating in the 4.8-37 GHz range. We discuss possible interpretations for the trends found in the remaining sources.
Determining the location of γ-ray emission in blazar jets is a challenging task. Pinpointing the exact location of γ-ray production within a relativistic jet can place strong constraints on our understanding of high-energy astrophysics and astroparticle physics. We present a study of the radio- and γ-bright flat-spectrum radio quasar (FSRQ) 4C +01.28 (PKS B1055+018) in which we try to pinpoint the emission site of several prominent GeV flares. This source shows prominent high-amplitude broadband variability on time scales ranging from days to years. We combine high-resolution VLBI observations provided by the Boston University (BU) Blazar Monitoring Program and multi-band radio light curves over a period of around nine years. We can associate two bright and compact newly ejected jet components with bright flares observed by the Fermi/LAT γ-ray telescope and at various radio frequencies. A cross-correlation analysis reveals the radio light curves systematically lag behind the γ-rays. In combination with the jet kinematics as measured by the VLBA, we use these cross-correlations to constrain a model in which the flares become observable at a given frequency when a plasma component passes through the region at which the bulk energy dissipation takes place at that frequency. We derive a lower limit of the distance of the γ-ray emitting region in 4C +01.28 of several parsecs from the jet base, well beyond the expected extent of the broad-line region. This observational limit challenges blazar-emission models that rely on the broad-line region as a source of seed photons for inverse-Compton scattering.
Context. Cygnus X-3 is a unique microquasar in the Galaxy hosting a Wolf-Rayet companion orbiting a compact object that most likely is a low-mass black hole. The unique source properties are likely due to the interaction of the compact object with the heavy stellar wind of the companion. Aim. In this paper, we concentrate on a very specific period of time prior to the massive outbursts observed from the source. During this period, Cygnus X-3 is in a so-called hypersoft state, in which the radio and hard X-ray fluxes are found to be at their lowest values (or non-detected), the soft X-ray flux is at its highest values, and sporadic γ-ray emission is observed. We use multiwavelength observations to study the nature of the hypersoft state. Methods. We observed Cygnus X-3 during the hypersoft state with Swift and NuSTAR in X-rays and SMA, AMI-LA, and RATAN-600 in the radio. We also considered X-ray monitoring data from MAXI and γ-ray monitoring data from AGILE and Fermi. Results. We found that the spectra and timing properties of the multiwavelength observations can be explained by a scenario in which the jet production is turned off or highly diminished in the hypersoft state and the missing jet pressure allows the wind to refill the region close to the black hole. The results provide proof of actual jet quenching in soft states of X-ray binaries.
Total intensity variability light curves offer a unique insight into the ongoing debate about the launching mechanism of jets. For this work, we utilised the availability of radio and γ-ray light curves over a few decades of the radio source 3C 84 (NGC 1275). We calculated the multi-band time-lags between the flares identified in the light curves via discrete cross-correlation and Gaussian process regression. We find that the jet particle and magnetic field energy densities are in equipartition (kr-=-1.08-±-0.18). The jet apex is located z91.5-GHz-=-22-645-Rs (2- -20 ×10-3-pc) upstream of the 3 mm radio core; at that position, the magnetic field amplitude is Bcore91.5-GHz-=-3-10 G. Our results are in good agreement with earlier studies that utilised very-long-baseline interferometry. Furthermore, we investigated the temporal relation between the ejection of radio and γ-ray flares. Our results are in favour of the γ-ray emission being associated with the radio emission. We are able to tentatively connect the ejection of features identified at 43 and 86 GHz to prominent γ-ray flares. Finally, we computed the multiplicity parameter λ and the Michel magnetisation Ï M, and find that they are consistent with a jet launched by the Blandford & Znajek (1977, MNRAS, 179, 433) mechanism.
BL Lac, the eponymous blazar, flared to historically high levels at millimeter, infrared, X-ray, and gamma-ray wavelengths in 2012. We present observations made with Herschel, Swift, NuSTAR, Fermi, the Submillimeter Array, CARMA, and the VLBA in 2012-2013, including three months with nearly daily sampling at several wavebands. We have also conducted an intensive campaign of 30 hr with every-orbit observations by Swift and NuSTAR, accompanied by Herschel, and Fermi observations. The source was highly variable at all bands. Time lags, correlations between bands, and the changing shapes of the spectral energy distributions can be explained by synchrotron radiation and inverse Compton emission from nonthermal seed photons originating from within the jet. The passage of four new superluminal very long baseline interferometry knots through the core and two stationary knots about 4 pc downstream accompanied the high flaring in 2012-2013. The seed photons for inverse Compton scattering may arise from the stationary knots and from a Mach disk near the core where relatively slow-moving plasma generates intense nonthermal radiation. The 95 spectral energy distributions obtained on consecutive days form the most densely sampled, broad wavelength coverage for any blazar. The observed spectral energy distributions and multi-waveband light curves are similar to simulated spectral energy distributions and light curves generated with a model in which turbulent plasma crosses a conical shock with a Mach disk.
Abstract Aims: The very high energy (VHE ≳100 GeV) γ-ray MAGIC observations of the blazar S4 0954+65, were triggered by an exceptionally high flux state of emission in the optical. This blazar has a disputed redshift of z = 0.368 or z ≥ 0.45 and an uncertain classification among blazar subclasses. The exceptional source state described here makes for an excellent opportunity to understand physical processes in the jet of S4 0954+65 and thus contribute to its classification. Methods: We investigated the multiwavelength (MWL) light curve and spectral energy distribution (SED) of the S4 0954+65 blazar during an enhanced state in February 2015 and have put it in context with possible emission scenarios. We collected photometric data in radio, optical, X-ray, and γ-ray. We studied both the optical polarization and the inner parsec-scale jet behavior with 43 GHz data. Results: Observations with the MAGIC telescopes led to the first detection of S4 0954+65 at VHE. Simultaneous data with Fermi-LAT at high energy γ-ray (HE, 100 MeV < E < 100 GeV) also show a period of increased activity. Imaging at 43 GHz reveals the emergence of a new feature in the radio jet in coincidence with the VHE flare. Simultaneous monitoring of the optical polarization angle reveals a rotation of approximately 100°. Conclusions: The high emission state during the flare allows us to compile the simultaneous broadband SED and to characterize it in the scope of blazar jet emission models. The broadband spectrum can be modeled with an emission mechanism commonly invoked for flat spectrum radio quasars (FSRQs), that is, inverse Compton scattering on an external soft photon fieldfrom the dust torus, also known as external Compton. The light curve and SED phenomenology is consistent with an interpretation of a blob propagating through a helical structured magnetic field and eventually crossing a standing shock in the jet, a scenario typically applied to FSRQs and low-frequency peaked BL Lac objects (LBL).
3C 84 (NGC 1275) is one of the brightest radio sources in the millimetre radio bands, which led to a plethora of very-long-baseline interferometry (VLBI) observations at numerous frequencies over the years. They reveal a two-sided jet structure, with an expanding but not well-collimated parsec-scale jet, pointing southward. High-resolution millimetre-VLBI observations allow the study and imaging of the jet base on a sub-parsec scale. This could facilitate the investigation of the nature of the jet origin, also in view of the previously detected two-railed jet structure and east-west oriented core region seen with RadioAstron at 22 GHz. We produced VLBI images of this core and inner jet region, observed over the past twenty years at 15, 43, and 86 GHz. We determined the kinematics of the inner jet and ejected features at 43 and 86 GHz and compared their ejection times with radio and -ray variability. For the moving jet features, we find an average velocity of βavg app = 0:055 - 0:22c (μavg = 0:04 - 0:18 mas yr-1). From the time-averaged VLBI images at the three frequencies, we measured the transverse jet width along the bulk flow. On the ≤1:5 parsec scale, we find a clear trend of the jet width being frequency dependent, with the jet being narrower at higher frequencies. This stratification is discussed in the context of a spine-sheath scenario, and we compare it to other possible interpretations. From quasi-simultaneous observations at 43 and 86 GHz, we obtain spectral index maps, revealing a time-variable orientation of the spectral index gradient due to structural variability of the inner jet.
We report the discovery of a rare new form of long-term radio variability in the light curves of active galaxies (AG)-symmetric achromatic variability (SAV)-a pair of opposed and strongly skewed peaks in the radio flux density observed over a broad frequency range. We propose that SAV arises through gravitational milli-lensing when relativistically moving features in AG jets move through gravitational lensing caustics created by 103 106 M subhalo condensates or black holes located within intervening galaxies. The lower end of this mass range has been inaccessible with previous gravitational lensing techniques. This new interpretation of some AG variability can easily be tested and if it passes these tests, will enable a new and powerful probe of cosmological matter distribution on these intermediate-mass scales, as well as provide, for the first time, micro-arcsecond resolution of the nuclei of AG-a factor of 30100 greater resolution than is possible with ground-based millimeter very-longbaseline interferometry.
Context. The γ-ray BL Lac object OJ 287 is known to exhibit inner-parsec "jet-wobbling", high degrees of variability at all wavelengths and quasi-stationary features, including an apparent (≈100°) position-angle change in projection on the sky plane. Aims: Sub-50 micro-arcsecond resolution 86 GHz observations with the global mm-VLBI array (GMVA) supplement ongoing multi-frequency VLBI blazar monitoring at lower frequencies. Using these maps, together with cm/mm total intensity and γ-ray observations from Fermi-LAT from 2008-2014, we aim to determine the location of γ-ray emission and to explain the inner-mas structural changes. Methods: Observations with the GMVA offer approximately double the angular resolution compared with 43 GHz VLBA observations and enable us to observe above the synchrotron self-absorption peak frequency. Fermi-LAT γ-ray data were reduced and analysed. The jet was spectrally decomposed at multiple locations along the jet. From this, we could derive estimates of the magnetic field using equipartition and synchrotron self-absorption arguments. How the field decreases down the jet provided an estimate of the distance to the jet apex and an estimate of the magnetic field strength at the jet apex and in the broad line region. Combined with accurate kinematics, we attempt to locate the site of γ-ray activity, radio flares, and spectral changes. Results: Strong γ-ray flares appeared to originate from either the so-called core region, a downstream stationary feature, or both, with γ-ray activity significantly correlated with radio flaring in the downstream quasi-stationary feature. Magnetic field estimates were determined at multiple locations along the jet, with the magnetic field found to be ≥1.6 G in the core and ≤0.4 G in the downstream quasi-stationary feature. We therefore found upper limits on the location of the VLBI core as ≲6.0 pc from the jet apex and determined an upper limit on the magnetic field near the jet base of the order of thousands of Gauss. The 3 mm GMVA data are only available at the CDS via anonymous ftp to http://cdsarc.u-strasbg.fr (http://130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/597/A80
The radio light curve of J1415+1320 (PKS 1413+135) shows time-symmetric and recurring U-shaped features across the centimeter-wave and millimeter-wave bands. The symmetry of these features points to lensing by an intervening object as the cause. U-shaped events in radio light curves in the centimeter-wave band have previously been attributed to Extreme scattering events (ESE). ESEs are thought to be the result of lensing by compact plasma structures in the Galactic interstellar medium, but the precise nature of these plasma structures remains unknown. Since the strength of a plasma lens evolves with wavelength λ as , the presence of correlated variations at over a wide wavelength range casts doubt on the canonical ESE interpretation for J1415+1320. In this paper, we critically examine the evidence for plasma lensing in J1415+1320. We compute limits on the lensing strength and the associated free-free opacity of the putative plasma lenses. We compare the observed and model ESE light curves, and also derive a lower limit on the lens distance based on the effects of parallax due to the Earth's orbit around the Sun. We conclude that plasma lensing is not a viable interpretation for J1415+1320's light curves and that symmetric U-shaped features in the radio light curves of extragalactic sources do not present prima facie evidence for ESEs. The methodology presented here is generic enough to be applicable to any plasma-lensing candidate.
Symmetric achromatic variability (SAV) is a rare form of radio variability in blazars that has been attributed to gravitational milli-lensing by a ∼102-105 M ⊙ mass condensate. Four SAVs have been identified between 1980 and 2020 in the long-term radio monitoring data of the blazar PKS 1413 + 135. We show that all four can be fitted with the same, unchanging, gravitational lens model. If SAV is due to gravitational milli-lensing, PKS 1413 + 135 provides a unique system for studying active galactic nuclei with unprecedented microarcsecond resolution, as well as for studying the nature of the milli-lens itself. We discuss two possible candidates for the putative milli-lens: a giant molecular cloud hosted in the intervening edge-on spiral galaxy, and an undetected dwarf galaxy with a massive black hole. We find a significant dependence of SAV crossing time on frequency, which could indicate a fast shock moving in a slower underlying flow. We also find tentative evidence for a 989 day periodicity in the SAVs, which, if real, makes possible the prediction of future SAVs: the next three windows for possible SAVs begin in 2022 August, 2025 May, and 2028 February.
We explore the connection between the γ-ray and radio emission in the jet of the blazar 0716+714 by using 15, 37, and 230 GHz radio and 0.1-200 GeV γ-ray light curves spanning 10.5 yr (2008-2019). We find significant positive and negative correlations between radio and γ-ray fluxes in different time ranges. The time delays between radio and γ-ray emission suggest that the observed γ-ray flares originated from multiple regions upstream of the radio core, within a few parsecs from the central engine. Using time-resolved 43 GHz Very Long Baseline Array maps we identified 14 jet components moving downstream along the jet. Their apparent speeds range from 6c to 26c, and they show notable variations in their position angles upstream from the stationary component (∼0.53 mas from the core). The brightness temperature declines as a function of distance from the core according to a power law that becomes shallower at the location of the stationary component. We also find that the periods at which significant correlations between radio and γ-ray emission occur overlap with the times when the jet was oriented to the north. Our results indicate that the passage of a propagating disturbance (or shock) through the radio core and the orientation of the jet might be responsible for the observed correlation between the radio and γ-ray variability. We present a scenario that connects the positive correlation and the unusual anticorrelation by combining the production of a flare and a dip at γ-rays by a strong moving shock at different distances from the jet apex.
Context. We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1, which includes, among other instruments, MAGIC, VERITAS, Whipple 10 m, and Fermi-LAT to cover the γ-ray range from 0.1 GeV to 20 TeV; RXTE and Swift to cover wavelengths from UV tohard X-rays; and GASP-WEBT, which provides coverage of radio and optical wavelengths. Optical polarization measurements were provided for a fraction of the campaign by the Steward and St. Petersburg observatories. We evaluate the variability of the source and interband correlations, the γ-ray flaring activity occurring in May 2009, and interpret the results within two synchrotron self-Compton (SSC) scenarios. Aims. The multiband variability observed during the full campaign is addressed in terms of the fractional variability, and the possible correlations are studied by calculating the discrete correlation function for each pair of energy bands where the significance was evaluated with dedicated Monte Carlo simulations. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters, hence providing a less biased interpretation than the "single-curve SSC model adjustment" typically reported in the literature. Methods. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found on the basis of the acquired data set. The SSC model grid-scan shows that the flaring activity around May 22 cannot be modeled adequately with a one-zone SSC scenario (using an electron energy distribution with two breaks), while it can be suitably described within a two (independent) zone SSC scenario. Here, one zone is responsible for the quiescent emission from the averaged 4.5-month observing period, while the other one, which is spatially separated from the first, dominates the flaring emission occurring at X-rays and very-high-energy (>100 GeV, VHE) γ rays. The flaring activity from May 1, which coincides with a rotation of the electric vector polarization angle (EVPA), cannot be satisfactorily reproduced by either a one-zone or a two-independent-zone SSC model, yet this is partially affected by the lack of strictly simultaneous observations and the presence of large flux changes on sub-hour timescales (detected at VHE γ rays). Results. The higher variability in the VHE emission and lack of correlation with the X-ray emission indicate that, at least during the 4.5-month observing campaign in 2009, the highest energy (and most variable) electrons that are responsible for the VHE γ rays do not make a dominant contribution to the ~1 keV emission. Alternatively, there could be a very variable component contributing to the VHE γ-ray emission in addition to that coming from the SSC scenario. The studies with our dedicated SSC grid-scan show that there is some degeneracy in both the one-zone and the two-zone SSC scenarios probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The observed γ-ray flaring activity, with the EVPA rotation coincident with the first γ-ray flare, resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties.
We report the results of decade-long (2008-2018) γ -ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, Fermi and Swift data, as well as polarimetric and spectroscopic data. The X-ray and γ -ray light curves correlate well, with no delay ≳ 3 h, implying general cospatiality of the emission regions. The γ -ray-optical flux-flux relation changes with activity state, ranging from a linear to amore complex dependence. The behaviour of the Stokes parameters at optical and radio wavelengths, including 43 GHz Very Long Baseline Array images, supports either a predominantly helical magnetic field or motion of the radiating plasma along a spiral path. Apparent speeds of emission knots range from 10 to 37c, with the highest values requiring bulk Lorentz factors close to those needed to explain γ -ray variability on very short time-scales. The MgII emission line flux in the 'blue' and 'red' wings correlates with the optical synchrotron continuum flux density, possibly providing a variable source of seed photons for inverse Compton scattering. In the radio bands, we find progressive delays of the most prominent light-curve maxima with decreasing frequency, as expected from the frequency dependence of the τ= 1 surface of synchrotron self-absorption. The global maximum in the 86 GHz light curve becomes less prominent at lower frequencies, while a local maximum, appearing in 2014, strengthens toward decreasing frequencies, becoming pronounced at ∼5 GHz. These tendencies suggest different Doppler boosting of stratified radio-emitting zones in the jet.
We present multi-wavelength studies of the radio galaxy 3C 120 and the blazar CTA 102 during unprecedented γ-ray flares for both sources. In both studies the analysis of γ-ray data has been compared with a series of 43 GHz VLBA images from the VLBA-BU-BLAZAR program, providing the necessary spatial resolution to probe the parsec scale jet evolution during the high energy events. To extend the radio dataset for 3C 120 we also used 15 GHz VLBA data from the MOJAVE sample. These two objects which represent very different classes of AGN, have similar properties during the γ-ray events. The γ-ray flares are associated with the passage of a new superluminal component through the mm VLBI core, but not all ejections of new components lead to γ-ray events. In both sources γ-ray events occurred only when the new components are moving in a direction closer to our line of sight. We locate the γ-ray dissipation zone a short distance from the radio core but outside of the broad line region, suggesting synchrotron self-Compton scattering as the probable mechanism for the γ-ray production.
The Event Horizon Telescope (EHT) observed the compact radio source, Sagittarius A* (Sgr A*), in the Galactic Center on 2017 April 5-11 in the 1.3 mm wavelength band. At the same time, interferometric array data from the Atacama Large Millimeter/submillimeter Array and the Submillimeter Array were collected, providing Sgr A* light curves simultaneous with the EHT observations. These data sets, complementing the EHT very long baseline interferometry, are characterized by a cadence and signal-to-noise ratio previously unattainable for Sgr A* at millimeter wavelengths, and they allow for the investigation of source variability on timescales as short as a minute. While most of the light curves correspond to a low variability state of Sgr A*, the April 11 observations follow an X-ray flare and exhibit strongly enhanced variability. All of the light curves are consistent with a red-noise process, with a power spectral density (PSD) slope measured to be between -2 and -3 on timescales between 1 minute and several hours. Our results indicate a steepening of the PSD slope for timescales shorter than 0.3 hr. The spectral energy distribution is flat at 220 GHz, and there are no time lags between the 213 and 229 GHz frequency bands, suggesting low optical depth for the event horizon scale source. We characterize Sgr A*'s variability, highlighting the different behavior observed just after the X-ray flare, and use Gaussian process modeling to extract a decorrelation timescale and a PSD slope. We also investigate the systematic calibration uncertainties by analyzing data from independent data reduction pipelines.
Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov and conservation of the Walker-Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87∗, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A∗ with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87∗. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images.
The black hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6 day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure-phase measurements on all six linearly independent nontrivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of ∼3°-5°. The only triangles that exhibit substantially higher variability (∼90°-180°) are the ones with baselines that cross the visibility amplitude minima on the u-v plane, as expected from theoretical modeling. We used two sets of general relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas.
The Event Horizon Telescope (EHT) provides the unprecedented ability to directly resolve the structure and dynamics of black hole emission regions on scales smaller than their horizons. This has the potential to critically probe the mechanisms by which black holes accrete and launch outflows, and the structure of supermassive black hole spacetimes. However, accessing this information is a formidable analysis challenge for two reasons. First, the EHT natively produces a variety of data types that encode information about the image structure in nontrivial ways; these are subject to a variety of systematic effects associated with very long baseline interferometry and are supplemented by a wide variety of auxiliary data on the primary EHT targets from decades of other observations. Second, models of the emission regions and their interaction with the black hole are complex, highly uncertain, and computationally expensive to construct. As a result, the scientific utilization of EHT observations requires a flexible, extensible, and powerful analysis framework. We present such a framework, Themis, which defines a set of interfaces between models, data, and sampling algorithms that facilitates future development. We describe the design and currently existing components of Themis, how Themis has been validated thus far, and present additional analyses made possible by Themis that illustrate its capabilities. Importantly, we demonstrate that Themis is able to reproduce prior EHT analyses, extend these, and do so in a computationally efficient manner that can efficiently exploit modern high-performance computing facilities. Themis has already been used extensively in the scientific analysis and interpretation of the first EHT observations of M87.
3C 279 is an archetypal blazar with a prominent radio jet that show broadband flux density variability across the entire electromagnetic spectrum. We use an ultra-high angular resolution technique - global Very Long Baseline Interferometry (VLBI) at 1.3mm (230 GHz) - to resolve the innermost jet of 3C 279 in order to study its fine-scale morphology close to the jet base where highly variable-ray emission is thought to originate, according to various models. The source was observed during four days in April 2017 with the Event Horizon Telescope at 230 GHz, including the phased Atacama Large Millimeter/submillimeter Array, at an angular resolution of ∼20 μas (at a redshift of z = 0:536 this corresponds to ∼0:13 pc ∼ 1700 Schwarzschild radii with a black hole mass MBH = 8 × 108 M⊙). Imaging and model-fitting techniques were applied to the data to parameterize the fine-scale source structure and its variation.We find a multicomponent inner jet morphology with the northernmost component elongated perpendicular to the direction of the jet, as imaged at longer wavelengths. The elongated nuclear structure is consistent on all four observing days and across diffierent imaging methods and model-fitting techniques, and therefore appears robust. Owing to its compactness and brightness, we associate the northern nuclear structure as the VLBI "core". This morphology can be interpreted as either a broad resolved jet base or a spatially bent jet.We also find significant day-to-day variations in the closure phases, which appear most pronounced on the triangles with the longest baselines. Our analysis shows that this variation is related to a systematic change of the source structure. Two inner jet components move non-radially at apparent speeds of ∼15 c and ∼20 c (∼1:3 and ∼1:7 μas day-1, respectively), which more strongly supports the scenario of traveling shocks or instabilities in a bent, possibly rotating jet. The observed apparent speeds are also coincident with the 3C 279 large-scale jet kinematics observed at longer (cm) wavelengths, suggesting no significant jet acceleration between the 1.3mm core and the outer jet. The intrinsic brightness temperature of the jet components are ≤1010 K, a magnitude or more lower than typical values seen at ≥7mm wavelengths. The low brightness temperature and morphological complexity suggest that the core region of 3C 279 becomes optically thin at short (mm) wavelengths.
The Event Horizon Telescope (EHT) is a millimeter very long baseline interferometry (VLBI) array that has imaged the apparent shadows of the supermassive black holes M87* and Sagittarius A*. Polarimetric data from these observations contain a wealth of information on the black hole and accretion flow properties. In this work, we develop polarimetric geometric modeling methods for mm-VLBI data, focusing on approaches that fit data products with differing degrees of invariance to broad classes of calibration errors. We establish a fitting procedure using a polarimetric “m-ring” model to approximate the image structure near a black hole. By fitting this model to synthetic EHT data from general relativistic magnetohydrodynamic models, we show that the linear and circular polarization structure can be successfully approximated with relatively few model parameters. We then fit this model to EHT observations of M87* taken in 2017. In total intensity and linear polarization, the m-ring fits are consistent with previous results from imaging methods. In circular polarization, the m-ring fits indicate the presence of event-horizon-scale circular polarization structure, with a persistent dipolar asymmetry and orientation across several days. The same structure was recovered independently of observing band, used data products, and model assumptions. Despite this broad agreement, imaging methods do not produce similarly consistent results. Our circular polarization results, which imposed additional assumptions on the source structure, should thus be interpreted with some caution. Polarimetric geometric modeling provides a useful and powerful method to constrain the properties of horizon-scale polarized emission, particularly for sparse arrays like the EHT.
In 2017 the Event Horizon Telescope (EHT) observed the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz (λ = 1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT data sets. The high observing frequency means that pulsars—which typically exhibit steep emission spectra—are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope, and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the fast folding algorithm, and single-pulse searches targeting both pulsars and burst-like transient emission. We use the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction (≲2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region.