VLA Snapshot Continuum Survey of FRI Quasar Candidates Selected from the LOFAR Two-Metre Sky Survey (LoTSS)

The radiative and jet power in active galactic nuclei is generated by accretion of material on to supermassive galactic-centre black holes. For quasars, where the radiative power is by definition very high, objects with high radio luminosities form around 10 per cent of the population, although it is not clear whether this is a stable phase. Traditionally, quasars with high radio luminosities have been thought to present jets with edge-brightened morphology (Fanaroff-Riley II - FR II) due to the limitations of previous radio surveys (i.e., FRIs were not observed as part of the quasar population). The LOw Frequency ARray (LOFAR) Two-metre Sky Survey (LoTSS) with its unprecedented sensitivity and resolution covering wide sky areas has enabled the first systematic selection and investigation of quasars with core-brightened morphology (Fanaroff-Riley I - FR). We carried out a Very Large Array (VLA) snapshot survey to reveal inner structures of jets in selected quasar candidates; 15 (25 per cent) out of 60 sources show clear inner jet structures that are diagnostic of FRI jets and 13 quasars (around 22 per cent) show extended structures similar to those of FRI jets. Black hole masses and Eddington ratios do not show a clear difference between FRI and FRII quasars. FRII quasars tend to have higher jet powers than FRI quasars. Our results show that the occurrence of FRI jets in powerful radiatively efficient systems is not common, probably mainly due to two factors: galaxy environment and jet power.


