WISDOM project – VIII. Multi-scale feedback cycles in the brightest cluster galaxy NGC 0708

We present high-resolution (synthesised beam size 0 . (cid:48)(cid:48) 088 × 0 . (cid:48)(cid:48) 083 or 25 × 23 pc 2 ) Atacama Large Millimetre/submillimetre Array (ALMA) 12 CO(2–1) line and 236 GHz continuum observations, as well as 5 GHz enhanced Multi-Element Radio Linked Interferometer Network (e-MERLIN) continuum observations, of NGC 0708; the brightest galaxy in the low-mass galaxy cluster Abell 262. The line observations reveal a turbulent, rotating disc of molecular gas in the core of the galaxy, and a high-velocity, blue-shifted feature ≈ 0 . (cid:48)(cid:48) 4 ( ≈ 113 pc) from its centre. The sub-millimetre continuum emission peaks at the nucleus, but extends towards this anomalous CO emission feature. No corresponding elongation is found on the same spatial scales at 5 GHz with e-MERLIN. We discuss potential causes for the anomalous blue-shifted emission detected in this source, and conclude that it is most likely to be a low-mass in-falling ﬁlament of material condensing from the hot intra-cluster medium via chaotic cold accretion, but it is also possible that it is a jet-driven molecular outﬂow. We estimate the physical properties this structure has in these two scenarios, and show that either explanation is viable. We suggest future observations with integral ﬁeld spectrographs will be able to determine the true cause of this anomalous emission, and provide further evidence for interaction between quenched cooling ﬂows and mechanical feedback on both small and large scales in this source.


