JWST has discovered a large population of active galactic nuclei AGN at high redshift which are weak in the Xrays. Here we present the NIRSpec spectrum of the most extreme of these objects GN28074 an AGN at inlineformulatexmath idTM0001 notationLaTeXz2.26texmathinlineformula with prominent hydrogen and He I broad lines and with the highest limit on the bolometric to Xray luminosity ratio among all spectroscopically confirmed AGN in GOODS. This source is also characterized by a midIR excess likely associated with the AGN torus hot dust. The high bolometric luminosity and moderate redshift of this AGN allow us to explore its properties more in depth relative to other JWSTdiscovered AGN. The NIRSpec spectrum reveals prominent slightly blueshifted absorption of Hinlineformulatexmath idTM0002 notationLaTeXalphatexmathinlineformula Hinlineformulatexmath idTM0003 notationLaTeXbetatexmathinlineformula and He I inlineformulatexmath idTM0004 notationLaTeXlambdatexmathinlineformula10830. The Balmer absorption lines require gas with densities of inlineformulatexmath idTM0005 notationLaTeXnrm Hgt 108rm cm3texmathinlineformula consistent with clouds in the broadline region BLR. This finding suggests that part of the Xray weakness is due to high Comptonthick Xray absorption by clouds in the BLR or in its outer regions. GN28074 is also extremely radioweak. The radio weakness can also be explained in terms of absorption as the inferred density of the BLR clouds makes them attenuate the radio emission through freefree absorption. Alternatively the nuclear magnetic field may be underdeveloped resulting both in intrinsically weak radio emission and lack of hot corona hence intrinsic Xray weakness. Finally we show that recently proposed scenarios invoking hyperdense outflows or Raman scattering to explain the broad Hinlineformulatexmath idTM0006 notationLaTeXalphatexmathinlineformula are ruled out.
['Astrophysics - Astrophysics of Galaxies']
JADES  the Rosetta stone of JWSTdiscovered AGN deciphering the intriguing nature of early AGN
{'JADES - The Rosetta Stone of JWST-discovered AGN: deciphering the intriguing nature of early AGN': "Ignas Juodžbalis, 1 , 2 ★ Xihan Ji, 1 , 2 Roberto Maiolino, 1 , 2 , 3 Francesco D'Eugenio, 1 , 2 Jan Scholtz, 1 , 2 Guido Risaliti, 4 , 5 Andrew C. Fabian, 6 Giovanni Mazzolari, 1 , 7 , 8 Roberto Gilli, 8 Isabella Prandoni, 9 Santiago Arribas, 10 Andrew J. Bunker, 11 Stefano Carniani, 12 Stéphane Charlot, 13 Emma Curtis-Lake, 14 Anna de Graaff, 15 Kevin Hainline, 16 Eleonora Parlanti, 12 Michele Perna, 10 Pablo G. Pérez-González, 10 Brant Robertson, 17 Sandro Tacchella, 1 , 2 Hannah Übler, 1 , 2 Christina C. Williams, 18 Chris Willott, 19 Joris Witstok 1 , 2  \n- 1 Kavli Institute for Cosmology, University of Cambridge, Madingley Road, Cambridge, CB3 OHA, UK.\n- 2 Cavendish Laboratory - Astrophysics Group, University of Cambridge, 19 JJ Thomson Avenue, Cambridge, CB3 OHE, UK.\n- 3 Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, UK\n- 4 Dipartimento di Fisica e Astronomia, Università di Firenze, via G. Sansone 1, 50019 Sesto Fiorentino, Firenze, Italy\n- 5 INAF - Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5, I-50125 Firenze, Italy\n- 6 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK\n- 7 Dipartimento di Fisica e Astronomia, Università di Bologna, Via Gobetti 93/2, I-40129 Bologna, Italy\n- 8 INAF - Osservatorio di Astrofisica e Scienza dello Spazio di Bologna, Via Gobetti 93/3, I-40129 Bologna, Italy\n- 9 INAF - Istituto di Radioastronomia, Via Gobetti 101, I-40129 Bologna, Italy\n- 10 Centro de Astrobiología (CAB), CSIC-INTA, Cra. de Ajalvir Km. 4, 28850- Torrejón de Ardoz, Madrid, Spain\n- 11 Department of Physics, University of Oxford, Denys Wilkinson Building, Keble Road, Oxford OX1 3RH, UK\n- 12 Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy\n- 13 Sorbonne Université, CNRS, UMR 7095, Institut d'Astrophysique de Paris, 98 bis bd Arago, 75014 Paris, France\n- 14 Centre for Astrophysics Research, Department of Physics, Astronomy and Mathematics, University of Hertfordshire, Hatfield AL10 9AB, UK\n- 15 Max-Planck-Institut für Astronomie, Königstuhl 17, D-69117, Heidelberg, Germany\n- 16 Steward Observatory, University of Arizona, 933 N. Cherry Avenue, Tucson, AZ 85721, USA\n- 17 Department of Astronomy and Astrophysics University of California, Santa Cruz, 1156 High Street, Santa Cruz CA 96054, USA\n- 18 NSF's National Optical-Infrared Astronomy Research Laboratory, 950 North Cherry Avenue, Tucson, AZ 85719, USA\n- 19 NRC Herzberg, 5071 West Saanich Rd, Victoria, BC V9E 2E7, Canada  \nAccepted XXX. Received YYY; in original form ZZZ", 'ABSTRACT': "JWST has discovered a large population of Active Galactic Nuclei (AGN) at high redshift, which are weak in the X-rays. Here we present the NIRSpec spectrum of the most extreme of these objects, GN-28074, an AGN at 𝑧 = 2 . 26 with prominent hydrogen and He /i.pc broad lines, and with the highest limit on the bolometric to X-ray luminosity ratio among all spectroscopically confirmed AGN in GOODS. This source is also characterized by a mid-IR excess, likely associated with the AGN torus' hot dust. The high bolometric luminosity and moderate redshift of this AGN allow us to explore its properties more in depth relative to other JWST-discovered AGN. The NIRSpec spectrum reveals prominent, slightly blueshifted absorption of H 𝛼 , H 𝛽 and He /i.pc 𝜆 10830. The Balmer absorption lines require gas with densities of 𝑛 H > 10 8 cm -3 , consistent with clouds in the Broad Line Region (BLR). This finding suggests that part of the X-ray weakness is due to high (Compton thick) X-ray absorption by clouds in the BLR, or in its outer regions. GN-28074 is also extremely radio-weak. The radio weakness can also be explained in terms of absorption, as the inferred density of the BLR clouds makes them attenuate the radio emission through free-free absorption. Alternatively, the nuclear magnetic field may be underdeveloped, resulting both in intrinsically weak radio emission and lack of hot corona, hence intrinsic X-ray weakness. Finally, we show that recently proposed scenarios, invoking hyper-dense outflows or Raman scattering to explain the broad H 𝛼 , are ruled out.  \nKey words: galaxies: active - quasars: supermassive black holes - galaxies: Seyfert", '1 INTRODUCTION': 'The first two years of operation of the James Webb Space Telescope (JWST) have revealed a large population of broad line active galactic  \nnuclei (AGN) at high redshift (Harikane et al. 2023; Kocevski et al. 2023; Matthee et al. 2024; Maiolino et al. 2024b, 2023a; Greene et al. 2024; Kokorev et al. 2023; Furtak et al. 2024; Wang et al. 2024b,c; Übler et al. 2024). These JWST-selected AGN span the redshift range from ∼ 3 to ∼ 11 and have bolometric luminosities ranging from 10 42 to 10 45 erg s -1 . Therefore, they cover the intermediate/lowluminosity regime of AGN activity inaccessible by pre-JWST highredshift quasar surveys.  \nThese objects display unambiguous evidence for the presence of an AGN, not only in terms of luminous, broad (FWHM ∼ 3000 km s -1 ) H 𝛼 and H 𝛽 lines (which lack a kinematically similar counterpart in the [O /i.pc/i.pc/i.pc] 𝜆𝜆 4959 , 5007 doublet, hence cannot be ascribed to outflow), but also some of them have other transitions typical of AGN, such as high ionization and even coronal lines (Übler et al. 2023; Ji et al. 2024a,b; Maiolino et al. 2024b).  \nHowever, these AGN discovered by JWST spectroscopy also exhibit unusual properties when compared to standard Type 1 AGN. A particularly notable property of JWST-selected low-luminosity Type 1 AGN is their X-ray weakness, which has been pointed out by various authors (Maiolino et al. 2024a; Yue et al. 2024a; Ananna et al. 2024). These works reported no X-ray detections for the vast majority of Type 1 AGN selected by JWST, even when stacked together. The same X-ray weakness has been found also for narrow line AGN candidates discovered by JWST (Scholtz et al. 2023; Chisholm et al. 2024) and for mid-IR selected AGN (Lyu et al. 2024). In addition, a fraction of JWST-selected Type 1 AGN show the presence of absorption in the Balmer lines (Matthee et al. 2024; Kocevski et al. 2024; Wang et al. 2024b). We should mention that examples of X-ray weak AGN with mild optical obscuration, and broad-lined AGN, were already found at high-z even before the advent of JWST data and using ground based surveys (e.g., Fujimoto et al. 2022; Ma et al. 2024).  \nIn the low-redshift Universe, X-ray weakness of type 1 AGN is not common, but it has been observed in a number of metal-poor and low-mass Type 1 AGN candidates that the X-ray emission is weak or lacking (e.g., Simmonds et al. 2016; Cann et al. 2020). Amongthelow-redshift and metal-poor AGN candidates, Burke et al. (2021) reported observations of an absorption feature in the broad H 𝛼 line in SDSS J102530.29+140207.3 at 𝑧 = 0 . 1, suggesting large column densities of neutral hydrogen and thus obscuration in the nuclear region. It is noteworthy that outside the metal-poor AGN sample, there have also been detections of strong broad Balmer-line absorptions in general Type 1 AGN up to a redshift of 𝑧 ∼ 2 in the pre-JWST era, but at a very low detection rate (e.g., Hutchings et al. 2002; Aoki et al. 2006; Aoki 2010; Hall 2007; Ji et al. 2012, 2013; Wang & Xu 2015; Zhang et al. 2015; Shi et al. 2016; Williams et al. 2017; Schulze et al. 2018; Hamann et al. 2019). With typical values of 10 5 -10 7 M ⊙ , the black hole masses of the local metal-poor AGN candidates are comparable to the Type 1 AGN at early times currently revealed by JWST.  \nSeveral explanations of this X-ray weakness in JWST-confirmed AGN have been proposed thus far. One possibility is that (hard) Xrays are absorbed by a Compton thick medium ( 𝑁 H > 10 24 cm -2 ). However, given that X-ray weakness is seen also for Type 1, broad line AGN, the absorbing medium should have comparatively little dust content as dust sufficient to absorb X-rays would also absorb the BLR emission, assuming typical dust-to gas ratios. One possibility that has been proposed is that the BLR clouds could be a potential absorbing medium, as they are expected to be dust free and with large column of gas. However, given that the X-ray weakness affects most JWSTselected AGN, the covering factor of the BLR clouds in these objects should be much higher than in classical AGN (which is typically of 𝐶 𝑓 ≡ Ω / 4 π ∼ 0 . 5 at low Eddington ratios; Ferland et al. 2020);  \nthis scenario seems to be confirmed by the larger equivalent width of the broad H 𝛼 in JWST-selected AGN relative to classical AGN (Maiolino et al. 2024a; Wang et al. 2024d). Alternatively, the X-ray weakness might be intrinsic to this population of AGN, possibly due to a lack of a corona. JWST-selected AGN tend to share properties similar to those of the so-called Narrow Line Sy1, characterized by a steep X-ray spectrum (Vasudevan & Fabian 2007; Tortosa et al. 2022, 2023) and some of them are characterised by high accretion rates, which typically results in a steep X-ray spectrum (Hagen et al. 2024). The steepened X-ray spectra would thus make the detection of hard X-ray emission in the highly accreting high𝑧 AGN very difficult.  \nHowever, further verification of these scenarios is difficult, because without X-ray detections there is no simple method of constraining the Compton thickness or the shape of the X-ray continuum (although the very few X-ray detected and X-ray weak JWST-selected AGN appear heavily absorbed; Maiolino et al. 2024a). However, one possible way of exploring heavy X-ray absorption and, possibly, Compton thickness would be via the detection of certain absorption lines that could be used to provide estimates on the total hydrogen column density obscuring the source.  \nOn the other hand, the radio properties of the X-ray weak AGN discovered by JWST have yet to be explored in detail. However, radio emission can provide crucial insights into the origin of the observed X-ray weakness, given the well-established connection between AGN X-ray and radio luminosities reported in several studies (e.g. Merloni et al. 2003; Panessa et al. 2019; Bariuan et al. 2022; Wang et al. 2024a, for a review). In the local universe, analyses of radio emission from X-ray weak sources have sometimes revealed corresponding radio weakness (Johnson et al. 2009; Arcodia et al. 2024), while in other cases, the radio emission is found at the expected levels (Berton et al. 2018; Paul et al. 2024), with detections of compact radio cores or small-scale jets (Meyer et al. 2024a; Singha et al. 2023; Yao et al. 2021). A more detailed discussion of the radio properties of JWST-detected BLAGN will be presented in a forthcoming paper (Mazzolari et al., in prep).  \nIn this paper we present the analysis of GN-28074 (RA = 189.06458, Dec = 62.27382; CANDELS ID GN-19549; 3D-HST ID 26945) a broad line AGN at 𝑧 = 2 . 26 detected by JWST in the GOODS-N field and exhibiting clear absorption features in the H 𝛼 , H 𝛽 and He /i.pc 𝜆 10830 lines. As we discuss in this paper, the presence of these lines reveals the existence of gas with moderately high temperature ( ∼ 10,000 K) and high density and high column densities, likely associated with the BLR or its outer parts. In addition, the continuum of the source exhibits similar colors to the recently identified red sub-population known as the Little Red Dots (LRDs); (Labbe et al. 2023; Furtak et al. 2023; Kokorev et al. 2024; Matthee et al. 2024), with it passing most the color cuts presented in Kokorev et al. (2024) if redshifted to 𝑧 = 7 (F277W - F356W = 0.96, F277W F444W = 1.64, F150W - F200W = 0.88). An additional component of dust emission at long wavelengths further confirms the presence of a reddened AGN or a dusty torus. However, despite the unambiguous presence of an AGN, this source is undetected both in the X-rays (Maiolino et al. 2024a) and in the radio (Mazzolari in prep.). We will show that, therefore, this AGN is a unique case study to understand the properties of the new class of AGN that is being discovered by JWST.  \nThe text is organized as follows. In section 2 we present a summary of the spectroscopic data reduction. section 3 contains a description of our line fitting procedure. The bulk of the analysis, pertaining to the detected absorption lines, properties of the central black hole and the continuum shape is presented in section 5. In the following section 6 we briefly explore alternative scenarios proposed to  \nexplain broad emission in JWST sources. Discussion on potential local analogues and BLR absorption in other spectroscopic JWST AGN is presented in section 8 and section 7 respectively. Finally, in section 9, we summarize our findings and put them in context of the previous observations. Throughout the paper, we assume flat Λ CDM cosmology with Ω 𝑚 = 0 . 315 and 𝐻 0 = 67 . 4 km s -1 Mpc -1 .', '2 DATA REDUCTION': "Data utilized in this paper has been obtained as part of the JWST Advanced Extragalactic Survey (JADES) program (Eisenstein et al. 2023). The spectroscopic observations are split into tiers depending on whether they were selected with HST or JWST instruments and depth of exposure. Spectroscopy of our source was obtained as part of the Medium/HST tier in the GOODS-N field and consisted of 1.7 hours of exposure per source in the prism and 1.7 hours per source in each of the medium-resolution gratings, G140M/F070LP, G235M/F170LP and G395M/F290LP. The observations used two dithered positions (grey rectangles in the inset figure of Figure 1), with three nods each.  \nData reduction procedures were largely similar to other JADES papers (Bunker et al. 2023; Carniani et al. 2023; D'Eugenio et al. 2024). In summary, the pipeline used was developed by the NIRSpec GTO team and the ESA NIRSpec Science Operations Team. As the absorption features appear to be unresolved, we employ 1D spectra extracted from the central 3 pixels (at 0.1 '' per pixel scale) of the 2D spectra.  \nAn additional post-processing procedure was carried out after the main data reduction stage as the sigma clipping employed by the GTOpipelinehaserroneously flagged strong emission lines and rapid changes in emission due to deep absorption features. This procedure consisted of utilizing the unclipped spectrum to visually identify incorrectly clipped regions and substitute them into the clipped data. The final spectrum obtained by combining the three grating spectra is shown in Figure 1.  \nWe also use the multi-band photometry from the NIRCam component of the JADES survey (Rieke et al. 2023) and from FRESCO (First Reionization Epoch Spectroscopically Complete Observations; Oesch et al. 2023), combined with HST and Spitzer photometry from the Rainbow database 1 .", '3 EMISSION AND ABSORPTION LINES FITTING': "The spectrum in Figure 1 shows clear broad components present in the bright permitted lines (hydrogen and helium lines), which lack kinematically comparable counterparts in the [O /i.pc/i.pc/i.pc] forbidden doublet. This unambiguously identifies their origin being the broad line region (BLR) of an active galactic nucleus (AGN). It should be noted that the profiles of [O /i.pc/i.pc/i.pc] 𝜆𝜆 4959,5007 and [S /i.pc/i.pc/i.pc] 𝜆 9531 lines do exhibit faint, broad wings, but these are likely associated with an ionized outflow as their kinematics differ strongly (much narrower) from the BLR and will be explored further in Appendix C.  \nIn addition, the spectrum of GN-28074 clearly shows deep blueshifted (by a few 100 km s -1 ) absorption features in the broad H 𝛼 , H 𝛽 and He /i.pc 𝜆 10830 lines in Figure 2. While the He /i.pc 𝜆 10830 is associated with a metastable level, hence often seen in absorption (provided that the conditions and column density are adequate), observing H 𝛼 and H 𝛽 in absorption is much more difficult, as it  \nrequires significant population of the 𝑛 = 2 level, which has a very short lifetime ( ∼ 10 -8 s) for radiative transition. As we will show, this implies large densities and large column densities of the medium along the line of sight (LOS).  \nThe absorption features were modeled using a standard attenuation model:  \n𝑓 𝜆 = 1 -𝐶 𝑓 + 𝐶 𝑓 𝑒 -𝜏 𝜆 , (1)  \nwhere 𝐶 𝑓 is the covering factor of the absorber and 𝜏 𝜆 - its optical depth profile. The optical depth profiles were assumed to be Gaussian:  \n𝜏 𝜆 = 𝜏 0 exp      -0 . 5 𝜆 -𝜆 0 𝑒 Δ 𝑣 / 𝑐 𝜎 ! 2      , (2)  \n  with 𝜏 0 corresponding to the optical depth at the core of the line, Δ 𝑣 - the velocity shift of the absorber and 𝜎 - its velocity dispersion. The attenuation was applied to the BLR and continuum emission. Indeed, the absorption is so deep that it cannot be absorbing only the continuum (as subtracting the broad line emission would result in negative absorption), it must therefore be absorbing, at least partially, also the broad lines. This implies that the absorbing medium must be located within or just outside the BLR.  \nAll emission lines in R1000 were modelled using Gaussian decomposition with separate Gaussians representing the narrow line, BLR and ionized outflow emission. Fitting was performed using a Bayesian method with linear uniform priors constraining the line widths to between 100 and 800 km s -1 for the narrow, 800 10,000 km s -1 for the broad and 600 - 2000 km s -1 for the ionized outflow components. The offset between the BLR and outflows was allowed to vary by up to ± 10,000 km s -1 with all BLR and ionized outflow components constrained to the same kinematics with the exception of He /i.pc 𝜆 10830 - due to its resonant nature its FWHM was treated separately. The widths of the broad Pa 𝛾 and Pa 𝛽 lines were also left independent, because their intrinsic widths could be different from those of the Balmer series (e.g., Landt et al. 2008). In addition, the Pa 𝛽 line is affected by wavelength calibration issues causing a ∼ 1 px offset from the other lines, which forces us to leave all its kinematics, aside from the width of the narrow line, independent of the remaining lines. All linewidths were convolved with the line spread function (LSF) for point sources from de Graaff et al. (2024) during the fitting procedure. Absorption in the Balmer lines was constrained to have the same kinematics, while the He /i.pc 𝜆 10830 absorption was fitted with independent Δ 𝑣 and 𝜎 . The prior for Δ 𝑣 was uniform, constrained to be between -1000 and 1000 km s -1 , while 𝜎 was allowed to vary between 0 and 400 km s -1 . The covering factor was constrained to be the same for all absorption lines. At this stage (i.e. for the purpose of the line fitting) the continuum was approximated empirically as a fourth-order polynomial.  \nThe posterior probability distribution of all model parameters was estimated using a Markov-Chain Monte-Carlo integrator (ForemanMackey et al. 2013); we initialised the chains in a small neighbourhood around the maximum-likelihood parameter values, estimated via least-squares minimisation. The moments of the posterior distribution for the neutral outflow model are presented in Table 1. We note that Figure 2 displays significant residuals around the peaks of bright lines. The likely cause of this is the line intensity saturating the detector in some exposures resulting in the pipeline failing to derive valid errors for their peaks. A summary of all fitted emission line fluxes is provided in Table 2 while the relevant line widths are summarized in Table 3. The fit uncertainties for the narrow lines include a ∼ 20% systematic uncertainty on the LSF, which is less significant for the well-resolved broad lines.  \nFigure 1. R1000 spectrum of GN-28074 with emission lines marked. The fluxes on the y axis are presented in logarithmic scale to better showcase weaker lines and the shape of the continuum. All bright permitted lines appear to have broad components while the continuum showcases the characteristic 'v' shape found in some objects of the LRD population. The false-colour NIRCam image (inset) highlights the point-source nature of the object, and the position of the MSA slitlets (grey rectangles). We note a pink region to the south-east, denoting a deficiency in F182M flux (see text for a discussion). The rest of the NIRCam single-band cutouts highlight the extended emission (only discernible at SW wavelengths); the angular size of these cutouts is the same as the RGB image.  \n<!-- image -->  \nTable 1. Asummary of absorption profile properties given by the fit. The first column gives the name of the fitted line, columns two and three - covering fraction and optical depth at the center of the absorption, columns four to five give the best-fit neutral outflow kinematics while the final column contains estimated column densities.", '4 CONTINUUM FITTING AND MID-IR EXCESS': "In the NIRCam images GN-28074 appears as a point source with a spatially extended component to the south-east (inset single-band cutouts in Figure 1). This extension is clearly seen in all available filters from the SW NIRCam module (F090W-F210M), but we do not rule out its presence at redder wavelengths (an in-depth morphological analysis is beyond the scope of this article). Shifting focus on the inset RGB image, we infer that this south-east extension is very likely to be at the same redshift as GN-28074. The extension  \nhas a strong line excess in F210M relative to F182M (and F210W). Notice that this excess feature (pink in the inset RGB) is stronger than what is seen at other position angles around the central source, which indicates it is an SED feature of the south-east extension, not of the point source itself. The most parsimonious explanation is that the line excess in F210M is tracing H 𝛼 in the south-east region, which, therefore, is at the same redshift as the point source. The pink colour in the RGB image also indicates that, compared to the point source, this extended component is stronger in F210M than in the rest-UV  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n(c) He /i.pc  \n+ Pa  \n𝛾  \n(b) H  \n𝛼  \n+ [N /i.pc/i.pc]  \nFigure 2. R1000 spectrum expanded around some of the permitted lines of GN-28074 and showing the associated spectral fits. The broad line region (BLR) and the (extended) outflow components are plotted as solid blue and red lines respectively. The narrow components are shown with dashed green lines. The combined fit, including the absorption absorption components, is shown in magenta. Top left: H 𝛽 and [O /i.pc/i.pc/i.pc] 𝜆𝜆 5007,4959 complex, showcasing the deep absorption in H 𝛽 and broadening of [O /i.pc/i.pc/i.pc] 𝜆𝜆 5007,4959 lines due to the presence an outflow. Top right: H 𝛼 , [N /i.pc/i.pc] 𝜆𝜆 6548,6583 and [S /i.pc/i.pc] 𝜆𝜆 6716,6731 . Notable is the relative weakness of the [N /i.pc/i.pc] 𝜆𝜆 6549,6585 and [S /i.pc/i.pc] 𝜆𝜆 6716,6731 doublets. Bottom left: Combined He/i.pc 𝜆 10830 and Pa 𝛾 line fit. The ionized outflow component is not required here by the data. Bottom right: Fit to the Pa 𝛽 line. Due to wavelength calibration issues, this line is slightly offset to the rest and thus requires separate kinematics to properly fit. The ionized outflow is absent just as in the Pa 𝛾 + He /i.pc 𝜆 10830 fit.  \n<!-- image -->  \nand rest-optical continuum. This however could result from at least two different scenarios: an SED dominated by star formation, or an SED dominated by nebular emission, such as the ionisation cone of an AGN (e.g., Tacchella et al. 2024). Interestingly, F150W seems to be more extended than F182M. Given that the NIRCam PSF is broader in F182M than F150W, this line of evidence may imply the presence of a diffuse region of [O /i.pc/i.pc/i.pc] emission.  \nIn a recent article, Wang et al. (2024b) report the discovery of a broad-line AGN with weak or even no dust emission, confirmed by MIRI. In addition, Matthee et al. (2024) find a tentative detection of broad H 𝛼 (4 𝜎 ) that was confirmed as MIR-weak by Williams et al. (2024), contrary to expectations that similar sources should exhibit strongly rising MIR continuum from AGN dust heating (e.g. Labbe et al. 2023). While a weak MIRI SED appears prevalent among photometrically identified samples (e.g. Williams et al. 2024; Akins et al. 2024; Pérez-González et al. 2024), unfortunately, we still lack  \nextensive MIRI coverage of spectroscopically confirmed type 1 AGN, therefore we cannot yet verify if the lack of MIR emission is a peculiarity of one-two sources, or if it is a more common feature of low-luminosity AGN, perhaps due to a lack of a dusty torus, or due to colder dust temperatures than in previously studied AGN. To investigate if GN-28074 also lacks MIR emission, we utilize Spitzer IRAC imaging across all four channels, covering 3.6 to 8 𝜇 m, together with MIPS 24 𝜇 m band. Note that MIPS may suffer from contamination due to a nearby source, and that we do not include Herschel non detections, because their high upper limits are not constraining. The rest-frame UV and optical range was covered using JWST and HST photometry, as discussed in section 2. Our object has S/N > 10 across all bands considered allowing for robust SED fitting of the rest-frame UV to near infrared continuum. This also marks GN-28074 as one of the few JWST AGN to have robust MIR constraints.  \nTable 2. Table summarizing all fitted emission line fluxes as measured without correcting for dust attenuation. The first column indicates the name of a given line, the second - its type (Narrow, Broad or Outflow). The final column gives the measured flux in units of 10 -18 erg s -1 cm -2 .  \nTable 3. Summary of widths for all fitted emission lines. The first three rows list the common widths fitted to the broad, narrow and outflow lines. The remaining rows give FWHM for the lines for which it was left independent.  \nWe carry out the fits using the CIGALE code (Yang et al. 2022). The AGN component of the templates was modeled using SKIRTOR continuum emission models (Stalevski et al. 2016) with AGN fraction allowed values of 0 to 0.99, the polar dust extinction of the continuum was modeled with the SMC extinction law, standard for dust reddened AGN (Reichard et al. 2003; Richards et al. 2003), and allowed E(B - V) values ranging from 0.03 to 4. Stellar populations were modeled using BC03 templates with Salpeter (1955) IMF.  \nStellar dust attenuation was modeled using the Calzetti attenuation law with E(B-V) allowed to vary between 0 and 4. Dust emission was modeled with Dale et al. (2014) models with emission slope 𝛼 = 2. Gas and stellar metallicities were set to either solar or 10% solar, however, the exact metallicity has little impact on the shape of the model SEDs. Parameters not mentioned were kept to CIGALE defaults.  \nThe resulting fit is shown in Figure 3 and shows that the bestfitting model prefers the presence of a strongly attenuated accretion disk component with E ( B -V ) = 1 . 04 ± 0 . 24, corresponding to an 𝐴 𝑉 = 3 . 02 ± 0 . 70, with accompanying infrared emission from hot dust forming the continuum emission displayed there in magenta. The stellar continuum component is attenuated considerably less, 𝐴 𝑉 = 0 . 68 ± 0 . 01.  \nTo test if the mid infrared photometry can be explained with stellar emission or associated dust reprocessing alone we refit the data by forcing the AGN fraction to 0. The resulting best fit without AGN is shown in Figure 4. As can be seen, the stellar-only models yield a considerably worse overall fit and can not explain the infrared photometry, necessitating the presence of a hot dust component, which can be readily identified as a dusty torus or hot polar dust. No matter the exact interpretation, our result indicates the presence of significant amounts of hot dust surrounding the core of GN-28074.  \nA notable feature, common to both stellar and AGN+stellar fits, is that both of them attempt to fit the discontinuity in flux between 1 and 2 𝜇 m with a Balmer break. However, both leave a significant residual in that location (Figure 3, Figure 4), thus a Balmer break fit does not explain the emission pattern well. If the feature is indeed a Balmer break, indicative of a more evolved stellar population, the stellar mass estimate gives log 𝑀 ∗ 𝑀 ⊙ = 10 . 78 ± 0 . 05. For comparison, we use the fits to the narrow lines to estimate an upper limit on the dynamical mass of the object using the virial estimator (calibrated on stellar kinematics) from Cappellari et al. (2006):  \n𝑀 dyn < 5 𝜎 2 𝑅 𝑒 𝐺 , (3)  \nwhere 𝜎 is the gas velocity dispersion and 𝑅 𝑒 the circularised half-light radius. The fits to the narrow lines give a FWHM = 150 + 36 -36 km s -1 , corresponding to 𝜎 ≈ 64 km s -1 . The value for 𝑅 𝑒 is trickier to estimate as the object is point-source dominated, thus we adopt 0.06 arcsec as an upper limit, using the value measured by van der Wel et al. (2014) from HST F160W imaging. This value is the semi-major axis, which translates into a circularised 𝑅 e = 0 . 035 arcsec after taking into account the measured shape (axis ratio 𝑞 = 0 . 65; van der Wel et al. 2014). The resulting mass is log 𝑀 dyn 𝑀 ⊙ ≲ 9 . 1, lower than the 'maximal' stellar mass inferred interpreting the spectral break as a Balmer break. The strong mismatch with 𝑀 ∗ does not depend on the specific virial calibrator we use; alternative formulae give log 𝑀 dyn 𝑀 ⊙ = 8 . 9-9 . 5 (Wisnioski et al. 2018; van der Wel et al. 2022). This inconsistency could be due to geometric effects, if the narrow-line gas kinematics are dominated by a disc seen close to face on. However, Wang et al. (2024c) reports three more examples of broad-line AGN with spectral breaks around rest-frame 3,600-4,000 Å, and where the virial mass is lower than the maximal stellar mass. Therefore, it is unlikely that the observed drop in the spectrum can be explained entirely by a Balmer break. We further note that in the R1000 spectrum (Figure 1) the break appears less steep and at redder wavelengths than the typical Balmer break. However, interpreting the break as a stellar 4000-Å break would further increase the stellar mass and, consequently, the tension with the dynamical-mass constraint.  \nWe do not discuss further the nature of the drop between 3,000  \nFigure 3. The best-fitting CIGALE model. It can be seen that the mid infrared emission is well explained by hot dust emission from an attenuated accretion disk emission (magenta) with comparatively little contribution from stellar dust (in red). The excess flux around 2.0 𝜇 mthat is not well reproduced can be attributed to the strong H 𝛼 emission.  \n<!-- image -->  \nand 4,000 Å rest frame, as we defer further analysis to another paper in which other similar sources are also analysed.", '5.1 Black hole properties': 'Weutilize our fits to the broad line region (BLR) emission to constrain the black hole mass using single-epoch virial relations. Recent results by Abuter et al. (2024) measured the BLR size via interferometric data, and have put into question the usage of BH virial relations using transitions such as CIV (Netzer et al. 2007; Negrete et al. 2012); however, their finding that the discrepancy is reduced to only a factor 2.5 when using the H 𝛼 line is reassuring, as this is within the calibration uncertainties. Furthermore, the deviation is associated with super-Eddington accretion influencing the size of the BLR, while the BH of 28074 is likely in a sub-Eddington regime.  \nOur fits to the BLR yield FWHM H 𝛼 = 3610 + 21 -22 km s -1 along with 𝐿 H 𝛼 = 1 . 21 + 0 . 01 -0 . 01 × 10 43 erg s -1 . Converting this to black hole  \nmass via the virial relation (Reines & Volonteri 2015):  \nlog 𝑀 BH 𝑀 ⊙ = 6 . 60 + 0 . 47 log GLYPH<18> 𝐿 𝐻𝛼 10 42 erg s -1 GLYPH<19> + 2 . 06 log GLYPH<18> FWHM 𝐻𝛼 1000kms -1 GLYPH<19> (4)  \nyields log 𝑀 BH 𝑀 ⊙ = 8 . 26 ± 0 . 30, the uncertainty on this value is completely dominated by the intrinsic scatter on Equation 4. The scaling relation between the luminosity of the broad H 𝛼 component and bolometric luminosity calibrated by Stern & Laor (2012) implies 𝐿 bol = 1 . 58 + 0 . 01 -0 . 01 × 10 45 erg s -1 , with a 0.3 dex intrinsic scatter, and an Eddington ratio of ∼ 7%.  \nEstimating the extinction correction to the BLR is not trivial even when a broad component in both H 𝛽 and H 𝛼 lines is seen, as the intrinsic Balmer decrement can be significantly steeper than the case Bscenario and can potentially reach values up to ten (Ilić et al. 2012). In the case of our object, we find H 𝛼 to H 𝛽 ratios to be 9 . 25 + 0 . 33 -0 . 31 and 4 . 14 + 0 . 16 -0 . 14 for the broad and narrow lines respectively. Assuming a standard case B decrement of 2.85 for the narrow lines along with the SMC extinction curve (Gordon et al. 2003), which might be  \nFigure 4. Same as Figure 3, except for models without AGN emission. It can readily be seen that stellar-only templates leave significant residuals from 5 𝜇 m onwards resulting in a poorer reduced 𝜒 2 value.  \n<!-- image -->  \nrepresentative in galaxies at 𝑧 ∼ 2 (e.g., Reddy et al. 2015), we obtain 𝐴 𝑉 = 1 . 13 + 0 . 11 -0 . 10 . Adopting this value for the dust correction increases the measured 𝐿 𝑏𝑜𝑙 by ∼ 0.5 dex and results in a corrected BH mass of log 𝑀 BH 𝑀 ⊙ = 8 . 47 ± 0 . 30 and an Eddington ratio of 0 . 12. If we instead take the broad lines Balmer ratio and conservatively assume an intrinsic case B ratio, we would get an upper limit on the extinction of 𝐴 𝑉 < 2 . 3, which would imply a maximum correction on the bolometric luminosity of 0.9 dex, a maximum black hole mass of log 𝑀 BH 𝑀 ⊙ < 8 . 67, and a maximum Eddington ratio of 0.2. One should also take into account the additional uncertainty on both black hole mass and Eddington ratio of 0.3 - 0.5 dex of the virial relations. Observed properties of the BLR are summarized in Table 4.', '5.2 Nuclear neutral absorption': '5.2.1 Modelling and evidence for very high gas densities  \nFrom Figure 2, it is clear that strong blueshifted absorption features are present on top of the broad components of H 𝛽 , H 𝛼 , and He/i.pc 𝜆 10830. Before JWST, strong blueshifted absorptions in hydrogen and/or helium lines have previously been observed in the spectra  \nTable 4. A summary of the observed properties of the BLR and the inferred mass of the supermassive black hole and its accretion rate. The first column denotes the name of the quantity, the second - its value. Listed luminosities are corrected for dust. The final column provides upper limits based on the conservative A 𝑉 estimation based on the observed broad Balmer line decrement.  \nof a few AGN (e.g., Hutchings et al. 2002; Aoki et al. 2006; Aoki 2010; Hall 2007; Ji et al. 2012, 2013; Wang & Xu 2015; Zhang et al. 2015; Shi et al. 2016; Williams et al. 2017; Schulze et al. 2018; Hamann et al. 2019), including the prototypical local Seyfert  \n1 NGC4151 (Hutchings et al. 2002). However, these were extremely rare cases, less than 0.1% of the quasar/AGN population. Yet, the JWST spectra are revealing an increasing number of these Balmer line absorption features in Type 1 AGN, possibly making about 10% of the AGN population discovered by JWST (Matthee et al. 2024; Wang et al. 2024b,c; Kocevski et al. 2024). Therefore, such absorption features are probably the signature of properties characterizing the new population of AGN discovered by JWST. The prominence of these features in GN-28074, along with the wealth of additional information, makes this object the perfect case study to explore these properties.  \nAs already discussed, absorption of H 𝛼 and H 𝛽 , and, to some extent, He /i.pc , requires peculiar conditions of the gas along the LOS. Specifically, these need substantial amount of (mostly) neutral gas along the LOS that has significant amounts of hydrogen and helium pumped to excited states, which are the 𝑛 = 2 level for hydrogen and the metastable triplet state 2 3 𝑆 for helium. Due to the short lifetime of H( 𝑛 = 2) compared to that of He(2 3 𝑆 ), it is much more difficult to have significant absorption in hydrogen lines without the presence of high-density gas, as we demonstrate in the following.  \nWe estimate the absorbing column densities of excited hydrogen and helium by integrating Equation 2 with the best-fit parameters using the apparent optical depth method (Wang et al. 2024b; Savage &Sembach 1991). The column density is then given by:  \n𝑁 = 𝑚 𝑒 𝑐 π 𝑒 2 𝑓 0 𝜆 0 𝜏, (5)  \nwhere 𝜏 is the integrated line optical depth and 𝑓 0 𝜆 0 is the product of the rest-frame wavelength and the oscillator strength of a given transition. Using this framework we find that the absorption in H 𝛼 and H 𝛽 give excited ( 𝑛 = 2) hydrogen column densities of log 𝑁 H ( n = 2 ) = 14 . 99 + 0 . 07 -0 . 07 cm -2 and 14 . 87 + 0 . 06 -0 . 06 cm -2 respectively, which are fully consistent between one another. The absorption profile of He /i.pc 𝜆 10830 yields log 𝑁 HeI ∗ = 14 . 56 + 0 . 01 -0 . 01 cm -2 . Given the relatively short lifetime of the excited hydrogen at the 𝑛 = 2 level, its large column density indicates a high total column density of hydrogen (Osterbrock & Ferland 2006).  \nThe absorption line kinematics of the fit gives 𝜎 = 126 + 6 -8 km s -1 and Δ 𝑣 = -351 + 14 -16 km s -1 for the Balmer lines and 𝜎 = 304 + 6 -6 km s -1 and Δ 𝑣 = -506 + 7 -7 km s -1 for the helium absorption lines. Estimating maximum outflow velocity as 𝑣 out = | Δ 𝑣 | + 2 𝜎 gives 𝑣 out = 604 + 11 -10 km s -1 and 𝑣 out = 1113 + 12 -12 km s -1 for the hydrogen and helium outflows respectively. We summarize our derived neutral gas outflow properties in Table 1.  \nWhile the outflow velocities inferred from the hydrogen and helium absorptions are significantly different, we note that the resonant nature of He /i.pc 𝜆 10830 could well lead to different kinematics and it is thus not a direct reflection of the outflow velocity. One might wonder whether the absorbing gas could have two components: one faster outflowing and low-density component primarily responsible for the absorption of He /i.pc 𝜆 10830, and the other slower outflowing and high-density component primarily responsible for the absorption of hydrogen Balmer lines. Regardless, we emphasize that a dense neutral outflow must exist, as we explain later in this section.  \nToinfer the physical conditions within the neutral outflow, we compute a series of simulations using the photoionization code C/l.pc/o.pc/u.pc/d.pc/y.pc (c17.03; Ferland et al. 2017) and compare them with observations. Specifically, we search for models with column densities in the 𝑛 = 2 level of the hydrogen atom and the first excited state 2 3 𝑆 of the helium atom that match our derivations.  \nTo start with, we need to assume what is the energy source that  \nShock only:  \nShock plus precursor:  \nFigure 5. Narrow emission-line diagnostic diagrams, showing our target (green circles with errorbars; the measurement uncertainties are often smaller than the size of the marker), and the distribution of galaxies at 𝑧 = 0 . 1 from SDSSDataRelease 7 (contour lines and scatter points). The demarcation lines in panels a-c are from Kewley et al. (2001); Kauffmann et al. (2003); Kewley et al. (2006); Schawinski et al. (2007). Our target lies in a scarcely populated region of the diagram, formally consistent with star-formation photoionisation (panels a and b) or AGN (panel c). Shock-driven emission in the NLR is also unlikely, because shock models (red and yellow lines) do not cover the location of our target in these diagnostic plots (see text for a description).  \n<!-- image -->  \npumps the hydrogen and helium in the absorber to their excited states. A natural candidate is the energetic photons from the AGN continuum emission in GN-28074. In principle, powerful shocks, as seen in some radio AGN, can also create excited hydrogen (e.g., Best et al. 2000). However, a shock in the NLR is heavily disfavoured by the BPT/VO diagrams (Baldwin et al. 1981; Veilleux & Osterbrock 1987) as shown in Figure 5: our target (green point with errorbars) is clearly inconsistent with a broad range of shock models indicated with yellow and red grid lines 2 .  \nOn the other hand, narrow line diagnostics can not be used to rule out a shock front in the BLR. However, shocks produced by radio jets (primary candidates for producing strong BLR shocks) are strongly disfavoured by the radio non detection, as discussed later on in Section 5.4.  \nWe use a standard AGN spectral energy distribution (SED) in C/l.pc/o.pc/u.pc/d.pc/y.pc as the source to pump hydrogen and helium to their excited  \nFigure 6. Photoionization models of pumped absorbing gas that has high column densities of hydrogen at the 𝑛 = 2 state and helium at the first excited state. Left: 𝑁 H ( n = 2 ) and 𝑁 He ( 2 3 S ) at a common depth of 𝑁 H ( Stop ) = 10 24 cm -2 for models spanning a range of hydrogen densities and ionization parameters. All models are plotted as circles color-coded by their hydrogen densities. The grey arrow indicates the increasing direction of the ionization parameter in this diagram. The star symbol represents our derived column densities for GN-28074. The open circles represent models that are capable of reproducing our derived column densities for GN-28074 at certain depths. The plausible models have high hydrogen densities and cover a narrow range in the ionization parameter. Right: hydrogen densities and column densities of plausible models colored coded by their ionization parameters. The dashed lines correspond to constant effective thickness of the absorbing gas. All the plausible models show small effective thickness of Δ 𝑙 < 10 -3 pc.  \n<!-- image -->  \nstates. The SED has a parametric form given by  \n𝐹 𝜈 = 𝜈 𝛼 uv exp (-ℎ𝜈 / 𝑘𝑇 𝐵𝐵 ) exp ( 0 . 01 Ryd / ℎ𝜈 ) + 𝑎𝜈 𝛼 x , (6)  \nwhere we set the big blue bump temperature to 𝑇 𝐵𝐵 = 10 6 K, the UVto-X-ray slope to 𝛼 ox = -1 . 4 (which is achieved by adjusting 𝑎 by the code), the UV slope to 𝛼 uv = -0 . 5, and the X-ray slope to 𝛼 x = -1 . 0 3 . We compute a total of 272 models spanning a range of hydrogen densities of 10 2 cm -3 ≤ 𝑛 H ≤ 10 10 cm -3 , and a range of ionization parameters of 10 -4 . 5 ≤ 𝑈 ≤ 10 0 . The dimensionless ionization parameter, 𝑈 , is defined as 𝑈 ≡ Φ 0 𝑛 H 𝑐 , where Φ 0 is the number flux of the hydrogen ionizing photons (impinging on the inner face of the cloud) and 𝑐 is the speed of light. For all of the models, we adopt a metallicity of 12 + log ( O / H ) = 7 . 85 (which is derived using narrow line fluxes of GN-28074, as described in Appendix B), and assuming the solar abundance pattern follows the values in Grevesse et al. (2010). While the metallicity in the absorber might be different from that in the forbidden-lines emitting region, overall the metallicity has negligible impact on our conclusions (as we are focusing on hydrogen and helium absorption). We do not include dust grains in our models. While dust can also absorb ionizing photons, at low metallicities, the amount of dust is significantly reduced in the ISM and its impact on the ionization structure of the gas becomes negligible. Furthermore, as we show later, the absorber has a density similar to the gas in the BLR, which is usually believed to be dustless (because it is inside the sublimation radius).  \nWe have included the effect of microturbulence in our simulations and varied the turbulent velocity over a range of 100 km s -1 ≤ v turb ≤ 300 kms -1 to match the velocity dispersions we fit for the hydrogen and helium absorptions. The microturbulence broadens the line profile and increases the chance of absorption. However, within  \nthe range of turbulent velocity we choose, the exact value of the turbulent velocity has little impact on our final results.  \nFinally, all models are truncated at a total hydrogen column density of 𝑁 H = 10 24 cm -2 . We search within each model whether there is a specific depth at which the cumulative column densities of the excited hydrogen and helium match observations within 2 𝜎 ( ∼ 0 . 14 dex).  \nThe left panel of Figure 6 shows the total column densities of excited hydrogen (at 𝑛 = 2) and helium (at 2 3 𝑆 ) in our models at a common depth of 𝑁 H = 10 24 cm -2 . Since this is the common stopping criterion of our models, the column densities shown are the highest values achievable by the models. The models are color-coded by their hydrogen volume densities and the direction of increasing ionization parameter is indicated by the arrow. The models that can reproduce the column densities we derive for GN-28074 (open star) within certain depths, quantified by column densities reached ( ≲ 𝑁 H = 10 24 cm -2 ), are highlighted by open circles (indeed the different combinations of 𝑁 H ( n = 2 ) and 𝑁 He ( 2 3 𝑆 ) are degenerate in reproducing the observed absorptions equally well). Clearly, the high 𝑁 H ( n = 2 ) we derive for GN-28074 strongly constrains the hydrogen density of the absorber to high values. This is because higher densities allow 𝑁 H ( n = 2 ) to build up at smaller depths. In addition, high densities boost collisional excitation of hydrogen, which helps to achieve high 𝑁 H ( n = 2 ) . Meanwhile, our derived 𝑁 He ( 2 3 S ) further narrows down the plausible range of densities and strongly constrains the ionization parameter.  \nWith these results, we go back to consider the possibility that the Balmer absorption and the He /i.pc absorption are produced by different absorbing clouds. From Figure 6, one can see the 𝑁 H ( n = 2 ) of GN-28074 alone would indicate a high 𝑁 He ( 2 3 S ) in the cloud that produces the Balmer absorption, unless the density is even higher. Therefore, the most likely picture would be either a single absorber making both absorptions, or two separate absorbers with the Balmerline absorber having a higher density. Either way, the hydrogen density is constrained to be 𝑛 H ≳ 10 8 cm -3 , indicative of a BLR origin. In what follows, we set our discussion in the context of a single absorber for simplicity, and we note that the more complex two-cloud possibility would not change our main conclusions.  \nIn the right panel of Figure 6, we show the densities and column densities of the plausible absorbing models for GN-28074. We note that the column density in this panel is not the total column density of the models (which has been fixed to 𝑁 H = 10 24 cm -2 ), but the cumulative hydrogen column densities in certain depths of the models where both column densities of He (2 3 𝑆 ) and H ( 𝑛 = 2) match the values of GN-28074. These are the minimum values of the actual column density of the cloud in order to match the observation, as the cloud can have an additional column of gas not playing a role in the absorption seen in the Balmer and He /i.pc lines (because the 𝑛 = 2 and 2 3 𝑆 levels are not properly excited in these regions). The models show an anti-correlation between the density and column density, which is expected since higher densities help to populate hydrogen to the 𝑛 = 2 level through collisional excitation. Overall, the plausible hydrogen (minimum) column densites are large, ranging from 10 21 . 5 cm -2 to 10 23 cm -2 given the range of densities we consider. The large column density facilitates Ly 𝛼 trapping, which is another important mechanism to produce excited hydrogen at 𝑛 = 2 (Ferland & Netzer 1979; Hall 2007). The physical scale of the absorber can be approximated as Δ 𝑙 = 𝑁 H 𝑛 H 𝜖 , where 𝜖 is the gas volume filling factor. If the absorber has 𝜖 ∼ 1, the physical scale implied by the matched models would be Δ 𝑙 ≲ 10 -3 pc, compatible with a single cloud within the BLR.  \nWe note that, for a Galactic gas-to-dust ratio, the very minimum column density of 𝑁 H > 10 21 . 5 cm -2 inferred for the absorbing clouds, would imply a minimum dust absorption of 𝐴 𝑉 > 2 . 2. This is inconsistent with the reddening inferred from the narrow lines, and also essentially inconsistent with the upper limit inferred from the broad lines. This is not a novel finding, as in AGN it has been known for a long time that dust absorption is much lower than what would be inferred from the absorbing gaseous column density (Maiolino et al. 2001). One of the primary interpretations of this discrepancy is that much of the gaseous column density (as typically inferred from the X-rays) is associated with (dust-free) BLR clouds (Risaliti et al. 2002, 2005; Maiolino et al. 2010). This perfectly fits the scenario that is emerging for GN-28074, in which we are finding evidence for absorbing clouds with density typical of the BLR clouds. As we will see below, these are also estimated to be within the sublimation radius.', '5.2.2 Location of the absorber': 'The inferred gas densities are consistent with those of the BLR clouds, or clouds close to the BLR (Netzer 1990).  \nIn principle, from the narrow range in log 𝑈 derived for the absorber, and from the bolometric luminosity derived in subsection 5.1 from the broad component of H 𝛼 , we can derive the distance of the absorber from the accretion disk. Assuming the absorbing clouds are embedded in tenuous and optically thin gas (i.e., a two-phase BLR model; Krolik et al. 1981), and assuming the ionizing fraction and average ionizing photon energy from the C/l.pc/o.pc/u.pc/d.pc/y.pc SED (which has a representative optical-to-X-ray slope of 𝛼 ox = -1 . 4), we use  \n𝑟 abs = √︄ 𝐿 bol 𝑓 0 4 π 𝑐𝑛 H 𝑈 ⟨ ℎ𝜈 0 ⟩ (7)  \nto estimate the location of the absorber. In this equation 𝑓 0 is the ionizing fraction and ⟨ ℎ𝜈 0 ⟩ is the average ionizing photon energy. We find the distance between the absorber and the source of the ionizing continuum to be from 1.9 pc to 0.19 pc, for 𝑛 H = 10 8 cm -3 to 𝑛 H = 10 10 cm -3 (this is without dust correction, but this and the following estimates of the sizes all scale as 𝐿 0 . 5 bol , and thus the relative  \nlocations are not affected by dust correction). These sizes (although extremely uncertain, given the several assumptions) are larger than the radius of the BLR estimated by using the scaling relation from Greene & Ho (2005), which would give 𝑟 BLR ∼ 0 . 036 pc. Inferring the BLR size from the H 𝛽 luminosity scaling relations of Kaspi et al. (2005) gives a consistent radius of 0.04 pc.  \nWe can also compare the estimated radial location of the absorber with the expected size of the dusty torus. Assuming a dust sublimation temperature of 𝑇 sub . ≈ 1 , 500 K, the thermal equilibrium condition  \ngives the inner radius of the torus to be 𝑟 torus = √︂ 𝐿 bol 4 π 𝜎𝑇 4 sub . ∼ 0 . 21  \npc without dust correction. Therefore, the location of the absorbing cloud can be between the BLR and the torus. Interestingly, this inferred location is consistent with previous interpretations of Balmer-line absorbing AGN at lower redshifts (e.g., Hutchings et al. 2002; Zhang et al. 2015), despite the fact that our derived column densities for GN-28074 are higher. However, we recall that these estimations, especially the location of the clouds based on the ionization parameter, are extremely uncertain. This is because we have assumed optically thin conditions for clouds inner to the absorber to translate the ionization parameter to the distance, which ignores the potential self-shielding of the continuum emission in the BLR (Gaskell 2009). It is still entirely possible that the absorbing clouds are actually associated with the outer parts of the BLR.  \nOne might wonder what is the physical state of the absorbing gas. In our C/l.pc/o.pc/u.pc/d.pc/y.pc simulations, the absorber has a warm temperature of 10 , 000 ∼ 20 , 000 K and is nearly isobaric, which places it in a thermally stable regime (Reynolds & Fabian 1995).  \nTo check the dynamical state of the absorber, we compare the integrated radiation pressure from the attenuated continuum emission with the gravitational pressure. The radiation pressure is from our C/l.pc/o.pc/u.pc/d.pc/y.pc models, integrated until the depths that match the derived column densities. The gravitational force induced solely by the central black hole is  \n𝐹 grav = 4 π 𝐶 𝑓 𝐺𝑀 BH 𝜇𝑚 p 𝑁 H , (8)  \nwhere 𝜇 = 1 . 4 is the mean atomic mass number and 𝑚 p is the proton mass. The ratio between the radiation pressure and the gravitational pressure is then given by  \n𝑃 rad 𝑃 grav = 𝑃 rad / 𝐺𝑀 BH 𝜇𝑚 p 𝑁 H 𝑟 2 abs ! , (9)  \nwhichlies in a range of 0 . 72 -25. Despite the large column density, the integrated radiation pressure upon the dense gas due to the attenuated continuum emission is comparable and even larger than the gravitational pressure, consistent with a radiatively driven outflow scenario. The radiation pressure is only lower than the gravitational pressure when the column density of the absorber reaches 𝑁 H ∼ 10 23 cm -2 . Still, we caution that this estimate does not take into account the possibility that the continuum emission is partially shielded by BLR clouds inner to the absorber, which would in principle reduce the radiation pressure.  \nAs we have discussed, the outflow velocity traced by He/i.pc absorption is larger than that traced by Balmer absorption. If we assume the difference is real despite the complexity from the resonant scattering of He /i.pc , this might imply a decelerating outflow as the Balmer absorption can preferentially occur in a more neutral and thus outer region. The deceleration could happen as the outflow sweeps up the surrounding gas. This is in contrast to the outflowing Balmer-line and He /i.pc -line absorber in a nearby AGN, NGC 4151, where the outflow is instead accelerating outwards (Hutchings et al. 2002). The decelerating scenario could indicate that the outflow will  \neventually fall back to the BLR and thus be self-replenished, maintaining the high covering factor. In fact, if one simply estimates the escape velocity for the absorber based on the distances derived from the ionization parameter, one obtains  \n𝑣 esc = √︄ 2 𝐺𝑀 BH 𝑟 abs ≈ 800 -3000 km s -1 , (10)  \nwhere the lowest 𝑣 esc corresponds to 𝑛 H = 10 8 cm -3 and the highest 𝑣 esc corresponds to 𝑛 H = 10 10 cm -3 . Comparing these values to the outflow velocities in Table 1 which lies within 600 -1000 km s -1 , it is indeed plausible that while the gas is outflowing at the time of observations, it will eventually stall and fall back. In addition, considering a two-phase or multi-phase BLR model (e.g., Krolik et al. 1981; Reynolds & Fabian 1995), as the dense absorber propagates outwards, the ram pressure from the volume-filling and highly ionized phase of the gas would help to slow and trap the outflow. The ram pressure would become important when the ambient gas has a number density of 𝑛 H ≳ 10 4 cm -3 , consistent with tenuous gas with temperatures close to the Compton temperature and in pressure equilibrium with the typical line-emitting clouds in the BLR. Furthermore, the absorber may continuously sweep up the gas in the outward direction, leading to increasing mass, which further works to slows the outflow. This would, in principle, also lead to strengthened absorption, but the effect might not be significant over a length scale of the size of the BLR, if the ambient density is 𝑛 H ∼ 10 4 cm -3 .', '5.2.3 Nuclear neutral outflow': "Finally, we investigate the properties of the nuclear outflow traced by the blueshifted absorption features. Estimating the total mass of the outflowing gas is not entirely straightforward as the absorbing clouds do not necessarily trace the entirety of outflowing gas.  \nWe then use the following equation derived for a spherically symmetric outflow with 𝐶 𝑓 = 1 to estimate the outflow rate:  \n/ 𝑀 out;BLR = 4 π 𝑅𝜇𝑚 p 𝑁 H 𝑣 out , (11)  \nwhere 𝜇 = 1 . 4 is the mean mass per proton and 𝑁 H is the hydrogen column density. We can take the size of the outflowing medium as estimated in the previous subsection, and a minimum column density of 𝑁 H > 10 21 . 5 cm -2 , as inferred by our models. We use both velocity values from Table 1 in conjunction. The inferred lower limit on the nuclear outflowing mass is 0 . 05 M ⊙ yr -1 .  \nThis is a very conservative lower limit, as we are assuming the lowest of the radius estimates above and the bare minimum of the possible column densities. However, it is interesting to note that, based on the luminosity of the broad H 𝛼 line and assuming a BLR gas density of 10 9 -10 10 cm -3 , the inferred mass of the BLR is about 10 -100 𝑀 ⊙ . If the nuclear outflow is tracing the outer part of the BLR, it would imply that the BLR would be wiped away in a few hundred/thousand years. These are very short timescales, suggesting that the BLR clouds are continuously replenished (as inferred by other independent evidence Maiolino et al. 2010), probably lifted from the accretion disc, or simply through the eventual falling back of the outflow gas. An alternative possibility is that they are the bloated atmospheres of nuclear giant stars (Scoville & Norman 1988; Alexander & Netzer 1994), in which case those stars would continuously replenish the 'clouds' on long timescales.  \nFor the decelerating nuclear outflow scenario, while the outflowing dense gas would be recycled at a later time, it is currently radiationpressure dominated and it might expand and become rarefied, creating a warm ambient medium, with a temperature below the Compton  \ntemperature of ∼ 10 7 K. It has been postulated that the warm gas that produces X-ray absorption features in AGN could be a product of a rarefied dense outflow (Reynolds & Fabian 1995). The warm clouds could have a temperature around 10 5 -10 6 K and lie at a location within or just outside the BLR (similar to what we find for the Balmer-line absorber), as inferred from modeling of previous observations of X-ray absorption in AGN (Reynolds & Fabian 1995; Reynolds 1997; Mathur et al. 1997; Komossa & Fink 1997; Nicastro et al. 1999; Chakravorty et al. 2009; Laha et al. 2021). Thus, under the isobaric condition, they would have a density of 𝑛 H ≳ 10 6 cm -3 . While not contributing to Balmer absorption, the warm clouds can contribute to the obscuration of the X-ray emission at the location of the absorber. Furthermore, the stall radius of the outflow can be calculated from the outflow velocity and the location of the absorber, which has a minimum value around 0.3 pc. If the warm medium extends to the stall radius, can easily build up a column density of 𝑁 H;warm ≳ 10 24 cm -2 and becomes Compton-thick. While more quantitative calculations are needed to verify such a scenario, it provides an interesting picture for maintaining the X-ray weakness.", '5.3 X-ray properties': 'GN-28074 is located in a region of the Chandra Deep Field North (CDF-N) with excellent sensitivity within the context of the 2Ms Chandra survey (Alexander et al. 2003). Its X-ray analysis is reported in detail in Maiolino et al. (2024a), where they infer an upper limit of < 2 . 9 × 10 41 erg s -1 on the 2-10 keV luminosity. When compared with the bolometric luminosity inferred from the broad lines, this implies a 𝐿 bol / 𝐿 𝑋 > 15400, placing the object far above the standard 𝐿 bol -𝐿 𝑋 relation from Duras et al. (2020), and actually making it the most X-ray weak among all AGN discovered by JWST and analysed in Maiolino et al. (2024a), with the measured upper limit being 3 dex lower than the expected X-ray luminosity of ∼ 10 44 erg s -1 . Fig.7 illustrates this more quantitatively by showing the distribution of the 𝐿 𝐵𝑜𝑙 / 𝐿 𝑋 for JWST-selected AGN (amber histogram) relative to the low-z relation, for optically selected AGN (black hollow histogram), and keeping in mind that the values for the JWST-selected AGN are mostly lower limits on 𝐿 𝐵𝑜𝑙 / 𝐿 𝑋 . The latter are clearly offset towards muchlarger 𝐿 𝐵𝑜𝑙 / 𝐿 𝑋 relative to low-z, optically-selected AGN, with GN-28074 being the most extreme of them all.  \nThe X-ray weakness of GN-28074 is comparable to what seen in SBS0335, a local metal poor dwarf galaxy, which was recently found to host an AGN, based on the detection of broad H 𝛼 and variable mid-IR emission (Hatano et al. 2023), but which is extremely weak in the hard X-rays (the weak detection is probably associated with the star forming host).  \nOf the scenarios discussed in Maiolino et al. (2024a) to explain the X-ray weaknesses of the JWST-discovered AGN, Compton thickness appears most likely, given the independent constraints on hydrogen column densities given by the absorption lines. However, we can not entirely discount intrinsic X-ray weakness, particularly given the extreme 𝐿 bol / 𝐿 𝑋 ratio, and additional information inferred from the radio emission discussed in the next section.', '5.4 Radio properties': "GN-28074 has been observed at 1.5 GHz, as part of the JVLA survey reported by Owen (2018) and also part of the e-MERGE Survey (Muxlow et al. 2020). While the radio constraints on this and other JWST-discovered AGN will be discussed in more detail in a separate paper (Mazzolari et al. in prep.), here we simply report that the source  \nFigure 7. Distribution of the bolometric to X-ray luminosity ratio, relative to the relation inferred for low-redshift optically-selected AGN (whose distribution is shown with a black, hollow histogram). The distribution of JWST-discovered AGN in GOODS (the vast majority of which are actually lower limits on 𝐿 𝐵𝑜𝑙 / 𝐿 𝑋 ) is shown with the solid, amber histogram. The lower limit for GN-280974 is indicated with the blue bar. We also show the lower limit observed in SBS 0335-052E, which is a nearby AGN with similar properties as the new population of the JWST-discovered AGN. We also show the value for the 'prototypical' Seyfert 1 NGC-4151 which shows H 𝛼 , H 𝛽 and HeI absorption, just like GN-28074, and strong (and variable) X-ray absorption.  \n<!-- image -->  \nis not detected at 1.5 GHz, with an upper limit of 1 . 6 × 10 39 erg s -1 on the luminosity at 1.5 GHz.  \nBased on the bolometric luminosity inferred from the broad lines, and using the bolometric-to-radio relationship for radio quiet AGN (Bariuan et al. 2022; Wang et al. 2024a) one would expect a radio luminosity that is at least a factor between 5 and 20 larger than observed.  \nOnce again, this radio weakness seems in line with the extreme radio weaknesses observed in SBS 0335 (Mazzolari et al. in prep.), which seems to be a local analogue of GN-28074 based on the X-ray and radio properties.  \nThe radio weakness may be connected with the X-ray weakness in two possible ways.  \nThe clouds responsible for the H 𝛼 absorption are also completely optically thick to free-free absorption. Indeed, the free-free absorption coefficient is given by  \n𝜅 𝑓 𝑓 = 3 . 3 × 10 -7 GLYPH<18> 𝑛 e cm -3 GLYPH<19> 2 GLYPH<18> 𝑇 e 10 4 K GLYPH<19> -1 . 35 GLYPH<16> 𝜈 GHz GLYPH<17> -2 . 1 pc -1  \nWhile the temperature of the clouds is always in the range between 5,000-25,000 K, the density of the clouds is so high that (despite their thickness of only 10 -6 -10 -3 pc) the expected optical thickness for free-free should be of the order of 𝜏 𝑓 𝑓 = 10 5 -10 6 . Therefore, the radio emission should be completely absorbed by the clouds that absorb H 𝛼 . However, obviously the question is whether the radio emission (possibly associated with small jets, in radio quiet AGN), is compact enough to be entirely covered by such clouds. This remains an open question, however we note that in local AGN evidence for nuclear, sub-pc scale jets being free-free absorbed has been seen (e.g. Taylor & Vermeulen 1996; Meyer et al. 2024b).  \nAnother possibility is that the magnetic field in the nuclear region of these early AGN has not yet developed properly. Some recent zoom-in simulations do expect a weaker nuclear magnetic field in the  \nearly phases of galaxy formation (S. Martin Alvarez, priv. comm.). A weak (or absent) nuclear magnetic field would obviously prevent significant radio emission as nuclear magnetic fields are believed to be a significant mechanism of launching AGN jets (Event Horizon Telescope Collaboration et al. 2019; Zamaninasab et al. 2014). At the same time, the hot corona responsible for the hard X-ray emission is thought to be produced by gas lifted from the accretion disc by the magnetic field (Haardt & Maraschi 1991; Cheng et al. 2020). Therefore, a weak magnetic field would also automatically result in to an intrinsic X-ray weakness.  \nLastly, it may be possible that this AGN has relatively recently launched a radio jet, which has weakened the corona resulting in Xray weakness. Such a jet may be visible at MHz frequencies, however, followup radio observations are required to test this scenario.", '6 RULING OUT ALTERNATIVE SCENARIOS FOR THE BROAD LINES': 'To explain the puzzling nature of JWST-confirmed AGN similar to GN-28074 that lack hard X-ray emission, it has been proposed recently that some of these sources are not real AGN but are SF-driven processes in disguise (e.g., Yue et al. 2024b; Kokubo & Harikane 2024). In this section, we discuss the possibility of alternative scenarios that create broad permitted lines.', '6.1 Ruling out hyperdense ultra-metal-poor galactic outflows': "The broad component of the permitted lines (H 𝛼 , H 𝛽 , HeI, Pa 𝛽 ) does not have a counterpart in the forbidden lines and, in particular, in bright [O /i.pc/i.pc/i.pc] 𝜆 5007. This has always been considered as a clear evidence for emission from the BLR, where the densities are much higher than the critical density of [O /i.pc/i.pc/i.pc] 𝜆 5007.  \nHowever, it has recently been suggested that the broad permitted lines (H 𝛼 , H 𝛽 , He /i.pc , Pa 𝛽 ) discovered by JWST might actually not be tracing the BLR of AGN but galactic outflows (Yue et al. 2024b; Kokubo & Harikane 2024). In this subsection we discuss that this scenario is totally untenable, based on multiple lines of evidence, specifically for GN-28074 and, more broadly, for the several other cases discovered by JWST.  \nAs already discussed, there is the indication of a galactic outflow component to all emission lines, with width intermediate between the narrow and the broad components of the permitted lines. However, here we are considering the more extreme scenario suggested by Yue et al. (2024b) and Kokubo & Harikane (2024) that also the very broad component of H 𝛼 , H 𝛽 , Pa 𝛽 , and He /i.pc is also tracing a (much faster) galactic outflow.  \nIn the outflow scenario there would be two possibilities for potentially suppressing the broad component of [O /i.pc/i.pc/i.pc] 𝜆 5007, which are an extremely low metallicity or a very high density. We discuss in the following that both scenarios are not possible.  \nAs it is a strong coolant of the ISM, the intensity of the [O/i.pc/i.pc/i.pc] 𝜆 5007, relative to the Balmer lines, actually increases as the metallicity decreases, down to about 0.15-0.2 Z ⊙ (Maiolino et al. 2008; Curti et al. 2017), which is about the gas metallicity of GN-28074, as measured from our spectra and discussed in Appendix B. Therefore, in the outflow scenario, the broad component of [O /i.pc/i.pc/i.pc] 𝜆 5007 should be even more prominent than H 𝛼 and H 𝛽 , contrary to what observed. As already discussed in Maiolino et al. (2024a), in order for the putative broad component of [O /i.pc/i.pc/i.pc] 𝜆 5007 to become significantly fainter than the corresponding broad component of the Balmer lines, the metallicity of the outflowing gas should  \nbe less than 0.01 Z ⊙ , which is more than one order of magnitude lower than the host galaxy generating the outflow. This is the opposite of what seen in all galactic outflows, in which the metallicity of the outflowing gas is higher than the host galaxy.  \nAlternatively, the density of the outflowing gas should be much higher than the critical density of [O /i.pc/i.pc/i.pc] 𝜆 5007. Specifically, in order not to be detectable the gas density should be higher than 10 7 cm -3 . With such high density, the implied column of gas in the outflow, assuming an extension of about 10-100 pc (i.e., about the size of the stellar component that should be driving the outflow) would result into a column of gas of about 10 27 cm -2 . For a dust-to-gas ratio expected at a metallicity of 0.15 Z ⊙ (see, e.g., D'Eugenio et al. 2024), would imply a dust attenuation of 𝐴 𝑉 ∼ 2 × 10 4 . This would totally absorb the receding side of the outflow, leaving visible only the blueshifted side, which is completely in contrast with the symmetric profile of all broad components. For the records, in galactic ouflows it only needs 𝐴 𝑉 ∼ 1 to result into a prominently blueshifted profile of the lines, especially in H 𝛽 . In order to comply with the hyper-dense outflow scenario, the symmetric profiles of the broad lines would require essentially no dust in the outflowing gas, or, equivalently, a metallicity of ∼ 10 -4 𝑍 ⊙ , coming out of the nucleus of a galaxy with a metallicity of 0 . 1 𝑍 ⊙ .  \nThe additional piece of evidence against the hyper-dense, ultrametal poor outflow scenario is the luminosity of the broad lines. If ascribed to an outflow, then the luminosity of the broad component of H 𝛼 should be associated with photoionization by star formation. Even in the most conservative extinction scenario, the H 𝛼 luminosity would imply an SFR ∼ 270 𝑀 ⊙ yr -1 (Pflamm-Altenburg et al. 2007). Such a huge SFR would have massive observable implications at other wavelengths. It should be a bright submm source (while this galaxy is undetected in deep SCUBA2 surveys of GOODS-N, Geach et al. 2017) and it should have a 24 𝜇 mflux at least one order of magnitude higher than observed. Additionally, the X-ray emission associated with SNe and X-ray binaries should have a luminosity higher than 1 . 4 × 10 42 erg s -1 , while the upper limit from the non-detection in the Chandra 2Ms observation is about an order of magnitude lower. Finally, given the small size of the source, the implied surface density of SFR would be about 3 × 10 4 M ⊙ yr -1 kpc -2 ; such humongous star formation rate density has never been observed in any system and is orders of mangnitude higher than the theoretical maximum startburst limit of 10 3 M ⊙ yr -1 kpc -2 , given by the Eddington limited starbursts (Thompson 2009).  \nSimilar considerations apply for many other galaxies with broad H 𝛼 discovered by JWST.", '6.2 Ruling out Raman scattering by neutral hydrogen': 'Another non-AGN scenario that can produce broad permitted emission in X-ray weak sources, as suggested recently by Kokubo & Harikane (2024), is the inelastic Raman scattering of continuum emission bluewards of Ly 𝛼 by neutral hydrogen at the ground state, where the continuum emission does not need to have an AGN origin. This phenomenon was initially discussed to explain the emission observed in symbiotic stars (Schmid & Schild 1990, 1994; Schmid 1992). Since the Raman scattering cross-sections are relatively high around the wavelengths of Lyman series with broad wings, it can lead to broad absorptions around Lyman lines and subsequently broad emission around Balmer lines and other hydrogen recombination lines remembling a BLR component.  \nSignificant Raman scattering requires the existence of neutral hydrogen at the 1 𝑠 state with 𝑁 H ≳ 10 21 cm -2 , which is possible for the case of GN-28074 given the presence of the partially neu-  \nFigure 8. Comparison between the FWHM ratios of the broad hydrogen lines measured in the NIRSpec R1000 spectrum of GN-28074 and those predicted by a single-cloud inelastic Raman scattering model by Kokubo (2024). Clearly, the profiles of the broad hydrogen lines measured in GN28074 cannot be explained by the Raman scattering of the UV continuum by neutral hydrogen.  \n<!-- image -->  \norber. However, as shown by Kokubo (2024), Raman scattering tends to create different line profiles for different hydrogen lines, with FWHM ( H 𝛼 ) ≈ 3 . 4 FWHM ( H 𝛽 ) ≈ 2 . 6 FWHM ( Pa 𝛽 ) for a single-cloud model (see Equations 32a-e in Kokubo 2024, assuming the FWHM scales with the velocity width of the line). In Figure 8 we compare the predicted line widths from a single-cloud Raman scattering model with those measured in the R1000 spectrum of GN-28074. While the FWHMs of hydrogen lines increase with increasing column densities under Raman scattering, the FWHM ratios are roughly constant, which means the comparison we make in Figure 8 can be used to test Raman scattering in a wide range of cloud depths. Based on our measurements for the broad emissionline component in GN-28074, FWHM ( H 𝛼 ) broad = 3610 ± 22 km s -1 and FWHM ( Pa 𝛽 ) broad = 2820 ± 31 km s -1 (see Table 3). These give FWHM ( H 𝛼 ) broad / FWHM ( Pa 𝛽 ) broad = 1 . 28 ± 0 . 02 ≪ 2 . 6, which is significantly lower than the prediction from the Raman scattering. We note that while the measured width of the broad H 𝛼 is indeed larger than that of the broad Pa 𝛽 , such a level of difference is not unexpected for a BLR origin based on previous observations of AGN (Landt et al. 2008; Lamperti et al. 2017). As suggested by Kim et al. (2010), the narrower Paschen lines could be explained by the scenario where they are preferentially emitted by clouds on the outskirts of the BLR compared to Balmer lines. In addition, the broad H 𝛼 and H 𝛽 lines are well fitted with the same width profiles, thus FWHM ( H 𝛼 ) broad / FWHM ( H 𝛽 ) broad ≈ 1 ≪ 3 . 4. When fitting their kinematics separately, FWHM ( H 𝛼 ) broad = 3600 + 22 -23 km s -1 while FWHM ( H 𝛽 ) broad = 4180 + 150 -160 km s -1 , giving a ratio of 0 . 86 ± 0 . 03, which is again inconsistent with Raman scattering. It should be noted here that, while the difference in line widths is significant ( ∼ 3.5 𝜎 ), separating the kinematics of H 𝛼 and H 𝛽 results in a worse fit in the outflow component for [O /i.pc/i.pc/i.pc] as well as creating significant residuals around the locations of Balmer absorptions. Thus, we prefer to keep the kinematics tied.  \nAlso, the broad component under Raman scattering would be tilted towards the red wing especially for Balmer lines (Kokubo 2024),  \nwhile the broad lines observed in GN-28074 are symmetric. The flux ratios between the broad lines would significantly deviate from Case B values, although this could be hidden by dust attenuation. Furthermore, such a high column density of gas to create significant Raman scattering ( 𝑁 H ≳ 10 21 cm -2 ) needs to be nearly dustless, otherwise the UV continuum needed for the inelastic scattering would be easily attenuated by dust. Clearly, there is dust in GN-28074 as indicated by the ratios of narrow lines.  \nFinally, the Raman scattering cannot explain the broad component in the He /i.pc 𝜆 10830 line as well as the broadening of the O /i.pc 𝜆 8446 line (see Figure A1 in Appendix A). Therefore, we conclude the broad components in the spectrum of GN-28074 cannot be produced by the Raman scattering and the BLR explanation is still preferred.  \nSimilar considerations apply for other galaxies with detections of broad H 𝛼 and other broad Balmer and Paschen lines discovered by JWST, although a careful comparison will be needed due to the weakness of broad lines other than H 𝛼 in these low-luminosity AGN.', '7 FRACTION OF JWST-DISCOVERED AGN WITH BALMER ABSORPTION': "As mentioned in the introduction, H 𝛼 absorption has been detected in many of the broad-line AGN discovered by JWST (Matthee et al. 2024; Kocevski et al. 2024; Wang et al. 2024c; Maiolino et al. 2023b). It is difficult to assess the fraction, due to the different selection methods and observing modes. However, as already pointed out by Maiolino et al. (2024a), the fraction of broad line AGN with H 𝛼 absorption in each of those surveys seems to be about 10%.  \nHowever, we argue that the true fraction is likely much higher still. Indeed, in order to see H 𝛼 in absorption with JWST/NIRSpec, lowresolution prism spectra (which are the most commonly used) are totally inadequate, and also the medium resolution grating (R ∼ 1000) might not have high enough resolution if the absorption features are less than a few hundred km s -1 in width. The high resolution gratings (R ∼ 2700) are in principle more appropriate, but this is the least used observing mode (primarily because spectra are truncated when this grating is used with the MSA); moreover, in many cases the resolution may come at expense of signal-to-noise ratio and only the brightest targets may allow the detection of absorption features. Finally, the absorption of H 𝛼 may be detectable only if associated with some significant blueshift (or redshift), while, if close to the systemic redshift, it would be difficult to disentangle from the narrow component of H 𝛼 . Indeed, the JADES survey is starting to reveal some cases of H 𝛼 in absorption with low velocity offset that is detectable simply because it is somewhat broader than the narrow component (D'Eugenio, priv. comm.). Therefore, it is likely that a much larger fraction of the type 1 AGN discovered by JWST have intrinsic H 𝛼 absorption, but it is not detected because of the inadequate observing mode used, insufficient S/N, or because the absorption is hidden by the emission features.  \nSummarising, the fraction of broad line AGN with H 𝛼 absorption is likely much larger than 10%, among the JWST-discovered AGN. This is much higher than found in previous surveys ( ∼ 0 . 1%), such as the SDSS at low redshift. This may fit the scenario suggested by Maiolino et al. (2024a) that the intermediate/low luminosity AGN discovered by JWST at high redshift may be characterised by a large covering factor of the BLR clouds, whose (dust-free/dust-poor) absorption in the X-rays may be responsible for their observed X-ray weakness.", '8 LOCAL ANALOGUES': "It has been pointed out that high-z AGN discovered by JWST are peculiar in being extremely X-ray weak (Maiolino et al. 2024a) and probably also radio weak (Mazzolari et al. in prep.), relative to the low redshift, optically selected (and, obviously, X-ray and radio selected AGN). In these terms, GN-28074 is among the most extreme JWST-discovered AGN, in the sense that it is clearly a luminous type 1 AGN with luminous broad permitted lines (H and He /i.pc ), strong mid-IR excess, but undetected in the X-ray and radio. Yet it is important to explore whether this prototype of JWST-discovered AGN has analogues in the local Universe, even if they may be rare (or possibly, overlooked).  \nTo begin with, it is interesting that what is considered to be the type 1 AGN 'prototype', NGC 4151, is actually characterised by blueshifted H 𝛼 , H 𝛽 and He /i.pc absorption, just as in GN-28074. In NGC 4151 there is clear evidence for time variable absorption in the X-ray; while the X-ray absorption does not reach the Compton thick values inferred for GN-28074, it does reach a few times 10 23 cm -2 , which can absorb substantially the hard X-rays (Zoghbi et al. 2019). More specifically, in the high-absorption regime, the observed 𝐿 bol / 𝐿 𝑋 ratio of NGC 4151 is in the regime of the JWSTdiscovered AGN at high-z (Fig.7), although not as extreme as GN28074.  \nMore compelling is the case of SBS0335-052E. This is a low metallicity (1/40 Z ⊙ ) dwarf galaxy that, until recently was thought not to host an AGN. However, it was recently found to host an AGN based on both the detection of broad permitted emission lines and variability in the mid-IR (Hatano et al. 2023). Interestingly, the AGN is undetected in the X-ray (the weak detection is ascribed to the host galaxy, hence used as an upper limit) and placing it at the same level of 𝐿 obs / 𝐿 𝑋 as GN-28074, as illustrated in Fig.7. In addition, SBS0335052E is also undetected in the radio at 1.5 GHz (once again the weak detection is likely associated with the host galaxy and used as upper limit), and with a radio luminosity much weaker than expected even for radio a quiet SED. This will be discussed further in Mazzolari et al. (in prep.). The similarity between GN-28074 and SBS0335052E in terms of broad emission lines, mid-IR excess, and X-ray and radio weakness, suggests similarities and connections between the AGN discovered by JWST with low metallicity dwarf galaxies. Indeed, more of these dwarf galaxies have been found to host AGN based on the detection of both high ionization lines and a broad component of the H 𝛼 (without forbidden counterpart), while being undetected in the X-rays (Hatano et al. 2024; Simmonds et al. 2016). Most interestingly, some of these show evidence for H 𝛼 absorption (Burke et al. 2021), just as in GN-28074, suggesting that also in (at least some of) these cases the X-ray weaknesses may be due to a dust-poor, Compton thick absorber.  \nMore recently, Arcodia et al. (2024) used eROSITA to explore the X-ray emission of AGN identified in dwarf galaxies identified via nuclear variability in the UV-optical-IR. Although they detect 10% of them, the remaining are undetected even in the stacking, just like the broad (and narrow) line JWST-discovered AGN at high redshift (Maiolino et al. 2024a). They suggest that these light AGN (hosting intermediate mass black holes) may be missing the canonical Xray corona and/or undergoing unusual accretion modes resulting in non-standard spectral energy distributions.", '9 DISCUSSION AND CONCLUSION': "GN-28074 can be considered an optimal case study to explore the nature and properties of the new population of AGN that is being  \ndiscovered by JWST at high redshift. Indeed, it shares many of their features and, being among the most luminous and at a lower redshift, it allows to characterize features that are more difficult to investigate in the more distant and fainter population.  \nGN-28074 has prominent permitted broad lines (H 𝛼 , H 𝛽 , Pa 𝛽 , He/i.pc ) without a counterpart in the forbidden lines (primarily [O /i.pc/i.pc/i.pc]), which clearly indicate that these originate from the Broad Line Region of a type 1 AGN with luminosity 𝐿 bol ∼ 6 × 10 45 erg s -1 . Its SED also has a clear mid-IR excess tracing hot dust emission typical dusty torii surrounding AGN.  \nAs with the vast majority of the AGN newly discovered by JWST, GN-28074 is not detected in the X-rays and, it is actually the most extreme of them all in terms of having the highest 𝐿 𝐵𝑜𝑙 / 𝐿 𝑋 with a lower limit of 15,000. We note that many of the other JWST-selected AGN may have similarly high bolometric to X-ray ratios, however in the case of GN-28074 we can set a more stringent lower limit because of its lower redshift and higher luminosity.  \nThe NIRSpec spectrum reveals additional insights into this source. In particular, the detection of absorption of H 𝛼 , H 𝛽 and He /i.pc 𝜆 10830, blueshifted by a few 100 km s -1 . Given that H 𝛼 and H 𝛽 are not resonant lines and that the 𝑛 = 2 level is very short lived, observing these two transitions in absorption requires very high densities in order to properly populate the 𝑛 = 2 level. More specifically, we have shown that, in order to properly reproduce the absorption seen in H 𝛼 , H 𝛽 and He /i.pc , the absorbing medium must have a density of at least 10 8 cm -3 , which is well above the values inferred for the ISM, while typical of the Broad Line Region (BLR) of AGN. The very small inferred sizes of the absorber (10 -6 -10 -4 pc) are also consistent with those expected for the BLR clouds (although the observed absorption is likely the superposition of multiple clouds). The inferred minimum column density of the absorbing gas is 10 21 . 5 -10 23 cm -2 , but the total column density of the multiphase gas is likely much larger. These findings support the scenario in which the extreme X-ray weakness is due (at least partly) to heavy absorption by Compton thick gas (dustpoor or dust-free), likely associated with the BLR clouds. The latter are known to be dust free and can therefore explain the widespread X-ray weakness of broad-lined, type 1 AGN discovered by JWST.  \nActually, although with a large uncertainty, we infer that the location of the absorbing medium is between the BLR and the dusty torus. Therefore, H 𝛼 , H 𝛽 and He /i.pc absorption may arise in dense neutral gas clouds in the outer part of the BLR and in the process of outflowing. However, we shall note that these outflow velocities are modest (a few hundred km s -1 ), much slower than the velocities seen in the nuclear region of quasars (e.g., in the case of BAL quasars), characterized by velocities of several thousands km s -1 in the optical lines and up to 10 4 -10 5 km s -1 in X-ray spectra.  \nWehave discussed that, because of observational difficulties in detecting H 𝛼 absorption, the fraction of broad-line AGN at high𝑧 characterized by this feature is likely much larger than the ∼ 10%currently estimated (which is already much larger than the fraction estimated for low𝑧 AGNobtained in previous surveys). This may indicate that the population of intermediate/low luminosity AGN discovered by JWST at high redshift is characterized by a larger covering factor of dense (and dust-poor/free) BLR clouds, than lower redshift counterpart. The associated (Compton-thick) absorption may explain, or contribute to their X-ray weakness.  \nAnother interesting feature of GN-28074 is its radio weakness. Once again, this seems to be a feature common to many other JWSTdiscovered AGN (Mazzolari et al., in prep.), but in GN-28074 it is most tightly constrained because of its high bolometric luminosity. The fact that it is radio weaker than radio quiet quasars (which are X-ray 'normal') suggest that the radio weakness and X-ray weakness  \nare somehow connected in this class of AGN. We have discussed that one possibility is free-free absorption. Indeed, the clouds responsible for H 𝛼 absorption are certainly optically thick in free-free absorption; in this case the absence of radio emission would be associated with the X-ray Compton thick absorption. The only concern is whether the extent of the distribution of absorbing clouds is large enough, and with a large enough covering factor, to absorb the radio emitting region (presumably a small jet) - however, direct observational evidence of jets absorbed by a free-free thick nuclear medium has been observed. The other possibility is that the nuclear magnetic field has not yet developed properly: this would both prevent the radio emission (absence of a radio jet) and would hamper the lifting of gas from the accretion disc to produce a hot corona, which is required for the hard X-ray emission.  \nWhile it has been suggested that the broad emission line components in these X-ray weak (and Radio weak) sources could originate from non-AGN activities, using our measurements for GN-28074, we rule out this possibility. Specifically, we explore alternative explanations for broad hydrogen lines including hyperdense or ultra metal-poor outflows and Raman scattering of the (stellar) UV continuum emission by a significant amount of neutral hydrogen. For the ultra metal-poor outflow scenario, the outflow needs to be a least an order of magnitude more metal-poor than the ISM, which is the opposite of what usually observed in galactic outflows. As for the hyperdense outflow scenario, the implied amount of dust attenuation would create significantly blueshifted profiles, in contrast to the symmetric broad profiles observed, unless the outflow is orders of magnitudes more metal-poor. In addition, the galactic outflow scenario would require an extreme SFR surface density above the theoretical Eddington limit and make GN-28074 very luminous in the submm, which is not detected in either SCUBA2 or NOEMA. The Raman scattering scenario, on the other hand, is also disfavored as the measured ratios of line widths of broad hydrogen lines are significantly different from the predictions. Detections of broad components in He/i.pc 𝜆 10830 and O /i.pc 𝜆 8446 further rule out the Raman scattering scenario. As a result, the AGN origin for the broad permitted lines is still favored over other scenarios.  \nWe finally note that local analogues of the AGN discovered by JWSTathigh𝑧 do exist. The fact that the, thoroughly studied, Seyfert 1 'prototype' NGC 4151 does have H 𝛼 , H 𝛽 and He /i.pc absorption, similar to GN-28074, and has substantial absorption in the X-rays, already indicates that GN-28074 (and likely many other JWSTdiscovered AGN) is not so peculiar. However, probably the closest and clearest example is SBS 0335-052E. This source has been traditionally considered to be an extremely metal-poor star-forming dwarf galaxy (Izotov et al. 1990). Only recently observations have revealed the clear presence of an AGN via the detection of a broad component of H 𝛼 (without a counterpart in the forbidden lines; Hatano et al. 2023), high ionization lines, and strong mid-IR variability (Hatano et al. 2023). While there are some differences (for instance, the SBS 0335-052E has a blue rest-optical continuum and a clear Balmer jump, whereas GN-28074 displays red optical continuum and a spectral break), there are also strong similarities. Like GN-28074 (and other JWST-discovered AGN), SBS 0335-052E is extremely X-ray weak (and the little X-ray emission detected is likely associated with the star-forming host) and also does not show evidence of nuclear radio emission at 1.5 GHz (Mazzolari et al, in prep.). Additionally, recent studies have shown that more local low-metallicity dwarf galaxies host similar AGN, which have been missed because of the lack of X-ray and radio emission.  \nIn summary, GN-28074 at intermediate redshift, and SBS 0335052E locally, are important 'Rosetta stones' to be investigated with  \nadditional observations and modelling to better understand the nature of the new, large population of AGN discovered by JWST.", 'ACKNOWLEDGEMENTS': "I.J. acknowledges support by the Huo Family Foundation through a P.C. Ho PhD Studentship. Ĳ, XJ, RM, FDE, JS, GM and JW acknowledge support by the Science and Technology Facilities Council (STFC), by the ERC through Advanced Grant 695671 'QUENCH', and by the UKRI Frontier Research grant RISEandFALL. RM also acknowledges funding from a research professorship from the Royal Society. SA acknowledges grant PID2021-127718NBI00 funded by the Spanish Ministry of Science and Innovation/State Agency of Research (MICIN/AEI/ 10.13039/501100011033). AJB acknowledges funding from the 'FirstGalaxies' Advanced Grant from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (Grant agreement No. 789056). ECL acknowledges support of an STFC Webb Fellowship (ST/W001438/1). M.P. acknowledges support from the research project PID2021-127718NB-I00 of the Spanish Ministry of Science and Innovation/State Agency of Research (MICIN/AEI/ 10.13039/501100011033), and the Programa Atracción de Talento de la Comunidad de Madrid via grant 2018T2/TIC-11715. PGP-G acknowledges support from grant PID2022139567NB-I00 funded by Spanish Ministerio de Ciencia e Innovación MCIN/AEI/10.13039/501100011033, FEDER, UE. BER acknowledges support from the NIRCam Science Team contract to the University of Arizona, NAS5-02015, and JWST Program 3215. ST acknowledges support by the Royal Society Research Grant G125142. HÜ gratefully acknowledges support by the Isaac Newton Trust and by the Kavli Foundation through a Newton-Kavli Junior Fellowship. JW acknowledges support from the Science and Technology Facilities Council (STFC), by the ERC through Advanced Grant 695671 'QUENCH', by the UKRI Frontier Research grant RISEandFALL. The work of CCW is supported by NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation.", 'DATA AVAILABILITY': 'All data used in this study are now public on the JWST DAWN archive.', 'APPENDIX A: CALCIUM TRIPLET, SODIUM ABSORPTION AND FLUORESCENT OXYGEN LINES': 'In addition to the clear absorption in the H 𝛼 , H 𝛽 and He /i.pc 𝜆 10830 lines 28074 also exhibits CaT and NaD absorption, shown in Figure A1.  \nAscanbeseeninthe figure, the sodium doublet absorption appears to coincide with the systematic redshift of the narrow lines, while the calcium triplet appears redshifted by about 300 km s -1 . As the triplet absorption is usually associated with stellar populations, a possible reason for the systematic redshift may be an underlying merger in the host galaxy. Following this interpretation, the reason for the lack of displacement of the sodium absorption may be resonant scattering in the underlying gas.  \nHowever, the kinematics of the calcium triplet region of the spectrum are difficult to precisely model due to underlying broad Ca /i.pc/i.pc emission as well as the nearby O /i.pc 𝜆 8446 line, while the sodium absorption is somewhat tentative (3-4 𝜎 significance). Therefore, we can not draw any strong conclusions about the origin of these features.  \nThe clearly detected O /i.pc 𝜆 8446 and O /i.pc 𝜆 11287 (Figure A2) lines have a photon flux ratio of around unity indicating that they are likely produced by Ly 𝛽 fluorescense, as is common in AGN (Rudy et al. 1989). However, constraints on O /i.pc 𝜆 1304 line are required to establish the Ly 𝛽 pumping scenario more robustly (Rodríguez-Ardila et al. 2002) while the O /i.pc 𝜆 1304 line is not covered by the wavelength range of the NIRSpec spectrum of GN-28074. Still, since there is no detection of O /i.pc 𝜆 13165 and O /i.pc 𝜆 7774 lines, the observed O/i.pc 𝜆 8446 is unlikely to have significant contributions from recombination, collisional excitation, and continuum fluorescence (Grandi 1980; Rodríguez-Ardila et al. 2002). The Ly 𝛽 pumping scenario, on other hand, is consistent with the presence of dense gas optically thick to Balmer lines. Recently, Choe et al. (2024) reported detection of O /i.pc 𝜆 8446 by JWST from an object in the gravitationally lensed Sunburst Arc at 𝑧 = 2 . 37. According to Choe et al. (2024), the O/i.pc 𝜆 8446 emitting object is likely an 𝜂 Carinae analog surrounded by dense gas, which produces O /i.pc 𝜆 8446 by Ly 𝛽 fluoresence. In the case of GN-28074 analysed in this paper, the dense gas that produces O/i.pc 𝜆 8446 is likely associated with regions surrounding the accreting black hole due to the presence of very broad components in permitted lines that are unlikely to be produced by stars.', 'APPENDIX B: GAS-PHASE METALLICITY': 'The clear presence of the [O /i.pc/i.pc/i.pc] 𝜆 4363 auroral line along with a tentative detection of the [O /i.pc/i.pc] 𝜆𝜆 7320,7331 doublet (Figure B1) allows for robust estimates of 𝑇 𝑒 in both O + and O ++ region. Utilizing the auroral to strong line ratios as well as an electron density of 𝑛 𝑒 = 170 + 232 -130 cm -3 , inferred from an observed [S /i.pc/i.pc] 𝜆 6716/[S /i.pc/i.pc] 𝜆 6731 ratio of 1 . 32 + 0 . 09 -0 . 13 , we estimate the 𝑇 𝑒 ( O ++ ) = 2 . 08 + 0 . 06 -0 . 07 × 10 4 K and 𝑇 𝑒 ( O + ) = 1 . 17 + 0 . 18 -0 . 15 × 10 4 K. Such a large difference in the temperatures of high-ionisation lines and low-ionisation lines are typical for the narrow line regions (NLRs) of AGN rather than starforming regions (e.g., Garnett 1992; Pérez-Montero 2014; Dors et al. 2020). With these temperatures, we derive a metallicity of 12 + log ( O / H ) = 7 . 85 + 0 . 17 -0 . 12 or about 15% Solar abundance when using the 𝐴 𝑉 derived in subsection 5.1.', 'APPENDIX C: LARGE SCALE IONIZED OUTFLOW': 'The fit of the [O /i.pc/i.pc/i.pc]doublet in the top left panel of Figure 2 clearly reveals the presence of weakly redshifted ionized outflow, seen also in the [S /i.pc/i.pc/i.pc] 𝜆 9531 line profile (Figure C1). A component with matching kinematics is required to fully explain the shape of the H 𝛼 and H 𝛽 lines, without which the absorption profiles yield inconsistent estimates for the column density of hydrogen in the n = 2 state. However, an outflow component is disfavoured by fitting in the [N /i.pc/i.pc] and [S /i.pc/i.pc] doublets as well as the He /i.pc 𝜆 10830 and Pa 𝛾 lines.  \nThe fit to the outflow component reveals a FWHM = 1310 + 19 -19 km s -1 and a slight, ∼ 130 km s -1 , redshift with respect to the narrow lines. We use the formalism of Carniani et al. (2015) to estimate the outflow mass traced by [O /i.pc/i.pc/i.pc] via:  \n𝑀 [ OIII ] = 1 . 7 × 10 3 𝑚 p 𝐶𝐿 [ OIII ] 10 [ O / H ]-[ O / H ]⊙ 𝑗 [ OIII ] 𝑛 e , (C1)  \nwhere 𝐶 = ⟨ 𝑛 e ⟩ 2 /⟨ 𝑛 2 e ⟩ , [ O / H ] ≡ log ( O / H ) , 𝑗 [ OIII ] is the emissivity of the [O /i.pc/i.pc/i.pc] doublet and 𝐿 [ OIII ] - the luminosity of the outflow component in [O /i.pc/i.pc/i.pc]. Our fits to the [S /i.pc/i.pc] doublet allow us to estimate 𝑛 e ∼ 100 cm -3 . As we do not detect the broad component in the [S /i.pc/i.pc] emission lines, we assume that the outflow has the same density as the ISM; this may not be entirely representative as AGN driven outflows tend to be somewhat denser than the host ISM (Davies et al. 2020; Santoro et al. 2020). Under these density assumptions, along with the previous 𝑇 𝑒 ( 𝑂 + +) = 2 × 10 4 K estimate to hold for the outflowing gas we estimate 𝑗 [ OIII ] = 1 . 2 × 10 -20 erg s -1 cm -3 using P/y.pcN/e.pc/b.pc (Luridiana et al. 2015). Using this emissivity along with the extinction corrected outflow luminosity yields 𝑀 [ OIII ] = 3 . 60 + 0 . 57 -0 . 46 × 10 7 M ⊙ . It should be noted that the mass given by [O/i.pc/i.pc/i.pc] emission should be treated as a lower limit as this emission generally does not trace the entire outflow (Carniani et al. 2015; Liu et al. 2013). We estimated outflow velocity to be 1240 + 18 -18 km s -1 by using the same formalism as before 𝑣 out = | Δ 𝑣 | + 2 𝜎 . Assuming a spherical geometry the mass outflow rate can then be estimated as (Maiolino et al. 2012):  \n/ 𝑀 𝑜𝑢𝑡 = 3 𝑣 𝑜𝑢𝑡 𝑀 𝑜𝑢𝑡 𝑅 𝑜𝑢𝑡 , (C2)  \nwhere 𝑣 𝑜𝑢𝑡 is the velocity, 𝑀 𝑜𝑢𝑡 - the mass and 𝑅 𝑜𝑢𝑡 - the radius of the outflow. Taking 𝑅 𝑜𝑢𝑡 = 1 kpc as an order of magnitude estimate based on the usual extent of AGN driven outflows (Harrison et al. 2014; Carniani et al. 2015), we obtain / 𝑀 [ 𝑂𝐼𝐼𝐼 ] = 137 + 21 -18 M ⊙ yr -1 .  \nIn addition to the [O /i.pc/i.pc/i.pc] doublet we also utilize the outflow emission detected in the H 𝛽 and H 𝛼 lines to estimate the mass and rate of outflowing hydrogen. We warn however, that these component are more uncertain than for the [O /i.pc/i.pc/i.pc] , as they are completely blended with the broad components of the same lines. Using the same derivation as in Carniani et al. (2015) we employ the following equation for the hydrogen outflow mass:  \n𝑀 = 0 . 8 𝑚 𝑝 𝐶𝐿 𝑖𝑘 𝑗 𝑖𝑘 𝑛 𝑒 , (C3)  \nwhere 𝐿 𝑖𝑘 and 𝑗 𝑖𝑘 are the luminosity and emissivity of a particular hydrogen transition respectively. Utilizing PyNeb as before we find 𝑗 𝐻𝛽 = 6 . 6 × 10 -26 erg s -1 cm -3 and 𝑗 𝐻𝛼 = 1 . 8 × 10 -25 erg s -1 cm -3 . Which lead to 𝑀 𝐻𝛽 = 13 . 0 + 2 . 7 -2 . 4 × 10 7 M ⊙ and 𝑀 𝐻𝛼 = 13 . 3 + 4 . 2 -3 . 2 × 10 7 M ⊙ - these estimates are consistent with each other and correspond to a total hydrogen outflow rate / 𝑀 𝐻 = 500 + 130 -110 M ⊙ yr -1 , with the [O /i.pc/i.pc/i.pc] derived outflow rate providing a lower limit. We note that the Paschen line emissivity in this regime is an order of magnitude below that of H 𝛽 thus the outflow component there is likely buried beneath the more luminous  \nFigure A1. Plots showcasing the positions of the calcium and sodium absorption features with respect to the OI 𝜆 8446 and He /i.pc 𝜆 5875 lines respectively. The observed spectrum is shown with the solid black line, errors are indicated by a shaded grey region.  \n<!-- image -->  \n˚ A]  \n<!-- image -->  \n˚  \n28074, z = 2.26  \nFigure A2. Region containing the O /i.pc 𝜆 11287 line, the position of which is indicated by a dashed green line. The R1000 grating data is shown by the black line, 1 𝜎 errors are indicated by a shaded green region. Sharp feature bluewards of the oxygen line is likely a data artefact as evidenced by the pipeline errors being considerably larger for that data.  \n<!-- image -->  \nBLR emission. On the other hand, the emissivity of He /i.pc 𝜆 10830 is estimated to be on the order of that of H 𝛼 , however, the outflow emission there can be strongly affected by resonant scattering thus its absence is not very constraining. All ionized outflow properties are summarized in Table C1.  \nIt is also possible to estimate the outflow kinetic power via the relation:  \n𝑃 𝑘 = 1 2 / 𝑀 𝑜𝑢𝑡 𝑣 2 𝑜𝑢𝑡 , (C4)  \nwhere / 𝑀 𝑜𝑢𝑡 and 𝑣 2 𝑜𝑢𝑡 are the outflow mass rate and velocity respectively. This equation yields 𝑃 𝑘 = 2 . 42 + 0 . 64 -0 . 52 × 10 44 erg s -1 when  \nTable C1. A summary of the observed ionized outflow properties. Column one gives the name of the line, second column - the 𝐴 𝑉 corrected outflow luminosity in units of 10 42 erg s -1 , third column - the outflowing mass in units of 10 7 M ⊙ , the fourth - outflow velocity in km s -1 . Column five shows the mass outflow rates in M ⊙ yr -1 , while the final column gives the outflow kinetic power in units of 10 43 erg s -1 .  \nusing the H 𝛼 emission and 𝑃 𝑘 = 6 . 85 + 1 . 14 -0 . 97 × 10 43 erg s -1 when using the [O /i.pc/i.pc/i.pc] doublet, corresponding to around 1-5% of total AGN luminosity. This ratio is about 1-2 dex above the expected value for quasars at similar redshifts, as illustrated in Figure C2, and is consistent with what found for molecular outflows. However, it should be noted that the [O /i.pc/i.pc/i.pc] outflow rate was calculated assuming the NLR metallicity for the outflowing gas, while outflows tend to be more metal rich than their hosts (Chisholm et al. 2018). Assuming Solar metallicity brings the outflow rate down by about 1 dex and results in a lower limit of 𝑃 𝑘 = 6 . 85 + 1 . 14 -0 . 97 × 10 42 erg s -1 . In addition, the comparisons in Figure C2 are somewhat complicated by differing outflow velocity estimations used across the literature. In the case of data from Harrison et al. (2014), the outflow velocities were taken to be FWHM/1.3 for single Gaussian profiles, which would lower our estimates on outflow power by about 50%. To reflect these systematic differences, 0.3 dex was added to the error budget on the outflow kinetic power.  \nThis paper has been typeset from a T E X/L A T E X file prepared by the author.  \n<!-- image -->  \nFigure B1. Fits to the [O /i.pc/i.pc] 𝜆𝜆 3726,3729 doublet and its auroral counterpart. The [O /i.pc/i.pc] 𝜆𝜆 7320,7331 doublet appears marginal in R1000, despite having a formal significance of ∼ 6 𝜎 .  \n<!-- image -->  \nFigure C1. Best fit model of the [S /i.pc/i.pc/i.pc] 𝜆 9532 line emission. The line profile exhibits clear broad wings that are fitted well by an outflow component.  \n<!-- image -->  \nFigure C2. A comparison between the AGN bolometric luminosity, on the x axis, and outflow kinetic power, on the y axis, for 28074 and QSOs from Carniani et al. (2015), Harrison et al. (2014) and Fluetsch et al. (2019). Molecular outflow data is marked with green squares, while ionized outflows - with circles. Our object is indicated by magenta and violet stars for estimates using H 𝛼 and [O /i.pc/i.pc/i.pc] emission at NLR metallicity respectively. The lower limit obtained from [O /i.pc/i.pc/i.pc] emission assuming solar metallicity is marked by the purple star. The error bars indicate ∼ 0.3 dex systematic uncertainties. The literature data has been adjusted for the factor of 3 discrepancy between assumed geometries. The lines denote constant ratios of 100%, 5% and 0.1%.  \n<!-- image -->'}