Introduction
Quasars are a class of Active Galactic Nuclei (AGN) that produce immense energy by the fast accretion onto supermassive black holes in the centre of massive galaxies. Based on the unification model (in which the main ingredient is the orientation effect) radio galaxies and quasars are the same objects seen at different angles. An object is called a quasar if it is viewed within a cone of half-angle approximately 45 • . They are common at moderate and high redshifts which make them excellent sites to probe the conditions of the early Universe, e.g., [1]. Investigations of quasars are not only crucial in their own right, but also essential for galaxy formation and AGN feedback studies, e.g., [2][3][4]. They allow us to constrain cosmological parameters, e.g., [5] and to investigate the radio emission in low-luminosity (or faint) AGN at moderate/high redshifts, owing to the strong lensing effect which is just one of the powerful tools for studying AGN at moderate/high redshifts, [6,7]. For a recent review on radio galaxies and radio-loud quasars we refer the reader to [8].
Fanaroff & Riley [9] classified radio sources based on their radio morphologies as FRI (core-brightened) and FRII (edge-brightened) and suggested that the division between FRI and FRII is not only in radio morphology but also in radio luminosity (i.e., the traditional FRI/FRII break is expected to be ∼ 5 × 10 25 W Hz −1 at 150 MHz). It is worth noting that this picture is changing with new sensitive radio surveys because many more FRII AGN below the traditional arXiv:2112.13160v1 [astro-ph.GA] 24 Dec 2021 FRI/FRII break have been discovered [Mingo et al. in prep] [10][11][12][13][14] but it is still the case that FRI sources tend not to be found at the highest radio powers.
Among AGN, optically powerful quasars have not traditionally been associated with FRI radio morphology. The reason behind this is not based on any motivated physical model, e.g., [15,16] but is due to a lack of these sources in the previous surveys (such as the 3C surveys, [17,18]). It should be mentioned that radiatively efficient AGN with a given radiative power can have jet powers that range continuously from zero up to as much as their radiative power ( [4] and references therein), but this does not rule out FRI quasars. Since previous radio surveys revealed quasars having extended emission and only associated with FRII morphology a quasar has broadly been expected to be either a radio-quiet source (i.e., with no significant radio emission except, perhaps core emission) or a radio-loud source with FRII morphology. There are certainly radio-loud sources showing compact structures like compact steep-spectrum sources and gigahertz peaked sources so naturally there also needs to be a continuum of sources between them. Refs. [19,20] discussed this bias and provided evidence that some quasars present FRI-like structures (six objects). However, radio surveys at that time were not sufficient to provide a large sample of FRI quasars only a few candidates from the dedicated efforts by [20] as the surveys were not sensitive enough to reveal any FRI structures. Recently, Ref. [21] showed that some of the X-ray selected quasars do not show any radio emission (not even a core emission is detected). These results also suggest that systems with varying jet power can exist independent of AGN power.
Hitherto, moderate to high-redshift FRIs have been poorly studied in comparison to FRII sources as FRIIs are typically brighter than FRI objects and can be detected out to high redshifts, e.g., [22]. In addition, the FRIs were found to be relatively low-power radio sources with jets that are relativistic on parsec scales, but decelerate to sub-relativistic bulk speeds on kpc scales. Quasars, on the other hand, are powerful sources amongst the AGN population. Understanding the occurrence and nature of quasars with FRI morphology on a scale of tens of kpc, how they compare to FRI radio galaxies on hundreds of kpc scales and how they fit in the interplay between AGN accretion, radio morphology and source environment will be only possible through studies based on statistically complete samples.
Systematic identification of large samples of FRI quasar candidates has only been possible with the LOw Frequency ARray (LOFAR) Two-metre Sky Survey (LoTSS, [23]) thanks to its unprecedented capabilities in terms of high sensitivity, angular resolution and survey speed at low radio frequencies. Tier 1 of LoTSS is designed to survey the northern sky using the LOFAR High Band Antennas (HBA) at an average frequency of 144 MHz. The survey has already collected data over large areas reaching the target sensitivity (∼70 µJy beam −1 ) and angular resolution of 6 arcsec [Data Release I/DR1] [23]. We [3] recently selected a sample of quasars from the Sloan Digital Sky Survey (SDSS) quasar catalogue fourteenth data release [DR14Q;] [24] over the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) region and evaluated their radio properties using the LoTSS data set. We found that optically selected quasars have a wide range of radio continuum properties for a given accretion power and a simple bimodality is not an adequate way to describe the radio properties of quasars. In a recent study we also presented a model of the radio luminosity distribution of the quasars that assumes that every quasar displays a superposition of two sources of radio emission: AGN (jets) and star formation [4]. The observed radio flux density distribution of the quasar sample that spans a wide range of redshift and optical luminosity is well described by the two-component model. This reinforces the idea that the jet-launching mechanism operates in all quasars but with different efficiencies, which also confirms the findings by [3] and is consistent with the lack of bimodality (see for studies suggesting the opposite [25,26]) observed for less powerful populations [27,28]. Therefore, quasars with low-power kinetic output (also quasars with moderate power but in particular in rich environments) are expected to form FRI structures.
With the motivation of investigating FRI quasars we carried out a pilot project using the Jansky Very Large Array (VLA) to perform a snapshot survey of a sample of 60 quasars showing extended structures that would be classified as FRI at 6 arcsec resolution. These 60 targets were homogeneously selected from LoTSS (selecting the sample based on low-frequency radio emission also enabled us to minimise contamination from Doppler boosting) in order to reveal the inner structures of the target objects, confirm their radio morphological classification and further evaluate their properties in a statistically meaningful manner. In this paper we present the results of the pilot study that used the recent VLA radio continuum observations. The layout of this paper is as follows. Descriptions of the sample and the data used in this work are given in Section 2. The primary analysis and key results are given in Section 3. We finish with a summary and discussion of future work in Section 4. Throughout the paper we use the most recent Planck cosmology [29]: H 0 = 67.7 km s −1 Mpc −1 , Ω m =0.308 and Ω Λ =0.692. Spectral index α is defined in the sense S ∼ ν −α .