INTRODUCTION
The lack of molecular gas in early-type galaxies (ETGs; ellipticals and S0s) has been a point of debate for some decades (e.g.Faber & Gallagher 1976;Lees et al. 1991; molecular gas reservoirs are an order of magnitude smaller than expected (based on cooling rates estimated from Xray emission; e.g.Lees et al. 1991;Fabian 1994;De Lucia & Blaizot 2007).
Searches for this cooled gas have persistently returned lower gas masses and fewer young stars than required by cooling flow observations (e.g.Johnstone et al. 1987;Heckman et al. 1989;McNamara & O'Connell 1989;Crawford et al. 1999;Donahue et al. 2000;Hoffer et al. 2012).For instance, in a range of brightest cluster galaxies the mass of molecular gas is found to be 5 − 10 per cent of that expected from the cooling flow (i.e.Edge 2001;Salomé & Combes 2004).A solution to this 'cooling flow problem', how gas leaves the hot phase but does not condense at the expected large rates (of 100-1000 M yr −1 ) on to the central galaxy, has been sought ever since.
High-angular resolution X-ray observations paved the way for answers, showing that despite appearing relaxed at low resolution, the centres of cooling flow clusters are in fact very dynamic.Active galactic nuclei (AGN) with powerful jets are found in essentially all cooling flow cluster central galaxies (Sun 2009), and they are the principal power source driving the ICM dynamics (e.g.Bîrzan et al. 2004Bîrzan et al. , 2012;;McNamara et al. 2005;Rafferty et al. 2006;McNamara & Nulsen 2007, 2012;Gaspari et al. 2013;Hlavacek-Larrondo et al. 2015).AGN jets appear able to inflate large bubbles in the hot ICM, that rise buoyantly and disrupt the cooling flows.Heat from the AGN is also distributed in the ICM through turbulent mixing and cocoon shocks (e.g.Gaspari et al. 2013;Yang et al. 2019;Wittor & Gaspari 2020).It has been shown that AGN jets have the mechanical power to balance the ICM's energy losses due to cooling, motivating the theory that mechanical (i.e.radio-mode) feedback is the principal regulator of ICM cooling, thus preventing a run away process (see reviews from e.g.McNamara & Nulsen 2007, 2012;Gaspari et al. 2020).Simulations also suggest AGN feedback is vital for the regulation of a galaxies gas reservoir and therefore its star formation rate.For instance, models including radio-mode feedback are in better agreement with e.g. the galaxy luminosity function (Bower et al. 2006(Bower et al. , 2008;;Croton et al. 2006;McCarthy et al. 2008;Davé et al. 2012).
The advent of high-resolution radio/sub-mm interferometry has begun to add to the growing picture of feedbackcontrolled galaxy evolution.In multiple cooling flow clusters significant (10 9 -10 10 M ) amounts of molecular gas have been detected in filaments which are sometimes coincident with buoyant X-ray bubbles rising through the ICM (e.g.McNamara et al. 2014;Russell et al. 2014Russell et al. , 2016Russell et al. , 2017b,a;,a;Vantyghem et al. 2016Vantyghem et al. , 2018)), and sometimes present throughout the inner ICM (Temi et al. 2018;Rose et al. 2019;Juráňová et al. 2019).It is not currently known if this cold gas has recently cooled from low-entropy gas lifted by the bubble, or is stimulated to cool in-situ by its passing.Many of these observed filaments also have star formation associated with them (e.g.Vantyghem et al. 2018).
A variety of works have used simulations to look at the formation of this multiphase ICM to ascertain how it is regulated.Gaspari et al. (2012), Sharma et al. (2012), Prasad et al. (2015) and Li et al. (2015) all see cycles within their simulations where dense, cold gas filaments condense out of the ICM and precipitate onto the central galaxy.This causes star formation and fuels the central SMBH (that becomes an AGN) via a process known as chaotic cold accretion (CCA; see Gaspari et al. 2019 for a recent overview).The AGN and supernova winds increase heating within the ICM, returning it to a high entropy state, stopping cooling and hence the fuel supply.When the heating stops, cooling resumes again.They find that cold gas filaments form when the instantaneous ratio of the thermal instability and free-fall (or eddy turnover) timescales is < ∼ 10 (e.g.Gaspari et al. 2012;Sharma et al. 2012, Li et al. 2015;Gaspari et al. 2018).
In all scenarios described above, the AGN is a crucial driver of the evolution of the gaseous material in brightest cluster galaxies, with jet-blown bubbles dominating on the scale of ten to hundreds of kilo-parsecs.The feedback produced by these AGN, however, interacts with multiphase condensation processes across a large range of physical scales.Tremblay et al. (2012aTremblay et al. ( ,b, 2016) ) found multiwavelength evidence of both large-and small-scale mechanical feedback in the BCG of galaxy cluster Abell 2597.They reported an extensive kpc-scale X-ray cavity network, with multiple rising buoyant bubbles, the largest of which coincides in both linear extent and position angle with the radio jet (Tremblay et al. 2012a).Hubble Space Telescope (HST) and Herschel observations reveal ongoing star formation co-spatial with knots in the X-ray emission (Tremblay et al. 2012b).Atacama Large Millimetre/submillimetre Array (ALMA) 12 CO(2-1) observations further add to this picture, exposing chaotic cold accretion onto the central SMBH (by revealing absorption features in the AGN continuum caused by clouds moving inwards towards the SMBH; Tremblay et al. 2016) and a molecular outflow associated with the jet itself.
Here we report on molecular gas observations of NGC 0708, the BCG in the low-mass galaxy cluster Abell 262, itself part of the Perseus-Pisces galaxy supercluster.NGC 0708 lies 58.3 ± 5.4 Mpc away (estimated using infrared surface brightness fluctuations; Jensen et al. 2003).It is a giant elliptical galaxy with a weak dust lane (Ebneter & Balick 1985;Wegner et al. 1996) and an effective radius of 33 (≈ 9.3 kpc;Wegner et al. 2012).See Fig. 1 for a HST image of NGC 0708.Abell 262 was identified as having an Xray emitting ICM by Jones &Forman (1984), andStewart et al. (1984) measured the cooling time to be 1.3 × 10 9 yr, smaller than the age of the universe so that the cluster is expected to form a cooling flow.The 20-cm observations of Parma et al. (1986) revealed a double-lobed, 'S'-shaped jet and led to the classification of NGC 0708 as a weak Fanaroff-Riley Class I radio source (Blanton et al. 2004).The top panel of Fig. 1 also has 330 MHz continuum observations from Clarke et al. (2009) overlaid (blue contours), to show the shape and orientation of the large-scale jet.Analysis of Chandra observations revealed a hole or bubble within the ICM, co-spatial with the eastern lobe of the jet (Blanton et al. 2004).Clarke et al. (2009) found additional 3-6 kpc radius cavities at differing position angles within the X-ray gas, and at a range of radial distances from the BCG (8 -29 kpc), indicating multiple episodes of AGN activity from a precessing SMBH jet.They concluded that the total AGN emission should be capable of counteracting the cooling flow over several outbursts.Using their multi-frequency observations of NGC 0708, Clarke et al. (2009) also calculated the radio spectral index (α) from 235 to 610 MHz, finding the spectrum to be flat in the core (α = −0.5),typical of new particles in a jet.They also estimated a lower limit on average outburst repetition timescales in Abell 262 to be τ rep ≥ 28 Myr.
NGC 0708 was thus observed to have large-scale feedback affecting the hot gas.In this work, we show that the cold interstellar medium (ISM) is also being affected on small scales by the feedback cycle in this object, as expected if AGN feeding/feedback is a multiscale self-regulated process.In Section 2, we present new ALMA and enhanced Multi-Element Radio Linked Interferometer Network (e-MERLIN) observations of NGC 0708.We present our analysis of these observations in Section 3. In Section 4 we discuss our results and compare NGC 0708 to other galaxies.We conclude in Section 5.
CO(2-1) in NGC 0708 was observed three times with ALMA (using different array configurations) as part of the mm-Wave Interferometric Survey of Dark Object Masses (WISDOM) project, aiming to measure its central SMBH mass.Previous work (e.g.Woo & Urry 2002;Donato et al. 2004) suggested an SMBH mass M BH ≈ 2.9 × 10 8 M , although this is highly uncertain.Olivares et al. (2019) published our intermediate-resolution observations from ALMA in a study of filaments in cool-core clusters.The observations at 0. 95 × 0. 61 (≈ 270 × 170 pc 2 ) show no filaments but instead a slightly warped rotating kiloparsec scale disc of molecular gas.Here we study this source in detail by including both higher and lower-resolution observations from WISDOM.

ALMA observations
Using ALMA, we observed the 12 CO(2-1) line in NGC 0708, first under programme 2015.1.00598.S at moderate angular resolution (0. 52 or ≈ 146 pc) on June 27th 2016 (published in Olivares et al. 2019), and then under programme 2017.1.00391.S at 0. 25 (≈ 70 pc) resolution on November 12th 2017 and 0. 03 (≈ 8.5 pc) resolution on September 19th 2018.The science target integration times for these were 11, 19 and 37 min, respectively.The baselines ranged from 15 m to 14 km, yielding a maximum recoverable scale of 7 (≈ 2 kpc), sufficient to cover the extent of the majority of the dust features in this source.For all observations a 1870 MHz (≈ 2500 km s −1 ) correlator window was centred at 226.8 GHz (the redshifted 12 CO(2-1) line frequency) with a raw channel width of 976.5 kHz (≈ 1.87 km s −1 ).To detect continuum emission, three additional low-spectral resolution correlator windows were included, each with a bandwidth of 2 GHz.
The raw data were calibrated using the standard ALMA pipeline, as provided by the European ALMA Regional Centre staff.The calibrators used for all observations were J0237+2848 for flux and bandpass calibration and J0205+3212 for phase calibration.The three observation tracks were combined and imaged using the Common Astronomy Software Applications (casa; McMullin et al. 2007).Continuum emission was detected, measured over the full line-free bandwidth, and then subtracted from the data in the uv-plane using the casa task uvcontsub.

e-MERLIN observations
We also observed NGC 0708 with e-MERLIN, to identify if any small-scale radio structures are present (e.g. a restarted radio jet from the recurrent, precessing AGN in the core of this galaxy).NGC 0708 was thus observed twice with e-MERLIN, the data providing sensitivity to 5 GHz emission distributed on the same angular scales as our 236 GHz data.The e-MERLIN data were processed through the standard e-MERLIN casa pipeline (eMCP 1 ) by the e-MERLIN fa-cility staff.The calibrators used were 0152+3616 for phase, 0319+4130 for pointing, 1331+3030 for flux and 1407+2827 for bandpass calibration.The total on-source integration time was 14.5 hours.
We additionally performed self-calibration to increase the sensitivity.The self-calibration involved 2 cycles, the first considering phase only, averaging over 240 s intervals, the second with both phase and amplitude.We imaged the data in casa using the tclean task, with a Briggs weighting robust parameter of 0.5 to balance sensitivity and resolution.This yielded a synthesised beam size of 0. 07 × 0. 03 (≈ 20 × 8 pc 2 ) and a RMS noise of 0.12 mJy beam −1 .

ALMA line emission
The moment maps, shown in Fig. 2, were created using the smooth-mask technique (e.g.Dame 2011).The mask was generated by taking a copy of the cleaned, primary beamcorrected cube and smoothing it, first spatially using a Gaussian of FWHM equal to that of the synthesised beam, and then spectrally using a Gaussian of FWHM of 4 channels.We then select pixels with an amplitude in the smoothed cube greater than 5 times the RMS noise in that cube.The mask was then applied to the un-smoothed cube to create the moment maps.Having said that, all quantitative analyses reported in this paper were performed using the unsmoothed, un-masked cube.The moment maps are shown only for illustrative purposes.
The zeroth (integrated intensity) and first (intensity-weighted mean line-of-sight velocity) moment maps reveal a rotating but warped molecular gas disc (see left and central panels Fig. 2).The moment zero is also shown in the bottom panel of Fig. 1 (blue contours), overplotted on a HST image, revealing that the molecular gas is coincident with dust features.The second moment (intensity weighted line-of-sight velocity dispersion; right panel Fig. 2) shows evidence of disturbance, with an off-centre peak significantly away (≈ 0. 4 or ≈ 113 pc) from the AGN position (that can also be independently measured from our data using the 236 GHz continuum emission; see Section 3.2).A major-axis position-velocity diagram (PVD; Fig. 3) was created by taking a 5-pixel wide pseudo-slit across the kinematic major axis of the cube, at a position angle of 349 • .This position angle was estimated by eye and agrees with that found by Pandya et al. (2017) for the ionised-gas disc traced by [O iii].
On the approaching side, the PVD has a sharp increase in velocity at a distance of ≈ 0. 4 (≈ 113 pc) from the centre, co-spatial with the aforementioned increase in velocity dispersion (see the right panel of Fig. 2 and the blue ellipse in Fig. 3).The properties of this blue-shifted emission will be considered further in Section 3.4.The global spectrum, shown in Fig. 4, was created by binning the data to 20 km s −1 channels and then integrating over the whole molecular gas disc, i.e. a 6 ×6 (1.7×1.7 kpc 2 ) area of the cube.This spectrum clearly shows the characteristic double-horned profile of a rotating disc, but with hints of an additional blue-shifted wing (highlighted by the magenta ellipse).The total 12 CO(2-1) flux detected in NGC 0708 is 57.9±0.1(stat)±5.8(sys)Jy km s −1 .The uncertainties quoted are the 1σ statistical uncertainty, and the ≈10 percent systematic flux calibration uncertainty, that typically dominates over the statistical uncertainty.
We estimate the H 2 gas mass present in this system, using the following standard equation: where m H is the mass of a hydrogen atom, λ is the rest wavelength of the observed molecular transition, is the line ratio (measured in beam temperature units) between the reference CO transition observed and the ground state CO(1-0) line, S v δ V is the integrated CO flux density and X CO is your CO-to-H 2 conversion factor of choice.This can be simplified to: and J upper is the upper state rotational quantum number for the transition observed (here J upper = 2).The other symbols are as defined above.
We here assume a typical Milky Way-like CO-to-H 2 conversion factor of 2×10 20 cm −2 (K km s −1 ) −1 (Dickman et al. 1986) (equivalent to α CO ≈ 4.6 M (K km s −1 ) −1 pc −2 ) and that the gas is has a line ratio T b,CO(2−1) /T b,CO(1−0) = 0.25 (Edge 2001), and thus R 2−1 = 4.The total flux therefore corresponds to a total molecular gas mass M tot = (1.5 ± 0.01(stat)±0.26(sys))×10 9M , which is slightly higher than that found by previous single dish measurements (Edge  2001; Salomé & Combes 2003;Ocaña Flaquer et al. 2010) suggesting that our combined ALMA observations do not resolve out a significant fraction of the extended flux in this object, despite high angular resolution.We note that the molecular gas mass derived here is different from that derived from a subset of these same data in Olivares et al. (2019).This difference arises because of different assumptions about the T b,CO(2−1) /T b,CO(1−0) ratio and X CO .

ALMA continuum emission
As mentioned in Section 2.1, NGC 0708 also has 1 mm continuum emission, detected by ALMA at a mean frequency of 236 GHz in the three low-resolution spectral windows and the line-free channels of the high-resolution spectral window.
The emission is resolved and has an extension to the South, clearly revealed in the left-hand panel of Fig. 6, showing the 236 GHz continuum emission (magenta contours) overlaid on the CO(2-1) velocity dispersion map.This highlights both the approximate coincidence between the extension of the continuum emission and the peak of the velocity dispersion, and the offset of that peak from the AGN position (i.e. the centre of the continuum emission).The magenta arrows in the left-hand panel of Fig. 6 show the direction of the largescale 330 MHz jet, to highlight the difference of orientation between that and the 236 GHz emission.The total 236 GHz continuum flux is 32.3 ± 0.2 mJy.Excluding the nucleus (by subtracting a point-source fit to the central peak from the total flux), the extended continuum emission component has a flux of 1.0±0.2mJy (where both figures are quoted with 1σ statistical uncertainties).

e-MERLIN 5 GHz continuum emission
The extension of the 236 GHz continuum emission in NGC 0708 is perpendicular to the large-scale jet (as traced by 330 MHz emission; see Figs. 1 and 6) and this prompted us to obtain additional 5 GHz continuum data, to ascertain if the 236 GHz continuum emission is from a small (potentially restarted) jet.These observations are described in Section 2.2.We detect a 5 GHz point source at the expected position of the SMBH in NGC 0708.We confirmed this source is spatially unresolved using the casa task imfit, that fits a Gaussian to the image, deconvolved from the synthesised beam.The integrated flux at 5 GHz is 5.25 ± 0.21 mJy (1σ statistical uncertainty).This is ≈ 5 times smaller than that measured by Clarke et al. (2009) at 5 GHz with the VLA on ≈ 4 (≈ 1.1 kpc) scales and 6 − 8 times smaller than singledish 5 GHz measurements (≈ 2. 6 or ≈ 44 kpc scales, Andernach et al. 1980; ≈ 10 or ≈ 170 kpc scales, Gregory et al. 1996).The disparity between these measurements and ours suggests that significant 5 GHz emission is associated with the large-scale radio jet, that we resolve out here.There is no obvious small-scale (restarted) jet visible at 5 GHz, although our observations are only deep enough to detect components with a 5 -236 GHz spectral index < −0.36 (at 3 σ ).If a young jet is present in this source this suggests that its 5 -236 GHz spectral index is shallow or inverted.

Blue-shifted component
Our ALMA data enable us to spatially and kinematically separate distinct components of the molecular gas distribution.The velocity dispersion map (right panel of Fig. 2) and PVD (Fig. 3) clearly indicate two kinematically-distinct components are present in NGC 0708, a disturbed but regularly rotating disc and a blue-shifted 'spike' feature.Here we investigate this anomalous blue-shifted feature.

Kinematic model
We begin by kinematically modelling the disc component at the centre of NGC 0708, to reveal the full extent of the anomalous emission.To construct this model we utilise the Python implementation of the KinMS2 kinematic modelling code (Davis et al. 2013(Davis et al. , 2020)), and fit the model to our data using the Bayesian Markov-Chain Monte Carlo (MCMC) code GAStimator3 .
We assume the mass distribution is axisymmetric (at least in the inner regions of interest here), and model the rotation curve of this system as a function of radius (R) with an arctangent function where V max is the asymptotic velocity, and R turn is the turnover radius of the rotation curve.We further assume that the molecular disc has an exponential surface brightness profile both as a function of radius and vertical extension (z), where I 0 is the central surface brightness, and R scale is the exponential scale length and z scale is the vertical scale height of the disc.These functional forms are very simple, and cannot reflect the full suite of rotation curves/surface brightness profiles present in nature.Despite this, here they allow us to construct a simple model which matches the data reasonably well.
KinMS takes these input functional forms, and allows us to construct a model datacube assuming the gas is in circular rotation around some unknown kinematic centre.We assume the gas has an intrinsic velocity dispersion (σ ) which is another free parameter in our fitting process.One can only estimate the gas dispersion from full kinematic modelling, because the observed moment two map is often dominated by beam smearing, especially in the central regions where the rotation curve rises quickly.We further assume here that all the gas is distributed in a single plane (i.e. the disc has the same kinematic position angle and inclination to our lightof-sight throughout the disc).While these assumptions are unlikely to be valid in the more disturbed outer disc of this galaxy, they do produce a reasonable fit in the central regions (where dynamical times are very short).We note that if we nevertheless allow radial variations of the inclination and position angle our best fit disc is consistent with being flat, suggesting that although it is morphologically lopsided, it may not be kinematically warped.Non-circular motions are almost certainly present within this gas disc, but here we aim to determine how much of the gas motion can be explained purely by regular rotation.
As mentioned above, we use a Bayesian MCMC analysis technique to identify the model which best matches our data, and simultaneously estimate uncertainties.Indeed, this allows us to obtain samples drawn from the posterior distribution of each model parameter.Our model, as described above, has a total of eleven free parameters (six key-parameters, and five nuisance parameters).The key parameters are the total flux of the system (which sets I 0 as described above), the two free parameters of the rotation curve, the scale radius of the molecular disc, its thickness in  the vertical direction, and its velocity dispersion.The nuisance parameters the position angle and inclination of the disc, and its kinematic centre (in RA, Dec, and velocity), Each parameter has a prior, that we typically set as a boxcar over a reasonable parameter range (an assumption of maximal ignorance).Details of all the priors are listed in Table 1.
Our MCMC procedure generates a model datacubes and compares them to our observed data using a simple log-likelihood (L ): where χ 2 is the standard χ 2 statistic.
As discussed in detail in Smith et al. (2019), because our ALMA data are noisy, the χ 2 statistic has an additional uncertainty associated with it, following the chi-squared distribution (Andrae 2010).Systematic effects can produce variations of χ 2 of the order of its variance (van den Bosch & van de Ven 2009), and ignoring this effect yields unrealistically small uncertainties.To mitigate this effect van den Bosch & van de Ven (2009) proposed to increase the 1σ confidence interval to ∆χ 2 = √ 2N, where N is the number of constraints.To achieve the same effect within our Bayesian MCMC approach we need to scale the log-likelihood, as done by Mitzkus et al. (2017).This is achieved here by increasing the input errors (i.e. the measured RMS noise in the cube) by (2N) 1/4 .This approach appears to yield physically credible formal uncertainties on the inferred parameters, whereas otherwise these uncertainties are unrealistically small.Within GAStimator we utilise a MCMC method with Gibbs sampling and adaptive stepping to explore the parameter space.The algorithm runs until convergence is achieved, and then the best chain is run (with a fixed step size) for 30,000 steps (with a 10% burn-in) to produce our final posterior probability distribution.We then marginalise over the probability surfaces for each model parameter to identify a best-fitting value (the median of the marginalised posterior samples) and associated 68% and 99% confidence levels (CLs).A quantitative description of the likelihood of each model parameter is presented in Table 1.
The best-fitting disc model requires a non-zero scale height (71±16 pc), and has a fairly large internal velocity dispersion (66 ± 5 km s −1 ).This suggests a reasonably turbulent disc, with a V max /σ ≈ 2.5 (compared to V/σ ≈ 10 found in typical relaxed low-redshift discs; Wisnioski et al. 2015).The major-and minor-axis PVDs extracted from the best fitting model are shown as black contours in the left and central panels of Figure 5, overlaid on the data (orange scale).Overall the model is a good fit to the data outside of the anomalous emission region.
In the right panel of Figure 5 we show the major-axis PVD created from the residual cube (data-model; orange scale), again overlaid with the best-fitting disc model (black contours).The major positive residuals are in the anomalous emission region, as expected, but residual emission is detected all the way to the systemic velocity.Additional, smaller (but still significant) residual features are detected at a similar velocity on the red-shifted side of the galaxy, and further out in the disc on the blue-shifted side.It is unclear if these features are in any way related to the dominant blueshifted component, or are simply brighter molecular knots within the disc of the galaxy.

Properties of the blue-shifted emission
As Figure 5 makes clear, the anomalous blue-shifted emission we detect in NGC 0708 cannot be explained by simple rotational motion of the ISM.To learn more we begin by isolating the gas in this anomalous component from that in the main gas disc.In Fig. 6 we constrain its spatial extent using the velocity dispersion peak seen in Fig. 2, adopting the region −0.32 < RA offset < 0 and −0.52 < Dec. offset < −0. 13 (indicated by the black box in Fig. 6) relative to the peak of the 236 GHz continuum (RA = 01 h 52 m 46.s 46 and Dec. = +36 • 09 06.46).In velocity, we impose v obs < 4550 km s −1 , indicated in the PVD (Fig. 3) by a grey dotdashed line.
We note that these cuts only include emission that is unambiguously outside the normal rotation of NGC 0708.One could estimate the properties of the anomalous emission using the residual cube presented in Figure 5, but it is unclear whether all the residual emission detected in the disc is related, or whether some might simply be due to brighter regions in the disc which are not well fitted by our axisymmetric model.Where appropriate below we discuss whether this choice would change our conclusions.
The right-hand panel of Fig. 6 shows the 12 CO(2-1) integrated spectrum of the spatial region containing the anomalous emission, with the channels satisfying the adopted velocity criterion (i.e. the blue-shifted wing) highlighted in magenta.We also show the result of fitting two Gaussians to this spectrum, one for the anomalous gas (shown in blue) and another for the gas in regular rotation in the galaxy disc (shown in orange).The sum of these two Gaussians is shown in red.The spectrum of this region has a 12 CO(2-1) integrated flux of 2.7 ± 0.1 Jy km s −1 , while the flux associated with the anomalous emission only (the area under the blue Gaussian in the right panel of Fig. 6) is 0.54 ± 0.05 Jy km s −1 (both 1σ statistical uncertainty).The anomalous emission is blue-shifted by 284 ± 6 km s −1 from the galaxy systemic velocity, and has a line-of-sight velocity width (RMS/dispersion of the blue Gaussian: σ v,los ) of 51 ± 5 km s −1 .The properties of the Gaussian fitted to the anomalous emission are listed in Table 2.We note that if we were to estimate the flux of all residual emission at this position, after subtraction of the disc model presented above, the flux of this anomalous emission would increase by a factor of ≈2.
Using our ALMA data we can estimate the extent of this anomalous emission in the plane-of-the-sky using the imfit task in CASA.In the channel where the fitted anomalous component reaches its peak intensity (4466 km s −1 ) the emission is marginally spatially resolved, with a diameter of (0. 28±0.07) × (0. 20±0.05), or ≈66±16 pc.
To estimate the amount of molecular gas associated with this feature, we need to assume a CO-to-H 2 conversion factor.However, the opacity and density of the gas in the anomalous feature is unknown, so we will conduct the analysis with three representative α CO values (following Morganti et al. 2015): a typical local/Milky Way factor (α CO, MW = 4.6 M (K km s −1 ) −1 pc −2 , as assumed for the bulk of the gas), a factor appropriate for the disturbed gas typically found in ultra-luminous infrared galaxies (ULIRGs; α CO, ULIRG = 0.8 M (K km s −1 ) −1 pc −2 ) and a factor appropriate for optically thin gas (α CO, thin = 0.34 M (K km s −1 ) −1 pc −2 ; see discussions of α CO in Bolatto et al. 2013 andGeach et al. 2014).As above, we assume a line ratio T b,CO(2−1) /T b,CO(1−0) = 0.25 (Edge 2001).The derived masses (tabulated in Table 2) vary from ≈ 10 6 to ≈ 10 7 M .This component thus contains only a small amount of the total molecular gas in this galaxy (<1%).
These derived masses enable us to derive plausible limits on the physical extent of the emitting component.We know that the anomalous gas is dense enough for CO molecules to survive and be excited enough to emit at their second energy level.Thus the volume density of this component is likely > ∼ 1000 cm −3 .If the emission we observe is coming from a spherical cloud of uniform density then its radius must be <41 pc (assuming α CO, MW ) or <17 pc (assuming α CO, thin ).Assuming a cylindrical geometry (appropriate e.g. if this material is a filament aligned with our line-of-sight) and a diameter of 66 pc (as measured above), the length of this filament must be <90 pc (assuming α CO, MW ) or <6 pc (assuming α CO, thin ).While this material is unlikely to have a constant density, these limits suggest that the anomalous emission is likely to arise from a physically-small component, that nonetheless has a large internal velocity gradient and/or dispersion.

DISCUSSION
In the above sections we have shown that NGC 0708 hosts a rotating disc of molecular gas at its core, coincident with the dust-lanes visible in optical extinction (Figure 1).It also contains an anomalous, kinematically-distinct blue-shifted component (shown clearly in Fig. 3), which has a large internal velocity gradient/dispersion.The velocity field and PVD of NGC 0708 indicate most of the gas is in regular rotation in the gravitational potential of the galaxy, although the disc is fairly turbulent.The blue-shifted feature, however, cannot arise from regularly-rotating material.The origin of this material is ambiguous, as both inflows and outflows can produce similar features.Below we therefore consider potential evidence for inflow or outflow, using our ALMA (and eMERLIN) observations and literature data from other wavelengths.

Inflow
Inflows are common in BCGs, especially those located inside cool-core clusters like Abell 262.Inflows can also be found in more normal galaxies, and can arise due to gas cooling or mergers.
In NGC 0708 we can rule out the type of inflows seen in isolated and satellite galaxies.In such systems gas cooling tends to be slow and to occur primarily along the disc plane (see e.g. the simulations of Agertz et al. 2009 andStewart et al. 2011).It would then most prominently appear in the moment 1 map as gas in the outer disc that is not following the expected rotation pattern.However, the blue-shifted feature we see in NGC 0708 is well-collimated, and essentially unresolved spatially, indicating that a significant component of the velocity is along the line of sight, i.e. out of the plane of the disc.The blue-shifted feature is also very close to the centre of the galaxy.
Mergers tend to be identified most readily in optical images, by the disturbed morphology of the galaxy and tails of gas and stars, as gas inflowing on to the more massive galaxy tends to form tidal tails extending over many kiloparsecs.In contrast, HST imaging of NGC 0708 (Fig. 1) shows no sign of a significant recent major merger, with no tidal tail or disturbance of the stellar body.The dust distribution is however messy and warped, suggesting that a minor, gasrich merger may have occurred.However, the PVD feature detected by ALMA is well collimated and has a small extent in the plane-of-the-sky.Such a small, localised and central feature is unlikely to arise from a tidal tail.Therefore, although we cannot definitively rule it out, the collimation, size and position of the blue-shifted feature indicate that secular or merger-driven inflow is an unlikely explanation.
This, however, leaves the most likely possibility.NGC0708 is a BCG which is expected to be fuelled by a quenched cooling flow4 .Molecular gas clumps and filaments are expected to precipitate out of the hot medium and rain down onto the galaxy.In large samples of observed BCGs these in-falling clumps/filaments are common, multiple typically being found around each BCG or brightest group galaxy (e.g.Temi et al. 2018;Gaspari et al. 2018;Rose et al. 2019;Olivares et al. 2019).These in-falling filaments typically contain a large fraction of the total molecular gas mass of the BCG.
Figure 7 shows the ratio of the molecular gas mass contained in in-falling filaments to the total BCG molecular gas mass for a range of clusters, plotted against the classical cooling rate of the hot ICM (McDonald et al. 2018), as presented in Olivares et al. (2019).If the anomalous emission in NGC 0708 comes from a single in-falling filament on the far side of the galaxy then we can place it on this diagram.The measurements from this paper for NGC 0708 are shown in Figure 7 with an orange star.We note here that we use the estimate of the molecular mass of this feature assuming a Galactic CO-to-H 2 conversion factor (to match the assumptions made in the literature sample), as listed in Table 2.If the true CO-to-H 2 conversion factor of the filament is lower, our measurement would move down.
In NGC 0708 we only see clear evidence of a single nonequilibrium structure, that only contains a very small fraction of the cold gas in this galaxy.Only a single BCG from the Olivares et al. ( 2019) sample (RXJ1539.5)has a similar low-mass in-falling structure (the west filament, one of three detected in this system).However, detecting low-mass structures like that observed here requires high-sensitivity and high spatial-resolution observations, which are not available for all BCGs.In addition, Abell 262 (the parent cluster of NGC 0708) has the lowest observed cooling rate in the Olivares et al. ( 2019) sample.It is thus possible that the low number of observed in-falling structures, and the low mass of the structure seen, could be typical of clusters with such low cooling rates.
To match the observed velocity structure of this blueshifted feature (which extends over > ∼ 200 km s −1 while having a physical extension perpendicular to our line of sight of < ∼ 66 pc), it is clear that any such clump must have significant internal velocity structure.This is because (for any realistic potential) purely ballistic infall would require an extremely long column of gas (in which molecular gas would not survive) to create such a velocity gradient, while keeping the elongated gas structure perfectly aligned with our line of sight.Gaspari et al. (2018) performed three-dimensional highresolution hydrodynamical simulations of CCA and the re-0.5 1.0 1.5 2.0 2.5 3.0 M cool (M yr 1 ) 10 2 10 1 10 0 M H2, filament /M H2, total Olivares+19 NGC 0708 (this work) Figure 7. Ratio of the molecular gas mass contained within inflowing filaments to the total molecular gas mass of the BCG, plotted as a function of the classical cooling rate of the hot intracluster medium in each galaxy cluster Olivares et al. (2019).NGC 0708 is shown as an orange star.It has a lower (classical) cooling rate than all other cluster objects (consistent with the lower total mass of Abell 262), and also contains a much lower fraction of its cold molecular gas in potentially inflowing filaments. .Velocity dispersion of molecular/ionised gas along individual lines of sight ('pencil beams') towards BCGs, plotted against the mean velocity shift of their line centroid from the host system's systemic velocity.The brown datapoints are reproduced from Table 2 of Gaspari et al. (2018), while our estimates from NGC 0708 are shown as an orange star.Green contours show the 1 to 3σ range in these quantities seen in the simulations of Gaspari et al. 2018.The anomalous emission in NGC 0708 has a significant internal velocity gradient, which is consistent at 1σ with the expectation from CCA simulations for gas which has condensed out of the turbulent ICM.
lated multiphase condensation cascade in group and cluster haloes, showing that the simulated clouds/filaments can have significant internal velocity dispersions, which match those observed in clusters and groups.In Figure 8 we plot the velocity dispersion of the anomalous emission detected in NGC 0708 against its mean velocity shift, and include for comparison the data (brown points; extracted from Table 2 of Gaspari et al. 2018 and references therein) and simulations (green contours; 1, 2 and 3σ confidence intervals of the emission detected along individual lines-of-sight) from Gaspari et al. (2018).The anomalous emission we detect in NGC 0708 is well within the scatter of the observational data, and deviates from the expectation of the simulations at only the 1σ level.
Further evidence exists supporting a cooling-flow interpretation for this blue-shifted feature.For instance, the optical image of NGC 0708 shown in Figure 1 shows that the dust on larger scales is disturbed, and thus there may be other filaments in the system, which are not (yet) traced by molecular gas (although dust is expected to be destroyed on very short timescales in the hot ICM, so its presence here would require fast dust formation mechanisms; Clemens et al. 2010).
The high velocity dispersion of the molecular gas disc of NGC 0708 could be caused by material, such as the putative clump detected here, condensing out of the hot ICM and raining down on the disc.The turbulent Taylor number of this gas (defined as V rot /σ ) is shown as a function of radius in Figure 9.We estimate V rot (R) using our best-fit parameters from Table 1 in Equation 3, and σ (R) by binning our observed moment two within elliptical apertures.We note that formally these estimates of the turbulent Taylor number are lower limits, as observational effects (e.g.beam smearing, channelisation and the line-of-sight integration through the disc) preferentially lead us to overestimate the velocity dispersion.We highlight in grey the central regions of the galaxy where this effect is expected to be dominant.Outside of this region the turbulent Taylor number in the inner 100 -300 pc of NGC 0708 is around 2.5, which suggests that one indeed expects the disc to be bombarded by condensing clouds and rare extended filaments (see e.g.Gaspari et al. 2015).The outer regions have lower velocity dispersions, and thus chaotic and filamentary condensation is less likely to be important there.We note, however, that this is not the only possible interpretation of the pervasive high dispersion of the molecular gas disc; see Section 4.2.
Overall we conclude that the blue-shifted material observed in NGC 0708 could be a low-mass filament/cloud of material condensing out of the hot ICM and falling onto the core of the galaxy.

Outflow
The other option for the origin of the blue-shifted feature is an outflow.Outflows can be caused by supernova-and/or AGN-driven winds as well as by jets directly impacting onto the ISM.
Winds driven by supernovae tend to be large scale and are expected to be roughly bipolar, depending on the gas distribution around the star-forming region.To lead to the feature shown in Fig. 3, the wind would have to be very localised, or currently only interacting with a single (or at most a few) giant molecular cloud(s).The total star formation rate (SFR) of NGC 0708 is also very low, making this an unlikely scenario.
The blue-shifted anomalous emission is also offset from the AGN position as traced in 236 and 5 GHz continuum emission (≈ 0. 4 or ≈ 113 pc offset; see the left-hand panel .Median turbulent Taylor number (Ta t ) of the molecular gas disc in NGC 0708 (estimated within elliptical apertures and shown as blue points with error bars showing the 16th to 84th percentiles of the distribution).Formally our measurements are lower limits, due to observational effects such as beam-smearing, and hence are plotted as upward facing arrows.The central regions, where these observational effects are expected to be dominant, are shaded in grey.Dashed black horizontal lines denote the expectation for gas arising from a strong CCA rain (Ta t < 1), a mild rain (which forms a clumpy disc; 1<Ta t <5) and a system little affected by CCA rain (Ta t > 5) , as described in Gaspari et al. (2015).The central parts of NGC 0708 have a low turbulent Taylor number, consistent with bombardment by a mild CCA rain.
of Fig. 6), making it unlikely that it is a quasar wind-driven outflow from the central AGN.The off-centre position of the feature could indicate a binary black hole system, with a dual AGN.However, neither radio nor X-ray observations detect accretion onto a second SMBH at this position, setting a stringent upper limit on the accretion power available to drive an outflow.
AGN jets, on the other hand, are well collimated, strongly directional and can do work significantly away from the centres of galaxies, corresponding closely to the characteristics of the blue-shifted feature we observe.Jet-driven outflows do require the chance alignment of the jet with the ISM, but this is known to happen in a variety of sources (e.g.Alatalo et al. 2011;Aalto et al. 2012;Morganti et al. 2015;Fernández-Ontiveros et al. 2020).
As mentioned previously, an old AGN-driven jet is present in NGC 0708 at 330 MHz, but it is too large and not at the correct orientation to drive the putative outflow associated with the blue-shifted feature.We do detect signs of a jet at millimetre frequencies, pointing towards the putative outflow location (the Southern extension of the 236 GHz emission discussed in Section 3.2; see also Fig. 6).However, we do not see signs of this jet at 5 GHz suggesting that, if present, it must have a flat ratio spectrum (and is thus likely to be very young).The effects of shocks from such a small jet would be difficult to discern, as e.g.X-ray telescopes do not have the angular resolution necessary to separate emission at the putative hotspot from that of the nuclear point source.NGC 0708 is known to have launched jets repeatedly within short timescales (with a repetition time of ≈28 Myr), Figure 10.As Figure 5, but adapting the KinMS model to include an expanding shell of material launched at an off centre position (see Section 4.2).In the left panel the full disc plus outflow model is displayed as contours overlaid on the observed data, while in the right panel only the outflowing material is displayed.Clearly, an expanding jet cocoon that drives an expanding shell of molecular gas can explain the anomalous, primarily blue-shifted emission in NGC 0708, as the redshifted outflow is not easily identifiable in projection against the disc material.
and at different position angles (likely due to precession of the central accretion disc).Indeed, cavities are detected in the hot ICM of Abell 262 along the north-south direction (Clarke et al. 2009), suggesting the AGN axis was orientated in this direction in the past, and thus that a jet-driven explanation is plausible.
The blue-shifted feature we find in NGC 0708 is very similar to the jet-driven molecular gas outflow found in the Seyfert 2 galaxy IC 5063 by Morganti et al. (2015), who observed a similarly jagged PVD with large deviations from the underlying large-scale rotation pattern at one specific off-centre position.Fernández-Ontiveros et al. (2020) also discovered a molecular gas outflow in the Seyfert 2 galaxy ESO 420 G13, by detecting a velocity dispersion peak 440 pc from the galaxy centre, again similar to the one detected here (see the right-hand panel of Fig. 2).In these systems a young radio jet is driving into the molecular gas disc, inflating cocoons that drive the observed outflows, and causing the host galaxy discs to have high velocity dispersions (as also observed here; Mukherjee et al. 2018).While CCA and inflows are ubiquitous in BCGs, this suggests that an outflow explanation for this specific feature cannot be ruled out.

Putative outflow geometry
Given the similarity between the anomalous blue-shifted emission we observe in NGC 0708 and known molecular outflows, we consider the potential geometry such an outflow would have.Our observations only reveal a blue-shifted feature, while we would typically expect outflows to be (bi-)symmetrical.However, given that the launch velocity of the putative outflow is likely to be significantly away from the galaxy systemic velocity, the red-shifted side of the outflow is likely to be (at least partially) obscured by emission from the bulk of the rotating disc, as in some other similar sources (e.g.Fernández-Ontiveros et al. 2020).To illustrate this, in Figure 10 we add an off-centre outflowing expanding sphere of gas (as in e.g.Morganti et al. 2015) to the kinematic model of a rotating disc presented in Section 3.4.1.
To model an expanding shell on top of this circular rotation model.we begin with the output 'cloudlets' from KinMS (using the return clouds option within the code).For every 'cloudlet' within a (3D) radius R shell of the assumed outflow launch point P launch (x, y, z) we add a fixed velocity (V shell ) directed radially away from P launch (x, y, z).This creates an expanding bubble or shell of material, superimposed on the large-scale galaxy rotation pattern.We include these five free parameters (R shell , V shell and P launch (x, y, z)), and run an MCMC exploration of the parameter space as discussed in Section 3.4.1.The best-fitting outflow position P launch (x, y, z) is offset from the galactic centre by 0. 00 +0.06 −0.05 in RA and -0.35±0.07 in Dec and 0. 04 +0.03 −0.04 in the plane of the sky.The best fitting shell radius R shell =0.14±0.05 (or ≈ 33 pc, fully consistent with the size estimated using imfit on the data directly in Section 3.4).The velocity of the outflowing shell is estimated to be V shell =160 +14 −13 km s −1 , similar to that observed in known jet-driven outflows and simulations (e.g.Mukherjee et al. 2018).As shown in Figure 10, such an expanding shell can well reproduce the kinematics of NGC0708, even the red/blue asymmetry.

Putative outflow properties
In this Section we calculate various properties of the putative outflow, and use these to discuss the viability of this explanation for the anomalous blue-shifted emission in NGC 0708.We assume that the the outflow is a shell, where the outflow rate can be calculated as in Lutz et al. (2020): .
where ∆R OF is the thickness of the shell, M OF is the total outflow mass, and V shell is the average velocity of the outflow.
We use here the masses estimated in Section 3.4 and listed in Table 2 which (if this scenario is correct) correspond to gas contained within the blue-shifted outflow (M OF ).However, we caution that the actual mass outflow rate may be underestimated, likely by a factor of ≈ 2 if the geometry presented in Section 4.2 is correct (to account for the hidden red-shifted outflow).We conservatively adopt the estimate of ∆R out presented in Section 4.2.1 (e.g.∆R OF = 33 ± 8 pc).In this case the outflow is a filled bubble, rather than a shell, and mass outflow rates estimated will be lower limits if a thinner shell is present.We assume V shell = 160 +14 −13 km s −1 as found in Section 4.2.1.
To calculate the kinetic power of the outflow (P kin,OF ), we use Equation 7of Holt et al. (2006) rescaled to CO(2-1) from [O iii] by Morganti et al. (2015, see their Eq. 1).Following both papers, we assume the relatively large line width of the outflowing gas reflects turbulent motion, so that the FWHM of the CO line represents the turbulent component of the outflow: where assume v turb ≈ FWHM ≈ 2.355 σ v,los , and other terms are as defined previously.We note that including the turbulent term in this equation is controversial, but in our case this term is sub-dominant.We include it here to allow direct comparisons to other works.
The final property we calculate is the depletion time for the gas reservoir in NGC 0708, assuming that the putative outflow continues at the same rate, via τ dep = M tot / .M OF .Table 3 lists the mass outflow rate, kinetic power and depletion times derived for each assumed α CO using the above equations.We derive mass outflow rates between ≈7 and 104 M yr −1 , and kinetic powers between ≈10 41 and 10 43 erg s −1 .The entire molecular medium would be depleted by such an outflow in between ≈22 and 309 Myr.The assumption of a given α CO is the dominant (systematic) uncertainty.We note that in other outflows the gas is typically (although not always) optically thin (Lutz et al. 2020), and thus we consider the lower power and outflow rate and longer depletion time estimates more likely.
To ascertain if the jets launched by NGC 0708 have enough energy to power the outflow, we calculate its jet power (Q jet ).Equation 11of Wu (2009) converts the radio luminosity at 151 MHz (L 151 ) to Q jet (we follow their analysis and use a normalisation factor f = 10).NGC 0708 was observed as part of the 6th Cambridge (6C) survey at 151 MHz at a resolution of ≈ 7. 2 × 7. 2 (≈ 120 × 120 kpc 2 ; Baldwin et al. 1985).The catalogue reports a continuum flux L 151 = 0.78 ± 0.075 Jy (Hales et al. 1993).The beam of these observations covers the whole of the old, large-scale jet in a spatially-unresolved manner, so we are forced to assume that an AGN restart would have produced a jet of similar power.
This assumption yields Q jet = (1.32 ± 0.01) × 10 43 erg s −1 , which would require a coupling factor with the ISM of 0.5 -6 percent (depending on α CO and the outflow geometry) to cause the outflow identified here.The efficiencies for an optically thin outflow are consistent with those seen in simulations (which yield jet-ISM energy transfer efficiencies of < ∼ 1 percent; e.g.Nesvadba et al. 2010;Wagner & Bicknell 2011) and in known jet-driven molecular outflows (e.g.Alatalo et al. 2011Alatalo et al. , 2015;;Davis et al. 2012;Aalto et al. 2012;Morganti et al. 2015;Fernández-Ontiveros et al. 2020).
Overall we conclude that, although inflow onto the BCG is the most likely explanation for the blue-shifted anomolous feature seen in NGC 0708, an jet-driver molecular outflow could also explain these observations.

SMBH mass
The original goal of the observations presented here was to estimate the mass of the SMBH in NGC 0708.As shown in Figs. 1 and 2, the gas in the galaxy is disturbed (especially in the outer regions), and it does not lie in the equatorial plane of the galaxy, making this difficult.Despite this, the kinematics of the gas in the very centre of NGC 0708 (around the AGN/continuum source seen; see Figs. 2 and 6) seem fairly regular.Given the short dynamical times in this region, it is possible that this gas is sufficiently relaxed to allow us to constrain the central potential reasonably accurately.
From the red-shifted side of the PVD, that appears fairly undisturbed by the outflow (see Fig. 3), we can make a crude estimate of the total mass enclosed within the inner resolution element of our data.At 0. 088 (≈24.8 pc; one synthesised beam major axis from the galaxy centre) the projected rotational velocity is ≈ 180 km s −1 , yielding an enclosed mass M enc = 2.15 × 10 8 M (assuming pure rotation in a spherical potential, so M enc = v 2 (r)r/G, where v(r) is the de-projected rotational velocity at radius r and G is the gravitational constant).The total molecular gas mass within this radius is ≈ 5.1 × 10 6 M , and from the HST F110W image we can estimate a stellar mass within the same radius of ≈ 3 × 10 6 M (assuming a very conservative F110W-band mass-to-light ratio of 2; see e.g.Fig. 11 of Balogh et al. 2001).This suggests a total dark mass of ≈ 2 × 10 8 M within 25 pc of the centre of NGC 0708.
While this estimate is very uncertain due to the unknown degree of kinematic disturbance in the gas (and the approximate stellar mass-to-light ratio, absence of dust correction and standard CO-to-H 2 conversion factor), it is consistent with the SMBH mass estimated from the black hole mass -central stellar velocity dispersion relations of Woo & Urry (2002) and Donato et al. (2004).
The Eddington luminosity of a black hole of this mass is L Edd ≈ 2.67 × 10 46 erg s −1 .In comparison, Clarke et al. (2009) estimated the total AGN kinetic luminosity to be L AGN, kin = 6.2 × 10 42 erg s −1 in NGC 0708.As a percentage of the Eddington luminosity, this suggests the SMBH in NGC 0708 is currently accreting at only ≈ 0.023 percent of its Eddington rate.

CONCLUSIONS
In this work we have presented 12 CO(2-1) line and 236 continuum ALMA observations (along with 5 GHz e-MERLIN continuum imaging), of the early-type galaxy NGC 0708, the BCG in the galaxy cluster Abell 262.Our observations reveal a turbulent, rotating disc of molecular gas in the core of this galaxy.A marginally spatially-resolved 'spike' of blueshifted anomalous emission with a large line width is also present, approximately 100 pc away from the galaxy nucleus.The (sub-)millimetre continuum emission peaks at the nucleus, but shows an extension towards this anomalous CO emission feature.No corresponding elongation is found on the same spatial scales at 5 GHz.
The central kinematics of NGC 0708 allow us to roughly estimate the non-luminous mass contained in the inner 25 pc as ≈ 2 × 10 8 M , although given the somewhat disturbed kinematic state of the gas this estimate must be treated with caution.This value is, however, consistent with the SMBH mass expected for this source based on the black hole mass -central stellar velocity dispersion relation.
We considered multiple potential causes for the anomalous blue-shifted emission detected in NGC 0708, and conclude that two explanations are viable.X-ray observations of Abell 262, the host cluster of NGC 0708, show that it is expected to have a turbulent and quenched cooling flow.Thus this feature could be a molecular filament that has condensed out of the hot ICM, and is falling onto the galaxy (consistent with predictions from CCA simulations).The fact that only a single such putative low-mass in-falling structure is detected in NGC0708 is consistent with it being located in a low-mass, low ICM cooling rate cluster.
While an inflowing filament explanation seems most likely, an alternative explanation is a molecular outflow, likely powered by a young, restarted radio jet.The kinematics of the anomalous emission in NGC 0708 can be explained by an expanding shell/bubble of molecular gas, as seen in a variety of other galaxies with known jet-driven molecular outflows.A jet-driven scenario is also consistent with the observed episodic nature of the AGN in NGC 0708, suggesting the outflow we are currently observing is young.This explanation would require the coincidence of the radio jet with the molecular disc, but the extension of the millimetre continuum provides tentative evidence of this, as long as the young jet has an inverted/flat synchrotron spectrum.
To fully resolve the nature of this anomalous blueshifted emission it would help to map the ionised gas at high spatial-resolution (e.g. with adaptive optics assisted integralfield spectroscopy), to see if the nebular emission in this system is consistent with cooling, or an outflow.Shock tracers (such as SiO) could also be observed with ALMA, to verify if a jet is impacting the molecular gas, and if so at what location.In either case, this source provides another intriguing place to study the interaction between cooling flows and AGN mechanical feedback.

Figure 1 .
Figure 1.Top: Large-scale (50 ×50 or ≈ 14 × 14 kpc 2 ) HST Wide Field Camera 3 F110W image of NGC 0708, with 330 MHz continuum Very Large Array (VLA) contours from Clarke et al. (2009) overlaid in blue.Bottom: Small-scale (12 ×12 or 3.2 × 3.2 kpc 2 ) HST combined Advanced Camera for Surveys and Wide Field Camera F435W image of NGC 0708, with our CO(2-1) integrated intensity (0. 088 or ≈25 pc resolution) contours overlaid in blue and 236 GHz continuum contours overlaid in magenta.The synthesised beam of the radio/millimetre data is shown in the bottom-left corner of each image.

Figure 3 .
Figure3. 12 CO(2-1) kinematic major-axis position-velocity digram of NGC 0708, taken at a position angle of 349 • .The grey dashed line denotes the systemic velocity, V sys = 4750 km s −1 .The grey dot-dashed line denotes the velocity cut used to isolate the blue-shifted feature (V obs = 4550 km s −1 ; see Section 3.4), itself indicated by the cyan polygon.

Figure 4 .
Figure 4. 12 CO(2-1) integrated spectrum of NGC 0708 smoothed to an effective resolution of 20 km s −1 , showing the characteristic double-horned shaped of a rotating disc.The anomalous blueshifted wing is indicated by the shaded grey region.

Figure 5 .
Figure5.Left and Centre: PVDs extracted from the NGC 0708 CO datacubes along the kinematic major (349 • ; left column) and minor axis (259 • ; central column) of the system.The observed CO data is shown as orange filled contours, while the best-fitting KinMS model of a rotating disc is overlaid in black contours.Right: Major axis PVD created from the residual (Data-model) cube (orange filled contours).Only positive residuals are shown due to the masking procedure used, but these are much larger than the negative residuals.The same disc model as present in the left panel is reproduced here, to guide the eye.This panel shows the full extent of the anomalous non-rotating emission in this system.

Figure 6 .
Figure 6.Left: 12 CO(2-1) moment 2 (intensity-weighted observed line-of-sight velocity dispersion) map of NGC 0708, with 236 GHz continuum emission isophotes overlaid in magenta.The magenta arrows indicate the direction of the large-scale jet traced by 330 MHz emission.The extension of the 236 GHz continuum emission matches well with the position of the velocity dispersion peak.Right: 12 CO(2-1) spectrum integrated over the spatial extent of the off-centred velocity dispersion peak (indicated by the black box in the left panel), with gas at anomalous velocities (V obs < 4550 km s −1 ) indicated in magenta.Gaussian fits are shown in blue for the anomalous component, orange for the disc emission, and red for the sum of the two components.The anomalous emission is significantly blue-shifted from the galaxy systemic.
Figure8.Velocity dispersion of molecular/ionised gas along individual lines of sight ('pencil beams') towards BCGs, plotted against the mean velocity shift of their line centroid from the host system's systemic velocity.The brown datapoints are reproduced from Table2ofGaspari et al. (2018), while our estimates from NGC 0708 are shown as an orange star.Green contours show the 1 to 3σ range in these quantities seen in the simulations ofGaspari et al. 2018.The anomalous emission in NGC 0708 has a significant internal velocity gradient, which is consistent at 1σ with the expectation from CCA simulations for gas which has condensed out of the turbulent ICM.
Figure9.Median turbulent Taylor number (Ta t ) of the molecular gas disc in NGC 0708 (estimated within elliptical apertures and shown as blue points with error bars showing the 16th to 84th percentiles of the distribution).Formally our measurements are lower limits, due to observational effects such as beam-smearing, and hence are plotted as upward facing arrows.The central regions, where these observational effects are expected to be dominant, are shaded in grey.Dashed black horizontal lines denote the expectation for gas arising from a strong CCA rain (Ta t < 1), a mild rain (which forms a clumpy disc; 1<Ta t <5) and a system little affected by CCA rain (Ta t > 5) , as described inGaspari et al. (2015).The central parts of NGC 0708 have a low turbulent Taylor number, consistent with bombardment by a mild CCA rain.

Table 1 .
Model parameters with their priors, best-fitting values and statistical uncertainties.

Table 2 .
Properties of the blue-shifted anomalous emission.CO , optically-thick ULIRG α CO and local/Milky Way α CO .Statistical uncertainties are quoted at 1σ .