No Arabic abstract
Controversial studies on the jet collimation profile of BL Lacertae (BL Lac), the eponymous blazar of BL Lac objects class, complicate the scenario in this already puzzling class of objects. Understanding the jet geometry, in connection with the jet kinematics and the physical conditions in the surrounding medium, is fundamental to better constrain the formation, acceleration and collimation mechanisms in extragalactic jets. With the aim of investigating the jet geometry in the innermost regions of the jet of BL Lac, and solving the controversy, we explore the radio jet in this source, using high resolution millimeter-wave VLBI data. We collect 86GHz GMVA and 43GHz VLBA data to obtain stacked images that we use to infer the jet collimation profile by means of two comparable methods. We analyze the kinematics at 86GHz, and we discuss it in the context of the jet expansion. Finally we consider a possible implication of the Bondi sphere in shaping the different expanding region observed along the jet. We found that the jet in BL Lac expands with an overall conical geometry. A higher expanding rate region is observed between ~5 and 10 pc (de-projected) from the black hole. Such a region is associated with the decrease in brightness usually observed in high-frequency VLBI images of BL Lac. The jet retrieves the original jet expansion around 17 pc, where the presence of a recollimation shock is supported by both the jet profile and the 15GHz kinematics (MOJAVE survey). The change in the jet expansion profile occurring at ~5 pc could be associated with a change in the external pressure profile in correspondence of the Bondi radius (~3.3X10$^5$$R_s$).
Parsec-scale VLBA images of BL Lac at 15 GHz show that the jet contains a permanent quasi-stationary emission feature 0.26 mas (0.34 pc projected) from the core, along with numerous moving features. In projection, the tracks of the moving features cluster around an axis at position angle -166.6 deg that connects the core with the standing feature. The moving features appear to emanate from the standing feature in a manner strikingly similar to the results of numerical 2-D relativistic magneto-hydrodynamic (RMHD) simulations in which moving shocks are generated at a recollimation shock. Because of this, and the close analogy to the jet feature HST-1 in M87, we identify the standing feature in BL Lac as a recollimation shock. We assume that the magnetic field dominates the dynamics in the jet, and that the field is predominantly toroidal. From this we suggest that the moving features are compressions established by slow and fast mode magneto-acoustic MHD waves. We illustrate the situation with a simple model in which the slowest moving feature is a slow-mode wave, and the fastest feature is a fast-mode wave. In the model the beam has Lorentz factor about 3.5 in the frame of the host galaxy, and the fast mode wave has Lorentz factor about 1.6 in the frame of the beam. This gives a maximum apparent speed for the moving features 10c. In this model the Lorentz factor of the pattern in the galaxy frame is approximately 3 times larger than that of the beam itself.
We study the kinematics of ridge lines on the pc-scale jet of the active galactic nucleus BL Lac. We show that the ridge lines display transverse patterns that move superluminally downstream, and that the moving patterns are analogous to waves on a whip. Their apparent speeds $beta_mathrm{app}$ (units of $c$) range from 3.9 to 13.5, corresponding to $beta_mathrm{wave}^mathrm{gal}= 0.981 - 0.998$ in the galaxy frame. We show that the magnetic field in the jet is well-ordered with a strong transverse component, and assume that it is helical and that the transverse patterns are Alfven waves propagating downstream on the longitudinal component of the magnetic field. The wave-induced transverse speed of the jet is non-relativistic ($beta_mathrm{tr}^mathrm{gal} lesssim 0.09$). In 2010 the wave activity subsided and the jet then displayed a mild wiggle that had a complex oscillatory behaviour. The Alfven waves appear to be excited by changes in the position angle of the recollimation shock, in analogy to exciting a wave on a whip by shaking the handle. A simple model of the system with plasma sound speed $beta_mathrm{s}=0.3$ and apparent speed of a slow MHD wave $beta_mathrm{app,S}=4$ yields Lorentz factor of the beam $Gamma_mathrm{beam} sim 4.5$, pitch angle of the helix (in the beam frame) $alphasim 67^circ$, Alfven speed $beta_mathrm{A}sim 0.64$, and magnetosonic Mach number $M_mathrm{ms}sim 4.7$. This describes a plasma in which the magnetic field is dominant and in a rather tight helix, and Alfven waves are responsible for the moving transverse patterns.
In an effort to locate the sites of emission at different frequencies and physical processes causing variability in blazar jets, we have obtained high time-resolution observations of BL Lacertae over a wide wavelength range: with the emph{Transiting Exoplanet Survey Satellite} (TESS) at 6,000-10,000 AA with 2-minute cadence; with the Neil Gehrels emph{Swift} satellite at optical, UV, and X-ray bands; with the Nuclear Spectroscopic Telescope Array at hard X-ray bands; with the emph{Fermi} Large Area Telescope at $gamma$-ray energies; and with the Whole Earth Blazar Telescope for measurement of the optical flux density and polarization. All light curves are correlated, with similar structure on timescales from hours to days. The shortest timescale of variability at optical frequencies observed with TESS is $sim 0.5$ hr. The most common timescale is $13pm1$~hr, comparable with the minimum timescale of X-ray variability, 14.5 hr. The multi-wavelength variability properties cannot be explained by a change solely in the Doppler factor of the emitting plasma. The polarization behavior implies that there are both ordered and turbulent components to the magnetic field in the jet. Correlation analysis indicates that the X-ray variations lag behind the $gamma$-ray and optical light curves by up to $sim 0.4$ days. The timescales of variability, cross-frequency lags, and polarization properties can be explained by turbulent plasma that is energized by a shock in the jet and subsequently loses energy to synchrotron and inverse Compton radiation in a magnetic field of strength $sim3$ G
One of the major challenges in studying the cosmic evolution of relativistic jets is the identification of the high-redshift ($z>3$) BL Lacertae objects, a class of jetted active galactic nuclei characterized by their quasi-featureless optical spectra. Here we report the identification of the first $gamma$-ray emitting BL Lac object, 4FGL~J1219.0+3653 (J1219), beyond $z=3$, i.e., within the first two billion years of the age of the Universe. The optical and near-infrared spectra of J1219 taken from 10.4 m Gran Telescopio Canarias exhibit no emission lines down to an equivalent width of $sim$3.5 A supporting its BL Lac nature. The detection of a strong Lyman-$alpha$ break at $sim$5570 A, on the other hand, confirms that J2119 is indeed a high-redshift ($zsim3.59$) quasar. Based on the prediction of a recent BL Lac evolution model, J1219 is one of the only two such objects expected to be present within the comoving volume at $z=3.5$. Future identifications of more $z>3$ $gamma$-ray emitting BL Lac sources, therefore, will be crucial to verify the theories of their cosmic evolution.
Since the launch of the Fermi satellite, BL Lacertae has been moderately active at gamma-rays and optical frequencies until May 2011, when the source started a series of strong flares. The exceptional optical sampling achieved by the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT) in collaboration with the Steward Observatory allows us to perform a detailed comparison with the daily gamma-ray observations by Fermi. Discrete correlation analysis between the optical and gamma-ray emission reveals correlation with a time lag of 0 +- 1 d, which suggests cospatiality of the corresponding jet emitting regions. A better definition of the time lag is hindered by the daily gaps in the sampling of the extremely fast flux variations. In general, optical flares present more structure and develop on longer time scales than corresponding gamma-ray flares. Observations at X-rays and at millimetre wavelengths reveal a common trend, which suggests that the region producing the mm and X-ray radiation is located downstream from the optical and gamma-ray-emitting zone in the jet. The mean optical degree of polarisation slightly decreases over the considered period and in general it is higher when the flux is lower. The optical electric vector polarisation angle (EVPA) shows a preferred orientation of about 15 deg, nearly aligned with the radio core EVPA and mean jet direction. Oscillations around it increase during the 2011-2012 outburst. We investigate the effects of a geometrical interpretation of the long-term flux variability on the polarisation. A helical magnetic field model predicts an evolution of the mean polarisation that is in reasonable agreement with the observations. These can be fully explained by introducing slight variations in the compression factor in a transverse shock waves model.