The warm molecular hydrogen of PKS B1718–649:

The warm molecular hydrogen of PKS B1718–649:

feeding a newly born radio AGN.
F. M. Maccagni Kapteyn Astronomical Institute, University of Groningen, Postbus 800, 9700 AV Groningen, The Netherlands ASTRON, Netherlands Institute for Radio Astronomy, Postbus 2, 7990 AA, Dwingeloo, The Netherlands    F. Santoro Kapteyn Astronomical Institute, University of Groningen, Postbus 800, 9700 AV Groningen, The Netherlands ASTRON, Netherlands Institute for Radio Astronomy, Postbus 2, 7990 AA, Dwingeloo, The Netherlands    R. Morganti Kapteyn Astronomical Institute, University of Groningen, Postbus 800, 9700 AV Groningen, The Netherlands ASTRON, Netherlands Institute for Radio Astronomy, Postbus 2, 7990 AA, Dwingeloo, The Netherlands    T. A. Oosterloo Kapteyn Astronomical Institute, University of Groningen, Postbus 800, 9700 AV Groningen, The Netherlands ASTRON, Netherlands Institute for Radio Astronomy, Postbus 2, 7990 AA, Dwingeloo, The Netherlands    J. B. R. Oonk ASTRON, Netherlands Institute for Radio Astronomy, Postbus 2, 7990 AA, Dwingeloo, The Netherlands Leiden Observatory, Leiden University, Postbus 9513, 2300 RA Leiden, the Netherlands    B. H. C. Emonts Centro de Astrobiología (INTA-CSIC), Ctra de Torrejón a Ajalvir, km 4, 28850 Torrejón de Ardoz, Madrid, Spain
Key Words.:
PKS B1718–649 – compact radio sources – active nuclei – neutral hydrogen – molecular hydrogen – ISM

We present new SINFONI VLT observations of molecular hydrogen (H) in the central regions ( kpc) of the youngest and closest radio source PKS B1718–649. We study the distribution of the H traced by the 1-0 S(1) ro-vibrational line, revealing a double disk structure with the kinematics of both disks characterised by rotation. An outer disk ( pc) is aligned with other components of the galaxy (atomic hydrogen, stars, dust), while the inner disk ( pc) is perpendicular to it and is polar with respect to the stellar distribution. However, in the innermost pc, the data show the presence of H gas redshifted with respect to the rotating inner disk (kms) which may trace gas falling into the super massive black hole associated with the central radio source. Along the same line of sight, earlier observations had shown the presence in the central regions of PKS B1718–649 of clouds of atomic hydrogen with similar unsettled kinematics. The range of velocities and mass of these unsettled clouds of H I and H suggest they may be actively contributing in fuelling the central newly-born radio source.

1 Introduction

Active Galactic Nuclei (AGN) are associated with the accretion of material onto the central super-massive black hole (SMBH) of galaxies. The gas surrounding the SMBH must lose angular momentum in order to fall into it so it can trigger and feed an active nucleus. Nevertheless, direct observational evidence of this process is still limited. Statistically, galaxies which have undergone a merger or an interaction event have higher probability to host an AGN (Ellison et al., 2008; Ramos Almeida et al., 2012; Hwang et al., 2012; Sabater et al., 2013). However, in several objects with signatures of past mergers or accretion, the time-scales associated with these phenomena can be much longer than the age of the AGN (Emonts et al., 2006; Tadhunter, 2008; Schawinski et al., 2010; Struve & Conway, 2012; Maccagni et al., 2014) suggesting that in these galaxies the link between mergers/accretion and AGN is, at most, indirect and other processes must occur to trigger the nuclear activity. Slow secular processes may help the gas lose angular momentum on short timescales ( years) and form a dense gas core in the central pc (Kormendy & Kennicutt, 2004; Wada, 2003; Combes, 2004, 2011). However, it is not clear whether these phenomena are efficient in the very innermost regions near the AGN (Athanassoula et al., 2005; Begelman & Shlosman, 2009). Thus, other processes taking place on small spatial and temporal scales are expected to be responsible for the direct fuelling onto the AGN (Wada & Tomisaka, 2004; King & Pringle, 2007; Hopkins & Quataert, 2010). Numerical simulations suggest that gravitational and thermal instabilities induce chaotic collisions in the interstellar medium ISM surrounding the SMBH (Soker et al., 2009; Gaspari et al., 2013; King et al., 2008; Nayakshin & Zubovas, 2012; King & Nixon, 2015). This causes small clouds or filaments of gas to lose angular momentum and begin a series of small-scale, randomly oriented accretion events, which then trigger the AGN. In this scenario, the gas deviating from regular rotation is responsible for the chaotic infall of clouds and, consequently, for the accretion onto of the AGN (Gaspari et al., 2015).

High spatial resolution observations tracing in particular the cold gas in the innermost regions of AGN are needed to investigate these hypothesis. Different type of AGN, i.e. Seyfert galaxies (Gallimore et al., 1999; Mundell et al., 2003; Hicks et al., 2009, 2013; Combes et al., 2014; Mezcua et al., 2015), low-ionisation nuclear emission region galaxies (LINER; García-Burillo et al., 2005; Müller-Sánchez et al., 2013), and radio galaxies (Neumayer et al., 2007; Dasyra & Combes, 2011; Guillard et al., 2012; Morganti et al., 2013) are rich in molecular (H) and atomic hydrogen (H I) which may represent the fuel reservoir for the nuclear activity. Indeed, the kinematics of at least part of this gas often appears to be unsettled with respect to the regular rotation of the galaxy, suggesting a strong interplay between the nuclear activity and the surrounding environment. On the one hand, it is likely that plasma ejected by the radio source perturbs the neutral and molecular hydrogen (Neumayer et al., 2007; Hicks et al., 2009; Dasyra & Combes, 2011; Guillard et al., 2012; Müller-Sánchez et al., 2013; Mezcua et al., 2015). On the other hand, it is also possible that this reflects the presence of processes as those described above that can cause the gas to stream towards the SMBH and trigger the nuclear activity (Hopkins & Quataert, 2010; Combes et al., 2014).