The Sample and Optical Data
Our quasar sample is drawn from the SDSS quasar catalogue DR14Q [24], which includes all SDSS-IV/the extended Baryon Oscillation Spectroscopic Survey [eBOSS;] [30] objects that were spectroscopically targeted as quasar candidates and that are confirmed as quasars via a new automated procedure combined with a partial visual inspection of spectra. The DRQ14 provides a compilation of all the spectroscopically-confirmed quasars identified in the course of any of the SDSS iterations and released as part of the SDSS. The SDSS, Wide-field Survey Explorer [WISE;] [31] and the Palomar Transient Factory [PTF;] [32,33] imaging data were used to select target quasars. The CORE sample is selected homogeneously to targets quasars at z > 0.9 based on the method that utilises optical and WISE-optical colours. There are also quasar candidates that are selected based on the photometric variability measured from the PTF. We refer the reader to [34] for the details of the SDSS-DR14Q sample selection. We cross-matched the value-added catalogue of the LoTSS-DR1 [35] over the HETDEX region with the SDSS DR14 quasar catalogue [34] using a 3 matching radius in order to select quasars for visual inspection. We then generated images for visual inspection, consisting of the Panoramic Survey Telescope and Rapid Response System [Pan-STARRS;] [36] image with LOFAR and the Faint Images of the Radio Sky at Twenty Centimeters [FIRST;] [37] contours overlaid. This process generated ∼6000 sources, all of which were inspected visually in order to classify sources as FRI and FRII. The requirement to classify a source as an FRI quasar is that it is brighter in the centre (not including the bright core component) and gets fainter as the distance from the centre increases. We identified 60 FRI quasar candidates of which the majority were convincing FRI quasars, together with a few more complex objects. It is important to emphasise that the combination of LOFAR and FIRST information is necessary for this efficient selection as FIRST picks out the central flat spectrum and compact components, whereas LOFAR is sensitive to the large-scale steep spectrum structures. Our classification has also been confirmed for some of our selected objects by a method which uses the surface brightness profile and the distance of extended structures of a radio source to its optical counterpart [LoMorph;] [11], which was developed for morphological classification of LoTSS sources. Only 519 quasars meet the angular size and flux requirements of LoMorph [11] and 109 of them were cleanly classified as FRI or FRII by this method (32 FRI quasars and 77 FRII quasars). Only 15 objects out of the selected 60 FRI quasar candidates met the criteria that was used by [11] (i.e., there are 15 common objects between our target sample and their cleanly classified sample). All of these 15 sources were classified as FRI or small FRI by LoMorph. In this paper we focus on the selected 60 FRI quasar candidates and compare their properties with quasars that are cleanly classified as FRI or FRII by [11] but not included in our VLA snapshot survey.
The DR14Q quasar catalogue contains i band absolute magnitudes (k-corrected to z = 2), SDSS magnitudes and fluxes at u,g,r,i,z bands together with various properties derived from these. These measurements were used in the analyses presented in this work. Additionally, we obtained extinction-corrected i band absolute magnitudes k-corrected to z = 0 by converting the i band absolute magnitudes k-corrected to z = 2 using the conversion given by [38] in the same way as done by [3]. Whenever available, we used black hole masses estimated using the MgII and CIV emission line widths (There are complications in the estimation of virial black hole masses using MgII and CIV lines that should be noted such as different line measurement techniques, the treatment of the narrow line component in both lines and the choice for the method of subtracting the continuum emission underneath the line. Addressing these caveats would go beyond the scope of this work), optical bolometric luminosities which are estimated using bolometric corrections by [38] using one of the quasar luminosities at 1350, 3000, 5100 Å at the source redshift, and Eddington ratios published by [39]. Otherwise, we used estimates of the same quasar properties from [40].

LoTSS DRI Image Data
The HETDEX Spring field (right ascension 10h 45 m 00 s to 15 h 30 m 00 s and declination 45 • 00 00 to 57 • 00 00 ) was observed with LOFAR as part of LoTSS DR1. This field was targeted as it is a large contiguous area at high elevation for LOFAR, whilst having a large overlap with the SDSS imaging and spectroscopic data. The creation of the radio images is described by [23]. Radio flux densities at 144 MHz for all SDSS quasars in our sample were directly taken from the value-added LOFAR/HETDEX catalogue [35]. These LoTSS radio flux densities at 144 MHz were used to derive the rest frame radio luminosities (using a spectral index of 0.7) of our sample sources which were used as an approximation of the kinetic energy output/jet power.
On the left panel of Figure 1 we show the distribution of 144-MHz luminosity of the selected quasars as a function of their redshifts. Different symbols are used to show quasars detected in the LoTSS DR1 (blue circles), sources showing extended (resolved in the LOFAR maps) emission (orange stars; these were obtained using the value-added catalogue) and our FRI quasar candidates (purple stars). The blue solid line shows the traditional FRI/FRII break at 144 MHz. The radio-loudness parameter has traditionally been used to separate RQ and RL AGN, although we are not using this parameter we show the distribution of this parameter as a reference for other studies. The radio loudness distribution of quasars (colour coded the same way as the left panel to show different classes) is seen in the right panel of Figure 1. Here we define the radio loudness parameter for our quasar sample using the ratio of L 144 to i band luminosity, as we did in [3]. As can be clearly seen in both panels of Figure 1 almost all of the FRI quasar candidates are above the FRI/FRII break (left panel) and they have a range of radio loudness parameter (covering a range of 0−5.5). The traditional radio-loud/radio-quiet dichotomy would predict that quasars with R < 2 (at 144 MHz, [3,41]) are radio-quiet but we know that all of the objects classified as FRI or FRII in the sample have radio emission dominated by AGN.

