Superluminal Motion and Jet Parameters in the Gamma-ray-Emitting Narrow-Line Seyfert 1 Galaxy TXS 1206+549

: Narrow-line Seyfert 1 (NLS1) galaxies are a peculiar subclass of active galactic nuclei (AGN). Among them, TXS 1206+549 belongs to a small group of radio-loud and γ -ray-emitting NLS1 galaxies. We focus on the radio properties of this galaxy by analysing archival, high-resolution, very long baseline interferometry (VLBI) imaging observations taken at 8 GHz frequency in six epochs between 1994 and 2018. Using the milliarcsecond-scale radio structure, we can resolve a core and a jet component whose angular separation increases by ( 0.055 ± 0.006 ) masyr − 1 . This corresponds to an apparent superluminal jet component motion of ( 3.5 ± 0.4 ) c . From the core brightness temperature and the jet component proper motion, we determine the characteristic Doppler-boosting factor, the bulk Lorentz factor, and the jet viewing angle. We find no compelling evidence for a very closely aligned blazar-type jet. The parameters for TXS 1206 + 549 resemble those of radio-loud quasar jets with a moderate Lorentz factor ( Γ ≈ 4) and ϑ ≈ 24 ◦ inclination to the line of sight.


Introduction
Seyfert galaxies [1] have bright, unresolved nuclei in optical images, with broad emission lines in their spectra.They typically reside in spiral host galaxies and belong to the general class of active galactic nuclei (AGN).These objects are divided into two main categories: Seyfert 1 (Sy1) and Seyfert 2 (Sy2) galaxies.According to the unified scheme of AGN, e.g., [2], Sy1 galaxies are seen nearly face-on and the broad-line region (BLR), composed of gas clouds orbiting around the central supermassive black hole (SMBH) with speeds of several thousand km s −1 , is visible from our line of sight.This gives rise to permitted lines in the HI Balmer series that are broader than the forbidden lines [3].In turn, Sy2 galaxies are viewed nearly edge-on and the dusty torus obscures the radiation from the BLR.Therefore, Sy2 spectra are characterised by narrower emission lines coming from the so-called narrow-line region (NLR), i.e., from ionised gas clouds orbiting further away from the SMBH.The full width at half maximum (FWHM) of the forbidden lines is typically 300-800 km s −1 for both Sy1 and Sy2, while the Balmer lines are much broader, reaching up to 6000 km s −1 in Sy1 galaxies [3].There are also Seyfert galaxies with intermediate properties [4].
There are also AGN with general Sy1 properties but a much narrower Hβ emission line in their spectra compared to what is typical for the Sy1 group.These are called narrow-line Sy1 (NLS1) galaxies [3].Their defining characteristics are the Hβ line with FWHM < 2000 km s −1 , and the relative weakness of the [OIII]λ5007 Å emission, [OIII]/Hβ < 3 [5].Because they usually exhibit strong FeII emissions, the presence of these lines in the spectrum is also often included in the definition.
The majority of AGN are radio-quiet or radio-silent; only a few of them (∼10%) are radio-loud [6].Because their radio emissions are dominated by the synchrotron emission of relativistic plasma jets ejected from the vicinity of the central SMBH, these are also called jetted AGN [7].Almost all extragalactic γ-ray sources identified with counterparts in other electromagnetic wavebands in the 4th Fermi Large-Area Telescope (LAT) catalogue are radio-loud jetted AGN [8].A similar fraction of radio-loud vs. radio-quiet sources is found for the NLS1 sub-class: only about ∼8% have a radio-loudness parameter, i.e., the radio to optical flux density ratio, R = f ν (5 GHz)/ f ν (4400 Å) > 10, and very few have R > 100 [9].The detection of γ-ray emissions from four radio-loud NLS1 (RL-NLS1) objects after the first year of operations of the Fermi LAT instrument [10,11] confirmed the existence of powerful relativistic jets in these AGN.
The observations suggest that radio-loud, jetted AGN are hosted in elliptical galaxies with high black hole (BH) mass.It seems that lower-mass (M ≲ 10 8 M ⊙ ) central BHs in late-type galaxies are generally not capable of launching powerful jets, establishing a strong connection between BH mass and radio loudness [12].However, with the discovery of not just NLS1s but other types of radio galaxies [13] exhibiting γ-ray emissions, this BH mass-radio loudness correlation has been under debate.NLS1s hosted in diverse types of galaxies, including spirals [14], are considered to have low-mass (∼10 6 -10 8 M ⊙ ) black holes and high accretion rates, close to or even above the Eddington limit [14,15].It is possible that the γ-ray-emitting RL-NLS1s occupy the upper end of the NLS1 mass spectrum [16].It has also been proposed that these are the low-luminosity counterparts of the flat-spectrum radio quasars (FSRQs), and there might be an evolutionary connection between them, as RL-NLS1 sources comprise the still-evolving, younger phases of FSRQs [17].There is no doubt that these peculiar γ-ray-emitting RL-NLS1s play an important role in understanding how jetted AGN differ from their radio-quiet counterparts and what contributes to the production of powerful jets [18,19].
One of the few γ-ray-detected RL-NLS1s that was recently found [20] is TXS 1206+549 (also known as J1208+5441) at a redshift of z = 1.344 [21].This is the most distant γray NLS1 found to date.In the Fermi-4LAC catalogue, it is classified as FSRQ [22].The discoverers obtained a near-infrared spectrum and measured the FMWH of the Hβ line as (1194 ± 77) km s −1 and the [OIII] to Hβ flux ratio as ∼0.7, so TXS 1206+549 satisfies the NLS1 criteria.Although an ultraviolet FeII line was found for the object, optical FeII lines were not seen in the spectrum [20].To estimate the BH mas in TXS 1206+549, different methods were employed, resulting in values of ∼10 7 -10 8 M ⊙ , consistent with what is found for NLS1 objects in general and γ-ray-detected NLS1s in particular.These values are one or two orders of magnitude lower than the SMBH masses of powerful jetted AGN [20].The high radio loudness parameter R ≈ 1300 indicates that TXS 1206+549 belongs to RL-NLS1s, although the variability known in both the optical and radio wavebands means that R is also variable.
Most recently, the presence of a quasiperiodic oscillation (QPO) with a period of about 134 d in the Zwicky Transient Facility optical light curve of TXS 1206+549 was proposed [23].The signal was present between 2018 and 2021.QPOs are usually found at multiple wavebands in AGN, including NLS1 sources, but their physical origin is not well understood a present.
TXS 1206+549 is a well-known, flat-spectrum radio source, with variable (∼120-400 mJy) total flux density across two orders of magnitude in frequency, from 150 MHz [24,25] to 15 GHz [26].Its extended (∼100 kpc scale) structure at low radio frequencies (144 and 642 MHz) is composed of two pairs of symmetric lobes that are perpendicular to each other (i.e., roughly east-west and north-south), hinting at a sudden realignment of the jet direction in the past [27].This extended relic radio structure is rarely seen in RL-NLS1 sources [28][29][30].To be compatible with the notion that γ-ray RL-NLS1 sources are "young" versions of blazars [19], the central source in TXS 1206+549 may indicate episodic jet activity.
When observed at a high resolution with very long baseline interferometry (VLBI), the source shows a compact, nearly unresolved milliarcsecond (mas)-scale structure at 5 GHz, slightly elongated to the south [31,32].This structure is similar to that of blazars with powerful relativistic jets, and, notably, of a couple of other γ-ray RL-NLS1 sources monitored with VLBI (see, e.g., [33] for a review).Its accurate coordinates in the VLBI-defined International Celestial Reference Frame are right ascension RA = 12 h 08 min (54.256432± 0.000014) s and declination Dec = 54 • 41 ′ (58.16329 ± 0.00014) ′′ [34].
In this paper, we analyse archival multi-epoch 8-GHz VLBI data to study the jet kinematics of TXS 1206+549 for the first time.Tracking the positional changes in a component located in the inner 1-mas section of the jet, identified across 5 different epochs covering a period of 7 yr, we determine its apparent proper motion.We also calculate core brightness temperatures, and estimate the bulk Lorentz factor and the inclination of the jet.In Section 2, we describe the observational data and the analysis method.In Section 3, we present the results of the imaging, brightness distribution modelling, and jet parameter estimation.In Section 4, we discuss our findings and the properties of TXS 1206+549.Finally, a brief summary is given in Section 5.
Throughout this paper, we assume a flat Λ Cold Dark Matter cosmological model with Hubble constant H 0 = 70 km s −1 Mpc −1 , vacuum energy density parameter Ω Λ = 0.7, and matter density parameter Ω m = 0.3.In this model, the luminosity distance of TXS 1206+549 at z = 1.344 is D L = 9524.4Mpc and the angular scale is 8.404 pc mas −1 [35].