Young radio sources in the first stages of their activity (Murgia, 2003; Fanti, 2009) are the best candidates to study the relation between the kinematics of the cold gas and the triggering of the AGN. They are often embedded in a dense gaseous environment where the fuelling of the AGN has just begun and is likely to be still on-going. Also, among all radio AGN, these sources show relatively often neutral and molecular gas with unsettled kinematics in proximity of the core (Emonts et al., 2010; Geréb et al., 2014, 2015; Curran et al., 2013; Guillard et al., 2014; Allison et al., 2015).

PKS B1718–649 is a compact radio source ( pc) with an optically classified LINER AGN (Filippenko, 1985) at a distance111; Mpc, 1 arcsec = 0.294 kpc; where CDM cosmology is assumed, kmsMpc, , . of 62 Mpc. The estimated age of the radio activity is 10 years (Tingay et al., 1997; Giroletti & Polatidis, 2009). PKS B1718–649 is morphologically classified as an S0-SABb early-type galaxy embedded in a disk of neutral hydrogen (see Fig. 1) which shows regular rotation out to large radii (23 kpc). Given the long time scale for such a regular disk to form, this excludes a merger or a disruptive event being directly responsible for the recent triggering of the central radio source. The accretion onto the SMBH, and the fuelling of the radio activity, could find its origin in a small-scale phenomenon. The detection, in H I absorption, of two separate clouds with kinematics deviating from the rotation of the disk, suggests that a population of clouds may be contributing to feed the AGN in the centre (Maccagni et al., 2014). The regions close to the radio source have been indirectly probed by the study of the variability of the radio continuum emission (Tingay et al., 2015), which has been attributed to changes in the free-free absorption due to a clumpy circum-nuclear medium around the radio source. The presence of such a clumpy medium has been also suggested by optical spectroscopic observations (Filippenko, 1985). The available information on PKS B1718–649 suggests that the origin of its newly born radio activity may be found in the kinematics of its circum-nuclear medium.

Figure 1: -band optical image of PKS B1718–649, overlaid with the column density contours (black) of the neutral gas. The contour levels range between cm and cm, with steps of cm . The unresolved continuum radio source is marked in white. The H I disk has the shape of an incomplete ring with asymmetries in the N-W and in the S of the disk (Maccagni et al., 2014).
Parameter Value
Field of view (2.372.37 kpc)
Pixel size (37 pc)
Spectral resolution (at 2.1 m) kms(=4000)
Spectral sampling (at 2.5 m) kms
Seeing (run 1, run 2, run 3) (131 pc); (735 pc); (154 pc)
Table 1: SINFONI observation specifications

Integral Field Unit (IFU) instruments allow us to analyse the spatial distribution and kinematics of the ISM in the innermost regions of low-redshift AGN. Here we present the results obtained for PKS B1718–649 using the Spectrograph for INtegral Field Observations in the Near Infrared (SINFONI) on the VLT. The detection of the molecular hydrogen traced by its ro-vibrational states ( K, H 1-0 S(0,1,2,3)) allows us to study its distribution and kinematics and to provide interesting insights on the role of the H in relation to the fuelling of the central radio source.

2 Observations and data reduction

We observed the inner kpc region of PKS B1718–649 in the -band (m), using SINFONI (Eisenhauer et al., 2003) mounted on the Very Large Telescope (VLT) UT4.

The observations were performed under seeing-limited conditions during three different nights (May 18-26-29, 2014) in period 93A. The spectral resolution is and the plate scale is pixel, yielding a field of view of . The Full Width Half Maximum (FWHM) of the sky lines is Å, with a spectral sampling of Å pixel. Due to bad seeing conditions, we exclude the observations performed on May 26th.

We perform the data reduction using the official ESO REFLEX workflow for the SINFONI pipeline (version ) and the standard calibration frames provided by ESO. The workflow allows us to derive and apply the corrections for dark subtraction, flat fielding, detector linearity, geometrical distortion, and wavelength calibration to each object and sky frame. Following Davies (2007), we subtract the sky from the data cubes of the two observing blocks. The typical error on the wavelength calibration is Å (kms).

Through IDL routines implemented by Piqueras López (2014), we calibrate the flux of each cube. First, we obtain the atmospheric transmission curve, extracting the spectra of the standard stars with an aperture of of the best 2D Gaussian fit of a collapsed cube. Then we normalize these spectra using a black-body profile at the temperature corresponding to the spectral type of the observed stars (as tabulated in the 2MASS catalogue, see Skrutskie et al. (2006)), using Table 5 in Pecaut & Mamajek (2013). We model the stellar hydrogen Br absorption line at m with a Lorentzian profile to determine the sensitivity function for the atmospheric transmission (Bedregal et al., 2009). We convert the star spectra from counts to physical units with a conversion factor extracted using the tabulated magnitudes in the 2MASS catalogue. We obtain the full-calibrated data cube dividing each spectrum by the sensitivity function and multiplying it by the conversion factor. The typical uncertainty for the flux calibration is . We combine the two cubes of the single observing runs into the final one by spatially matching the peaks of the emission of the galaxy. The spatial resolution of the final data cube is equal to the FWHM of a 2D Gaussian we fit to the central region of the collapsed cube: arcseconds.