VLA Observations
The 60 FRI quasar candidates were observed with the VLA on 29 September 2019 in the A-configuration (giving a resolution of around 1.3 arcsec). Observations were made in the broad-band L-band system in 8-bit mode (comprising 16 contiguous subbands of 64 MHz to cover the entire 1-2 GHz of the L-band receiver). Since we only required snapshot observations, exposure times were 10 min for each target (scheduled as two identical 5-h sessions, which allowed to improve the snapshot uv coverage by observing each source for 2 × 5 min ), reaching background rms levels of ∼30 µJy beam −1 in the final combined images.The observations were targeted at the SDSS optical counterpart position for each source.
The quasar 3C 286 was used as the flux and bandpass calibrator, and since the large area of sky covered by our sources (the HETDEX-Spring region is ∼424 sq deg) five phase calibrators were used in total to correct for ionospheric variations over the course of ten hour observation time. The two epochs of observations were reduced separately. Prior to calibration the AOFlagger software [42] was used on each data set to automatically flag radio-frequency interference (RFI). The data were then reduced using the CASA VLA pipeline version 5.3.1 using CASA version 5.1.10 [43]. Various calibration tables were inspected to check the quality of calibration, and baselines displaying residual RFI or erratic phase variations were flagged manually in CASA. Images were produced by combining both epochs of observations in the uv-plane using the CASA image reconstruction technique CLEAN [making use of the Clarke algorithm;] [44].
Flux densities of quasars not showing any extended structure (i.e., point sources) were estimated by Gaussian fitting in CASA. For sources with extended emission we applied a mask to select pixels at least three times the off-source noise and measured the total flux (for each source separately). Errors on these flux densities were calculated by taking into account the noise and calibration errors [45]. In Figure 2 we show the VLA luminosity distribution of our sample sources as a function of their redshifts.

Classification Statistics of the Targets Using VLA Observations
We initially used the final VLA images of the target sources to have better constraints on their radio morphology. The VLA images of our target sources along with other radio data (LOFAR, FIRST, and VLASS) are shown in Figure A1 in the Appendix A. We visually inspected these images to classify their radio morphology and these results can be found in Table 1. Due to the limitation of the snapshot observations along with the effect of the source orientation to the line of sight the inner structures of some target quasars were not clearly revealed. The decision on a classification was made based on collective agreement from individual classifications by six authors. If a source was classified in the same category by more than three classifiers it is considered to be a genuine FRI or FRII. If three classifiers agreed upon a category then this source was labelled as uncertain FRI/FRII class. There are sources with no resolved components in the VLA images which were classified as core dominated sources and the remaining form the uncertain class. Example images for each morphological class can be seen in Figure 3. A summary of our target sources along with VLA measurements and their classifications are given in Table A1

