VLBI Analysis of a Potential High-Energy Neutrino Emitter Blazar

Recent studies suggest that high-energy neutrinos can be produced in the jets of blazars, radio-loud active galactic nuclei (AGN) with jets pointing close to the line of sight. Due to the relatively poor angular resolution of current neutrino detectors, several sources can be regarded as the possible counterpart of a given neutrino event. Therefore, follow-up observations of counterpart candidates in the electromagnetic regime are essential. Since the Very Long Baseline Interferometry (VLBI) technique provides the highest angular resolution to study the radio jets of blazars, a growing number of investigations are being conducted to connect individual blazars to given high-energy neutrino events. We analyzed more than 20 years of available archival VLBI data of the blazar CTD 74, which has been listed as a possible counterpart of a neutrino event. Using cm-wavelength data, we investigated the jet structure, determined the apparent speed of jet components, and the core flux density before and after the neutrino event. Our results indicate stationary jet features and a significant brightening of the core after the neutrino event.


Introduction
Blazars belong to the group of radio-loud active galactic nuclei (AGN) with their jets pointing in a small angle to the line of sight [1].In recent studies [2], blazars were associated with high-energy neutrino events detected by the IceCube Neutrino Observatory located at the geographical South Pole [3].While as of now, except for the first blazar associated to a neutrino event, TXS 0506+056 [4], no definitive link has been established between individual high-energy IceCube events and radio-loud AGN, there have been several assertions of associations with AGN at various statistical levels [5][6][7][8][9].However, other authors found that the contribution of the blazar class as a whole to the IceCube neutrino signal is not dominant [10,11].Other studies using very similar blazar catalogs found no meaningful correlation between IceCube neutrino events and radio-bright AGN [12,13].In addition to the IceCube, similar studies were conducted with the ANTARES neutrino telescope [14], which hinted a possible connection between the ANTARES-detected neutrino candidates with blazars.Apart from blazars, other candidates for the source of high-energy neutrinos were proposed, such as gamma-ray bursts, supernova remnants, starburst galaxies [e.g., 15], tidal disruption events [16], white dwarf mergers [17], and other non-blazar AGN, e.g., NGC 1068 [18].Thus, since the first association, analysing possible neutrino-source blazars has become an important task [19][20][21] as these objects may help us understand the astrophysical processes that lead to the production of high-energy neutrinos.
In this paper, we present high-resolution Very Long Baseline Interferometry (VLBI) analysis of the blazar CTD 74 (alternatively: TXS 1123+264, J1125+2610) that has been statistically linked to an extremely high-energy track-like neutrino alert (EHEA2012-05-23) observed in 2012, due to its positional proximity.CTD 74 was in the sample of radio-loud AGN having VLBI observations at 8 GHz, and monitored at 22 GHz by the Russian RATAN-600 radio telescope [22] that were analyzed in connection with high-energy (> 200 TeV) track-like neutrino events detected by the IceCube instrument [6].The selection criteria of Extremely High Energy Alerts and alert-like events are described in [23].The directions of such events can be determined with an uncertainty of less than 1 • [24], and CTD 74 fell into the error region of such event.The positional association of the high-energy event combined with the available 25-yr long multi-frequency VLBI observations, which covered the time of the neutrino event, made this object a promising candidate to study individually in the context of possible neutrino-emitter blazars.However, we note that the consideration of this particular blazar as a high-energy neutrino source was based on positional coincidence in previous works.We emphasize the possibility that the detected neutrino event might originate from a different source or might not even be astrophysical.Individual AGN sources have been studied in the past in a similar manner [18,[25][26][27].However, we are not aware of any study concerning the high-resolution radio emission of CTD 74.
CTD 74 [28] is a blazar at a redshift of z = 2.3502 ± 0.0001 [29].Its coordinates in the International Celestial Reference Frame (ICRF) are 11 h 25 m 53.7119 s right ascension and +26 • 10 ′ 19.979 ′′ declination [30].The source has been observed in several epochs with VLBI, but to our knowledge, a comprehensive analysis of the data has not been published yet.It is not included in the latest Fermi Large Area Telescope catalog as a γ-ray source [31] and no information was found about its X-ray detection.
Assuming a flat Λ Cold Dark Matter cosmological model with Hubble constant H 0 = 70 km s −1 Mpc −1 , matter density parameter Ω m = 0.27, and vacuum energy density parameter Ω vac = 0.73, the angular diameter distance of the source is D A = 1734.2Mpc, and the angular scale is 8.407 pc mas −1 [32].

Observations and data reduction
We analyzed multi-frequency, multi-epoch pre-calibrated VLBI visibility data, obtained from the Astrogeo (accessed on 28 September 2023) (http://astrogeo.org,see also https: //astrogeo.smce.nasa.gov/vlbi_images/)website maintained by L. Petrov.Most of the observations of CTD 74 were conducted at the 8-GHz frequency band (central frequencies 7.62 − 8.67 GHz) in 23 epochs between 1996 and 2021.The majority of these were astrometric observations conducted with a dual-band receiver, therefore in 19 epochs, 2-GHz-band measurements are also available.There were additional 3 observations done at the 5-GHz band (Table 1).Since the central frequencies of the observations were slightly different from epoch to epoch, we refer to the three main bands of the observing frequencies as 8, 5, and 2 GHz in the text, while for calculations, we naturally use the central frequencies of the given data sets.
These observations were performed by heterogeneous arrays of radio telescopes.Primarily the Very Long Baseline Array (VLBA) of the U.S. National Radio Astronomy Observatory (NRAO) was used.It consists of 10 antennas at Brewster (BR), Fort Davis (FD), Hancock (HN), Kitt Peak (KP), Los Alamos (LA), Mauna Kea (MK), North Liberty (NL), Owens Valley (OV), Pie Town (PT), and St. Croix (SC).But occasionally, multiple antennas from the European VLBI Network (EVN), Medicina (MC, Italy), Sheshan (SH, China), Onsala (ON, Sweden), Hartebeesthoek 26-m (HH, South Africa), as well as other geodetic VLBI radio telescopes, Goddard (GG, USA), Gilmore Creek (GC, USA), Hartebeesthoek 15m (HT, South Africa), Kokee Park (KK, USA), Matera (MA, Italy), Wettzell (WZ, Germany), Tsukuba (TS, Japan), Ny-Ålesund (NY, Norway), Westford (WF, USA), and Yarragadee (YG, Australia) participated in various experiments.The details of the observations are summarized in Table 1.It is indicated with the corresponding two-letter station codes preceded by a minus sign if one or two antennas from the VLBA were missing from the array.It is important to note that the data set is not homogeneous, with vastly different on-source integration times (as seen in Table 1) resulting in different sampling of the visibility functions that can affect the quality of the obtained images.In Figure 1, we show the distributions of baseline vectors, commonly referred as (u, v) coverage, for two observations, one with the shortest integration time (14 s, performed on 2017 March 27), and one with a significantly longer integration time (12 892 s, performed on 2020 November 18).
Table 1.Details of the analysed VLBI observations.The observational epochs, central frequencies in GHz (ν), the codes of the participating antennas (Stations), the on-source integration times of the measurements, the number of intermediate frequency channels (IFs), their bandwidth in MHz (BW), and the project codes and references (where applicable) are given in the columns.We used the Difmap software (version number 2.5k) [41] to produce images of CTD 74 using the hybrid mapping method [e.g.42,43].First, we created the dirty map via Fouriertransformation of the visibilities, after setting the appropriate map and pixel size.We defined a small region around the brightest pixel value in the dirty map, and using the the CLEAN algorithm [44], we subtracted point-source responses with 5% of the brightest pixel value from the image in 50 or 100 iterations.This first CLEAN component model was then used to calibrate the phases of the visibilities in the phase self-calibration process [45].Then a new, residual map, from which the previous CLEAN component model had already been subtracted, was checked for the value and position of the brightest pixel, and the above described procedure was repeated.In order not to include spurious features in the point source model, i.e. to avoid cleaning the artifacts arising in the residual maps from the nonideal visibility sampling, the brightest pixel search was always guided by the definitions of small regions for the cleaning process.That way the point source model was iteratively refined in small steps, in parallel with the improvement of the phase self-calibration.

Epoch
When the model could not be improved further as judged from the noise level in the residual map, typically when the brightest pixel in the residual image was ≲ 5 times of the noise level, we performed iterations of amplitude and phase self-calibration, gradually decreasing the solution interval by a factor of two in each iteration starting with a few hours down to 1 − 2 minutes.We note here that because the noise in the residual image does not follow Gaussian statistics, the usage of root mean square value to estimate the total uncertainty is not strictly correct [46,47].The maps shown in Figure 2 are chosen as examples to represent the final images that were obtained for all epochs and frequencies.
In order to quantify the source brightness distribution, we fitted the self-calibrated visibility data with two-dimensional Gaussian model components [48].The parameters of these components can be found in Tables A1, A2, and A3 for the 8-, 2-, and 5-GHz data, respectively.To minimize the number of free parameters, we initially attempted to fit circular Gaussian components only.All four parameters of the components, positions, flux density (S [Jy]), and full width at half maximum (FWHM) size (W 1 [mas]), were set to freely variable.In seven cases, we could not obtain a stable fit with using only circular Gaussian components, therefore we included elliptical components also.For those, the additional two parameters, the minor axis size (W 2 [mas]) and the position angle of the major axis (Φ [ • ]) were free parameters as well.In most of these cases (5), the elliptical components were needed to describe the brightness distribution of the core, while in two other cases, a jet feature was fitted by elliptical Gaussians, as seen in Tables A1, A2, and A3.

Error calculation of model parameters
When calculating the errors of the model parameters, we used the formulae given in [49,50].These are based on the prescriptions of [51], however, they take into account the image artefacts arising from the sparse (u, v) coverage of VLBI observations via the inclusion of the restoring beam.
We determined the post-fit root mean square error of the image, σ S p within a rectangular region around the position of the fitted components in the residual map.The error of the component diameter (FWHM) is given as: where S p is the peak intensity of the components with an error of σ S p .θ denotes the size of the restoring beam as θ = θ 2 min + θ 2 maj , where θ min is the minor and θ maj is the major axis of the restoring beam.The index i can take a value of 1 and 2 for the two axes of the fitted component.
If the size of a component is smaller than the smallest resolvable size with the interferometer array of the given observation, only an upper limit of the source size, θ lim , can be given.To calculate θ lim , we used the equation in [52]: where SNR stands for the signal-to-noise ratio, calculated from the ratio of peak intensity to the noise level.
The error of the flux density (σ S ) includes the squared sum of the error given by the analytical formula in [51], and the contribution arising from the amplitude calibration error: where The amplitude calibration error is commonly taken as 5 % of the flux density in case of VLBA observations [53].The source brightness distribution could be fitted with a core X0 and one or more additional jet components.For most of the 23 epochs of 8-GHz measurements, we could consistently identify two jet components, J1 and J2 (in 20 epochs).In 7 epochs, due to the denser sampling of the visibility function and the better sensitivity, we could fit an additional component (J1.5) located in between J1 and J2.An image from one of these measurements, the one performed on 2014 Feb 12, is shown in Figure 3 with the model components labeled.Finally, in the earliest two epochs, we could not securely identify the fitted jet components with those after 2000.These components are named as X1 and X2 (Table A1).

Analysis
At 2 GHz, 11 out of the 19 epochs could be described with three Gaussian components.In the other cases, two Gaussian components were needed to adequately fit the visibility data (Table A2).The parameters of the fitted Gaussian model components at 5 GHz in 3 epochs are given in Table A3.

Results
In all the analysed data, CTD 74 consistently showed a core-jet structure, with a short jet oriented to the northwest.Even at 2 GHz, the jet emission cannot be traced farther than ∼ 10 mas from the core.This can be caused by the non-ideal sampling of the visibilities at the shortest baselines causing the observations to be less sensitive to the largest spatial scales of the emitting regions.To study the apparent motion of the jet components, we calculated their separation from the core at each epoch.Since 8 GHz is the highest frequency and thus these observations provide the best angular resolution, the 8-GHz data are the most suitable to detect structural changes and jet component motions in the compact environment around the core.The core separations of these jet components are shown in Figure 4. We fitted linear functions to the core-jet component separations versus time, to determine the apparent proper motion of the jet components.The slopes of the lines fitted for the components J1 and J2 are (−0.002± 0.006) mas yr −1 and (−0.009 ± 0.008) mas yr −1 , respectively.These values are consistent with no proper motion within the uncertainties (Figure 4).Thus, no ballistic movement could be detected at 8 GHz in either of the jet components.We reached the same conclusion for the 2-GHz data, after fitting linear functions to the core-jet components separations.

Jet morphology in various epochs
When studying core separation, we assume that the components are moving linearly, neglecting their two-dimensional motion in the plane of the sky.To investigate this motion, we examined how the position angles of the jet components change over time.In Figure 5, we show the position angles of the jet components that are the most distant from the core at the given frequency, and having a large number of detections, i.e., J2 at 8 GHz and S2 at 2 GHz.The average position angles are (−50 ± 8) • and (−69 ± 6) • for J2 and S2, respectively.In case of the components closer to the core, the average position angles are (−52 ± 8) • and (−71 ± 7) • for J1 and S1, respectively.Thus, the northwestern direction of the jet remains constant through the epochs, but there is a slight bend of ∼ 20 • between the inner ≲ 3 mas and outer regions of the jet.This could be traced with the observations taken at different frequencies.
In summary, the jet components appear stationary in the plane of the sky, not showing significant changes in the their core separation and position angle.

Flux density
We examined how the flux density of the core component changes over time.The average flux density value of the core is 0.62 Jy at 8 GHz, but it shows significant variability during the ∼ 15 yr of analysed VLBI observations.Most notably, it shows a brightening around 2014-2015 (see Figure 6).During its brightening in 2014-2015, the core was 38% brighter than its average value, reaching (0.86 ± 0.04) Jy on 2014 May 31.Then the core faded back quickly after 2017.As there was no 8-GHz VLBI observation taken of the object between 2003 and 2010, the core flux density changes cannot be discussed in this time range.
Since the flux density derived from the Gaussian components fitted to the visibilities may in principle be influenced by the heterogeneous resolution of the different observations, we also compared the peak intensity value of the images created with the same restoring beams, 0.98 mas × 2.47 mas with a position angle of −22.51 • .We chose the observation with the shortest on-source time (2017 March 27).The same trend as in the fitted core component flux densities can be seen in the time variability of the intensity values as well (Figure 6).
We also examined the change in core flux density at 2 GHz.Here, with less frequent time sampling, we also observe the brightening and fading that was seen at 8 GHz.Knowing the redshift of the object, from the parameters of the fitted Gaussian components, their brightness temperatures, T b,VLBI , can be calculated [55].For the core component:

Brightness temperature
where S core is the flux density expressed in Jy, ν the observing frequency in GHz.The component major and minor axes (FWHM), W 1 and W 2 , are given in mas.
Assuming that the intrinsic brightness temperature of the source is equal to the equipartition value, T int = T eq ≈ 5 • 10 10 K [54], the Doppler boosting factor, which quantifies the relativistic beaming effect, can be calculated as follows: As seen in Figure 7, every measured brightness temperature value is 1 − 2 orders of magnitude higher than the red horizontal line denoting the assumed intrinsic brightness temperature, T eq .Under the assumption of the intrinsic brightness temperature being close to the equipartition value, we can interpret these data as the result of Doppler boosting present at every epoch.The average brightness temperature is 1.87 • 10 12 K, and the average Doppler boosting factor is 37.The highest T b,VLBI value is actually a lower limit, since the fitted core size was smaller than the smallest resolvable angular size of the array in that particular epoch.Thus only an upper limit could be given for the size of the core component.

Spectral index
We calculated the spectral index of the core, α, using the convention S core ∝ ν α .For the estimation, we selected the 1996 measurements, since in this year, all measurements at 3 different frequencies were conducted close in time, within a few weeks.
The fitted power-law function has an exponent of α = −0.29 ± 0.13.Alternatively, we can also calculate a two-point spectral index, α CX , using the simultaneous 8.3 and 4.9 GHz measurements on 2016 July 17.The value obtained is α CX = 0.0 ± 0.1.
Both of these spectral indices indicate a flat radio spectrum, which is in line with the general expectations that the cores of blazars have either flat (α ≈ 0) or inverted (α ≳ 0.5) spectra, while the jet typically has a steep (α ≲ −0.7) spectrum [56].

Discussion
While CTD 74 shows brightness temperature values significantly exceeding the equipartition limit, no apparent motion is shown by the jet features even in the highest-resolution VLBI data at 8.3 GHz, therefore superluminal motion could not be derived in its jet.At the redshift of the object, a maximum proper motion of ∼ 0.1 mas yr −1 would already imply apparent superluminal motion of a jet feature.Assuming a moderate Lorentz factor of 10, this would give a critical inclination angle of 1/Γ ∼ 6 • .However, the core-jet component separations in the analyzed VLBI data do not allow for such proper motion value (see Figure 4).The undetectable proper motions might be attributed to the very small inclination angle of the jet to the line of sight, rendering the movement of the components projected on the sky undetectable.On the other hand, it was shown using probabilistic arguments that all of the BL Lac objects with slow moving components cannot have jets with low Lorentz factors and small viewing angles.It was suggested that in those sources, the pattern speed may not be equal to the speed of the jet beam [57].Nevertheless, as long as the proper motion of the jet components cannot be derived, the inclination angle of the jet cannot be safely deduced from the brightness temperature measurements only.
CTD 74 is part of the sample of blazars monitored at 15 GHz with the Owens Valley Radio Observatory (OVRO) as part of the OVRO 40-m telescope flux density monitoring program [58].For a large number of the monitored sources, the variability brightness temperatures were calculated [59].For CTD 74, the obtained value is log 10 54  −0.42 , and the inferred variability Doppler factor is δ var = 56.8+26.95 −20.26 .This is the highest Doppler factor among 20 neutrino-emitter candidate blazars studied using their OVRO 15-GHz variability light curves [7].Within the uncertainties, this value is consistent with the average Doppler factor, δ avg ≈ 37, obtained from the 8-GHz VLBI observations.The slightly lower VLBI value may be because the 8-GHz core can be further resolved at higher frequencies, thus the region responsible for the variability measured at 15 GHz (the 15-GHz core) is more compact.This could be checked with higher-frequency VLBI imaging observations.
The core flux density measured at 8 GHz showed significant variability in the analysed VLBI observations (Figure 6).At the time and until about 4 yr after the high-energy neutrino event, EHEA2012-05-23, the flux density remained in an elevated state, above 0.8 Jy.Additionally, CTD 74 showed significantly higher T b,VLBI 14 months after the neutrino event.The more densely time-sampled long-term OVRO light curve of CTD 74 shows that the radio brightening seen in the flux density of the 8-GHz VLBI-detected core began before the neutrino event, in 2011 ( [7], Figure 8), which time interval was not sampled by VLBI observations.The neutrino event occurred during a relatively steeply rising, ∼ 1.5-yr long section of the 15-GHz total flux density light curve.According to the OVRO light curve, the initial maximum in late 2012 was followed by at least two smaller flares.After 2016, the flux density gradually decreased and reached its minimum value according to the published OVRO light curve in late 2019.Similar fading can be traced in the 8-GHz flux density of the VLBI-detected core component.More than 80% of the 8-GHz VLBI data were obtained after 2010.The OVRO light curve suggests that there was a potentially even brighter flare before 2008, during the time when no 8-GHz VLBI observation was conducted of CTD 74, highlighting the importance of the OVRO monitoring program.
In the standard shock-in-jet models, it is expected that the radio flares are delayed by a few tens to a few hundreds of days from the higher to the lower frequencies [60].The closest peak before our highest estimated T b,VLBI value (2013 Jul 24) in the OVRO light curve is in late November 2012.However, due to the sparse VLBI observations, we cannot securely connect the 8-GHz brightening in the VLBI core to any specific feature in the 15-GHz OVRO light curve.
An increase in the 10-GHz radio emission close in time to the high-energy neutrino events was reported in a sample of VLBI-detected blazars selected as potential counterparts due to their positions with respect to the high-energy neutrino events [6].PKS 1502+106 was identified as showing the highest temporal correlation between the radio flare and the high-energy neutrino event.The neutrino event associated with PKS 1502+106 was detected on 2019 Jul 30 and took place towards the end of a 4-yr long radio flare seen at various radio frequencies between 2016 and 2020 [5].In other two of the reported five highprobability high-energy neutrino candidate sources [5], PKS 1741−038 and PKS 0735+178, the corresponding neutrino events coincided with the rising part of major radio flares.Additional flare-neutrino coincidences include most notably TXS 0506+056, as well as PKS 0215+015 [61] and PKS 1424−418.In the latter case, an increase in γ-ray, X-ray, and optical emission has been reported in temporal and positional coincidence with a PeV energy neutrino event [62].In the context of radio flare-neutrino coincidences, CTD 74 showed similar behaviour.However, enhanced flux density variability (across the entire electromagnetic spectrum) is one of the major characteristics of blazars.Therefore a clear connection between the neutrino event and the observed flux density variability in the radio cannot be established for CTD 74 from the available data.
High-energy neutrinos can be produced in blazar jets via, e.g., the proton-photon process (see [63] and references therein).The ingredients of that, the high-energy photons and protons are principally available in the jets, the former via the synchrotron self-Compton mechanism, while the latter can be accelerated in standing shocks [63].In CTD 74, all detected VLBI features seem to be stationary (Figure 4) which can be explained as the jet having a very small inclination angle to line of sight, but also one (or more) could indeed be stationary, a standing shock formed in the jet.
In the proton-photon formation channel of high-energy neutrinos, photons at GeV energies can also be created in an electromagnetic cascade.However, these photons could hardly be detected due to the expected high optical depth for γ-rays at this region of the jet [63].Thus, a high-energy neutrino candidate source does not necessarily have to be a γ-ray loud blazar.CTD 74 does not appear in the most recent Fermi-LAT catalog of γ-ray-detected extragalactic sources [31].

Summary
We analyzed a series of VLBI measurements conducted between 1996 and 2021 of CTD 74, a blazar that was positionally linked to a high-energy neutrino event detected by IceCube.We reconstructed the brightness distribution of the source and subsequently fitted Gaussian model components to the calibrated visibilities.We were able to identify a couple of jet components both at 2 and 8 GHz in a short northwest-directed jet.During the ∼ 25 yr spanned by the VLBI observations, we did not detect significant motion in any of the jet components.
The 8-GHz brightness temperature of the core exceeded the equipartition brightness temperature value in every epoch, implying Doppler boosting of the emission, and Doppler factors of several tens.
The 8 GHz flux density of the core component revealed a brightening at around the end of 2011.This several year-long flare can be seen in the densely-sampled 15-GHz total flux density light curve from the OVRO monitoring program as well.The neutrino event took place during the rising phase of this radio brightening, preceding the maximum by about half a year.
We found that CTD 74, which was positionally associated [6] with a high-energy neutrino event in 2012, shows several properties that are in line with other blazars listed as possible neutrino sources.These are the high Doppler factor, the presence of stationary features in the jet, and enhanced flux density variability around the time of the neutrino event.However, since blazars generally have Doppler-boosted jet emission and are variable, a comprehensive theoretical framework of the neutrino emission from blazars, with observationally testable predictions would be essential for reaching firmer conclusion in cases like CTD 74.Moreover, rapid-response quasi-contemporaneous VLBI imaging of candidate high-energy neutrino source blazars, triggered by new neutrino detections, would help reveal if there is clear connection between the neutrino emission and the changes in the radio structure.Obviously, continuing high-cadence flux density monitoring of a large sample of potential neutrino-emitting blazars is needed for associating the neutrino events to outburst or any other light curve feature.

Figure 1 .
Figure 1.The (u, v) coverages of two observations performed at 8.6 GHz.The axes represent the baseline vectors in u and v coordinates in units of million wavelengths.Left: The (u, v) coverage of the observation performed on 2017 March 27 with VLBA-only baselines.Right: The (u, v) coverage of the observation performed on 2020 November 18 with VLBA and global baselines.

Figure 2 .
Figure 2. VLBI images of CTD 74 obtained through hybrid mapping for the measurements taken in 1996 at 2, 5, and 8 GHz in panels (a), (b), and (c), respectively.In all maps, the positions of the Gaussian components fitted to the visibility data are marked.Panel (a): The peak intensity is 1160 mJy beam −1 .The lowest contours are at ±3.44 mJy beam −1 .The restoring beam size is 7.54 mas × 3.34 mas.The position angle of the major axis is PA = −8.35• .Panel (b): The peak intensity is 1020 mJy beam −1 .The lowest contours are at ±5.44 mJy beam −1 .The restoring beam size is 3.92 mas × 1.44 mas, PA = −4.63• .Panel (c): The peak intensity is 746 mJy beam −1 .The lowest contours are at ±4.45 mJy beam −1 .The restoring beam size is 2.07 mas × 0.926 mas, PA = −8.59• .In each map, the lowest positive contour is drawn at 3σ image noise level, and subsequent contours increase by a factor of two.The restoring beams are shown in the lower left corners of each panel.

Figure 3 .
Figure 3. VLBI image of CTD 74 taken on 2014 Feb 12 at 8 GHz, created via Gaussian model fitting to the self-calibrated visibilities.The location and size (FWHM) of the model components are indicated with crosses and circles.The peak intensity is 867 mJy beam −1 .The lowest contour is at ±3.7 mJy beam −1 , at 3 times the image noise level.The restoring beam size is 1.0 mas × 0.77 mas, PA = −78.8• .It is shown in the lower left corner of the image.

Figure 4 .
Figure 4. Separation of the jet components from the core as a function of time at 8 GHz.The apparent proper motion of J1 and J2 components that could consistently be identified through most of the epochs was modeled with a linear function.These functions are indicated by the dashed red and solid yellow lines, respectively.The blue vertical line marks the time of the neutrino event EHEA2012-05-23 (2012 May 23).

Figure 5 .
Figure 5.The position angles of components J2 (green) and S2 (magenta) as a function of time.The angles are measured from north through east.The average and standard deviation of the position angles are indicated.The blue vertical line marks the time of the neutrino event EHEA2012-05-23 (2012 May 23).

Figure 6 .
Figure 6.Top: Flux density of the core component as a function of time at 8 GHz (black squares).Bottom: Peak intensity (purple circles) in each image restored with the same beam (see details in the text).The errors of the peak intensities are comparable to the size of the symbol.The blue vertical line marks the time of the neutrino event EHEA2012-05-23 (2012 May 23).

Figure 7 .
Figure 7. Brightness temperature of the core component as a function of time, as measured at 8 GHz.The blue vertical line marks the time of the neutrino event EHEA2012-05-23 (2012 May 23).The red horizontal line indicates the value of the equipartition brightness temperature, 5 • 10 10 K [54].

Figure 8 .
Figure 8.The radio light curve of CTD 74 at 15 GHz as measured by the OVRO monitoring program (red crosses, [7]) along with the flux densities of the core component obtained from the 8-GHz VLBI measurements (black squares).The time of the neutrino event EHEA2012-05-23 is indicated by the blue vertical line.