In the final cube, we determine whether a line of the ro-vibrational states of the H is detected considering the spectra extracted over regions of size equal to the spatial resolution of the observations. In Fig. 3 we show the spectra extracted in five different regions of the H distribution, approximately in the north (), the south (), the west (), the east () and in the centre () (see Fig. 2). In all five regions, we detect the H 1-0 S(1) line, the brightest in the spectrum, at m, as well as the H 1-0 S(3) line, while we do not detect the H 1-0 S(0,2), and the higher excited states of molecular hydrogen, H 2-1 S(1,3). In this case, we determine the upper limits assuming a FWHM of the line equal to the one of the H 1-0 S(1) line. Br and [Si VI], tracers of high-excitation ionized gas, are also expected in the same spectral range, but lie below the detection limit of these observations. In Table 2, we list the fluxes and upper limits of the lines for the five different regions.

We focus on the H 1-0 S(1) line to determine the distribution and the kinematics of the molecular hydrogen and we use the integrated fluxes or the upper limits of the H 1-0 S(0,1,2,3) lines to determine the temperature of the molecular hydrogen and its mass. We derive the distribution and kinematics of the H 1-0 S(1) emission line using two independent methods, which provide consistent results. In the first method, we spatially smooth the cube with a Gaussian with FWHM 0.5 arcsec and we fit a single Gaussian component to the H 1-0 S(1) line in each pixel of the field of view. We build a mask of the regions of pure line emission, selecting the pixels where the H 1-0 S(1) line is detected with signal-to-noise ratio of . We extend the mask to also consider regions neighbouring the ones where the fit is successful. Following this, for each pixel within the mask, we extract the spectrum and we fit the H  S(1,3) lines with a single gaussian component. We derive the intensity and velocity fields of the H 1-0 S(1) line only in the regions where the two lines are detected.

Line [] R R R R R
H 1-0 S(3)   [m] 20.21.29 6.881.17 11.41.52 14.11.62 8.353.08
[Si VI]   [m] <3.88 <3.52 <4.55 <4.86 <8.88
H 1-0 S(2)   [m] <5.85 <6.59 <5.51 <7.33 <10.3
H 2-1 S(3)   [m] <6.41 <3.64 <4.89 <5.90 <8.04
H 1-0 S(1)   [m] 24.20.737 11.61.62 19.10.659 16.10.835 44.66.15
Br   [m] <4.01 <4.01 <4.56 <6.19 <4.55
H 1-0 S(0)   [m] <6.39 <3.61 <5.63 <7.54 <10.2
H 2-1 S(1)   [m] <4.73 <6.81 <4.59 <9.09 <10.4
222 The fluxes are given in units of ergscm. The upper limits indicate the 3- noise level, measured in the wavelength ranges where we expect to detect the lines. The spectra are extracted in five regions, R, R, R, R, R of the field of view shown in Fig. 2. Their sizes correspond to the spatial resolution of the observations ().
Table 2: Line fluxes of the molecular hydrogen gas in five regions of the SINFONI field of view.
Figure 2: Distribution of molecular hydrogen (black contours) overlaid on to the Hubble Space Telescope WFPC2 image. A dust lane is visible, oriented in the north-south direction. The grey crosses mark the regions where we extracted the spectra to measure the temperature of the H. The grey dashed square marks the SINFONI field of view.

In the second method, we build a cube free of emission-line signal masking out channel-by-channel the regions where the H 1-0 S(1) line is detected above the 2.5- level. Next, we smooth the edges of the masked regions to completely exclude any residual emission-line signal. From this cube, we determine the template of the stellar continuum spectrum. We subtract this stellar spectrum from every pixel where the H 1-0 S(1) line is detected and we obtain a data cube of pure emission-line spectra. We determine the distribution of the H as the zeroth moment map of this cube, summing along the velocity axis all emission above the 3- level in at least two consecutive pixels. The velocity field corresponds to the first moment map and is centred on the systemic velocity of the H I disk, kms (Maccagni et al., 2014).

This method is less conservative than the first, but does not rely on the quality of the Gaussian fitting, which may bias the characterisation of the morphology of the H 1-0 S(1) emission. Interestingly, the total intensity and velocity field determined from the two methods are very similar. In the next section, we use the results of the second method to analyse the kinematics of the molecular hydrogen in PKS B1718–649.

The Gaussian fitting of the H 1-0 S(1) line describes the line across the field of view well, except for in the central , in proximity of the AGN. This is the only region where the fit with a single Gaussian component leaves substantial residuals above of the noise (see Fig. 3 right panels). The profile is also clearly more asymmetric and the velocity dispersion of the profile is higher than in the other regions of the field of view. This suggests that in the centre, more than one component is needed to fully characterise the kinematics of the H. A further analysis of the morphology of the H 1-0 S(1) emission in the central pc of PKS B1718–649 is given in Section 3.2.