Nuclear Properties of Quasars Associated with Different Radio Morphologies
Although FRIs make up almost 5 percent of the LoTSS clean AGN sample [11], among the extended quasars the fraction of FRIs is much lower: we found 32 FRIs (30 per cent) and 77 FRIIs (70 per cent). In this study out of the 60 quasars (we selected using SDSS DR14 and LoTSS DR1) the follow-up VLA snapshot observations allowed us to confirm 15 of these to be genuine FRIs. There are 13 quasars that show extended kpc structures similar to FRIs but there are bright locations in these jets which might be a hot spot or a jet knot so these were classified as uncertain FRIs. There are also three clear FRIIs and six uncertain FRIIs. Quasars which were not resolved by the VLA (A-config) beam were classified as core dominated. There are 14 objects which do not fit in any of these classifications so these are labelled as uncertain. The sample presented here form the largest sample of FRI quasars (and that of the possible candidates) to date. Moreover, our sample is not limited to powerful, bright objects in the sky as they are selected from the most sensitive low-frequency radio survey (going down to L 150 MHz >10 25 W Hz −1 ). Taking into account the size of our sample and its selection these statistics suggest that quasars with FRI radio morphology are not common in the universe. These results reinforces the idea that it is possible for a radiatively efficient system to form FRI jets but this probably requires certain conditions.
It is important to compare the properties of quasars associated with different radio morphology (FRI and FRII). We know that by definition quasars have radiatively efficient accretion so the accretion mode is not a parameter that plays a role here. The properties of the central engine in AGN might be important for understanding the rarity of FRI quasars. In Figure 4 we show the histogram of quasar black hole masses, where available. In the left panel we show the black hole mass distribution of our sample observed as part of the VLA snapshot survey. From this panel only, we do not see any apparent difference in the range of black hole masses of quasars associated with different radio morphology, though the sample size is small. In the right panel of Figure 4 we only show FRI, FRII and uncertain FRI quasars as well as FRI and FRII quasars from LoTSS sample classified by LoMorph. Here we see a similar result to the left panel: there is no clear difference between the black hole masses of FRI and FRII quasars.  The black hole mass histogram of quasars that were classified as FRI, FRII and uncertain FRI as part of our pilot study along with FRI and FRII quasars classified by [11] using LoMorph. For a clean representation, we show FRII quasars and uncertain FRIIs using the same symbol.
Although our sample size is small we used (k-sample) Anderson-Darling (AD) test to assess whether underlying probability distributions of black hole masses of FRI and FRII quasars differ from each other ( Table 2). As can be seen from the results we cannot reject the null hypothesis (that two samples are drawn from the same population) at 5% level for the black hole masses of quasars classified as FRI or FRII using our VLA observations. In terms of comparing FRI quasars (including uncertain FRIs) with FRII quasars similarly we fail to reject the null hypothesis at 5% level. The underlying probability distribution of the black hole masses of all FRI quasars (FRIs, uncertain FRIs and LoMorph FRIs) and that of all FRII quasars (FRIIs and LoMorph FRIIs) do not differ from each other at 5% level. AGN jet power is another key parameter (particularly the jet interaction with AGN environment) that defines the structural difference between FRI and FRII sources. The low-frequency radio luminosity is not only related to the lifetime-averaged jet kinetic power but also to the age of the source and to the properties of the external environment, e.g., [46]. As it is not straightforward to estimate the lifetime-averaged jet power we initially evaluate the low-frequency radio luminosity (which can be used to approximate the instantaneous jet power with the caveat of uncertainties) of our sources. The low-frequency radio luminosity distribution of our sample sources is seen in the left panel of Figure 5. As can been seen in this panel, core-dominated quasars have a range of radio luminosities (10 25 <L 150 / W Hz −1 <10 27 ) whereas FRI quasars tend to have lower radio luminosities and FRII quasars have relatively higher radio luminosities. This is more apparent in the right panel of Figure 5 where we have better statistics in terms of classification, albeit not all quasars in this panel have follow-up observations with the VLA. The results of the Anderson-Darling hypothesis testing also support this picture (Table 3). We can reject the null hypothesis at 5 per cent level that the radio luminosities of FRI and FRII quasars are drawn from the same distribution. We obtain similar results when we assess the radio luminosities of FRIs (including uncertain FRIs) and FRIIs (at 5 per cent significance level). Taking into account all FRIs (FRIs, uncertain FRIs and LoMorph FRI quasars) and FRIIs (FRIIs and LoMorph FRIIs) do not indicate any different result: the low-frequency radio luminosities (an approximation for AGN jet power) of FRI and FRII quasars are not drawn from the same distribution (at 5% significance level). These results do not change when we use the jet powers of our sources estimated using a relation given by [47] instead of their low-frequency radio luminosity (also see below). The low-frequency radio luminosity histogram of quasars that were classified as FRI, FRII and uncertain FRI as part of our pilot study along with FRI and FRII quasars classified by LoMorph. For a clean representation, we show FRII quasars and uncertain FRIIs using the same symbol. Table 3. Results of the Anderson-Darling test that has been used to assess the low-frequency radio luminosity (as an indicator for AGN jet power) of FRI and FRII quasars. Radio-loud quasars, as radiatively efficient systems, produce both radiative and kinetic output as the matter is channelled onto the supermassive black hole via the accretion disk. In order to compare the radiative power of FRI and FRII quasars with their jet powers we estimate the total bulk kinetic power of quasars using the low-frequency radio luminosity following the relation given by [47]. It is worthwhile to remind the reader that jets are capable of producing radio emission at a very wide range of luminosities, with the radio luminosity depending on the jet power, the environment and the source age [48] and the jet power−radio luminosity relation might be expected to be different for radio populations classified as FRI and FRII [49]. In Figure 6 we compare the radiative (i.e., the bolometric luminosity) and the kinetic (i.e., jet power) products of quasars, colour coded by the Eddington ratios (Eddington ratio here is the ratio of the AGN bolometric luminosity to the Eddington luminosity.). Qualitatively the Eddington ratios of quasars with different radio morphology do not show any significant difference. The results of AD test (Table 4) support this qualitative result: Eddington ratios of FRI and FRII quasars are not different from each other.