Observations and Data Analysis
We analysed archival observations conducted with the U.S. Very Long Baseline Array (VLBA) [36] over 6 different epochs covering more than 24 yr between 1994 and 2018.
The calibrated visibility data were acquired from the Astrogeo VLBI data base (http:// astrogeo.org/cgi-bin/imdb_get_source.csh?source=J1208%2B5441, accessed on 20 January 2024) maintained by L. Petrov.From the three different frequency bands (around 2, 5, and 8 GHz) available in the archive, we selected X-band data at 8.34-8.67GHz frequencies because their angular resolution is the highest and the time coverage is the best.This allows us to detect possible changes in the inner jet structure of TXS 1206+549.The short snapshot VLBA data were taken with all 10 radio telescopes of the network with right circular polarisation.Further details (mean epochs, central observing frequencies within the 8-GHz band, on-source integration times, bandwidths) are listed in Table 1.To produce radio images, we imported the calibrated visibility data into the Caltech DIFMAP program [39].We followed the standard hybrid mapping procedure prescribed in the DIFMAP Cookbook, with several iterations of CLEAN deconvolution [40], phase-only, and finally phase and amplitude self-calibration, until we reached the lowest achievable noise level in the residual image.Then, we used the self-calibrated visibility data to model the source brightness distribution with two-dimensional circular Gaussian components [41].The parameters of these model components describe the source structure.

Source Structure
The mas-scale radio structure of TXS 1206+549 is compact at 8 GHz, as indicated in Figure 1, shown as an example.In all six epochs, from 1994 to 2018, the brightness distribution modelling reveals a compact core (C), which corresponds to the synchrotron, self-absorbed base of the relativistic jet, and a more extended and weaker jet component can be seen south-southeast of the core.The parameters of the fitted circular Gaussian model components are listed in Table 2.The sizes of the components always exceeded the minimum resolvable angular size with the interferometer array [42].We generally adopted the method for parameter error estimates from [43].However, the formal errors of flux density values were conservatively set at 10% to account for the possible VLBI absolute amplitude calibration uncertainties [44].The uncertainty of the jet component position angle is taken as 1 • .For an error estimation of the size of the fitted Gaussian components, we used the approach in [45], which provides more realistic errors for fitted components that are smaller than the restoring beam size.In the first epoch, the core-jet component separation is 1.56 mas, while in later epochs, starting from 2011, the separation is within 1 mas.We interpret these as two different components in the jet structure of TXS 1206+549, denoted by J1 (in 1994) and J2 (from 2011 to 2018).Apparently the J1 component moved away from the core and became too faint and extended for detection with VLBI by the 2010s.During the 17 yr long gap in the archival observations, a new component, J2, must have emerged, apparently following the track of the earlier J1.Since the time sampling of the data in the 2010s is much better, we could use the last five epochs between 2011 and 2018 for the jet kinematic study of J2.

Jet Component Proper Motion
Based on our fitted brightness distribution models (Table 2), the angular separation between the core and J2 components monotonically increased from 2011 to 2018, consistently with an outward jet motion.The values are plotted in Figure 2. Assuming there is a linear change in the core-jet separation, the apparent proper motion of J2 is µ = (0.055 ± 0.006) mas yr −1 .Extrapolating back in time, the component ejection could have occurred around the year 2002.This proper motion can be expressed as the speed of light (c) as follows: where µ is measured in mas yr −1 and the luminosity distance D L is expressed in Mpc, e.g., [46].By substituting z = 1.344 and D L = 9524.4Mpc, we can obtain β app = 3.5 ± 0.4 apparent transverse speed.As β app clearly exceeds unity, the J2 component is characterised by apparent superluminal motion.

Jet Parameters
The knowledge of the apparent transverse speed in the inner jet, together with the measurement of the brightness temperature in the compact core, allow us to estimate the basic physical and geometric parameters of the relativistic jet in TXS 1206+549.
The core brightness temperature values listed in the last column of Table 2 for each observed epoch are calculated as follows: [47] Here, the core flux density S is measured in Jy, the observing frequency ν in GHz, and the fitted circular Gaussian component diameter (FWHM) d in mas.The measured brightness temperatures are close to or exceed the equipartition value, T b,eq ≈ 5 × 10 10 K [48].The latter indicate Doppler-boosted radio emissions from the jet.Since brightness temperatures vary considerably (by more than an order of magnitude) from epoch to epoch, we followed the approach of [49] and took the median value, T b,med = 9.1 × 10 10 K as characteristic to the jet.The advantage of using the median brightness temperature is that this value is insensitive to outliers due to source variability.Note that by far the highest measured brightness temperature, T b = (5.6 ± 0.2) × 10 11 K, is caused by a major outburst, as indicated by the high core flux density in 2017.According to Equation (2), the small angular size of the core component also contributes to the high T b (Table 2).The characteristic Doppler-boosting factor is calculated as follows: where we adopt T b,int = 4.1 × 10 10 K as the intrinsic brightness temperature.This is a typical value found for a sample of radio quasars that are not in the outburst [49].Notably, T b,int ≲ T b,eq .For the jet in TXS 1206+549, we obtained δ = 2.2.
From the characteristic Doppler-boosting factor and the apparent transverse speed of the inner sub-mas jet component J2, the bulk Lorentz factor (Γ) of the plasma and the viewing angle of the jet (ϑ) with respect to the line of sight can be calculated, e.g., [6]: After substituting the numerical parameters derived above, β app = 3.5 ± 0.4 and δ = 2.2, we obtain Γ = 4.1 ± 0.6 and ϑ = (23.5 ± 0.9) • .
To investigate the uncertainties inherent in our method of using T b,med for estimating the characteristic Doppler factor (Equation (3)), we calculated possible lower and upper limits for the Doppler factor, based on individual T b measurements taken from all observing epochs but the one in 2017 (Table 2).The exclusion of the 2017 epoch when the source was in outburst is justified because equipartition conditions are known to be violated in the maximum brightness state [50].Therefore, using the intrinsic brightness temperature characteristic of the low state in the denominator of Equation ( 3) would lead to an unrealistically high Doppler factor value.Assuming the range 3.4 × 10 10 K ≤ T b ≤ 1.75 × 10 11 K, and again adopting T b,int = 4.1 × 10 10 K as the intrinsic brightness temperature [49], the allowed range of Doppler factors is 0.8 ≤ δ ≤ 4.3.While the lower limit implies no Doppler boosting, the upper limit, together with the apparent superluminal speed found for the jet component J2, would result in a 3.7 Lorentz factor and 13.2 • jet inclination (Equations ( 4) and ( 5)).

Discussion
According to our analysis of 8-GHz VLBA data, the jet in TXS 1206+549 is characterised by a moderate Lorentz factor, Γ ≈ 4. Studies of powerful jetted AGN samples e.g., [51] found that the distribution of Lorentz factors peaks between 5 and 15, with high values occasionally reaching Γ max ≈ 50.The inclination of the jet (ϑ ≈ 24 • ) is larger than expected from a blazar-type source, where the jet is supposed to closely align with the line of sight within about 10 • .Note, however, that even a moderate jet bending well inside the region resolvable with VLBI at 8 GHz (i.e., within <1 mas, corresponding to <8 pc in projected linear size) might cause the innermost section of the jet to be better aligned to the line of sight.This could possibly be tested with higher-resolution VLBI imaging at higher frequencies.Nevertheless, the jet inclination is generally consistent with the Sy1 classification of the galaxy (i.e., seen face-on) if we naturally assume that the jet direction is perpendicular to the plane of the accretion disk and the obscuring "torus".
With the apparent speed measured for the jet component J2 (β app = 3.5), a Doppler factor of δ ≳ 5.3 would be required for a blazar-like jet inclination.Note that we characterised the jet using the median core brightness temperature T b,med , and there were individual T b measurements that could, in principle, lead to higher Doppler factor estimates according to Equation (3).Indeed, the range of brightness temperatures measured at epochs when the source was in a low state would allow for a δ of as high as 4.3 and a jet inclination as small as 13.2 • .However, the extremely high T b = (5.6 ± 0.2) × 10 11 K measured in 2017 is clearly associated with an outburst when the intrinsic brightness temperature can significantly exceed the equipartition value, with characteristic value estimated as T b,int ∼ 2 × 10 11 K [50].If so, the Doppler factor is comparable to the value of δ = 2.2 we determined from the median brightness temperature that represents the quiescent state of the jet.Therefore, we believe our characteristic Doppler factor estimate is robust.TXS 1206+549 is known to have rapid, high-amplitude flux density variations in the radio.It is part of the AGN sample monitored at 15 GHz with the Owens Valley Radio Observatory (OVRO) 40-m radio telescope [52].For example, the published light curve [26] (available at https://vizier.cds.unistra.fr/viz-bin/VizieR?-source=J/MNRAS/438/3058, accessed on 20 January 2024) indicates almost a factor of 4 change in total flux density during the relatively short interval between April 2009 and December 2011, with a typical peak-to-peak time difference of just over ∼1 yr.Unfortunately, this monitoring interval has little overlap with our archival VLBA data (Table 1).Close to our only epoch falling in the available OVRO monitoring time range, 2011.808yr, the OVRO 15-GHz flux density was at around its local minimum, with ∼130 mJy.With an OVRO light curve continuing into the 2010s, it would be possible to independently verify the outburst seen in 2017 in our sparsely time-sampled VLBA data.
The OVRO flux density observations were applied to estimate variability brightness temperatures for the sources in the monitored sample, including TXS 1206+549 [53].The method is based on fitting flares identified in the radio light curve with exponential rise and decay functions.The value T var = 5 × 10 14 K found for TXS 1206+549 led to the estimate of the variability Doppler factor δ var = 29.09+10.28 −10.02 .This is significantly, by an order of magnitude, higher than our estimate based on the measured median brightness temperature and the assumed intrinsic brightness temperature, which is close to the equipartition value.A comparison of TXS 1206+549 with other γ-ray-emitting RL-NLS1 sources indicates that the characteristic brightness temperature we determined is similar to those found for other sources, in the order of T b ∼ 10 10−11 K [32,[54][55][56].
Under the name of the Fermi γ-ray source 4FGL J1208.9+5441,our RL-NLS1 object has various estimates of the Doppler factor found in the literature.The value of 15.24 was determined for the so-called luminosity Doppler factor derived from the observed γ-ray luminosity and the photon index [57].Another publication dealing with Fermi blazar jet properties [58] arrived at δ > 8.12 by comparing the jet radiation power and the γray luminosity, assuming that the γ-ray luminosity is representative of the jet bolometric luminosity and δ = Γ for blazars.(The latter assumption does not hold according to our VLBI results.)This paper [58] also provides an alternative value of 15.12 calculated based on a correlation between the jet radiative power and the luminosity of the BLR [59].However, another method, based on a fitting model for the broad-band spectral energy distribution (SED), again assuming equality between the Doppler factor and the Lorentz factor, led to a different value: δ = 7.8 [60].Most recently, δ = 25.69 was derived for 4FGL J1208.9+5441, using SED fitting to quasi-simultaneous multi-wavelength data and assuming a one-zone leptonic model [61].
As we can see from the values collected above, a wide range of Doppler factors can be estimated for TXS 1206+549, with our δ = 2.2 being the lowest.Clearly, there are several differences in the methods used in the literature, for example, in the frequency bands that are observed, and therefore the regions of the jet that are probed, in the observed epoch (i.e., in quiescence or outburst), and in the physical model assumptions used.It is known that while different determinations of the Doppler factors are generally statistically consistent with each other for populations of sources, there might indeed be large discrepancies for individual objects [62].
The cited studies estimating the Doppler factor from high-energy data consider TXS 1206+549 as a blazar.However, there are no generally accepted criteria for the blazar classification of AGN.The working definition of blazars could vary from waveband to waveband, depending of the types of observing data that are available and the purpose of the study.The basic concept is that blazars are jetted AGN, whose relativistic jet (more specifically, one of the intrinsically symmetric bipolar jets on opposite sides of the central BH) is directed very closely to the observer.Exactly how closely is a matter of definition.The ϑ < 10 • criterion, e.g., [63][64][65], that we adopt in this paper is probably the simplest, although some authors cite ϑ < 20 • , e.g., [66].Alternatively, for statistical studies, it can be practical to apply ϑ ≤ 1/Γ, e.g., [67,68].Note that, according to this latter definition and based on our Lorentz factor estimate, Γ = 4.1, the jet inclination angle upper limit in our case would be ϑ lim = 0.2439 rad ≈ 14 • , still lower than the ϑ ≈ 24 • we derived.In any case, the jet inclination angle can only be measured indirectly from VLBI kinematic studies similar to ours (Section 3) and using certain assumptions for the jet physics.Moreover, in large part because such measurements require considerable observational resources at multiple epochs, this method is only applicable to a very limited sample of the sources compared to all blazar candidates.Thus, viewing angles are not too practical to use for a definition.
As a consequence of the closely aligned relativistic jets with respect to the line of sight in blazars, VLBI observations often indicate the apparent superluminal motion of the jet components.However, on the one hand, very small inclinations (ϑ ≈ 0 • ) do not lead to detectable transverse motion in the plane of the sky.On the other hand, β app > 1 values are possible even at relatively high jet inclinations (see, e.g., Figure 21 in [6]).According to our results, this is the case for TXS 1206+549 as well.Therefore, the detection of superluminal motion cannot be used as a defining criterion for blazars.
Small-inclination relativistic jets in AGN produce highly enhanced synchrotron radiation via Doppler boosting, high-energy emissions via the inverse Compton process, and may show rapid and high-amplitude variability in the broadest possible range of electromagnetic wavebands.Therefore, these criteria are often applied for a "loose" definition of blazars.For example, the important and widely used multi-frequency catalogue of blazars, the Roma-BZCAT [69,70] collection contains AGN with optical identification and blazar-like spectrum, high X-ray luminosity, and radio detection at ∼1 GHz frequency with "compact" radio morphology.However, the radio emission confined to within ∼10 ′′ angular scales does not guarantee the kind of compactness one would expect from a closely-aligned jet.Indeed, it was found recently that a significant fraction, 3% of the so-called confirmed Roma-BZCAT entries show an extended radio structure at ∼1 ′′ scale and are likely misclassified as blazars [66], not to mention VLBI imaging, which provides a higher resolution by orders of magnitude to potentially refute the blazar nature of certain candidates, e.g., [71].Finally, to demonstrate the confusion about the blazar classification in the literature, we shall briefly mention the so-called "radio-weak" BL Lac objects, e.g., [72][73][74].The term itself is contradictory since the presence of a small-inclination relativistic jet is an essential ingredient of a blazar, and should be a source of Doppler-beamed radio emission via the synchrotron process.
As was found from multi-waveband studies of a large sample of 42 flat-spectrum RL-NLS1 sources [75], their central BH masses are lower than those for blazars, while the accretion rates are generally comparable.The jet power per unit mass is very similar for these different types of jetted AGN, suggesting that the observational differences are mainly caused by scaling effects.In other words, the jet launching mechanism in the RL-NLS1 sources is similar to that in blazars, e.g., [76], the only major difference being the BH mass, which is lower in RL-NLS1 sources.For a sample of 16 γ-ray RL-NLS1 sources, it was found that the jet power in almost all of them is consistent with a hybrid jet formation model involving a mixture of both the Blandford-Payne and Blandford-Znajek mechanisms [77].If the general picture regarding the similar jet launching required for blazars and RL-NLS1 sources is true, then there is no obvious evolutionary connection between RL-NLS1 sources and FSRQs that harbour one or two orders of magnitude more massive central BHs.Because of the low-mass BH, variability time scales in RL-NLS1 sources are shorter than those expected from Doppler boosting only [75].In the case of TXS 1206+549, this might explain the much higher variability Doppler factor compared to our jet kinematic estimate.
Concerning the dichotomy between the radio-loud and radio-quiet AGN, the properties of the SMBHs, as well as their environments, should play a role.Simulations suggest that the jet power is a strong function of the BH spin in the presence of a thick accretion disk [78].However, observations indicate that most SMBHs rotate rapidly [79]; therefore, a high black hole spin alone cannot be responsible for the jet production, and disk magnetization could also be an important factor (see the discussion and references in [80]).In any case, as the radio-loud fraction of NLS1 sources is similar to that of the general AGN population, it is reasonable to assume that the same effects govern these objects, which have lower-mass BHs in their central engine.

Summary and Conclusions
We analyzed the mas-scale radio structure of the γ-ray-detected RL-NLS1 galaxy TXS 1206+549 at z = 1.344.This is the highest-redshift γ-ray NLS1 known to date [20].We used the archival VLBA data available in the Astrogeo data base in 6 epochs from 1994 to 2018.Imaging and brightness distribution modelling allowed for us to track the motion of a jet component within 1 mas from the core from 2011 to 2018, and to determine its superluminal apparent transverse speed β app = 3.5 ± 0.4.We characterised the brightness temperature of the compact core with the median of the measured values, and derived the Doppler factor δ = 2.2 by assuming an intrinsic brightness temperature typical for jetted AGN that are not in a flux density outburst.From the apparent transverse speed and the characteristic Doppler factor, we obtained Γ = 4.1 ± 0.6 for the bulk Lorentz factor for the jet plasma and ϑ = (23.5 ± 0.9) • for the jet inclination angle with respect to the line of sight.The latter is consistent with the Sy1 classification of the source.Some properties of TXS 1206+549-the γ-ray emission [20], the superluminal jet motion, and the rapid flux density variability in multiple wavebands [20] in particular-are similar to those of powerful blazars.Because the jet inclination angle is found to exceed (10-20) • , strictly speaking, we did not find strong evidence for its blazar-like nature from the 8-GHz VLBI imaging data.However, blazar definitions are usually more liberal in the literature, and ϑ is not particularly large.A moderate bending of the innermost jet section could help to reconcile TXS 1206+549 with a blazar-type object in terms of jet inclination angle.Dedicated higher-resolution VLBI imaging at higher observing frequencies would be beneficial to test this possibility and probably also to constrain the core brightness temperature more tightly.

Figure 1 .
Figure 1.Naturally weighted 8.7-GHz VLBA images of TXS 1206+549 from 30 July 2018, deconvolved with CLEAN components (left) and a model composed of two Gaussian components whose positions and FWHM diameters are indicated with cross-hairs in yellow circles (right).In both images, the peak intensity is 167 mJy beam −1 , the contour levels start at ±0.6 mJy beam −1 (∼3σ image noise), and further positive contours increase by a factor of 2. The elliptical Gaussian restoring beam size is 1.69 mas × 0.92 mas (FWHM), the major axis position angle is 7.2 • (measured from north through east), as indicated in the bottom left corner of the images.

Figure 2 .
Figure 2. The core-jet separation for the J2 component in TXS 1206+549 as a function of time, based on VLBI observations and model fitting at 8 GHz in five epochs covering a period of about 7 yr.Linear regression indicates µ = (0.055 ± 0.006) mas yr −1 apparent proper motion.

Table 1 .
Details of the archival 8-GHz VLBA observations obtained from the Astrogeo data base Notes: Col. 1: mean observing epoch; Col. 2: central observing frequency; Col. 3: on-source integration time; Col. 4: number of intermediate frequency channels times bandwidth; Col. 5: project code and the literature reference (if available).