Figure 3: Spectra extracted from the five regions (R, R, R, R, R) illustrated in Fig. 2. The right panels show is a zoom-in on the H 1-0 S(1) line, centred at the systemic velocity of PKS B1718–649. The dashed line shows the fit with a single Gaussian component, while the dotted line shows the residuals. In the bottom panel, solid lines mark the locations of the H 1-0 S(0,1,2,3) lines, while dashed and dotted lines show the H 2-1 S(1,3) and [Si VI] and Br lines, respectively.
Figure 4: (a): Intensity map of the H 1-0 S(1) line in the inner kpc of PKS B1718–649. The position angle of the radio source is shown in dashed green (PA). Intensity contours at 1.5, 3, 5, 7, 9, 12, 15- are shown in black, starting from 1.5ergscm. The isophotes in grey show the distribution of the stellar component. (b): Velocity field of the H 1-0 S(1) line with the contours of the intensity map overlaid, the position angle of the radio source is shown in black. Velocities are given relative to the systemic value, kms. The regions marked by crosses are those where we extracted the profiles of Fig. 5, (see the text for further details).

3 Results

3.1 Distribution and kinematics of the molecular hydrogen

The intensity map and velocity field of the H 1-0 S(1) line emission in the central regions of PKS B1718–649 are shown in Figs 4 (a) and (b). At radii pc, the H is assembled in a disk aligned in the N-S direction, the same direction as the H I disk, H and the dust lanes at larger radii ( kpc) (Keel & Windhorst, 1991; Maccagni et al., 2014). The H disk reaches velocities of kms, similar to the rotational velocities of the large scale H I disk (Maccagni et al., 2014). From now on, we will refer to it as the ‘outer disk’ of H. At radii pc, the major axis of the H disk abruptly changes orientation, from approximately north-south (PA) to east-west (PA). Hence, we refer to this as the ‘inner disk’ of H. The outer disk has asymmetries extending towards the inner disk, possibly suggesting that inner and outer disk are and part of a single, strongly warped structure. From the stellar continuum, we determine the distribution of the stars in the field of view; see the grey isophotes in Fig. 4 (a). The outer disk is aligned with the stellar component in the N-S direction. Conversely, the inner disk is polar.

In the inner disk, the major axis is aligned in the E-W direction perpendicular to the outer disk. As we move towards the centre, the velocity field suggests that the kinematic minor axis of the disk (green velocities in Fig. 4 (b)) may change its orientation within from the radio source. There, its axis of rotation appears to be aligned with the direction of propagation of the radio jets (dashed line in the Figures). However, given the quality of the data, the presence of this warp should be considered only suggestive.

3.2 The H 1-0 S(1) line in the innermost pc

The velocity field in Fig. 4 (a) suggests that, overall, the disk is dominated by rotation. However, as mentioned above, the very central region () shows a broader profile compared to the neighbouring regions, suggesting a much larger velocity dispersion in the innermost pc. In Fig. 5, we show the H 1-0 S(1) line profile extracted from the nucleus () and from two adjacent regions on either side of the nucleus (marked by crosses in Fig. 4 (b)). The line extracted from the central region (C) is shown in black. The dotted red and dashed blue lines, instead, show the H 1-0 S(1) line from two regions adjacent to the centre, on the east (E) and on the west (W), respectively. All spectra are centred on the systemic velocity of the galaxy. From the figure it is clear that in the centre the H line has an asymmetric profile, with a second component peaking at velocities kms. This component lies outside the range of velocities of the rotation, limited by the flanks of the blue and red lines. The centre of the galaxy ( pc) is the only region of the galaxy where the H 1-0 S(1) line has this feature.

This can also be illustrated by the position-velocity diagram extracted along the major axis of the inner disk (). Fig. 6 shows that the kinematics are characterised by rotation overall, as suggested by the smooth gradient in velocity along the x-axis, symmetric with respect to the centre of the galaxy and with respect to its systemic velocity. Nevertheless, in the centre ( pc), the profile appears to be broader and more asymmetric towards redshifted velocities (kms) than in the rest of the disk.

Some considerations of the rotation curve of the inner disk of PKS B1718–649 allow us to explore this kinematics in more detail. Willett et al. (2010) estimate that the mass of the SMBH is 4. Assuming that the velocity dispersion of the stars is 200 kms, this means that the SMBH dominates the kinematics of the galaxy out to pc; while beyond that radius the stellar mass distribution, which is well described by a de Vaucouleurs profile  (Veron-Cetty et al., 1995), also contributes. Fig. 6 shows a rotation curve based on such a model where we have assumed a total mass of PKS B1718–649 of , an effective radius of 9.7 kpc (Veron-Cetty et al., 1995), circular orbits of rotation, and corrected for the inclination of the inner disk333The ratio between the minor and major axis of the inner disk, assuming a finite thickness, indicates the disk is oriented approximately edge-on ().. Looking at the central pc, part of the broad profile can be described by the effect of the SMBH on the gas rotation. However, at redshifted velocities (kms) there is gas extending beyond the velocity range of the rotation curve expected for the mass model: kms. Although our model is fairly qualitative, it suggests gas with anomalous velocities () exists very close to the SMBH which may be directly involved in its fuelling (see Section 4 for further details) .

Figure 5: Spectra of the H 1-0 S(1) line, centred at the systemic velocity of PKS B1718–649. Spectra are extracted along the line of sight to the radio source (black), and on a region on the east (red) and on the west (blue) of the inner disk, at from the centre. Only the spectrum in front of the radio source appears broader and redshifted.
Figure 6: Position velocity plot of the H 1-0 S(1) line extracted along the major axis of the inner disk. Contour levels are –3, –2, 2, 3, 5, 7,9, and 12-. The black dashed line shows the rotation curve predicted from the stellar photometry, while the fine dashed line shows the contribution of the SMBH to the rotation. The solid line is the total rotation curve derived from the two. In the centre, at velocities kms, we identify a component of H deviating from the predicted rotation curve (see Section 3.2 for further details).