Compared
In the left panel we show only FRI, uncertain FRI and core-dominated quasars from our sample with VLA follow-up observations along with FRIs quasars from LoTSS. All cleanly classified FRI quasars, the majority of uncertain FRIs and almost all of FRIs from the LoTSS exhibit a narrow range of jet powers whereas their accretion powers have a wider range of values (37< log 10 L bol <40.5). This picture emerges more clearly in the right panel of Figure 6 where we include FRII quasars from both our sample and the LoTSS sample.  Figure 6. The distribution of total bulk kinetic power as function of radiative power of our target sources: FRI, uncertain FRI, core-dominated quasars, respectively. We also show FRI quasars classified by LoMorph left) [11]. The right panel shows the same distribution as the right-hand side but we further include quasars with uncertain morphological classification as well as LoMorph-FRII quasars. For a clean representation, we show FRII quasars and uncertain FRIIs using the same symbol. Pink shaded regions just highlight the narrow range of jet powers of FRIs across a wider range of their radiative powers. Colour coding shows the Eddington ratios (i.e., Bolometric luminosity/Eddington luminosity).

Host Galaxy Colours
Recent studies focusing on SDSS quasars by [50] using FIRST data and [51] using LoTSS data investigated quasars with red and blue optical colours and found that quasars with red optical colours present enhanced radio emission from AGN. In their morphological classification most of the red quasars show compact features and some have FRI-like morphology. It is therefore interesting to evaluate the optical colours of our sample sources to see whether there is any trend in the optical colours of FRI quasars and how these compare to FRII quasars. Figure  7 shows the optical colour distribution of quasars from the SDSS DR16 (this is most recent version which was not available at the time a sample for this project was selected), the quasar sample classified by LoMorph and our target sources from the VLA survey. Here the optical colours have been defined in the same way as done by [51] using SDSS g and i magnitudes (corrected for Galactic extinction using the extinction values provided in the DR16Q catalogue [52]). Here in Figure 7 red quasars are shown as red points and the control sample as defined by [51] are shown as blue and the remaining quasars are shown as gray points (both quasars that are shown as gray and blue points form the blue quasar sample by the aforementioned work). It is clear from this figure that: (i) The majority of FRI quasars do not have red optical colours and (ii) there is no clear difference between the optical colours of FRI and FRII quasars. Figure 7. The optical colour (SDSS g − i) distribution of quasars from the DR16Q catalogue, our target sample and the quasars sample with clean classification by LoMorph. Red quasar are quasars with red optical colours, control quasar indicate the control sample (the 10% of quasars around the median of the g − i distribution) and pure blue colour quasars are given as blue quasar. The optical colour classification is based on [50,51]. The dashed lines show the range of redshift within which the optical colours are reliable as described by [50,51]. Symbols of our quasars sample are as above. For a clean representation, we show FRII quasars and uncertain FRIIs using the same symbol.

