Context. Giant gt 100 kpc nebulae associated with active galaxies provide rich information about the circumgalactic medium around galaxies its link with the interstellar medium of the hosts and the mechanisms involved in their evolution. Aims. We have studied the giant nebula associated with the Teacup z  0.085 quasar based on VLT MUSE integral field spectroscopy to investigate whether the wellknown giant 10 kpc active galactic nucleus AGN induced outflow has an impact on the distribution of heavy elements in and outside the host galaxy. Methods. We have mapped the oxygen and nitrogen gas relative abundances OH and NO in two spatial dimensions across the giant nebula and within the galaxy by means of comparing emission line ratios with photoionisation model predictions. Results. The widely studied AGNdriven outflow responsible for the 10 kpc ionised bubble is enhancing the gas metal abundance up to 10 kpc from the AGN. OH is solar or slightly higher at the edges of the bubble in comparison with the subsolar abundances across the rest of the nebula median OH 0.63 OHSUBSUB. Conclusions. AGN feedback can produce metal enrichment at large extranuclear distances in galaxies 10 kpc.
['galaxies: abundances', 'galaxies: active', 'galaxies: evolution', 'quasars: individual: Teacup', 'Astrophysics - Astrophysics of Galaxies']
AGN feedback can produce metal enrichment on galaxy scales
{'AGN feedback can produce metal enrichment on galaxy scales': 'M. Villar Martín 1 , C. López Cobá 2 , S. Cazzoli 3 , E. Pérez Montero 3 , A. Cabrera Lavers 4 , 5 ⋆  \n- 1 Centro de Astrobiología (CAB), CSIC-INTA, Ctra. de Ajalvir, km 4, 28850 Torrejón de Ardoz, Madrid, Spain\n- 2 Institute of Astronomy and Astrophysics, Academia Sinica, No. 1, Section 4, Roosevelt Road, Taipei 10617, Taiwan\n- 3 Instituto de Astrofísica de Andalucía, CSIC, Apartado de correos 3004, E-18080 Granada, Spain\n- 4 GRANTECAN, Cuesta de San José s/n, 38712 Breña Baja, La Palma, Spain\n- 5\n- Instituto de Astrofísica de Canarias, Vía Láctea s/n, 38200 La Laguna, Tenerife, Spain e-mail: villarmm@cab.inta-csic.es', 'ABSTRACT': 'Context. Giant ( > 100 kpc) nebulae associated with active galaxies provide rich information about the circumgalactic medium (CGM) around galaxies, its link with the interstellar medium (ISM) of the hosts and the mechanisms involved in their evolution.  \nAims. We have studied the giant nebula associated with the Teacup ( z =0.085) quasar based on VLT MUSE integral field spectroscopy to investigate whether the well known giant ( ∼ 10 kpc) active galactic nucleus (AGN) induced outflow has an impact on the distribution of heavy elements in and outside the host galaxy.  \nMethods. We have mapped the oxygen and nitrogen gas relative abundances (O/H and N/O) in two spatial dimensions across the giant nebula and within the galaxy by means of comparing emission line ratios with photoionisation model predictions.  \nResults. The widely studied AGN driven outflow responsible for the ∼ 10 kpc ionised bubble is enhancing the gas metal abundance up to ∼ 10 kpc from the AGN. O/H is solar or slightly higher in the bubble edges, in comparison with the subsolar abundances across the rest of the nebula median (O/H ∼ 0.63 (O/H) ⊙ ).  \nConclusions. AGN feedback can produce metal enrichment at large extranuclear distances in galaxies ( ≥ 10 kpc).  \nKey words. galaxies: evolution - galaxies: active - galaxies: abundances - quasars: individual: Teacup', '1. Introduction': "The circumgalactic medium (CGM) is the bound gas halo surrounding galaxies outside their interstellar medium (ISM) and inside their virial radius, extending out to a few hundred kpc. The CGM may be the key regulator of the galactic gas supply. Gas flows occurring between the CGM and the interstellar medium (ISM) are thought to shape galaxies and drive their evolution, via feedback, accretion and recycling of gas. Thus, investigating the metallicity, structure and kinematics of the different gas phases will help us to understand how galaxies gain, eject and recycle the gas during their existence (see Tumlinson et al. 2017 for a review).  \nBecause the CGM is very diffuse and therefore almost invisible in emission, its physical properties remain largely unconstrained. Its understanding has come so far mostly from studies based on absorption lines produced by the CGM around galaxies, imprinted on the spectrum of background objects such as quasars and radio galaxies. The CGM of a few non-active nearby galaxies has also been studied in emission (e.g. Hayes et al. 2016), but the procedure is extremely challenging. The presence of a powerful active galactic nucleus (AGN) can render the CGM gas observable in emission around galaxies up to many 10s of kpc, well into the CGM. Giant emission line nebulae (size > ∼ 60 kpc and sometimes > 100 kpc) associated with quasars  \nand radio galaxies at different redshifts have been studied since the 80's (e.g.Baum et al. 1988; McCarthy et al. 1990; McCarthy 1993; van Ojik et al. 1996; Villar-Martín et al. 2003; Borisova et al. 2016; Villar-Martín et al. 2018; Fossati et al. 2021; Balmaverde et al. 2022; Wang et al. 2023). The discovery of giant ( ≥ 70 kpc) widely spread reservoirs of molecular gas (e.g. Emonts et al. 2016; Falkendal et al. 2021) associated with several high z ( z > ∼ 2) radio galaxies reveals a multiphase CGM that was chemically enriched when the Universe was 3 Gyr old and supplies gas from which galaxies grow.  \nWide field integral field spectroscopic instruments on 810 meter telescopes, such as the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT), have opened excellent opportunities to detect and study in great detail the elusive material from the warm ( T < 10 5 K) ionised CGM around powerful active galaxies at different redshifts thanks to the illumination by the powerful AGN. If it were not for the excitation by the AGN continuum, this gas might not be detected or would only be observable through absorption line studies.  \nThe object of this work is the well known Teacup radio quiet type 2 quasar at z = 0.085. It shows a ∼ 10 kpc loop of ionised gas resembling a handle of a teacup (hereof its nickname), which was discovered by volunteers of the Galaxy Zoo project (Keel et al. 2012) and has been widely studied in the context of AGN feedback and its potential impact on galaxy evolution. The system has been proposed to be the scenario of a giant outflow generated either by an  \nAGN wind or induced by a 1 kpc radio jet whose effects are noticed up to at least ∼ 10 kpc from the AGN and might be responsible for the bubble-like morphology (Gagne et al. 2014; Harrison et al. 2015; Keel et al. 2015; Ramos Almeida et al. 2017; Villar-Martín et al. 2018; Moiseev & Ikhsanova 2023; Venturi et al. (2023, hereafter V23).  \nWe present in this paper a detailed optical spectroscopic study of the Teacup nebula based on VLT-MUSE archival data with the main goal of mapping the gas abundances in two spatial dimensions. The ultimate science goal is to establish whether the giant outflow has an impact on the distribution of heavy elements from the nucleus on large spatial scales. We will discuss the results in the context of other studies of the CGM and its role in the evolution of galaxies.  \nVillar-Martín et al. (2018) discovered a > 100 kpc ionised nebula associated with this object (see also Villar Martín et al. 2021, which could be a product of a merger that occurred 1-2 Gyr ago (Keel et al. 2015). This rich gas reservoir, which extends into the CGM, has been rendered visible due to the activity of the quasar nucleus. AGN photoionisation dominates the excitation of the spatially resolved gas emission up to its boundary (Gagne et al. 2014; Villar-Martín et al. 2018; V23; Moiseev & Ikhsanova 2023), except at some locations outside the putative quasar ionisation cones where evidence for shock excitation has been found (V23). Stellar photoionisation could also contribute to ionise the gas locally in some tidal features. The large scale kinematics are strongly reminiscent of rotation (Villar-Martín et al. 2018; Moiseev & Ikhsanova 2023) and tentative results suggest subsolar nebular abundances ( ∼ 0.5 Z ⊙ , Villar-Martín et al. 2018). The well known bubble appears to be expanding from the nucleus and out into the nebula.  \nThe paper is organised as follows: the data and analysis method are presented in Sect. 2 and 3 respectively. The latter includes the description of the methods used to derive the gas chemical abundances and physical properties of the gas. The results are presented in Sect. 4 and discussed in Sect. 5. The main conclusions are summarised in Sect. 6.  \nThroughout this paper, we assume flat Λ CDM cosmology following Planck Collaboration et al. (2020), with H 0 = 67.4 km s -1 Mpc -1 , Ω m = 0.31. This gives a spatial scale of 1.65 kpc arcsec -1 at z = 0 . 085 .", '2. Data': "The data were collected for the 0102.B-0107 program (principal investigator, PI: L. Sartori; see V23 for details) with the European Southern Observatory (ESO) Very Large Telescope (VLT) and the Multi Unit Spectroscopic Explorer (MUSE, Bacon et al. 2010). This instrument covers a 1 ' × 1 ' field of view (FoV) in the Wide Field Mode (WFM), with a spatial sampling of 0.2 '' pix -1 . The wavelength coverage is ∼ 4650-9300 Å with a 1.25 Å pix -1 spectral sampling and a resolving power R = λ / ∆ λ ∼ 1700-3400 ( ∆ V ∼ 176-88 km s -1 ).  \nThe observations were performed in March 2019. The processed archive data cube was used for this study. The Teacup ionised nebula fills a large fraction of the MUSE FoV. Separate sky cubes are not available from these observations. Sky oversubtraction due to contamination of the sky spectrum by object emission was identified at a few spatial locations of very low surface brightness which appeared as artefacts that mimic absorption features adja-  \ncent to the strongest emission lines (specially [OIII] λ 5007). Because our attempts to improve the sky subtraction did not achieve significant improvements, we finally used the archive cube. To further evaluate the potential impact of this effect, a comparison with the processed archive datacube from program 0103.B-0071 (PI: C. Harrison) was also performed. Although separate sky cubes were obtained for this program, the resulting sky subtraction was not significantly better. The comparison was in any case valuable to characterise the impact of the sky subtraction on the data quality. The artefacts are in general faint in comparison with the emission lines and are shifted in wavelength so that they could be efficiently identified and isolated. We confirm that imperfections on the sky subtraction do not affect our results and conclusions. These are all based on the 0102.B-0107 datacube because of the significantly better seeing (FWHM=0.74 '' vs. 1.3 '' for the 0103.B-0071 program) and of somewhat better signal-to-noise ratio (SN).", '3.1. Spatially resolved emission line flux measurements': "Our main goal is to map the electron temperature, T e , and the oxygen abundance O/H in two spatial dimensions to investigate whether the giant outflow has an impact on the distribution of heavy elements across the galaxy and out into the CGM. For the T e determination, it is essential to measure the flux ratio of a nebular to an auroral emission line, such as the [OIII] λ 5007/ λ 4363 (or [NII] λ 6583/ λ 5755). For this, we defined spatial apertures at different locations through the nebula and extracted the integrated spectrum from each one.  \nThe apertures were selected based on the visual inspection of the [OIII] λ 5007 morphology at different wavelengths (i.e. velocities) scanned through the line profile. These scans reveal striking morphological changes with velocity and are specially useful to identify some faint structures, such as tidal tails, arcs, knots, filaments, which have similar velocities (Fig. 1). By restricting the central λ and the width, ∆ λ , of the narrow band images, the contrast of the lowest surface brightness features is enhanced and the spatial apertures covering them can be defined with greater accuracy to maximise the SN of the integrated spectra. This was useful to detect faint emission lines through the largest possible extension, while preserving at the same time the morphological information on the diversity of nebular structural elements.  \nThe map of the positions of the resulting 64 selected apertures is shown if Fig. 2 (left panel). We also show in adjacent panels of the same figure the W 80 and V s maps based on [OIII] λ 5007. W 80 is the velocity width that encloses 80% of the total line flux and V s is the velocity shift relative to the narrow core of the nuclear [OIII] line, considered here as indicator of the systemic velocity (Greene & Ho 2005). These maps can be directly compared with those of V23, that were based on a spaxel by spaxel analysis. The huge dimensions of the nebula both along and perpendicular to the radio axis are apparent. Although some prominent features stand out (knots, filaments, etc), gas emission seems to fill the entire area within the nebular outer boundaries. The large scale rotation pattern (Gagne et al. 2014; Villar-Martín et al. 2018; Moiseev & Ikhsanova 2023) and  \n-1  \nFig. 1. [OIII] continuum subtracted images covering different spectral windows that were selected to highlight the diversity of nebular morphological features. Each image covers a different velocity (i.e. spectral) range relative to the nuclear systemic velocity as indicated on top. The nebular morphology strongly varies with velocity. The left panel shows the total [OIII] flux narrow band image. The well known ∼ 10 kpc ionised bubble is marked with a tiny yellow star in the left panel. The green lines indicate the position angles of the radio axis to the NE and to the W (Harrison et al. 2015). To guide the reader, the letters 'A' to 'F' mark some emission line features that can also be identified in the mask map of Fig. 2.  \n<!-- image -->  \nFig. 2. Map of the masks used in our analysis (left) and [OIII] kinematic maps (middle and right panels). The colours in the first map have no particular meaning but help differentiating the apertures. A 1D spectrum was extracted from each one, so that a single W 80 and V s (middle and left panels) values are associated to each aperture. V s is the velocity shift relative to the narrow core of the nuclear [OIII] line. The maps cover the total MUSE FoV ( ∼ 1 ' × 1 ' ). W 80 and V s are in km s -1 . Letters 'A' to 'F' , yellow star and solid lines have the same meaning as in Fig. 1.  \n<!-- image -->  \nthe enhanced width, W 80 , of the lines in the direction perpendicular to the radio axis are clear (V23).  \nThe 1D spectra extracted from the individual apertures were used to measure the fluxes of the lines involved in our analysis: [OIII] λλ 4959,5007 and λ 4363, H γ , H β , H α , [NII] λλ 6548,6583, and λ 5755 and [SII] λ 6716,6731.  \nFor a subset of spectra the stellar continuum was relatively strong compared with the emission line fluxes, that is, the equivalent width of the lines was low, so that it was first necessary to fit it and subtract it. The effects of stellar absorption can affect the Balmer emission lines and need to be accounted for to obtain accurate reddening estimates. It was also necessary in several spectra to reconstruct the H γ and [OIII] λ 4363 baseline and recover both line fluxes. The continuum was fitted for each spectrum with the pypipe3d tool (e.g., Lacerda et al. 2022), that performs a decomposition of the observed stellar spectra into multiple simple stellar populations, each of different age and metallicity. As a test, we also subtracted the continuum with the method described in Cazzoli et al. (2022). The results were consistent.  \nWe also checked the reddening inferred from H γ /H β and H α /H β in the final, corrected spectra. Case B H α /H β =2.87 and H γ /H β =0.47 were assumed. Reddening estimations from both ratios are in general consistent within the errors. Changing these values within the range H α /H β =2.76-3.05 and H γ /H β =0.45-0.47 (Osterbrock 1989) has no significant impact. The E B -V values implied by H α /H β were used for extinction correction since they have smaller errors. Using instead the E B -V implied by H γ /H β has also a negligible impact on the results.", '3.2. Derivation of temperatures and gas chemical abundances': 'We used the python code PyNeb (Luridiana et al. 2015) to estimate T e of [OIII] for spectra with detected [OIII] λ 4363 using [OIII] λ 5007/ λ 4363, and, when possible, of [NII] using [NII] λ 6583/ λ 5755. This code was also used to infer the electron density n e using [SII] λλ 6717,6731. We calculated the corresponding errors by applying to the calculations a  \nMonte Carlo iteration using the nominal fluxes perturbed with the observational errors.  \nFor the derivation of abundances, the direct method (i.e. the derivation of total chemical abundances from their ionic fractions and the measured electron temperature) cannot be satisfactorily applied for the case of the Narrow Line Regions (NLR) in Active Galactic Nuclei (AGN) (e.g. Dors et al. 2015). Instead, we used the code HII-CHI-mistry (hereafter HCm , Pérez-Montero 2014). The approach of this code consists of a bayesian-like comparison between certain reddening-corrected ratios of emission line fluxes relative to a recombination H line (H β in our case), that are sensitive to total oxygen abundance (O/H), nitrogento-oxygen ratio (N/O), and ionisation parameter ( U ), with the predictions from a large grid of photoionisation models. In particular, we used version 5.3 of HCm , which considers models calculated under the most usual conditions in the NLR of active galaxies. As a result, the most probable values and their uncertainties for O/H, N/O and log ( U ) are obtained.  \nThe grid of models is explained in detail in PérezMontero et al. (2019). The gas is assumed to be distributed homogeneously with a filling factor of 0.1 and a constant electron density of n e =500 cm -3 . The densities across the Teacup nebula are however significantly lower ( n e < ∼ 150 cm -3 , see Table .1 and V23). Since collisional de-excitation effects are not relevant for the line ratios under consideration in this density regime, the inferred abundances are not affected by the comparatively high density used in the models (Pérez-Montero et al. 2019).  \nThe Spectral Energy Distribution (SED) was considered to be composed by two components: one representing the Big Blue Bump peaking at 1 Ryd, and the other a power law with spectral index α OX =-1.2 representing the nonthermal X-rays radiation (see Pérez-Montero et al. 2019 for additional details). The stopping criterion to measure the resulting emergent spectrum is that the proportion of free electrons in the ionised gas is lower than 2%. The models consider a Mathis et al. (1977) grain size distribution and a dust-to-gas mass ratio of 7.5 × 10 -3 (Rémy-Ruyer et al. 2014).  \nThe grid of models vary the oxygen and nitrogen abundances within 12+log(O/H)=6.9-9.1 in bins of 0.1 dex and log(N/O)=-2.0-0.0 in bins of 0.125 dex, respectively. The rest of chemical abundances are scaled in the models with respect to oxygen following the solar proportions, 12+log(O/H)=8.69 ± 0.04, given by (Asplund et al. 2021). In accordance with these authors, we also assumed 12+log(N/H)=7.83 ± 0.07 for the Sun.  \nThe original grid of models in Pérez-Montero et al. (2019) covers values of log ( U ) from -4.0 to -0.5 in bins of 0.25 dex, although we only assumed log ( U ) > -2.5 to break the degeneracy with U for the ratios of high- to lowexcitation lines in the optical range predicted by the models. This restriction does not imply any significant change in the derivation of the chemical abundances (see PérezMontero et al. 2019 for a better clarification). The code also provides uncertainties to the derived quantities calculated as the quadratical addition of the standard deviation of the obtained bayesian distribution with the uncertainty obtained from a Monte Carlo iteration through the given input nominal flux of each line perturbed with its corresponding observational error.', '4. Results': 'We show in Table .1 the O/H and N/O abundances predicted by photoionisation models for different apertures, as well as the [OIII] and [NII] electron temperatures inferred with pyNeb for the densities estimated with [SII] λλ 6716,6731 (see Sect. 3.2). The observed extinction corrected line ratios relative to H β and measured with the 1dim spectra are also shown. This analysis was performed for a set of 23 (including the nucleus) apertures where [OIII] λ 4363 has been detected (see left panel in Fig. 3). The top (lower) part of the table comprises the apertures in which [NII] λ 5755 is detected (undetected). This does not affect the estimated chemical abundances, but simply indicates whether T [NII] could be inferred or not. The O/H, N/O and T [OIII] information in Table .1 is also shown as 2dim maps in Fig. 3 and 4.  \nThe main result which stands out in Fig. 3 is that the giant ionised bubble shows significantly higher O/H and lower T [OIII] than the gas at most locations across the giant nebula. This is confirmed in all apertures along the bubble edge, where 12 + log (O/H) is in the range 8.68 ± 0.098.77 ± 0.12. As in the nucleus (8.66 ± 0.11), these abundances are consistent with the solar value or slightly higher.  \nFor the rest of the nebula (this is, the extended gas outside and beyond the bubble), except in Ap. 13, the oxygen abundances are subsolar everywhere (as tentatively found by Villar-Martín et al. 2018). Considering all apertures across the nebula, the log of the median is 12+ log (O/H)=8.49 (or 63% solar), with values as low as 8.37 ± 0.08-8.44 ± 0.13 (48%-56% solar) in Ap. 15, 16, 19, 21, 22. Interestingly, the gas encircled by the bubble edges (Ap. 11), with 12+ log (O/H)=8.53 ± 0.17 (69% solar), shows an intermediate abundance between the most enriched gas (the bubble edges) and the more metal poor gas in the rest of the nebula.  \nIt is important to highlight that the abundance predictions are based on the assumption that the gas is photoionised by the AGN everywhere. This is reasonable for most apertures except, possibly, at a few positions located in the direction perpendicular to the main axis of the nebula, where shocks might be contributing to the excitation of the gas (V23). These apertures are labelled Ap. 12, 18, 19 in Fig. 3, left panel. In spite of this uncertainty, the fact that (except for Ap. 13) this gas follows the general behaviour of the O/H and T [OIII] maps suggests that the derived abundances are reliable.  \nIn general, the T [OIII] map (right panel in Fig. 3) mimics the behaviour of the O/H maps in an anticorrelation. The nucleus shows the minimum temperature, T 4 = T [OIII] /10 4 =1.27 ± 0.08 K, in comparison with the rest of the gas, which shows T 4 > 1.5 everywhere. As before, the bubble stands out in this map, with the edges being significantly colder ( T 4 ∼ 1.5-1.6 in Ap. 3 to 9) than the giant nebula, which is very hot with T 4 ∼ 1.7-1.9 almost everywhere. V23 presented a 2dim T [OIII] map tracing the nucleus and the ionised bubble. They inferred a narrow range of significantly lower T [OIII] ∼ (1.3-1.4) ∼ 10 4 K compared to our estimations. The reason for this discrepancy remains unknown.  \nIn regard to the N/O abundance (Fig. 4 and Table .1), the most obvious result is the much higher nuclear value ( log (N/O)=-0.56 ± 0.08, i.e. roughly twice the solar abundance ratio of -0.86 ± 0.07), compared with the so-  \nFig. 3. Maps of (O/H)/(O/H) ⊙ ratio and T [OIII] . The masks used in this analysis are shown in the left panel. They correspond to the apertures for which [OIII] λ 4363 is detected. Notice that the FoV is smaller than in Fig. 1 and 2. The exact values of O/H and T [OIII] are in Table .1. T [OIII] is in units of 10 4 K.  \n<!-- image -->  \nFig. 4. Map of (N/O)/(N/O) ⊙ . The FoV is the same as in Fig. 3. The nuclear value is 1.99, which is outside the colour bar range to enhance the contrast for visualisation purposes. The (N/O)/(N/O) ⊙ values are in Table .1.  \n<!-- image -->  \nr or somewhat below solar N/O everywhere else (-1.03 ± 0.09 ≤ log (N / O) ≤ -0.84 ± 0.20), including the bubble edge. There is not such a striking trend of N/O with location as found for O/H. The bubble edge shows among the lowest N/O ratio ( ∼ 68-85% N/O ⊙ ), although not unique in comparison with other regions of the giant extended nebula. The bubble edge N/O values are lower than the gas encircled by it (ap. 11), the nucleus and the gas in between (ap. 1 and 2).  \nT [NII] could be measured only for the nucleus and up to the bubble edge (Table .1). The temperature difference between these two regions is less pronounced than for T [OIII] . T 4 [NII] ∼ 0.983 ± 0.066 K in the nucleus and somewhat higher in the bubble ( ∼ 1.01-1.12 K).', '5. Discussion': 'With a maximum extension as traced by the MUSE data of ∼ 126 kpc (this is a lower limit since the gas fills the FoV in this direction), the Teacup giant nebula traces part of the CGM. Its properties are strongly conditioned by the AGN, but it still provides valuable information about the CGM. If it was not for the nuclear activity, most (if not all) of this gas would remain invisible.  \nThe giant nebula shows subsolar abundances, with O/H ∼ (48%-84%) × (O/H) ⊙ almost everywhere, and median 63%, well below the nuclear, roughly solar abundance. For comparison, different works based on absorption line studies have shown that the dense gas in the CGM of z < ∼ 1 galaxies has a bimodal metallicity distribution function, with an equal number of absorbers in the low-metallicity ( Z < ∼ 0.03 Z ⊙ ) and high-metallicity ( Z ∼ 0.4 Z ⊙ ) branches (e.g. Lehner et al. 2013; Wotta et al. 2016). The abundance of the Teacup ionised nebula falls in the latter group (clearly, this does not discard the presence of lower metallicity gas). The high-metallicity branch has been proposed to trace galactic winds, recycled outflows, and tidally stripped gas (Lehner et al. 2013).  \nLower CGM metallicities compared with the ISM have been found for different galaxy types, based on studies of absorption line systems at < ∼ 200 kpc from their host galaxies (e.g. Kacprzak et al. 2019). These authors found a large range of abundances in the CGM of isolated star forming galaxies, 0.01 < ∼ Z/Z ⊙ < 1 and an offset of log ( dZ ) =-1.17 ± 0.11 between the CGM and ISM, which shows no dependence with stellar mass. The relation has a large scatter of 1 σ =0.72. The offset for the Teacup is smaller for O/H ( log ( dZ ) ∼ -0.17 using the median nebular O/H), although still within the scatter. Different processes may be at work in this system, related to the nuclear activity and/or to its merger history.  \nWe have shown that the ∼ 10 kpc Teacup bubble to the East of the nucleus, that is known to be driven by an AGN wind or the small nuclear radio jet, is associated with obvious changes in the gas abundance. The bubble edge shows significantly higher O/H (solar or slightly solar, similar to the nucleus) in comparison with the subsolar O/H across the rest of the nebula. Most likely as a consequence, the  \nbubble edge is also significantly less hot ( T 4 ∼ 1.5-1.6) than the rest of the nebula (1.7 < ∼ T 4 < ∼ 1.9). Therefore, the outflow appears to be causing a change in the gas metal content from the nucleus up to ∼ 10 kpc.  \nThis mechanism may also explain the nuclear deficit of O/H. For the Teacup values of log( M ∗ /M ⊙ )=11.15 ± 0.05 and star forming rate SFR ∼ 10 M ⊙ yr -1 (Jarvis et al. 2020; Ramos Almeida et al. 2022), 12+log(O/H) ∼ 8.85 is expected, according to the mass-metallicity-SFR ( M ∗ -ZSFR) relation by Pérez-Montero et al. (2013) or ∼ 8.8 for the extrapolation of the M ∗ -Z-SFR relation by Andrews & Martini (2013) to high M ∗ . Seyfert 2 galaxies with similar M ∗ tend to show O/H consistent with these predictions (Pérez-Díaz et al. 2021). The Teacup nucleus, on the contrary, shows lower than expected O/H (8.66 ± 0.05), while the nuclear N/O (-0.56 ± 0.08, Table .1) is consistent with that expected for its M ∗ ( log (N/O)=-0.54, Pérez-Montero et al. 2013; Andrews & Martini 2013).  \nThe dilution produced by inward flows of low-metallicity gas (for instance, from the outskirts of the two merging galaxies) could explain the nuclear O/H and N/O (Edmunds 1990; Rupke et al. 2010). The outflow is also an interesting possibility. This scenario is supported by the similar O/H of the bubble edge and the nuclear gas, which is moreover enhanced in comparison with the rest of the nebula. On one hand, radiative outflows can couple more efficiently with metals via resonance line scattering (e.g. Pauldrach et al. 1994; Arav et al. 1994; Higginbottom et al. 2024; see also Edmunds 1990). On the other, galaxies show metallicity gradients such that the inner regions have higher abundances (Searle 1971; Kewley et al. 2019; Maiolino & Mannucci 2019). Based on this, a gradient is expected to exist within the ∼ 4 kpc diameter nuclear aperture used in the Teacup analysis. The outflow has been generated in the inner regions ( < ∼ 1 kpc) close to the AGN (Harrison et al. 2015; Ramos Almeida et al. 2017; V23), where the gas is expected to be more metal rich. If it drags gas out to large distances, as proposed in this scenario, the global metallicity of the residual gas within the nuclear aperture would be lower as a consequence.  \nAn implication of this scenario is that the outflow has been capable of ejecting gas from the galaxy center and dragged it up to ∼ 10 kpc. If the outflow expanded without displacing significant amounts of gas to large distances, the bubble abundance would be similar to the rest of the nebula. The implications are important. It supports that metal-enriched galactic outflows (driven by an AGN in this particular case) shape the mass-metallicity relationship, by removing metals from galaxy potential wells and ejecting them to large distances, possibly out into the CGM (Tremonti et al. 2004; Peeples et al. 2014; Chisholm et al. 2018; Tortora et al. 2022).  \nThe behaviour of N/O remains to be explained. An outflow could preserve N/O (Edmunds 1990 (specific modeling would be valuable for AGN generated winds). On the contrary, this ratio is depleted in the bubble edge in comparison with the nucleus and its gradient shows no obvious spatial correlation with the bubble, but just a tentative trend to show among the lowest N/O. This is not necessarily a discrepancy. Given the complexity of the N/O behaviour in terms of secondary and primary stellar production processes, and the fact that the ejected gas would mix with gas across the nebula with a non-uniform N/O distribution, it  \nis difficult to predict how this ratio would behave as the bubble expands and mixes with the pre-existing reservoir.  \nYet another possible scenario to explain the abundance values in the bubble relates to the depletion by dust of metals from the gas-phase ISM. If shocks destroyed dust as the bubble expands (e.g. Dopita et al. 2018), metals could be released to the gas. Since oxygen is more sensitive to depletion than nitrogen (Jones & Ysard 2019), this could explain the higher O/H of the bubble in comparison with the rest of the nebula, while having a tentatively lower N/O.  \nAn alternative scenario to gas ejection from the centre is that local chemical enrichment has been produced by young stars. This is supported by the detection of blue-coloured continuum emission co-spatial with the bubble edge due to a population of stars that are younger ( < ∼ 100-150 Myr) than in the rest of the galaxy ( > ∼ 0.5-1 Gyr, V23). According to these authors, widespread star formation has been triggered at the edge of the bubble due to the compressing action of the jet and outflow (positive feedback). The time scale could be long enough to enrich the local gas with O, but not with secondary N (Mollá et al. 2006; Kumari et al. 2018). This would explain the enhanced O/H in comparison with the rest of the nebula, while N/O is not clearly different, and being at the same time significantly lower than in the nucleus.  \nWhether due to nuclear gas ejection, to local star formation or dust destruction, in all three scenarios the AGN induced outflow is responsible for the metal enrichment of the gas at distances as large as ∼ 10 kpc. The implications are different. In the ejection scenario, the behaviour of O/H provides observational evidence of how AGN induced outflows can deprive the central regions of galaxies from metals and transfer them up to very large distances, possibly out of the galaxy and into the CGM. In the second scenario, the behaviour of O/H provides evidence of how AGN induced outflows can produce local metal enrichment (this is also the case in the third, the dust depletion scenario) at very large distances from the nucleus, with a delay between the quasar onset and the induced metal enrichment of ∼ 100-150 Myr.', '6. Conclusions': 'The giant ( > ∼ 126 kpc) nebula associated with the Teacup QSO2 at z = 0.085 traces part of its CGM. Its properties are strongly influenced by the nuclear activity up to the outer detected emission line regions, where it still provides valuable information about the CGM surrounding the quasar host galaxy. If it was not for the nuclear activity, most (if not all) of this gas would remain invisible. This study is an example of the great potential of studying giant nebulae to investigate in emission the CGM around active galaxies at all redshifts.  \nThis could be a consequence of the ejection of metal rich gas from the nucleus. In such scenario, the Teacup provides observational evidence for how AGN feedback can deprive the central regions of galaxies from gas and displace metals out to very large distances, possibly out of the galaxy. It  \nThe widely studied AGN driven outflow responsible for the well known ionised bubble is enhancing the gas metal content (O/H) up to ∼ 10 kpc from the AGN. The giant nebula shows subsolar metallicity almost everywhere, except the bubble, which has about solar or slightly supersolar metallicity.  \nsupports that metal-enriched AGN outflows can shape the mass-metallicity relationship of galaxies. Alternatively, the O/H enrichment could have been produced locally by the young stellar population formed in the bubble edge, possibly formed as a consequence of positive feedback ∼ 100-150 Myr (V23). A third possibility is the release of oxygen to the gas phase as a consequence of dust destruction in the bubble by shocks triggered by the expanding outflow.  \nIn any of the scenarios considered, the nuclear activity is the ultimate mechanism responsible for the metal enrichment of the gas at large extranuclear distances ( ∼ 10 kpc).  \nAcknowledgements. We are grateful to Luc Binette, Bjorn Emonts and Bruno Rodríguez for useful feedback. We also thank the referee for the careful revision of the manuscript and valuable suggestions. MVM and ACL research has been funded by grant Nr. PID2021-124665NBI00 by the Spanish Ministry of Science and Innovation/State Agency of Research MCIN/AEI/ 10.13039/501100011033 and by "ERDF A way of making Europe". SC acknowledges financial support from the Severo Ochoa grant CEX2021-001131-S and the Ministry of Science, Innovation and Universities (MCIU) under grants PID2019106027GB-C41. EPM acknowledges financial support by project Estallidos8 PID2022-136598NB-C32 (Spanish Ministry of Science and Innovation). The Cosmology calculator by Wright (2006) has been used.', 'References': "Andrews, B. H. & Martini, P. 2013, ApJ, 765, 140  \nArav, N., Li, Z.-Y., & Begelman, M. C. 1994, ApJ, 432, 62 Asplund, M., Amarsi, A. M., & Grevesse, N. 2021, A&A, 653, A141  \nBacon, R., Accardo, M., Adjali, L., et al. 2010, Proc. SPIE, 7735, 773508  \nBalmaverde, B., Capetti, A., Baldi, R. D., et al. 2022, A&A, 662, A23 Baum, S. A., Heckman, T. M., Bridle, A., et al. 1988, ApJS, 68, 643 Borisova, E., Cantalupo, S., Lilly, S. J., et al. 2016, ApJ, 831, 39  \n- Cazzoli, S., Hermosa Muñoz, L., Márquez, I., et al. 2022, A&A, 664, A135  \nChisholm, J., Tremonti, C., & Leitherer, C. 2018, MNRAS, 481, 1690 Dopita, M. A., Vogt, F. P. A., Sutherland, R. S., et al. 2018, ApJS, 237, 10  \n- Dors, O. L., Cardaci, M. V., Hägele, G. F., et al. 2015, MNRAS, 453, 4102  \nEdmunds, M. G. 1990, MNRAS, 246, 678  \n- Emonts, B. H. C., Lehnert, M. D., Villar-Martín, M., et al. 2016, Science, 354, 1128\n- Falkendal, T., Lehnert, M. D., Vernet, J., et al. 2021, A&A, 645, A120 Fossati, M., Fumagalli, M., Lofthouse, E. K., et al. 2021, MNRAS, 503, 3044\n- Gagne, J. P., Crenshaw, D. M., Kraemer, S. B., et al. 2014, ApJ, 792, 72\n- Greene, J. E. & Ho, L. C. 2005, ApJ, 627, 721\n- Harrison C. M., Thomson A.P., Alexander D. M., Bauer F.E., Edge A.C., Hogan M.T., Mullaney J. R., Swinbank A. M., 2015, ApJ, 800, 45\n- Hayes, M., Melinder, J., Östlin, G., et al. 2016, ApJ, 828, 49 Higginbottom, N., Scepi, N., Knigge, C., et al. 2024, MNRAS, 527,\n- 9236\n- Jarvis, M. E., Harrison, C. M., Mainieri, V., et al. 2020, MNRAS, 498, 1560\n- Jones, A. P. & Ysard, N. 2019, A&A, 627, A38.\n- Kacprzak, G. G., Pointon, S. K., Nielsen, N. M., et al. 2019, ApJ, 886, 91\n- Keel, W. C., Chojnowski, S. D., Bennert, V. N., et al. 2012, MNRAS, 420, 878\n- Keel, W. C., Maksym, W. P., Bennert, V. N., et al. 2015, AJ, 149, 155\n- Kewley, L. J., Nicholls, D. C., & Sutherland, R. S. 2019, ARA&A, 57, 511\n- Kumari, N., James, B. L., Irwin, M. J., et al. 2018, MNRAS, 476, 3793\n- Lacerda, E. A. D., Sánchez, S. F., Mejía-Narváez, A., et al. 2022, New A, 97, 101895\n- Lehner, N., Howk, J. C., Tripp, T. M., et al. 2013, ApJ, 770, 138 Luridiana, V., Morisset, C., & Shaw, R. A. 2015, A&A, 573, A42  \nMaiolino, R. & Mannucci, F. 2019, A&A Rev., 27, 3  \nMathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425  \n- McCarthy, P. J., Spinrad, H., van Breugel, W., et al. 1990, ApJ, 365, 487\n- McCarthy, P. J. 1993, ARA&A, 31, 639\n- Mannucci, F., Cresci, G., Maiolino, R., et al. 2010, MNRAS, 408, 2115 Martin, D. C., Matuszewski, M., Morrissey, P., et al. 2015, Nature, 524, 192\n- Martin, D. C., O'Sullivan, D., Matuszewski, M., et al. 2019, Nature Astronomy, 3, 822\n- Mollá, M., Vílchez, J. M., Gavilán, M., et al. 2006, MNRAS, 372, 1069 Moiseev, A. V. & Ikhsanova, A. I. 2023, Universe, 9, 66\n- Osterbrock, D. E. 1989, Astrophysics of Gaseous Nebulae and Active Galactic Nuclei, by Donald E. Osterbrock. Published by University Science Books, ISBN 0-935702-22-9\n- Pauldrach, A. W. A., Kudritzki, R. P., Puls, J., et al. 1994, A&A, 283, 525\n- Peeples, M. S., Werk, J. K., Tumlinson, J., et al. 2014, ApJ, 786, 54 Pérez-Díaz, B., Masegosa, J., Márquez, I., et al. 2021, MNRAS, 505, 4289\n- Pérez-Montero, E., Contini, T., Lamareille, F., et al. 2013, A&A, 549, A25\n- Pérez-Montero, E. 2014, MNRAS, 441, 2663\n- Pérez-Montero, E., Dors, O. L., Vílchez, J. M., et al. 2019, MNRAS, 489, 2652\n- Planck Collaboration, Aghanim, N., Akrami, Y., et al. 2020, A&A, 641, A6.\n- Prescott, M. K. M., Martin, C. L., & Dey, A. 2015, ApJ, 799, 62\n- Ramos Almeida, C., Piqueras López, J., Villar-Martín, M., et al. 2017, MNRAS, 470, 964\n- Ramos Almeida, C., Bischetti, M., García-Burillo, S., et al. 2022, A&A, 658, A155.\n- Rémy-Ruyer, A., Madden, S. C., Galliano, F., et al. 2014, A&A, 563, A31\n- Rupke, D. S. N., Kewley, L. J., & Barnes, J. E. 2010, ApJ, 710, L156. Searle, L. 1971, ApJ, 168, 327\n- Tortora, C., Hunt, L. K., & Ginolfi, M. 2022, A&A, 657, A19\n- Tremonti, C. A., Heckman, T. M., Kauffmann, G., et al. 2004, ApJ, 613, 898\n- Tumlinson, J., Peeples, M. S., & Werk, J. K. 2017, ARA&A, 55, 389 van Ojik, R., Roettgering, H. J. A., Carilli, C. L., et al. 1996, A&A, 313, 25\n- Venturi, G., Treister, E., Finlez, C., et al. 2023, A&A, 678, A127 Villar-Martín, M., Vernet, J., di Serego Alighieri, S., et al. 2003, MNRAS, 346, 273\n- Villar Martín, M., Emonts, B., Humphrey, A., et al. 2014, MNRAS, 440, 3202\n- Villar-Martín, M., Cabrera-Lavers, A., Humphrey, A., et al. 2018, MNRAS, 474, 2302\n- Villar Martín, M., Emonts, B. H. C., Cabrera Lavers, A., et al. 2021, A&A, 650, A84  \nWang, W., Wylezalek, D., Vernet, J., et al. 2023, A&A, 680, A70 Wotta, C. B., Lehner, N., Howk, J. C., et al. 2016, ApJ, 831, 95 Wright, E.L. 2006 PASP, 118, 1711", 'T ab le .1. Photoionisation mo del predictions': '<!-- image -->  \nNotes. Mo del predictions of the O/H and N/O relativ e abundances, th e electron temp eratures, T [OI I I] and T [NI I] and electron densities n e inferred from the [SI I] λ 6716 λ 6731 ratio. The observ ed, extinction corrected line fluxes used in the mo dels are quoted relativ e to H β . E(B-V) v alues w ere inferred from H α /H β . The spatial lo cation of the ap ertures can b e seen in Fig. 3. Ap ertures 3 to 9 co v e r the edge of the ionised bubble. The top (lo w er) part of the table comprises the ap ertures in whic h [NI I] λ 5755 is detected (undetected). N/A in the [SI I] λ 6716 λ 6731 column means that the doublet ratio could not b e measured accurately due to residual sky artefacts. The solar abundances are 12+ l og (O/H)=8.69 ± 0.04 and l og (N/O)=-0.86 ± 0.07 for the Sun (Asplund et al. 2021) .  \nArticle number, page 8 of 8'}