3.3 The temperature and mass of the H

The relative intensity of the H emission lines can be used to infer the temperature and mass of the molecular gas. We estimate the temperature in five different regions within the field of view from the fluxes shown in Table 2. As shown in Fig. 2, we choose two regions in the outer disk (R and R), two in the inner disk (R and R) and one in the centre (R). Following Jaffe et al. (2001); Wilman et al. (2005); Oonk et al. (2010) (and references therein), assuming that the gas is in local thermal equilibrium, the logarithm of the ratio between the flux of a H line and the flux of the H 1-0 S(1) line depends linearly on the excitation temperature () of the gas itself. The flux ratios of the H 1-0 lines suggest that the inner and outer disk have a temperature between K and K, while in the centre the temperature is lower K. In the R and R regions we measure upper limits for the H 1-0 S(0,2) line fluxes which are inconsistent with the local thermal equilibrium (LTE) temperatures derived from the 1-0 S(3) over 1-0 S(1) ratio. Deeper observations are necessary to further investigate these possible deviations from LTE.

The X-ray emission is localized in the innermost of the galaxy. Considering that X-rays in the energy range could be strongly absorbed, we estimated the lower limit on the luminosity ergs (Maccagni et al., 2014). This is only one order of magnitude higher than the H luminosity of the inner disk, ergs(where the flux is ergscm  at Mpc). This agrees with the thermal excitation scenario, and suggests that only a small fraction of high-energy photons is required to produce molecular hydrogen emission. The small size of the radio source ( pc) hints that shocks, if present, may excite the warm molecular gas only in the regions right next to the radio jet ( pc). Hence, this cannot be the main excitation process of the H, suggesting that thermal excitation is likely the main responsible mechanism for the warm H emission.

Given the temperature of the molecular hydrogen, we determine the mass of the Hof the inner disk, , from the flux of the H 1-0 S(1) line; as shown in Turner et al. (1977); Scoville et al. (1982) and Dale et al. (2005). From the data cube, we also measure the flux of the unsettled H component in the innermost pc (see Section 3.2): ergscm. This corresponds to .

The amount of warm gas found in the inner disk of PKS B1718–649 is in the same range of masses found in the innermost hundreds of parsecs of other LINER galaxies (Müller-Sánchez et al., 2013). Since the H is mainly thermally excited, the H 1-0 S(1) line may reflect the total mass of the cold molecular component, i.e. the H in its ground state ( K) commonly traced by the CO lines and the H 0-0 S(0,…,7) rotational lines. Within one order of magnitude, the mass of the cold H can be estimated from the mass of the warm H. We find  for the inner disk, and  for the unsettled H in the central pc. Where we use the relation found for a sample of galaxies with similar morphological classification as PKS B1718–649, (Mueller Sánchez et al., 2006; Dale et al., 2005; Mazzalay et al., 2013; Emonts et al., 2014),.

Parameter Inner disk Deviating component
Radius [pc] <650 < 75
Flux [ergscm]
Luminosity [ergs]
Temperature [K]
Mass H (warm) []
Mass H (cold ) []
Tab. 3: Main properties of the molecular hydrogen in the innermost regions of PKS B1718–649.

4 Relating the kinematics of the gas to the radio nuclear activity

The SINFONI observations in the innermost kilo-parsec of PKS B1718–649 reveal two disks of molecular hydrogen. The outer disk ( pc), oriented in the N-S direction, follows the rotation of the stars and of the other gaseous components of the galaxy. The inner disk ( pc) is oriented E-W with kinematics overall characterised by rotation. In Section 3.2, we showed that in the innermost pc of PKS B1718–649, the H 1-0 S(1) line is brightest and asymmetric, suggesting the presence of a second component of H with unsettled kinematics deviating from the rotation with redshifted velocities kms.

The H is not the only gaseous component with unsettled kinematics near the radio source. Along the same line of sight (and in particular only in front of the central pc of the radio source) the H I also shows kinematics deviating from regular rotation (Maccagni et al., 2014). Two separate absorption lines, with opposite velocities with respect to the systemic value, suggest the presence of small clouds of cold gas close to the AGN that deviate from the rotation of the other components of the galaxy. From the separation between the two lines, we estimate that these clouds have unsettled velocities of kms. Given that the distribution of the H in the inner disk is not homogeneous, it is reasonable to assume that the H I clouds are located in the same region of the H with unsettled kinematics, i.e. in the innermost pc of the galaxy. The presence of a clumpy multiphase environment around the radio source is also suggested by the variability of its radio-continuum (Tingay et al., 2015), which has been attributed to changing conditions in the free-free absorption in a surrounding clumpy cold medium. Moreover, optical spectroscopic observations also suggest the presence of a clumpy circum-nuclear medium (Filippenko, 1985). A similar distribution of H, with increasing velocity dispersion in the central 100 pc, has been detected in a number of different AGN and Seyfert galaxies  (Hicks et al., 2009, 2013; Davies et al., 2014; Guillard et al., 2012; Müller-Sánchez et al., 2013; Mazzalay et al., 2013; Mezcua et al., 2015).