Discussion
Based on the available data, we have evaluated various parameters with respect to the nuclear engine in the interest of testing whether there are differences between FRI and FRII quasars. As shown above, FRI and FRII quasars have a similar range of black hole masses ( Figure 4 and Table 2) and Eddington ratios ( Figure 6 and Table 4). There is a tendency for FRII quasars to have relatively higher jet powers than FRI quasars ( Figure 5). Both FRI and FRII quasars can produce a variety of accretion power. For a given range of accretion powers FRII quasars can produce a broad range of jet powers (∼10 36.6 −10 40 W). FRI quasars, on the other hand, exhibit a quite narrow range of jet powers (∼10 37.5 −10 38.5 W). It is also possible that based on our selection criteria we are limited with less luminous objects (with the full survey we will be able to cover a wider range of luminosities).
The formation of different radio morphologies has been argued to be due to either the influence of external environmental factors on the jet structure, e.g., [53,54] or by parameters associated with the jet production mechanism [55]. It is thought that both FRI and FRII sources initially expand supersonically and FRII objects then might evolve to FRIs once their jet disrupted, e.g., [56]. There is no reason not to think that both effects can contribute to this dichotomy. In addition this, we know that AGN bolometric luminosity is a function of accretion rate, which is a function of black hole mass and radiative efficiency. For a standard accretion disk the radiative efficiency is expected to depend on black hole spin. Taking into account black holes masses of FRI and FRII quasars are not different for our sample and they are capable of producing a similar range of accretion powers quasars showing jet structures in the sample (independent of their radio morphology) should have spin values above a certain threshold. On the other hand, the radio luminosity of AGN is a complex function of AGN jet power (along with radiative losses, time, cosmic epoch and finally AGN environment) which depends on black hole spin and magnetic flux. As indicated above, quasars both with FRI and FRII morphology should have a black hole spin above a certain threshold in order to produce the relativistic jets but the intrinsic scatter on black hole spins and the strength of magnetic fields might be the key factors in determination of the low jet powers in FRI quasars, e.g., [57]. Quasars with jets have been found be clustered more strongly than their radio-quiet counterparts across a wide redshift range, e.g., [58]. It has not been possible to compare the clustering/environments of FRI and FRII quasars to date but this is likely to change with the current and upcoming surveys. Nevertheless, previous studies on the radio galaxy environments found that FRI radio galaxies reside in denser (host scale) environments than FRII radio galaxies [53,61,66], but see [59]. Therefore, we think that a combination of the scatter on black hole spin and magnetic flux, and the influence of external environmental factors define the radio morphology. Additionally, it is worth noting that based on the unification model quasar jets are seen at small angles to the line of sight [Θ < 45 deg;] [60] so it is possible that the projection effects coupled with the sensitivity limited observations restrict the number of sources we can classify.
We do not see any particular difference in the host galaxy colours probed using g − i (extinction corrected) magnitudes from the SDSS data (Figure 7): the majority of FRI and FRII quasars have mainly blue optical colours. All these results indicate that another key parameter is important for quasars to produce FRI type of jets. A potential parameter for this is the galaxy environment. The current paradigm suggests that [61][62][63] jets decelerate by entrainment of matter from relativistic speeds on pc scales to sub-relativistic speeds on kpc scales to form FRI jets. It is still debated whether entrainment is defined by the power of the jets and the galaxy external environment or it is due to a more fundamental difference. With the current data available to us we are unable to comment on the environment properties of quasars as our target sources are mainly at moderate and high redshifts.

Conclusions
We homogeneously selected 60 FRI quasar candidates using the LoTSS DR1 data and further observed these target sources at higher angular resolution using the VLA with the aim of revealing the inner structures of the AGN jets. Our sensitivity-limited snapshot survey allowed us to validate the FRI classification of these quasars and increase the sample of known FRI quasars to date. We utilised the available SDSS, LoTSS and value-added data for understanding this population further and our results are as follows: • 25 per cent of the selected target quasars have been confirmed to be FRIs and around 25 per cent show extended structures similar to FRI jets but with the current observations it is not clear whether the bright jet parts are jet knots of an FRI jet or hot spots in FRII jets, • Although our sample size is small we do not see any difference in the range of black hole masses and Eddington ratios of FRI and FRII quasars. On the other hand, FRI quasars have relatively low jet powers compared to FRII quasars in the samples explored in this study, consistent with the expectation from low-z radio galaxies, • Investigation of optical galaxy colours of quasars showed that FRI and FRII quasars mainly have blue optical colours. Based on the sample studied here we conclude that quasars, as the most powerful, radiatively efficient systems among the AGN population are more likely to produce FRII type jets. For a quasar to form FRI jets should have low jet power (which is expected to be characterized by black hole spin and magnetic flux) and reside in a significantly dense environment.
Future studies using the VLA polarisation data of our target sources as well as e-MERLIN observations of selected quasars will provide more insights into the nature of these objects. As the ongoing/near future radio surveys progress the sample of FRI quasars will increase and eventually with the SKA (in conjunction with the optical surveys) we will advance our understanding of this rare population.