PKS B1718–649 is thus a newly born compact radio AGN surrounded by a rotating clumpy multi-phase circum-nuclear disk, where we measure deviations from rotation in the H I and the H only in the innermost pc. While the H I, because it is detected in absorption, must be located in front of the radio source, the H is detected in emission and can be located either in front and/or behind the radio source. Hence, in principle, the redshifted velocities of the unsettled gas could correspont to either an infall or an outflow. It is difficult to disentangle this from the available data. However, given the properties of this AGN we note it is unlikely that it is driving an outflow. The small scale of the radio source ( pc) and its low jet-power (ergs) would exclude a jet-driven outflow. Since PKS B1718–649 is a LINER galaxy, the radiation from the optical AGN is also limited (ergs) and an outflow is not likely to occur on energetic grounds. These considerations make it plausible to assume that the redshifted unsettled velocities of the H are connected to gas falling into the AGN and perhaps being responsible for its fuelling. PKS B1718–649 does not show traces of previous periods of radio activity which could have perturbed the gas444since the 1.4 GHz continuum flux over pc (beam of the VLBI observations,Tingay et al. (2002)) is the same as over kpc (resolution of the ATCA observations, Maccagni et al. (2014)).. Hence, it is likely that, when the radio source was triggered, these clouds with unsettled kinematics (kms) were already present in the innermost pc of the circum-nuclear disk. The double disk structure of the H, and the large-scale strongly warped H I disk, suggest that the gas in PKS B1718–649 is still settling in the gravitational potential and that stellar torques are acting on the gas to align into a stable configuration. These torques may strip gas clouds from the inner two-disk configuration so that the clouds subsequently become unsettled and fall toward the SMBH.

Simulations of black hole accretion in rotating environments (King et al., 2008; Nayakshin & Zubovas, 2012; Gaspari et al., 2013, 2015) have suggested that, because of the local instabilities of the medium, chaotic collisions between clouds, cold filaments and the clumpy circum-nuclear disk may unsettle the kinematics of the gas and promote the cancellation of angular momentum. This may lead to the triggering of accretion into the SMBH. If the velocity dispersion does not exceed the rotational velocity, the accretion rates onto the AGN are predicted to be yr. This scenario could be an alternative explanation for the disturbed kinematics observed in the H I and the H of PKS B1718–649, where the deviations from rotation (kms) are on the same order of magnitude as the rotational velocity (kms). This scenario predicts inefficient accretion onto the AGN, which is also suggested by the radio power and by the LINER nature of PKS B1718–649.

In low-efficiency radio AGN, the radio power may set a constraint on the accretion rate onto the SMBH (Allen et al., 2006; Balmaverde et al., 2008). In PKS B1718–649  this is equal to yr. In Maccagni et al. (2014), we derived a limit for the contribution of the H I on the accretion using a very uncertain distance of the clouds from the nucleus due to the large beam of the H I observations. The H emission allows us to constrain this distance to pc, hence we determine the accretion rate of the H I clouds and of the H with unsettled kinematics, and investigate whether this could sustain the radio activity. Assuming the velocities deviating from rotation are equal to the in-fall velocity into the black hole ( kms) and assuming a distance of the clouds from the SMBH pc, we determine a typical timescale of accretion of these components to be years. The mass of the H I clouds is constrained by the column density of the absorption lines; assuming these are located within pc from the radio source, we determine . From this, it follows yr, which is insufficient, alone, to sustain the radio activity. In the innermost pc, the warm molecular hydrogen with unsettled kinematics has a mass of , which also gives an accretion rate of yr. If some of the cold H ( K, see Section 3.3) is also involved in feeding the AGN, we may obtain an accretion rate sufficient to power such radio source.

PKS B1718–649 has some interesting features in common with the nearest radio galaxy Centaurus A. Like PKS B1718–649, Centaurus A has a young radio core surrounded by a circum-nuclear rotating disk of H that is embedded in a large-scale H I disk (Struve et al., 2010). Centaurus A also shows a brighter and asymmetric H line profile, in the innermost pc. This can be explained by gas streaming down into the AGN (Neumayer et al., 2007). PKS B1718–649 appears to be another example where we witness the fuelling of a radio-loud AGN. We plan to furhter investigate this with future observations.

5 Conclusions

Our SINFONI m observations of the innermost of PKS B1718–649 have shown the presence of molecular hydrogen assembled into two orthogonal disks. The outer ( pc) disk of H is oriented along the north-south direction aligned with the stellar distribution and of which the kinematics connects smoothly to that of the large-scale H I disk. At radii pc, the H is assembled in an inner circum-nuclear disk, aligned in the east-west direction and polar with respect to the stars. The kinematics of the disks is characterised by rotation with velocities of about kms. Assuming thermal equilibrium within the disk, we determine the temperature of the H to be K and its mass , which may trace up to  of cold molecular hydrogen.

The kinematics of the inner disk of H is characterised by rotation due to the combination of the stellar distribution and the SMBH (see Section 3.2). In proximity of the radio source, at radii pc, we detect H deviating from such rotation. In the innermost pc the H has unsettled kinematics in the range kms. This component of warm H has a mass  , which may trace  of cold molecular hydrogen. The H I clouds detected in absorption against the compact radio core by Maccagni et al. (2014) have similar velocities deviating from rotation and could be located in the same region close to the radio source. These observations, along with the information collected from the variability of the radio continuum (Tingay et al., 2015) and the line ratios of the optical forbidden lines (Filippenko, 1985), suggest that the circum-nuclear ISM is clumpy and may represent the fuel reservoir of the radio source. The mass traced by the H I clouds and by the warm H alone is insufficient to fuel the AGN to power the radio jets. Instead, the mass of total cold H ( K) traced by the warm unsettled H in the innermost pc could fuel the radio source at the required accretion rate. Given the low power of the AGN, inefficient accretion is most likely to occur in PKS B1718–649. Given the double disk structure of the H, which is part of the larger ( kpc) H I disk, the gas configuration could be caused by the stellar torques acting on the gas to align into a stable configuration and which may give rise to small clouds with unsettled kinematics. The small clouds of H I and H with unsettled velocities of kms, that we detect, could be falling into the AGN, contributing to the fuelling of the radio source.


The research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Programme (FP/2007-2013) / ERC Advanced Grant RADIOLIFE-320745. BE acknowledges funding by the European Union 7th Framework Programme (FP7-PEOPLE-2013-IEF) grant 624351. The authors wish to thank J. Piqueras López for the help in the data reduction and the development of the IDL routines. SINFONI is an adaptive optics assisted near-infrared integral field spectrometer for the ESO VLT. The observations presented in this paper have been taken at the La Silla-Paranal Observatory under programme 093.B-0458(A).


  • Allen et al. (2006) Allen, S. W., Dunn, R. J. H., Fabian, A. C., Taylor, G. B., & Reynolds, C. S. 2006, MNRAS, 372, 21
  • Allison et al. (2015) Allison, J. R., Sadler, E. M., Moss, V. A., et al. 2015, MNRAS, 453, 1249
  • Athanassoula et al. (2005) Athanassoula, E., Lambert, J. C., & Dehnen, W. 2005, MNRAS, 363, 496
  • Balmaverde et al. (2008) Balmaverde, B., Baldi, R. D., & Capetti, A. 2008, A&A, 486, 119
  • Bedregal et al. (2009) Bedregal, A. G., Colina, L., Alonso-Herrero, A., & Arribas, S. 2009, ApJ, 698, 1852
  • Begelman & Shlosman (2009) Begelman, M. C. & Shlosman, I. 2009, ApJ, 702, L5
  • Combes (2004) Combes, F. 2004, ASSL, 319, 57
  • Combes (2011) Combes, F. 2011, IAUS, 6, 119
  • Combes et al. (2014) Combes, F., García-Burillo, S., Casasola, V., et al. 2014, A&A, 565, A97
  • Curran et al. (2013) Curran, S. J., Allison, J. R., Glowacki, M., Whiting, M. T., & Sadler, E. M. 2013, MNRAS, 431, 3408
  • Dale et al. (2005) Dale, D., Sheth, K., Helou, G., Regan, M., & Hüttemeister, S. 2005, AJ, 129, 2197
  • Dasyra & Combes (2011) Dasyra, K. M. & Combes, F. 2011, A&A, 533, L10
  • Davies (2007) Davies, R. I. 2007, MNRAS, 375, 1099
  • Davies et al. (2014) Davies, R. I., Maciejewski, W., Hicks, E. K. S., et al. 2014, ApJ, 792, 101
  • Eisenhauer et al. (2003) Eisenhauer, F., Bonnet, H., Abuter, R., et al. 2003, 4841, 1548
  • Ellison et al. (2008) Ellison, S. L., Patton, D. R., Simard, L., & McConnachie, A. W. 2008, AJ, 135, 1877
  • Emonts et al. (2010) Emonts, B. H. C., Morganti, R., Struve, C., et al. 2010, MNRAS, 406, 987
  • Emonts et al. (2006) Emonts, B. H. C., Morganti, R., Tadhunter, C. N., et al. 2006, A&A, 454, 125
  • Emonts et al. (2014) Emonts, B. H. C., Piqueras-López, J., Colina, L., et al. 2014, A&A, 572, A40
  • Fanti (2009) Fanti, C. 2009, Astron. Nachr., 330, 120
  • Filippenko (1985) Filippenko, A. 1985, ApJ, 289, 475
  • Gallimore et al. (1999) Gallimore, J. F., Baum, S. a., O’Dea, C. P., Pedlar, A., & Brinks, E. 1999, ApJ, 20, 48
  • García-Burillo et al. (2005) García-Burillo, S., Combes, F., Schinnerer, E., Boone, F., & Hunt, L. K. 2005, A&A, 441, 1011
  • Gaspari et al. (2013) Gaspari, M., Ruszkowski, M., & Oh, S. P. 2013, MNRAS, 432, 3401
  • Gaspari et al. (2015) Gaspari, M., Ruszkowski, M., & Oh, S. P. 2015, A&A, 432, 3401
  • Geréb et al. (2015) Geréb, K., Maccagni, F. M., Morganti, R., & Oosterloo, T. A. 2015, A&A, 575, A44
  • Geréb et al. (2014) Geréb, K., Morganti, R., & Oosterloo, T. A. 2014, A&A, 569, A35
  • Giroletti & Polatidis (2009) Giroletti, M. & Polatidis, A. 2009, Astron. Nachr., 330, 193
  • Guillard et al. (2014) Guillard, P., Boulanger, F., Lehnert, M. D., et al. 2014, A&A, 574, A32
  • Guillard et al. (2012) Guillard, P., Ogle, P., Emonts, B., et al. 2012, ApJ, 747, 95
  • Hicks et al. (2013) Hicks, E., Davies, R., Maciejewski, W., et al. 2013, ApJ, 768, 107
  • Hicks et al. (2009) Hicks, E., Davies, R., Malkan, M., et al. 2009, ApJ, 696, 448
  • Hopkins & Quataert (2010) Hopkins, P. F. & Quataert, E. 2010, MNRAS, 407, 1529
  • Hwang et al. (2012) Hwang, H. S., Park, C., Elbaz, D., & Choi, Y.-Y. 2012, A&A, 538, A15
  • Jaffe et al. (2001) Jaffe, W., Bremer, M. N., & van der Werf, P. P. 2001, MNRAS, 324, 443
  • Keel & Windhorst (1991) Keel, W. C. & Windhorst, R. A. 1991, ApJ, 383, 135
  • King & Nixon (2015) King, A. & Nixon, C. 2015, MNRAS, 453, L46
  • King & Pringle (2007) King, A. R. & Pringle, J. E. 2007, MNRAS, 377, L25
  • King et al. (2008) King, A. R., Pringle, J. E., & Hofmann, J. A. 2008, MNRAS, 385, 1621
  • Kormendy & Kennicutt (2004) Kormendy, J. & Kennicutt, R. 2004, ARA&A, 42, 603
  • Maccagni et al. (2014) Maccagni, F. M., Morganti, R., Oosterloo, T. A., & Mahony, E. K. 2014, A&As, 571, A67
  • Mazzalay et al. (2013) Mazzalay, X., Saglia, R. P., Erwin, P., et al. 2013, MNRAS, 428, 2389
  • Mezcua et al. (2015) Mezcua, M., Prieto, M. a., Fernández-Ontiveros, J. a., et al. 2015, MNRAS, 452, 4128
  • Morganti et al. (2013) Morganti, R., Fogasy, J., Paragi, Z., Oosterloo, T., & Orienti, M. 2013, Science, 341, 1082
  • Mueller Sánchez et al. (2006) Mueller Sánchez, F., Davies, R. I., Eisenhauer, F., et al. 2006, A&A, 454, 481
  • Müller-Sánchez et al. (2013) Müller-Sánchez, F., Prieto, M., Mezcua, M., et al. 2013, ApJL, 763, L1
  • Mundell et al. (2003) Mundell, C. G., Wrobel, J. M., Pedlar, A., & Gallimore, J. F. 2003, ApJ, 583, 192
  • Murgia (2003) Murgia, M. 2003, PASA, 20, 19
  • Nayakshin & Zubovas (2012) Nayakshin, S. & Zubovas, K. 2012, MNRAS, 427, 372
  • Neumayer et al. (2007) Neumayer, N., Cappellari, M., Reunanen, J., et al. 2007, ApJ, 671, 1329
  • Oonk et al. (2010) Oonk, J. B. R., Jaffe, W., Bremer, M. N., & Van Weeren, R. J. 2010, MNRAS, 932, no
  • Pecaut & Mamajek (2013) Pecaut, M. J. & Mamajek, E. E. 2013, ApJS, 208, 9
  • Piqueras López (2014) Piqueras López, J. 2014, PhD thesis, Universidad Complutense de Madrid
  • Ramos Almeida et al. (2012) Ramos Almeida, C., Bessiere, P. S., Tadhunter, C. N., et al. 2012, MNRAS, 419, 687
  • Sabater et al. (2013) Sabater, J., Best, P. N., & Argudo-Fernandez, M. 2013, MNRAS, 430, 638
  • Schawinski et al. (2010) Schawinski, K., Dowlin, N., Thomas, D., Urry, C. M., & Edmondson, E. 2010, ApJ, 714, L108
  • Scoville et al. (1982) Scoville, N. Z., Hall, D. N. B., Ridgway, S. T., & Kleinmann, S. G. 1982, ApJ, 253, 136
  • Skrutskie et al. (2006) Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163
  • Soker et al. (2009) Soker, N., Sternberg, A., Pizzolato, F., Heinz, S., & Wilcots, E. 2009, in American Institute of Physics Conference Series, Vol. 1201, AIPC, ed. S. Heinz & E. Wilcots, 321–325
  • Struve & Conway (2012) Struve, C. & Conway, J. 2012, A&A, 546, A22
  • Struve et al. (2010) Struve, C., Oosterloo, T. A., Morganti, R., & Saripalli, L. 2010, A&A, 515, A67
  • Tadhunter (2008) Tadhunter, C. 2008, NAR, 52, 227
  • Tingay et al. (1997) Tingay, S., Jauncey, D., Reynolds, J., et al. 1997, AJ, 113, 2025
  • Tingay et al. (2015) Tingay, S., Macquart, J.-P., Collier, J., et al. 2015, AJ, 149, 74
  • Tingay et al. (2002) Tingay, S. J., Reynolds, J. E., Tzioumis, A. K., et al. 2002, ApJS, 141, 311
  • Turner et al. (1977) Turner, J., Kirby-Docken, K., & Dalgarno, A. 1977, ApJS, 35, 281
  • Veron-Cetty et al. (1995) Veron-Cetty, M.-P., Woltjer, L., Ekers, R., et al. 1995, A&A, 297, L79
  • Wada (2003) Wada, K. 2003, Carnegie Observatories Astrophysics Series, 1, 16
  • Wada & Tomisaka (2004) Wada, K. & Tomisaka, K. 2004, ASPC, 320
  • Willett et al. (2010) Willett, K., Stocke, J., Darling, J., & Perlman, E. 2010, ApJ, 713, 1393
  • Wilman et al. (2005) Wilman, R. J., Edge, A. C., & Johnstone, R. M. 2005, MNRAS, 359, 755
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description