Here we present high spectral resolution inlineformulatexmath idTM0001 notationLaTeXit JWSTtexmathinlineformula NIRSpec observations of GN 42437 a lowmass logMinlineformulatexmath idTM0002 notationLaTeXast rm Modot7.9texmathinlineformula compact inlineformulatexmath idTM0003 notationLaTeXre lt 500texmathinlineformulapc extreme starburst galaxy at inlineformulatexmath idTM0004 notationLaTeXz5.59texmathinlineformula with 13 emissionline detections. GN 42437 has a low metallicity 510 per cent Zinlineformulatexmath idTM0005 notationLaTeXodottexmathinlineformula and its restframe Hinlineformulatexmath idTM0006 notationLaTeXalphatexmathinlineformula equivalent width suggests nearly all of the observed stellar mass formed within the last 3 Myr. GN 42437 has an extraordinary 7inlineformulatexmath idTM0007 notationLaTeXsigmatexmathinlineformula significant Ne V 3427  detection. The Ne V line has a restframe equivalent width of inlineformulatexmath idTM0008 notationLaTeX11pm 2texmathinlineformula Ne VHinlineformulatexmath idTM0009 notationLaTeXalpha 0.04pm 0.007texmathinlineformula Ne VNe III 3870 inlineformulatexmath idTM0010 notationLaTeX 0.26pm 0.04texmathinlineformula and Ne VHe II 4687 inlineformulatexmath idTM0011 notationLaTeX  1.2pm 0.5texmathinlineformula. Ionization from massive stars shocks or highmass Xray binaries cannot simultaneously produce these Ne V and other lowionization line ratios. Reproducing the complete nebular structure requires both massive stars and accretion on to a black hole. We do not detect broad lines nor do the traditional diagnostics indicate that GN 42437 has an accreting black hole. Thus the very highionization emission lines powerfully diagnose faint narrowline black holes at high redshift. We approximate the black hole mass in a variety of ways as logMinlineformulatexmath idTM0012 notationLaTeXrm BHrm Modot sim 57texmathinlineformula. This black hole mass is consistent with local relations between the black hole mass and the observed velocity dispersion but significantly more massive than the stellar mass would predict. Very highionization emission lines may reveal samples to probe the formation and growth of the first black holes in the universe.
['Astrophysics - Astrophysics of Galaxies']
Ne v emission from a faint epoch of reionizationera galaxy evidence for a narrowline intermediatemass black hole
{'[Ne v] emission from a faint epoch of reionization-era galaxy: evidence for a narrow-line intermediate mass black hole': '- J. Chisholm 1 ★ , D. A. Berg 1 , R. Endsley 1 , S. Gazagnes 1 , C. T. Richardson 2 , E. Lambrides 3 , 4 , J. Greene 5 ,\n- S. Finkelstein 1 , S. Flury 6 , 7 , N. G. Guseva 8 , A. Henry 9 , T. A. Hutchison 3 , 4 , Y. I. Izotov 8 ,', 'R. Marques-Chaves 10 , P. Oesch 10 , C. Papovich 11 , A. Saldana-Lopez 12 , D. Schaerer 10 , 13 , M. G. Stephenson 1': '1 Department of Astronomy, University of Texas at Austin, 2515 Speedway, Austin, Texas 78712, USA  \n- 2 Elon University, 100 Campus Drive, Elon, NC 27244  \n3 Astrophysics Science Division, Code 662, NASA Goddard Space Flight Center, 8800 Greenbelt Rd, Greenbelt, MD 20771, USA  \n4 NASA Postdoctoral Fellow  \n- 5 Department of Astrophysical Sciences, Princeton University, 4 Ivy Lane, Princeton, NJ08544, USA\n- 6 Department of Astronomy, University of Massachusetts Amherst, Amherst, MA 01002, United States\n- 7 NASA FINESST Fellow\n- 8 Bogolyubov Institute for Theoretical Physics, National Academy of Sciences of Ukraine, 14-b Metrolohichna str., Kyiv, 03143, Ukraine\n- 9 Space Telescope Science Institute, 3700 San Martin Drive Baltimore, MD 21218, United States\n- 10 Geneva Observatory, Department of Astronomy, University of Geneva, Chemin Pegasi 51, CH-1290 Versoix, Switzerland\n- 11 Department of Physics and Astronomy, Texas A&M University, College Station, TX, 77843-4242 USA\n- 12 Department of Astronomy, Stockholm University, Oscar Klein Centre, AlbaNova University Centre, 106 91 Stockholm, Sweden\n- 13 CNRS, IRAP, 14 Avenue E. Belin, 31400 Toulouse, France  \nAccepted XXX. Received YYY; in original form ZZZ', 'ABSTRACT': 'Here we present high spectral resolution JWST NIRSpec observations of GN 42437, a low-mass (log(M ∗ / 𝑀 ⊙ ) = 7 . 9), compact ( 𝑟 𝑒 < 500pc), extreme starburst galaxy at 𝑧 = 5 . 59 with 13 emission line detections. GN 42437 has a low-metallicity (5-10% Z ⊙ ) and its rest-frame H 𝛼 equivalent width suggests nearly all of the observed stellar mass formed within the last 3 Myr. GN 42437 has an extraordinary 7 𝜎 significant [Ne /v.pc] 3427 Å detection. The [Ne /v.pc] line has a rest-frame equivalent width of 11 ± 2Å, [Ne /v.pc]/H 𝛼 = 0 . 04 ± 0 . 007, [Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] 3870Å = 0 . 26 ± 0 . 04, and [Ne /v.pc]/He /i.pc/i.pc 4687Å = 1 . 2 ± 0 . 5. Ionization from massive stars, shocks, or high-mass X-ray binaries cannot simultaneously produce these [Ne /v.pc] and other low-ionization line ratios. Reproducing the complete nebular structure requires both massive stars and accretion onto a black hole. We do not detect broad lines nor do the traditional diagnostics indicate that GN42437 has an accreting black hole. Thus, the very-high-ionization emission lines powerfully diagnose faint narrow-line black holes at high-redshift. We approximate the black hole mass in a variety of ways as log(M BH / 𝑀 ⊙ ) ∼ 5 -7. This black hole mass is consistent with local relations between the black hole mass and the observed velocity dispersion, but significantly more massive than the stellar mass would predict. Very-high-ionization emission lines may reveal samples to probe the formation and growth of the first black holes in the universe.  \nKey words: galaxies: high-redshift - galaxies: evolution - galaxies: formation', '1 INTRODUCTION': "JWST observations have discovered a rich and ubiquitous population of super massive black holes (SMBH) in fainter galaxies at higher redshifts than previously known (Brinchmann 2023; Furtak et al. 2023; Goulding et al. 2023; Greene et al. 2023; Harikane et al. 2023; Kocevski et al. 2023; Kokorev et al. 2023; Labbé et al. 2023; Larson et al. 2023; Maiolino et al. 2023; Matthee et al. 2023; Onoue et al. 2023; Scholtz et al. 2023; Übler et al. 2023; Kokorev et al. 2024; Lambrides et al. 2024). These recently detected black holes have masses that range up to ∼ 10 8 M ⊙ only a few 100 Myr after the Big Bang. Energy and momentum from these intrinsically luminous  \nsources may impact the formation and evolution of their host galaxies (Gebhardt et al. 2000; Ferrarese & Merritt 2000). Fully accounting for the energy and momentum that accretion onto early black holes injects into the gas within galaxies requires revealing the population demographics of black holes in the early universe.  \nCombined with observations of UV-bright quasars (Becker et al. 2001; Fan et al. 2006; Willott et al. 2010; Jiang et al. 2016; Bañados et al. 2018; Matsuoka et al. 2018, 2019; Eilers et al. 2023; Yang et al. 2023), these observations have firmly established the existence of SMBHs at 𝑧 > 5, but theory still debates how these SMBHs formed. The age of the universe is just 1 Gyr at 𝑧 = 5 . 5. If black holes start as stellar mass black holes ( ∼ 10 -100 M ⊙ ), Eddingtonlimited accretion just barely grows the black hole to 10 8 M ⊙ in the age of the universe (Haiman & Loeb 2001; Bromm & Larson 2004).  \nThis leaves little time for the universe to expand and sufficiently cool to form the first black holes without super-Eddington accretion prescriptions (Madau et al. 2014; Schneider et al. 2023).  \nAlternatively, the first black holes in the early universe may have formed as intermediate mass black holes (IMBHs; 10 3 -6 M ⊙ ; Greene & Ho 2004; Mezcua 2017; Greene et al. 2020). These larger seeds start closer to the observed mass of SMBHs, partially relieving the pressure to accrete so rapidly to become a SMBH. Heroic effort has attempted to uncover IMBHs in the local universe, revealing limited robust detections (e.g. Gebhardt et al. 2005; Izotov et al. 2007; Greene et al. 2008; Noyola et al. 2010; Lützgendorf et al. 2013; Moran et al. 2014; Reines & Volonteri 2015; Nguyen et al. 2017; Lin et al. 2018; Nguyen et al. 2019; Pechetti et al. 2024). As such, observations still have not provided the crucial link between stellar mass black holes and SMBHs.  \nIn the local universe, SMBHs are distinguishable from starforming galaxies using the ratio of strong optical emission lines (Baldwin et al. 1981; Veilleux & Osterbrock 1987; Kewley et al. 2001; Kauffmann et al. 2003) or broad ( > 1000 km s -1 ) Balmer emission lines (Osterbrock et al. 1976). Galaxies with broad Balmer lines are often called broad-line Active Galactic Nuclei (AGN), or Type /i.pc AGN. The broad-lines provide direct kinematic evidence for the sphere of influence of the black hole and enable robust black hole mass estimates (M BH ; Vestergaard & Peterson 2006; Greene & Ho 2007; Shen & Kelly 2012; Padovani et al. 2017). The other type of AGN, Type /i.pc/i.pc, are galaxies where ratios of their optical emission lines indicate the presence of a black hole, but the emission lines are narrow ( < 1000 km s -1 ). Emission line ratios have uncovered significant populations of Type /i.pc/i.pc AGN in the local universe (Heckman et al. 2004; Kauffmann & Heckman 2009; Liu et al. 2009; Heckman &Best 2014; Moran et al. 2014; Hickox & Alexander 2018).  \nLow-redshift methods to distinguish AGN have challenges at highredshift, where galaxies have lower masses and metallicities (Groves et al. 2004; Kewley et al. 2013; Dors et al. 2024). Most high-redshift SMBHs have been classified based on the presence of broad Balmer emission lines (Greene et al. 2023; Harikane et al. 2023; Kocevski et al. 2023; Larson et al. 2023; Maiolino et al. 2023; Matthee et al. 2023). However, high-redshift black holes are likely significantly lower mass. The expected broad-line Gaussian velocity widths ( 𝜎 ) reduce from > 1000 km s -1 for SMBHs to 20-250 km s -1 for IMBHs. These narrow line widths are similar to H /i.pc/i.pc regions within low-mass star-forming galaxies, making IMBHs challenging to spectroscopically tease out. Extremely high spectral resolution observations of these intrinsically faint sources must detect the relatively narrow emission lines. As cosmological dimming makes more distant galaxies increasingly fainter, determining the kinematic properties of IMBHs becomes increasingly challenging. While most of the high-redshift SMBHs have been discovered using broad-lines, only about one-quarter to a half of the low-redshift AGN are Type /i.pc (Ho et al. 1997; Ho 2008). To fully characterize the high-redshift AGN population a complete census of Type /i.pc/i.pc AGN are required.  \nUnfortunately, metal-poor stars have extremely hard ionizing spectra (Leitherer et al. 1999; Schaerer 2002; Eldridge et al. 2017; Chisholm et al. 2019; Telford et al. 2023), which can mimic the ionizing spectral shapes of accretion disks onto black holes. Lowmetallicity massive stars and IMBHs have similar ionizing slopes and shapes up to ≈ 54 eV (Feltre et al. 2016). These harder spectra ensure that low-metallicity star-forming galaxies have similar [O /i.pc/i.pc/i.pc]/H 𝛽 ratios as high-metallicity AGN. Additionally, the hard-ionizing spectra of low-metallicity stars can fully ionize low-ionization gas ([N /i.pc/i.pc] 6585 Å, [S /i.pc/i.pc] 6718,6733 Å, [O /i.pc] 6300 Å). Low-metallicity AGN and star-forming galaxies reside in nearly the same locations  \nin the traditional diagnostic diagrams (Shirazi & Brinchmann 2012; Moran et al. 2014; Reines & Volonteri 2015). A different diagnostic is required to reveal narrow-line AGN in low-metallicity systems.  \nStellar ionizing spectra do not extend much below the second helium ionization edge (54 eV), except at zero metallicity (Tumlinson & Shull 2000; Schaerer 2002). Metal poor stars cannot produce very-high-ionization gas phases, which are defined as having ionization potentials greater than 54 eV (Leitherer et al. 1999; Izotov et al. 2004; Thuan & Izotov 2005; Eldridge et al. 2017; Stanway & Eldridge 2019; Berg et al. 2021; Olivier et al. 2022). Meanwhile, the ionizing spectrum of accretion disks around black holes extends to these ionization energies (Cann et al. 2018, 2019; Satyapal et al. 2021; Richardson et al. 2022; Hatano et al. 2023), and very-highionization emission lines are commonly observed in the spectra of AGNatlow-redshift (Osterbrock 1989; Zakamska et al. 2003; Dasyra et al. 2008; Gilli et al. 2010). Therefore, very-high-ionization emission lines promise to robustly separate low-metallicity star-forming galaxies and IMBHs. Powerful diagnostics using C /i.pc/v.pc 1550 Å (with an ionization energy of 64 eV), He /i.pc/i.pc 1640 Å (with an ionization energy of 54 eV), [Ne /i.pc/v.pc] 2423 Å (with an ionization energy of 97 eV), and N /i.pc/v.pc 1720 (with an ionization energy of 77 eV) have been proposed to separate low-metallicity star-forming galaxies from AGN (Shirazi & Brinchmann 2012; Feltre et al. 2016; Berg et al. 2018; Nakajima et al. 2018; Senchyna et al. 2020; Schaerer et al. 2022; Olivier et al. 2022).  \nHere, we present JWST NIRSpec high-spectral resolution observations of a low-mass, extreme starburst galaxy, hereafter called GN 42437, at 𝑧 = 5 . 59 that has a strong [Ne /v.pc] 3427 Å detection. Ionizing gas from the [Ne /i.pc/v.pc] to the [Ne /v.pc] state requires photons with energy above 97 eV, energies that even the most massive stars cannot produce. In section 2 and section 3 we introduce the observations and data reduction. We use these observations in section 4 and section 5 to determine the galaxy properties and emission line characteristics. We discuss the very-high-ionization emission lines (section 6), the full nebular structure (subsection 7.1), and how a black hole is the most likely source of the observed very-high-ionization emission lines (subsection 7.2). In section 8 we discuss the properties of the black hole in GN 42437. IMBHs during the early stages of the universe are discussed in section 9 and we conclude in section 10.  \nThroughout this paper, we use a cosmology with H 0 = 67 . 4, and Ω M = 0 . 315 (Planck Collaboration et al. 2020), such that 0.1' corresponds to 600 pc at 𝑧 = 5 . 59. All magnitudes are AB magnitudes. Solar metallicity is defined as 12+log(O/H ) = 8 . 69 (Asplund et al. 2009). All rest-frame wavelengths of emission lines are given as Angstroms (Å) in the vacuum-frame using the National Institute of Standards and Technology (NIST) database (Kramida et al. 2018).", '2 OBSERVATIONS': 'We obtained JWST (McElwain et al. 2023; Rigby et al. 2023) NIRSpec (Böker et al. 2023) G235H/F170LP and G395H/F290LP observations of 20 galaxies in the GOODS-North (GN) field ( JWST Project ID: 1871, PI: Chisholm) using the Multi-Slit Array (MSA; Ferruit et al. 2022) on February 10, 2023. The main goal of this project is to measure both the production and escape of ionizing photons from a sample of galaxies during the epoch of reionization. We selected the high-resolution grating configuration to obtain velocity-resolved profiles of the Mg /i.pc/i.pc emission lines. Mg /i.pc/i.pc profiles offer crucial clues for the neutral gas column density (Henry et al. 2018; Chisholm et al. 2020; Xu et al. 2023). Detecting faint Mg /i.pc/i.pc emission requires deep observations. Additionally, weak Balmer emission lines mea-  \nTable 1. Galaxy properties of GN 42437 derived from the HST+ JWST photometry and spectra. The right ascension and declination are the coordinates of GN 42437 in the MSA configuration file. The spectroscopic redshift, 𝑧 spec , was measured as the median redshift of the H 𝛼 , [O /i.pc/i.pc/i.pc] 5008 Å, [O /i.pc/i.pc/i.pc] 4960 Å, and [Ne /i.pc/i.pc/i.pc] 3970 Å emission lines. The stellar mass (log( 𝑀 ∗ )) was calculated from the /b.pc/a.pc/g.pc/p.pc/i.pc/p.pc/e.pc/s.pc SED model fitting using both the imaging and spectra (section 4) and the star formation rate (SFR) was determined using the H 𝛼 emission line and the Kennicutt & Evans (2012) prescription. The dynamical mass (M dyn ) was calculated using the intrinsic H 𝛼 line width and the size of the galaxy in the F444W image that traces the H 𝛼 (see section 8). The absolute UV magnitude ( 𝑀 UV ) was measured from the HST imaging (Finkelstein et al. 2015).  \nure the dust attenuation and nebular properties (Stephenson et al. in preparation) of the individual galaxies. The Mg /i.pc/i.pc results will be featured in Gazagnes et al. (in preparation), while here we focus on the [Ne /v.pc] 3427 Å emission from a galaxy in the MSA footprint, GN 42437.  \nWe centered the MSA footprint on a well-known, bright, 𝑧 = 7 . 5 Ly 𝛼 emitter in the GN field (Finkelstein et al. 2013; Hutchison et al. 2019; Jung et al. 2020). We used the deep Hubble Space Telescope ( HST ) H-band observations to provide possible targets with spectroscopically confirmed (Jung et al. 2020) or photometrically suggested redshifts greater than 5. We created a catalog of 1,036 possible sources from Finkelstein et al. (2015) and Bouwens et al. (2015) and used the 100 mas pixel-scale Complete Hubble Archive for Galaxy Evolution (CHArGE; Fujimoto et al. 2022) 12 HST F160W imaging to define the source locations. CHArGE aligns all of the HST images onto Gaia DR2 reference frames and, preJWST imaging, provided the highest spatial resolution and most accurate source locations. From the catalog of 1,036 candidates, we selected 20 sources within one MSA footprint of the bright Ly 𝛼 emitter. We required these targets to have F160W magnitudes less than 28 mag and > 3 𝜎 Spitzer 4.5 𝜇 mdetections. This ensured that the targeted galaxies resembled star-forming galaxies within the epoch of reionization that were sufficiently bright to detect Mg /i.pc/i.pc and other nebular emission lines. The full galaxy sample will be discussed in an upcoming paper (SaldanaLopez et al. in preparation). Here we focus on GN 42437, which was designated as source "63" in the MSA planning tool because it was the 63rd source drawn from the Finkelstein et al. (2015) catalog (see Table 1 for the Right Ascension, Declination, and properties of GN 42437).  \nImportantly, GN 42437 was selected as one of the faintest galaxies in the sample (F160W magnitude of 27.6 mag; Table 1) with a photometric redshift of 5.813 (Finkelstein et al. 2015). This galaxy was not selected specifically based on any of its rest-frame FUV or optical properties. We selected it only because it was within one NIRSpec MSA footprint of the main target, its photometric redshift  \nsuggested it was in the tail-end of the epoch of reionization, and the HST and Spitzer observations suggested it was sufficiently bright for NIRSpec to detect rest-frame optical emission lines.  \nWe used the Skelton et al. (2014) catalog to define five reference stars for MSA alignment. These five sources were distributed in three of the MSA quadrants. We estimated the expected JWST colors of these reference stars using the HST photometry (Documentation 2016). We included two of these reference stars on the MSA slit design to aid in flux and wavelength calibration. For the target acquisition, we used the MSATA method with the F140X filter and the NRSRAPIDD6 readout pattern. The acquisition used three groups per integration and four total integrations for a total exposure time of 687 s. We then took a final confirmation image with one exposure and 30 groups per integration with the NRSIRS2RAPID readout pattern for a total time of 452 s.  \nThe science observations were split between the G235H and G395H grating configurations. The G235H observations typically cover rest-frame wavelengths between 2550-3500 Å for galaxies within the epoch of reionization (5 . 5 < 𝑧 < 9 . 5). This wavelength range contains the intrinsically weak Mg /i.pc/i.pc emission lines and the [Ne /v.pc] 3427 Å line. Thus, we prioritized the G235H grating with longer total exposure times. For the G235H observations, we exposed using the NRSIRS2 readout mode for 20 groups per integration with 6 integrations per exposure and 36 total exposures for a total integration time of 53,044 s, or 14.7 hours. The G395H observations cover rest-frame optical emission lines that are brighter than Mg /i.pc/i.pc and require significantly shorter integration times. For the G395H grating we observed using the NRSIRS2 readout mode with 22 groups per integration and two integrations per exposure with 6 total integrations for a total exposure time of 9716 s, or 2.7 hours. We observed with the standard three shutter nod pattern for background subtraction. These three NIRSpec MSA shutters are overlaid on the NIRCam image of GN 42437 in Figure 1. The background shutters (upper and lower rectangles) do not contain any discernible galaxy emission. We intentionally designed the MSA configurations to include all of the crucial rest-frame optical and NUV lines of GN 42437, including H 𝛽 , H 𝛾 , and [O /i.pc/i.pc/i.pc] 4363 Å. However, the slight difference between the photometric and spectroscopic redshift ( 𝑧 phot = 5 . 81 versus 𝑧 spec = 5 . 58) meant that some of these lines fell in the detector gaps or between the two gratings.  \nNIRCam(Riekeetal.2023)imagingoftheGNfieldwastakenwith the First Reionization Epoch Spectroscopically Complete Observations (FRESCO; Oesch et al. 2023). These observations included medium-band observations with the F182M and F210M bands in the short-wavelength channel, and F444W observations in the longwavelength channel. The FRESCO imaging was designed to reach the 5 𝜎 significance for a 28.2 mag source in each filter with 4456, 3522, and 934 s of total exposure time for the F182M, F210M, and F444W filters, respectively. GN 42437 is significantly detected in all FRESCO images (see subsection 3.1). The GN field has a rich data set in multi-wavelength observations, including deep HST observations. We also used the Chandra Deep Field-North Survey (Alexander et al. 2005) to provide X-ray information on GN 42437. We use these NIRCam and Chandra images to verify the NIRSpec data reduction and to infer the properties of GN 42437.', '3 REDUCTION': 'Herewedescribethereductionofthe JWST NIRCam(subsection3.1) and JWST NIRSpec (subsection 3.2) data.  \nTable 2. HST /ACS, HST /WFC3, NIRCam, and NIRSpec photometry of GN 42437. The NIRSpec flux densities are estimated from the NIRSpec data in synthetic bandpasses corresponding to the NIRCam filters. The filters above the horizontal line are either HST /ACS, HST /WFC3, or have both NIRCam(Oeschetal. 2023) and NIRSpec observations. HST /ACSor HST /WFC3 filters are labeled with an ACS or WFC3 before their filter name. The flux densities listed below the last horizontal line only have NIRSpec data and are used to determine equivalent widths of individual emission lines (see section 4). The last three columns provide the fitted reference size (r ref ), Sérsic index ( 𝑛 ), and ellipticity ( 𝑞 ) fitted from the NIRCam data. /m.pc/s.pc/a.pc/f.pc/i.pc/t.pc (de Graaff et al. 2023) requires these last three columns to estimate the NIRSpec spectral resolution. The dynamical mass estimate in section 8 also uses the NIRCam profile fit.', '3.1 NIRCam Data Reduction': 'We reduce the FRESCO NIRCam imaging following the methods outlined in Endsley et al. (2023), utilizing both custom scripts as well as the JWST Science Calibration Pipeline 3 (v1.11.3). We use custom snowball and wisp subtraction scripts, that incorporate sky flats and wisp templates constructed from public data, and adopt the photometric zeropoints from Boyer et al. (2022) as part of the jwst\\_1106.pmap reference file. We subtract 1 / 𝑓 noise and the 2D background from the *\\_cal.fits files on an amplifier-by-amplifier basis using the /s.pc/e.pc/p.pc package (Barbary 2016). Each *\\_cal.fits file is individually aligned to the Gaia astrometric frame using the CHArGE HST WFC3/F160Wreductions of the GN field. All NIRCam mosaics are resampled onto the same World Coordinate System with a scale of 30 mas pixel -1 during the final stage of the pipeline. The individual images are convolved to the point-spread function (PSF) of the F444W filter (the longest wavelength filter) using empirically defined PSFs from within the FRESCO mosaics. Figure 1 shows the FRESCO three-color PSF-convolved NIRCam image of GN 42437 with the 3-shutter NIRSpec MSA pattern overlaid. GN 42437 has an obvious red color from a F444W excess.', '3.2 NIRSpec Data Reduction': 'As in-flight reference files have been added to the NIRSpec reduction pipeline, significant improvements have been made to the data quality. We used the reference files designated as jwst\\_1174.pmap , which were released on the CDRS website 4 on December 14, 2023.  \n0  \nFigure 1. Three color NIRCam image (blue is F182M, green is F210M, and red is F444W) of GN 42437 with north up and east to the left. GN 42437 has a distinct red color in these filters due to strong H 𝛼 emission. The location of the three NIRSpec shutters are overlaid in orange. GN 42437 is fully within the slit and shows a compact morphology. There exists a secondary component to the North East of GN 42437 that is not within the slit. Photometric fitting does not definitively distinguish whether the nearby galaxy is a true companion of GN 42437 or a foreground interloper.  \n<!-- image -->  \nThese updated reference files include in-orbit dark, mask, and bias files that were taken between July 2022 and December 2022. The improvement from the reference files has dramatically changed the delivered data quality since the observations were initial taken. These latest reference files significantly impact both the relative and absolute fluxing of the NIRSpec data as well as the outlier rejection. Previous reductions using older reference files produced significantly non-physical Balmer decrement values and contained a large number of flux spikes that compromised the interpretation of the faint emission lines presented here.  \nWe reduced the NIRSpec data using /m.pc/s.pc/a.pc/e.pc/x.pc/p.pc v06.17 (Brammer 2023). /m.pc/s.pc/a.pc/e.pc/x.pc/p.pc is a publicly available 5 python code that acts as a wrapper for the default Space Telescope Science Institute data reduction pipeline, while providing a few additional steps. In particular, /m.pc/s.pc/a.pc/e.pc/x.pc/p.pc uses the /g.pc/r.pc/i.pc/z.pc/l.pc/i.pc code (Brammer 2019) to drizzle and combine the individual exposures. We further augment the /m.pc/s.pc/a.pc/e.pc/x.pc/p.pc pipeline using the residual 1/f-noise correction NSC/l.pc/e.pc/a.pc/n.pc which was recently released by Space Telescope Science Institute (Rauscher 2023). We used the Space Telescope Science Institute default pathloss correction algorithm, assuming a point source. The NIRCam imaging of GN 42437 suggests that it is near the diffraction limit of the NIRSpec system (Table 2). Figure 1 also suggests that the object is well-captured by the MSA slit but slightly off-center. Comparisons between the NIRCam and NIRSpec fluxes (see below) suggest that this NIRSpec reduction process accurately reproduces the NIRCam fluxes of GN 42437.  \nFor each of the 20 sources that we opened MSA slits for, we manually set the initial location, size, and shape of the spectral extraction trace in /m.pc/s.pc/a.pc/e.pc/x.pc/p.pc based on a by-eye assessment of the source location and then fit for the optimal extraction from the rectified spectra. The final spectrum for GN 42437 was extracted from a location 3.5 pixels below the center of the slit using an optimal extraction with an initial guess at the width of 1 pixel for both grating configurations. We visually inspected all of the individual exposures for contamination and shutter malfunctions, which we did not find, and included all exposures of GN 42437 in the final co-addition of the data.  \nThe final data product of /m.pc/s.pc/a.pc/e.pc/x.pc/p.pc is a two-dimensional and onedimensional extracted flux and error spectra of each configuration. We compared the continuum standard deviation to running medians of the √︁ < Error 2 > estimated from the /m.pc/s.pc/a.pc/e.pc/x.pc/p.pc pipeline. We found that the 1 𝜎 noise estimates accurately reproduced the variability in the continuum. Therefore, we did not scale the uncertainties further.  \nPrevious NIRSpec programs have normalized the NIRSpec data to the observed NIRCam observations to finalize the NIRSpec absolute flux calibration. In the past this has been done to correct for uncertain flux calibration or for aperture losses. The FRESCO NIRCam imaging provides three separate wavelength regimes to compare both the relative and absolute fluxing of the NIRSpec data. We computed synthetic magnitudes using /s.pc/y.pc/n.pc/p.pc/h.pc/o.pc/t.pc (STScI Development Team 2018, 2020) in the respective FRESCO filters. The three upper panels of Table 2 tabulates the synthetic NIRSpec flux values in the second column and compares them to the NIRCam aperture photometry. We found excellent (within 2 𝜎 significance) agreement at all wavelengths for the NIRSpec and NIRCam flux densities (see the comparison of the gold and blue points in Figure 2). The F182M synthetic NIRSpec photometry has a larger uncertainty because the NIRSpec F182M synthetic magnitude contains the NRS detector gap. This reduces the total number of NIRSpec data points. The NIRSpec and NIRCam observations have consistent fluxing both within individual gratings (a comparison of the F182M and F210M flux densities tests the G235H relative fluxing) and between gratings (comparing the F210M and F444W flux densities tests the G235H and the G395H). Thus, we do not observe the tendency for NIRSpec grating observations to have systematically higher flux density values than NIRCam observations (e.g., Bunker et al. 2023; Cameron et al. 2023). When we use older reference files, we did measure significant flux density discrepancies between NIRSpec and NIRCam, suggesting that the updated reference files produce more reliable NIRSpec absolute and relative flux calibration.', '4 GALAXY PROPERTIES': "Using the NIRCam and NIRSpec data, we constrain the properties of GN 42437. We detect both the continuum (Table 2) and emission lines (Table 3) in the NIRSpec observations. This is emphasized in Figure 2 where some of the gold points include emission lines from the underlying SED fit in gray, but some of the NIRCam and NIRSpec points only probe the galaxy continuum. This disentangles the contribution of the high equivalent width emission lines to the integrated photometry. In Figure 2 we highlight that the observations probe light as red as 4.98 𝜇 m (rest-frame 7500 Å), suitable to measure stellar masses (M ∗ ) and star formation histories (SFHs). We infer M ∗ and the SFH of GN 42437 by fitting both the NIRCam imaging and the NIRSpec spectra with the Bayesian Analysis of Galaxies for Physical Inference and Parameter EStimation (/b.pc/a.pc/g.pc/p.pc/i.pc/p.pc/e.pc/s.pc; Carnall et al. 2018) code. /b.pc/a.pc/g.pc/p.pc/i.pc/p.pc/e.pc/s.pc uses the updated Bruzual & Charlot (2003) stellar population synthesis templates, the Kroupa & Boily (2002) stellar  \nFigure 2. The observed spectral energy distribution (SED) of GN 42437. The NIRCam and HST imaging are shown as dark-blue squares and the gold circles are the NIRSpec synthetic photometry. All data are given in Table 3. The /b.pc/a.pc/g.pc/p.pc/i.pc/p.pc/e.pc/s.pc fit to the photometry and the spectra is given as the gray model. Using both the photometry and the spectra, we more finely sample the SED and probe emission-line-free regions that estimate the star formation history and stellar mass of GN 42437.  \n<!-- image -->  \ninitial mass function, and includes nebular continuum and emission by processing the stellar emission through /c.pc/l.pc/o.pc/u.pc/d.pc/y.pc v17.00 (Ferland et al. 2017). In fitting the data, we adopt the ''bursty continuity'' SFH prior of Tacchella et al. (2022) motivated by the large equivalent widths of GN 42437. We allow for a wide range of stellar masses, metallicities, and ionization parameters, adopting log-uniform priors for all three physical properties. To measure the M ∗ and SFH, we focus on data that constrain the underlying rest-frame UV through optical continuum, as well as the strength of strong optical lines ([O /i.pc/i.pc/i.pc] and H 𝛼 ). An AGN may contribute to a fraction of the strong optical nebular lines (see subsection 7.2), but we are unable to provide strong constraints on the AGN contamination from the photometry. Indeed, the SED is well-fit by a template that only includes massive stars, as is common for Type /i.pc/i.pc AGN. We thus use the three FRESCO photometric measurements, as well as the NIRSpec spectroscopic redshift (z spec ), NIRSpec narrowband synthetic photometry, and observed nebular emission lines to estimate the stellar continuum and nebular emission lines.  \nWith /b.pc/a.pc/g.pc/p.pc/i.pc/p.pc/e.pc/s.pc, we infer log( 𝑀 ∗ / 𝑀 ⊙ ) = 7.9 ± 0.2. This is consistent with the typical stellar mass of 𝑀 UV ≈ -19 Lyman-break galaxies at 𝑧 ∼ 6inEndsleyetal. (2023) inferred using similar SFH assumptions. The SFR over the past 3 Myr inferred from the photometry plus spectra /b.pc/a.pc/g.pc/p.pc/i.pc/p.pc/e.pc/s.pc fitting is ≈ 30 𝑀 ⊙ yr -1 , while the average SFR over the past 3-10 Myr is < 0.2 𝑀 ⊙ yr -1 . We do not observe significant stellar mass formed beyond 10 Myr. Over 3 Myrs, 30 M ⊙ yr -1 is roughly consistent with the SFR of 11 M ⊙ yr -1 averaged over 10 Myr estimated from the attenuation-corrected H 𝛼 flux and the Kennicutt & Evans (2012) conversion (Table 1). This suggests that nearly all of the observed stellar light within GN 42437 has formed in the past 3 Myr. The extremely high [O /i.pc/i.pc/i.pc] 5008 Å and H 𝛼 equivalent widths (1644 ± 260 and 901 ± 145 Å, respectively) can only be explained by very young bursts of star formation (Leitherer et al. 1999). These  \nextreme [O /i.pc/i.pc/i.pc] equivalent widths are ≈ 3 times larger than the typical 𝑀 UV ≈ -19 𝑧 ∼ 6 Lyman-break galaxy (Endsley et al. 2023).  \nThe Chandra imaging tests whether there is a strong AGN within GN42437, but we do not detect X-ray emission from the galaxy. We calculated the upper-limit on the 2-10 keV flux in the 2 Ms combined event image using the CIAO tools function aprates Fruscione et al. (2006). The flux upper-limit is found to be 𝐹 𝑋, 2 -10 keV < 1 . 9 × 10 -16 erg s -1 cm -2 . We then calculated the X-ray 2-10 keV luminosity (L 𝑋 ) using the X-ray flux upper-limit and the spectroscopic redshift, and found 𝐿 𝑋 [ 2 -10keV ] < 4 . 5 × 10 43 erg s -1 (Table 1).  \nThe NIRCam imaging in Figure 1 shows an apparent companion to the North East of GN 42437. The SED fitting with the full HST + JWST observations puts a nearly equal probability that the companion source is at 𝑧 = 1 and 𝑧 = 5. This is largely because there is a 3 𝜎 detection in the HST/ACS F606W filter (0 . 012 ± 0 . 04 nJy). F606W probes rest-frame ∼ 900 Å at 𝑧 = 5 . 58, an unlikely wavelength to detect a 𝑧 = 5 . 5 galaxy. Unlike GN 42437, the apparent companion does not have a FRESCO grism detection (Oesch et al. 2023). The /b.pc/a.pc/g.pc/p.pc/i.pc/p.pc/e.pc/s.pc posterior suggests that the H 𝛼 flux of the companion would be between 1-8 × 10 -18 erg s cm -2 . There is a significant probability that the companion is too faint to be detected with FRESCO grism observations. The companion has a NIRSpec slit placed on it as part of the JWST Program 1211 (PI: Luetzgendorf). This data will stringently test whether this apparent companion galaxy is at a similar redshift as GN 42437.", '5 EMISSION LINE MEASUREMENTS': "In the top panel of Figure 3 we show the fully extracted onedimensional spectrum of GN 42437. There are many strong emission lines at 4.32, 3.3, 3.2, and 2.55 𝜇 m. We fit these lines using the A/s.pc/t.pc/r.pc/o.pc/p.pc/y.pc python package S/p.pc/e.pc/c.pcU/t.pc/i.pc/l.pc/s.pc (Astropy Collaboration et al. 2013, 2018, 2022). We assume a single constant local continuum (using an average continuum band-pass of 0.02 𝜇 m, or 200 Å) for each line and fit for the Gaussian line properties. For most spectral features we use a single Gaussian, but the high-resolution G235H+G395H observations resolve many interesting lines (for instance the [O /i.pc/i.pc] 3727, 3730 Å doublet in the middle left panel of Figure 3). We set upper limits to non-detected lines (such as [O /i.pc] 6300 Å and [N /i.pc/i.pc] 6583 Å) by integrating the error spectrum over ± 250 km s -1 around the expected line center. This is a conservative upper limit because the full extent of the lines are typically 300-400 km s -1 .  \nWe fit the Gaussian line width ( 𝜎 ) of the strong H 𝛼 emission line observed at 4.325 𝜇 m and find that it is 50 km s -1 . The publicly available python package /m.pc/s.pc/a.pc/f.pc/i.pc/t.pc 6 (de Graaff et al. 2023) estimates the spectral resolution at the observed wavelength of the H 𝛼 emission line using the measured size of the source, the location of the source in the MSA array (quadrant 1, 𝑖 = 350, 𝑗 = 141), and the positioning of the object within the slit to be 35 km s -1 . This suggests that the H 𝛼 emission is spectroscopically resolved with an intrinsic Gaussian width ( 𝜎 int = √︃ 𝜎 2 obs -𝜎 2 res ) of 36 km s -1 . The H 𝛼 line is well-fit by a single Gaussian profile (see discussion in section 8) and we do not detect a fainter broad emission component in any emission line (Saldana-Lopez et al. in preparation).  \nWeassumethat all of the observed nebular emission lines originate in the same nebular region and fix the velocity width of each line to the observed velocity width (intrinsic plus instrumental broadening) of H 𝛼 for each line using /m.pc/s.pc/a.pc/f.pc/i.pc/t.pc (de Graaff et al. 2023).  \nTable 3. Measured integrated fluxes and intensities for GN 42437 for the 13 emission lines detected at the > 3 𝜎 significance and upper-limits of nondetected lines. Upper-limits are indicated with ' < '. The first column gives the ionic species and rest-frame vacuum wavelength. The second column gives the observed integrated flux ratio of the emission line ( 𝐹 𝜆 ) relative to the integrated H 𝛼 flux ( 𝐹 H 𝛼 ). The third column gives the integrated intensity (attenuation-corrected) ratio relative to H 𝛼 , and the fourth column gives the intensity ratio relative to H 𝛽 assuming a temperature of 20,000 K and a density of 2033 cm -3 that we measure from the [O /i.pc/i.pc] doublet ratio (Stephenson et al. in preparation). We use a model to give the fluxes relative to H 𝛽 because H 𝛽 was between the NRS chips. Below the horizontal line we provide the observed H 𝛼 integrated flux ( 𝐹 H 𝛼 ; in units of 10 -20 erg s -1 cm -2 ); the H 𝛼 , [O /i.pc/i.pc/i.pc] 5008 Å, and [Ne /v.pc] 3427 Å rest-frame equivalent widths ( 𝐸𝑊 H 𝛼 , 𝐸𝑊 [ O/i.pc/i.pc/i.pc ] , 𝐸𝑊 [ Ne/v.pc ] ; in units of Å); [O /i.pc/i.pc/i.pc] 5008 Å attenuation-corrected luminosity ( 𝐿 [ O/i.pc/i.pc/i.pc ] in units of 10 42 erg s -1 ); and the nebular reddening ( 𝐸 ( 𝐵 -𝑉 ) ) measured from the Balmer line ratios of H 𝛿 /H 𝛼 and H7/H 𝛼 using the Cardelli et al. (1989) attenuation law.  \nThis improves the detection reliability of extremely faint lines. The lower right panel of Figure 3 shows the [S /i.pc/i.pc] 6718,6733 Å spectral region (considered a non-detection). The redder 6733 Å [S /i.pc/i.pc] line appears to have a significantly broader profile than the H 𝛼 emission line (compare the gold fit to the blue observations). If we did not fix the [S /i.pc/i.pc] 6733 Å line width, we would find a 3.6 𝜎 significant [S /i.pc/i.pc] 6733 Å total integrated line flux, but a line width that is more than twice the H 𝛼 line width. Combined with the fact that we do not detect significant [S /i.pc/i.pc] 6718 Å flux, we consider [S /i.pc/i.pc] 6733 Å to be heavily contaminated by noise. We fix the [S /i.pc/i.pc] 6733 Å line velocity width to that of H 𝛼 and conclude that the [S /i.pc/i.pc] 6733 Å lines are not significantly detected. This highlights the benefit of constraining the velocity widths of weaker lines. We tested whether this assumption holds for other strong lines in the spectrum of GN 42437. We find that every line is well-fit with a Gaussian that has an intrinsic velocity width of 36 km s -1 , except for [Ne /v.pc] 3427 Å. We return to the velocity width of [Ne /v.pc] in section 8.  \nWe detect 13 nebular emission lines at the > 3 𝜎 significance (Table 3). These lines range from [Ne /v.pc] 3427 Å in the blue portion of the spectrum (and highest ionization state) to a strong H 𝛼 in the red  \nTable 4. Attenuation-corrected emission line ratios of GN 42437 derived from the JWST NIRSpec spectra. The first column shows how we define the ratio. The second column gives the linear value of the ratio and the third column gives the logarithm of the ratio.  \nportion. Table 3 also lists a few upper-limits for lines that are below our detection threshold. The observed fluxes in the second column of Table 3 are given relative to the strong H 𝛼 emission line ( 𝐹 𝜆 / 𝐹 𝐻𝛼 ), along with the errors on the line ratios by propagating the uncertainty of the H 𝛼 integrated flux. We give the integrated fluxes relative to H 𝛼 instead of the canonical H 𝛽 line because the H 𝛽 line falls in the gap between the two NRS detectors. We estimate the equivalent widths (EW) of certain lines using the integrated line flux and the continuum estimated from the spectrum in mock JWST medium band filters using /s.pc/t.pc/s.pc/y.pc/n.pc/p.pc/h.pc/o.pc/t.pc (STScI Development Team 2018, 2020) in adjacent regions to the emission lines of interest: F210M (for [Ne /v.pc] 3427 Å), F360M (for [O /i.pc/i.pc/i.pc] 5008 Å), and F410M (for H 𝛼 ; Table 2). The integrated H 𝛼 fluxes and EWs are placed in the rest-frame using the estimated redshift and tabulated in the lower portion of Table 3.  \nWe detect 5 Balmer emission lines with a signal-to-noise ratio greater than 5 (Table 3). We used /p.pc/y.pc/n.pc/e.pc/b.pc (Luridiana et al. 2015) to determine the nebular reddening (E(B-V)) to be 0.031 mag using the observed H 𝛿 /H 𝛼 ratio, the Cardelli et al. (1989) attenuation law, the electron density measured from the [O /i.pc/i.pc] doublet of 2033 cm -3 (Stephenson et al. in preparation), and an assumed electron temperature of 2 × 10 4 K. The largest source of uncertainty in this measurement is not the error on the H 𝛿 /H 𝛼 ratio (which has a 10% error), but rather the unknown electron temperature. A robust electron temperature could be estimated using the temperature-sensitive [O/i.pc/i.pc/i.pc] 4363 Å auroral line, but [O /i.pc/i.pc/i.pc] 4363 Å is in the gap between the two gratings. We estimate the uncertainty on 𝐸 ( 𝐵 -𝑉 ) by taking the variation of possible 𝐸 ( 𝐵 -𝑉 ) values using temperatures of 1 . 5 × 10 4 (0.009 mag) and 2 . 25 × 10 4 (0.033 mag). This estimates an E(B-V) uncertainty of 0.013 mag. The H7/H 𝛼 ratio is more uncertain than the H 𝛿 /H 𝛼 ratio (a 15% error), but both the H 𝛿 /H 𝛼 and H7/H 𝛼 ratios provide consistent E(B-V) values using 2 × 10 4 K and the measured density. In the third column of the upper panel of Table 3 we provide the attenuation-corrected line ratios relative to H 𝛼 ( 𝐼 𝜆 / 𝐼 𝐻𝛼 ), and the forth column gives these attenuation-corrected line ratios relative to the unobserved, but more common, H 𝛽 line (assuming an electron temperature of 2 × 10 4 K and density of 2033 cm -3 ). Table 4 lists the attenuation-corrected line ratios used throughout the paper with both linear and logarithmic values.  \nThe observed-frame wavelength of the four strongest emission lines in the spectrum ([O /i.pc/i.pc/i.pc] 5008 Å, H 𝛼 , [O/i.pc/i.pc/i.pc] 4960 Å, and [Ne /i.pc/i.pc/i.pc] 3870 Å) define the spectroscopic redshift ( 𝑧 spec ) of GN 42437. These four emission lines are all observed on the NRS2 detector, but are found in both the G235H and G395H gratings. Regardless, the median of these four lines gives extremely consistent 𝑧 spec estimates for GN 42437 of 5 . 58724 ± 0 . 00005, where the error is the standard deviation of the four individual 𝑧 spec estimates. In fact,  \nTable 5. Kinematics of selected emission lines from GN 42437. The first column gives the line, the second and third columns give the observed Gaussian velocity dispersion ( 𝜎 obs ) in 𝜇 m and velocity (km s -1 ). The forth column gives the NIRSpec instrumental resolution using /m.pc/s.pc/a.pc/f.pc/i.pc/t.pc (de Graaff et al. 2023) using the observed NIRCam size of GN 42437 listed in Table 2. The last column gives the intrinsic line width by subtracting the instrumental resolution in quadrature from the observed width.  \n𝑧 spec of all four lines agrees to within 0.2% (5 km s -1 ). This velocity accuracy remains consistent both within the G395H grating (where there are three lines) and between the G235H and G395H grating. This 5 km s -1 is about one-third the quoted high-resolution grating accuracy (15 km s -1 ; Böker et al. 2023). Most of the discussion in section 8 centers around the very-high-ionization [Ne /v.pc] 3427 Å line which has a velocity offset of -16 ± 15 km s -1 from the median of the four strong emission lines. This high-quality wavelength consistency illustrates the exquisite delivered velocity accuracy of the high-resolution JWST gratings.  \nIn subsection 3.2 we found excellent agreement between the synthetic NIRSpec magnitudes and the observed NIRCam magnitudes. The fluxing can also be tested using emission doublets that have emissivity ratios set by atomic physics. We find that the attenuationcorrected [O /i.pc/i.pc/i.pc] 5008/4960 intensity ratio (observed in the G395H grating) of 2 . 95 ± 0 . 11 is easily within 1 𝜎 of the emissivity ratio (calculated using /p.pc/y.pc/n.pc/e.pc/b.pc to be 2.98). Similarly, the attenuation corrected [Ne /i.pc/i.pc/i.pc] 3870/3969 ratio (observed in the G235H grating and spectroscopically deblended with the high-resolution) is 4 . 3 ± 1 . 0, consistent within 1 𝜎 of the emissivity ratio (3.31). Thus, the relative and absolute fluxing of both NIRSpec gratings agree within 1 𝜎 of the expectations from atomic physics.", '6 UNPRECEDENTED VERY-HIGH-IONIZATION EMISSION LINE STRENGTHS': "The very-high-ionization emission lines (transitions arising from gas with ionization potentials greater than 54 eV) uniquely distinguish GN 42437's spectrum. The [Ne /v.pc] 3427 Å and He /i.pc/i.pc 4687 Å are detected at the 6.8 𝜎 and 3 𝜎 significance, respectively. The higher  \nFigure 3. The extracted spectrum of GN 42437. Upper Row: The full G235H (light blue) and G395H (dark blue) observed flux density. The lower x-axis gives the observed wavelength, while the upper x-axis gives the rest-frame wavelength using the redshift derived from the strong [O /i.pc/i.pc/i.pc]4960+5008 Å, H 𝛼 , and [Ne /i.pc/i.pc/i.pc] 3870 Å emission lines. Lower 2 Rows: Zoom-ins on individual emission lines from both gratings in ascending wavelength order. The rest-frame wavelengths of each feature is marked by a gray vertical line and the fitted Gaussian emission line is overplotted in gold. We detect 13 nebular emission lines at the > 3 𝜎 significance from the spectrum of GN 42437. The [S /i.pc/i.pc] 6718,6733 Å doublet is not detected at > 3 𝜎 , but shown to emphasize the upper-limits on the lines. H 𝛼 and [O /i.pc/i.pc/i.pc] 5008 have large rest-frame equivalent widths of 901 and 1644 Å, respectively.  \n<!-- image -->  \nRest Wavelength [Å]  \nsignificance of the [Ne /v.pc] emission is largely accounted for by the fact that the G235H observations have 5.4 times longer integrations than the G395H observations. The 2-dimensional [Ne /v.pc] 3427 Å line profile has a slightly extended spatial and spectral morphology (upper panel of Figure 4) that is at the same spatial (y-axis in Figure 4) location with a similar spatial morphology as other bright lines on the G235H spectrum. The [Ne /v.pc] line profile is spectroscopically resolved by the G235H observations with more than 8 spectroscopic pixels above the local noise level. The line peaks -16 ± 15 km s -1  \nfrom the rest-frame wavelength set by the other strong emission lines (see cyan line in Figure 4). Figure 5 shows both of these very-highionization lines. He /i.pc/i.pc 4687 Å has a lower significance, but the line peaks exactly at the expected rest-frame wavelength and has multiple pixels significantly above the noise level.  \nGN 42437 has unprecedented [Ne /v.pc] strength compared to local star-forming galaxies of similar stellar mass ( ∼ 10 8 M ⊙ ). The [Ne /v.pc]/H 𝛼 (or [Ne /v.pc]/H 𝛽 ) is 0 . 04 ± 0 . 006 (0 . 12 ± 0 . 02). This [Ne /v.pc]/H 𝛼 is 14 times larger than the median value from Izotov et al.  \nFigure 4. Upper Panel: The two-dimensional spectrum near the [Ne /v.pc] 3427 Å emission line. Lower Panel: Collapsed 1-dimensional spectrum of the same spectral region. The rest wavelength, defined by the other strong optical emission lines, is indicated by the cyan line in the upper panel. The [Ne /v.pc] 3427 Å emission is spatially distributed and centered in the middle of the extracted trace. The [Ne /v.pc] 3427 Å line is spectroscopically resolved and has an integrated 7 𝜎 significance (the error on the flux density is given as the gray band in the lower panel).  \n<!-- image -->  \n(2021) and seven times larger than the largest in that sample. The Izotov et al. (2021) sample is one of the largest 𝑧 ∼ 0 samples of [Ne /v.pc] emitting galaxies in low-metallicity star-forming galaxies in the local universe. GN 42437 has extreme rest-frame equivalent widths of many lines (H 𝛼 is 901 Å), but the [Ne /v.pc] rest-frame equivalent width of 11 ± 2 Å is 4 times larger than in Tol 1214-277. Tol 1214-277 is often used as the local example of a low-mass, star-forming galaxy with extreme [Ne /v.pc] emission (Izotov et al. 2004). The [Ne /v.pc] equivalent width is more consistent with values of 6-21 Å from z ∼ 1 AGN (e.g., Gilli et al. 2010; Mignoli et al. 2013; Cleri et al. 2023a). Using these 𝑧 ∼ 1 AGN, Ne5Ne3 = [Ne /v.pc]3427/[Ne /i.pc/i.pc/i.pc] 3870 has been suggested as a relatively metallicity-insensitive AGN diagnostic that tests the shape of the very-high ionizing spectrum (see subsection 7.2 below; Abel & Satyapal 2008; Cleri et al. 2023b). Ne5Ne3 values greater than 0.5, the gray dashed line in Figure 6, have been suggested to indicate the presence of pure AGN ionization (Cleri et al. 2023b). In Figure 6, we show the Ne5Ne3 value for GN 42437 as the light-blue square and the Izotov et al. (2021) sample as the gold circles. GN 42437 has a Ne5Ne3 = 0 . 26 ± 0 . 04 that is 13 times larger than the median of the low-redshift low-mass star-forming galaxies (Izotov et al. 2021). While GN 42437 has [Ne /v.pc] 3427 Å flux ratios and an equivalent width similar to local AGN, the [Ne /v.pc] line width is only 73 km s -1 , a factor of 5 smaller than the 400 km s -1 typically observed from 𝑧 ∼ 1 Type /i.pc/i.pc AGN (see section 8; Mignoli et al. 2013). This [Ne /v.pc] line width is narrower than the typical AGN at low-redshifts (e.g., Greene & Ho 2004; Mullaney et al. 2013; Kormendy & Ho 2013).  \nThe He/i.pc/i.pc/H 𝛽 ratio of 0 . 1 ± 0 . 04 is similarly extreme, although with larger uncertainties, compared to local star-forming galaxies with similar stellar mass. This He /i.pc/i.pc/H 𝛽 ratio puts GN 42437 in the AGN portion of the Shirazi & Brinchmann (2012) diagram (horizontal gray dashed line in Figure 6). While these He /i.pc/i.pc/H 𝛽 values typically indicate AGN, the Shirazi & Brinchmann (2012) demarcation does not differentiate between other hard-ionizing sources such as high-mass X-ray binaries (HMXBs) or Ultra-luminous X-ray  \nSources ULXs (see subsection 7.2). The He /i.pc/i.pc/H 𝛽 from GN 42437 is five times stronger than the largest value in the low-mass low-redshift sample of Berg et al. (2019), and 1.4 times larger than the largest local value recorded in Thuan & Izotov (2005). The [Ne /v.pc]/He /i.pc/i.pc ratio of 1 . 2 ± 0 . 5 is higher than all but one galaxy in the Izotov et al. (2021) [Ne /v.pc] sample, J1222 + 3602 (also discussed in Izotov et al. 2007). J1222 + 3602 is the only low-mass galaxy that Izotov et al. (2021) suggest harbors an AGN.  \nA final component of the high ionization structure of GN 42437 is the X-ray non-detection. X-ray detections are strong indicators of HMXBs and Compton-thin AGNs. Using the archival Chandra observations, we find that GN 42437 has a 2-10 keV X-ray Luminosity ( 𝐿 𝑋 ) upper-limit of < 4 . 5 × 10 43 erg s -1 , an 𝐿 𝑋 /SFR upper limit of < 4 × 10 42 erg s -1 M ⊙ yr -1 , and the upper limit of 𝐿 𝑋 / 𝑀 ∗ < 5 . 6 × 10 35 . The X-ray non-detection rules out the presence of a luminous AGN, but it does not rule out a lower luminosity or a Compton thick AGN.  \nWhile local galaxies with similar stellar mass have strong [Ne /v.pc] and He /i.pc/i.pc emission, the strength of the very-high-ionization emission in GN 42437 exceeds all of these examples. The very-high-ionization emission lines of GN 42437 are most comparable to local AGN.", '7 A UNIFIED IONIZATION STRUCTURE FOR GN 42437': 'GN 42437 has an impressive array of observed and stringent upperlimits of a vast array of nebular emission lines (Table 3). These observations probe the classical low and mid-ionization zones (e.g. [S /i.pc/i.pc] 6733 Å and [O /i.pc/i.pc] 3727,3730), to the high-ionization (typically defined by [O /i.pc/i.pc/i.pc] 5008 Å). However, what sets GN 42437 apart from most observations of either low or high-redshift star-forming galaxies is the detection of the very-high-ionization zone that is defined by the He /i.pc/i.pc 4687 Å and [Ne /v.pc] 3427 Å detections described above. Here we walk through the full ionization structure of the galaxy, looking for a single simplified model that explains all of the nebular emission within GN 42437. The most robust way to do this is to compare the ratios of nebular emission lines to photoionization models. The goal is to reproduce as many of the emission lines as possible with a single theoretical model. We start in subsection 7.1 with the low to high-ionization zones ([S /i.pc/i.pc] to [O /i.pc/i.pc/i.pc]) because these are commonly observed in star-forming galaxies across all redshifts and provide a baseline to compare to GN 42437. We then move on to the very-highionization zone and explore whether any single model explains the unprecedentedly strong [Ne /v.pc] 3427 Å emission (subsection 7.2).', '7.1 Low to high-ionization zones resemble a low-metallicity star-forming galaxy': 'Many of the nebular emission line strengths and ratios of GN 42437 look like a typical low-metallicity star-forming galaxy. In Figure 7 we plot the [O /i.pc/i.pc/i.pc]/H 𝛽 ratio versus the [N /i.pc/i.pc] 6585 Å/H 𝛼 (left panel) and the [S /i.pc/i.pc] 6718+6733/H 𝛼 (right panel). As a reminder, the H 𝛽 flux used for the [O /i.pc/i.pc/i.pc]/H 𝛽 ratio here, and elsewhere in the text, has been extrapolated from the attenuation-corrected H 𝛼 emission because H 𝛽 fell between the NRS1 and NRS2 detectors (see section 5). Both [N/i.pc/i.pc] and [S /i.pc/i.pc] have detection significance less than 3 𝜎 (Table 4) and we only mark their upper-limits on Figure 7. Regardless, the upper-limits place GN 42437 within the "star-forming" locus of the classical Baldwin, Phillips and Terlevich (BPT; Baldwin et al. 1981) diagrams (Kewley et al. 2001; Kauffmann et al. 2003). GN 42437 resides in a similar location as other 𝑧 > 4 star-forming galaxies in the BPT diagrams (dark-blue points; Cameron et al. 2023; Sanders  \nFigure 5. The very-high-ionization emission lines (ionization potentials greater than 54 eV) in GN 42437. The left panel shows [Ne /v.pc] 3427 Å (ionization potential greater than 97 eV) and the right panel shows the He /i.pc/i.pc 4687 Å (ionization potential greater than 54 eV) features. The profiles are placed into velocity space using the redshift of the strong optical emission lines. Zero-velocity is marked by a gray dashed line. The gray ribbons show the 1 𝜎 error on the flux density. The [Ne /v.pc] and He /i.pc/i.pc lines are detected at the 6.8 𝜎 and 3.0 𝜎 significance, respectively. The gold lines show a fit using a velocity width ( 𝜎 ) that is fixed to the H 𝛼 velocity width. The light-blue line in the [Ne /v.pc] 3427 Å panel shows a fit where the velocity width is free to vary. At the depth of the observations, the observed [Ne /v.pc] profile is 2.5 𝜎 broader than the H 𝛼 profile.  \n<!-- image -->  \net al. 2023). The upper limits of GN 42437 also place it in a similar location in the BPT as local galaxies that have very-high-ionization emission lines such as [Ne /v.pc] and He /i.pc/i.pc (Izotov et al. 2004; Thuan & Izotov 2005; Izotov et al. 2007; Shirazi & Brinchmann 2012; Kehrig et al. 2015; Feltre et al. 2016; Senchyna et al. 2017; Berg et al. 2018, 2019; Izotov et al. 2021). These upper limits place strong physical constraints on the low- and high-ionization emission within GN 42437.  \nThe comparison in section 6 of [Ne /v.pc] and He /i.pc/i.pc in GN 42437 to local galaxies suggested that GN 42437 likely hosts an AGN, but the low ionization lines traditionally used to separate AGN from starforming galaxies indicate that GN 42437 could be a star-forming galaxy (GN 42437 is well-within the star-forming locus in Figure 7). Asimilar conclusion is found when considering the [O /i.pc] 6300 Å nondetection. However, Figure 7 also includes the broad-line selected (Type /i.pc) AGN from Maiolino et al. (2023) as light-blue diamonds. These AGN have strong-line ratios indistinguishable from GN 42437 or star-forming galaxies at either low- or high-redshift. As discussed further in subsection 7.2, this strongly cautions against using the traditional strong-line diagnostics to separate AGN and star-forming systems in low-metallicity, low-mass, strongly star-forming systems with either high- or low-redshift observations (Groves et al. 2004, 2006; Kewley et al. 2013; Reines & Volonteri 2015; Scholtz et al. 2023; Dors et al. 2024).  \nThe observed strong-lines provide constraints on the metallicity and ionization of GN 42437. While the temperature-sensitive [O/i.pc/i.pc/i.pc] 4363 auroral line falls between the G235H and G395H grating configuration, the R23, Ne3O2, and O 32 strong-line ratios can approximate the 12+log(O/H). Using the recent calibrations with  \nJWST data from Sanders et al. (2023), we estimate 12+log(O/H) to be 7.36, 7.19, and 7.31 ( ∼ 5% Z ⊙ ) using these three different line ratios, respectively. A recent calibration of local high-ionization emission-line galaxies that uses O 32 from Izotov et al. (2024) predicts a slightly higher 12+log(O/H) = 7 . 54, or 7% Z ⊙ . This implies a low, but non-zero, metallicity of 5-7% Z ⊙ for GN 42437.  \nHowever, strong-line metallicity calibrations typically fail to accurately estimate the metallicity of AGNs (Groves et al. 2006). This is largely because the traditional strong-line calibrations do not include very-high-ionization states that contain a large fraction of the total metal emission. In these situations models of the nebular structure provide a better handle on the full ionization structure and metallicity of the gas. We use photoionization models from Richardson et al. (2022) 7 to simultaneously test the ionization and metallicity of GN 42437. These models are largely used because they also include contributions from an intermediate mass black hole (IMBHs), which we will use in subsection 7.2 to explain the [Ne /v.pc]. These models use CLOUDYv17.00 (Ferland et al. 2017) and stellar population synthesis models from BPASS v2.0 (Eldridge & Stanway 2009; Stanway et al. 2016; Eldridge et al. 2017) that include binary synthesis stellar evolution, a 20 Myr instantaneous burst stellar population with a Kroupa initial mass function (Kroupa 2001), and a high-mass cutoff at 300 M ⊙ . We tested a larger grid of these models tailored to GN 42437, but that expanded grid does not significantly alter the results discussed here. The models include metallicities in 11 steps that have the Nicholls et al. (2017) solar abundance pattern. The stellar and nebular abundances are forced to match each other. We  \nIMBH Fraction  \nFigure 6. The He /i.pc/i.pc 4687/H 𝛽 versus [Ne /v.pc] 3427/[Ne /i.pc/i.pc/i.pc] 3870 emission line diagnostic diagram. This probes the strength of the very-high-ionization emission lines found in GN 42437 (light blue square) in a relatively metallicity independent way. The gold circles show a local sample of low-metallicity low-mass star-forming galaxies with [Ne /v.pc] detections (Izotov et al. 2021). GN42437is above the He /i.pc/i.pc/H 𝛽 demarcation proposed to separate AGN from star-forming galaxies (gray horizontal dashed-dotted line; Shirazi & Brinchmann 2012) and near the [Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] ratio proposed to separate AGN (gray vertical dashed line; Cleri et al. 2023b). Stellar population only models (darkblue star) produce ∼ 6 orders of magnitude lower [Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] than observed in GN 42437. In comparison, Intermediate Mass Black Holes (IMBH) plus massive star models from Richardson et al. (2022) are also shown. The best-fit model to the full ionization structure of GN 42437, with a 10 5 M ⊙ black hole, log 𝑈 = -2, and Z = 0.1Z ⊙ , is shown as a solid blue line. Each model point on the blue line is color-coded by the fraction of the hydrogen ionizing photons provided by the IMBH, from 0% (blue) to 100% (yellow). GN 42437 most closely matches the ∼ 30% IMBH model. Variations in ionization parameter are shown by the gold-shaded region, from log 𝑈 = -3 (top) to log 𝑈 = -1 (bottom), while variations in metallicity are shown by the blue-shaded region, from Z = 0.05 Z ⊙ (bottom) to Z = 0.25 Z ⊙ (top). The ionization parameter has a large impact on the very-high-ionization emission lines, while metallicity is a sub-dominate effect. We also compare to a 10 3 M ⊙ IMBH model (blue dashed line). The lower mass IMBH model shifts the line ratios by ∼ 0.4 dex to higher ionization. Finally, a mixture of 70% massive stars and 30% shocks with varying shock velocities is shown as a gold line (Allen et al. 2008). This mixture is required to approach the non-detection of the low-ionization lines ([S /i.pc/i.pc], [N /i.pc/i.pc], and [O /i.pc/i.pc]), but fails to reproduce the observed [Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] or He/i.pc/i.pc/H 𝛽 ratio. The IMBH plus massive star model simultaneously explains the very-high-ionization and low-ionization emission lines in Figure 7 and Figure 9.  \n<!-- image -->  \nuse the intrinsic fluxes, a closed geometry, and a constant density of 100 cm -3 . In this section, all models have no contributions from intermediate mass black holes. We discuss models with IMBHs in the next section.  \nFigure 8 shows the [Ne /i.pc/i.pc/i.pc]/[O /i.pc/i.pc] line ratio versus the [O /i.pc/i.pc/i.pc]/H 𝛽 line ratio for the photoionization models using three different metallicities (0.05, 0.1, 0.25 Z ⊙ ). Each point corresponds to a changing ionization parameter (log 𝑈 ), as highlighted by the 0.1 Z ⊙ track. This diagram is a function of both ionization and metallicity, with four independent emission lines. The observations of GN 42437 (light blue  \nsquare) are consistent with the 0.1 Z ⊙ and log 𝑈 = -2 to -1 . 5 model. This model also satisfies the upper limits on the low-ionization emission lines in the optical BPT (Figure 7). The inferred log 𝑈 ≈ -2 is very similar to the value found for other low-redshift very-highionization emission line galaxies (Berg et al. 2021). For the rest of this paper we use these parameters as our \'default\' nebular properties because it robustly matches the low-ionization structure. While the low-ionization and high-ionization lines can be reproduced with a stellar population-only model, in the next section we show that these models under produce the very-high-ionization lines, like [Ne /v.pc], by six orders of magnitude. Reproducing the very-high-ionization emission lines in low-metallicity star-forming galaxies has historically been theoretically challenging (Izotov et al. 2004; Thuan & Izotov 2005; Shirazi & Brinchmann 2012; Senchyna et al. 2017; Schaerer et al. 2019; Berg et al. 2021; Izotov et al. 2021; Olivier et al. 2022).', '7.2 Source of full ionization structure': "The nebular structure of GN 42437 is characterized by very-highionization [Ne /v.pc] and He /i.pc/i.pc, but also with relatively typical lower ionization emission lines. Can we explain both with a single selfconsistent model? subsection 7.1 showed that stellar-population-only photoionization models can reproduce the low and high-ionization emission lines with a metallicity between 5-10% Z ⊙ and log 𝑈 between -2 and -1.5. To test whether massive stars produce the observed very-high-ionization emission lines, we show the very-highionization emission line diagnostic of [Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] versus He /i.pc/i.pc/H 𝛽 in Figure 6. The blue star shows a photoionization model that only uses a massive star population as the ionizing source of the gas. The stellar population-only model fails to generate the observed [Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] in GN 42437 (light blue square) by over six orders of magnitude. Metallicity (blue shaded region) and log 𝑈 (gold shaded region) changes the predicted very-high-ionization line ratios, but no massive star only model gets closer than 5 orders of magnitude below the observations. If massive stars are the only ionizing source within this galaxy, we should not detect He /i.pc/i.pc or [Ne /v.pc] 3427 Å at any significance (Izotov et al. 2004; Thuan & Izotov 2005; Shirazi & Brinchmann 2012; Feltre et al. 2016; Senchyna et al. 2017; Schaerer et al. 2019; Berg et al. 2021; Izotov et al. 2021; Olivier et al. 2022; Feltre et al. 2023; Garofali et al. 2024). Yet GN 42437 has a 7 𝜎 [Ne /v.pc] detection that is 4% of the flux of the H 𝛼 emission. What could possibly produce these extreme very-high-ionization emission lines?  \nHigh-mass X-ray Binaries (HMXBs) are one of the most compelling very-high-ionization photon sources (HMXBs; Fragos et al. 2013; Schaerer et al. 2019; Senchyna et al. 2020; Simmonds et al. 2021; Garofali et al. 2024; Lecroq et al. 2024). HMXBs are approximately stellar mass (e.g., 10-100 M ⊙ ) black holes that are in a binary pair with a star. Mass-transfer from the star to the black hole heats gas to extreme temperatures and produces prodigious amounts of hard-ionizing and X-ray photons. These hard-ionizing photons are usually, but not always, produced with a temporal lag after the onset of star formation because massive stars must first evolve through their lives and explode as supernovae to create the black hole in the binary pair. This is why there is usually a correlation of the HMXB L 𝑋 and stellar age or H 𝛽 equivalent width (Schaerer et al. 2019; Senchyna et al. 2020; Garofali et al. 2024). The extreme H 𝛼 equivalent width (901 ± 145 Å; Table 4) and young SED fit to the JWST + HST observations suggest that GN 42437 is too young to have a significant HMXB population. Finally, the [Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] and He /i.pc/i.pc/H 𝛽 values are more extreme than HMXB models can produce even at older ages  \nFigure 7. Classical ionization diagnostics of the [O /i.pc/i.pc/i.pc]/H 𝛽 and [N /i.pc/i.pc]6585/H 𝛼 ( left panel ) and [S /i.pc/i.pc]6718+6733/H 𝛼 ( right panel ). The H 𝛽 is inferred from the attenuation corrected H 𝛼 because H 𝛽 is not covered by our observations (see section 5). [S /i.pc/i.pc] and [N /i.pc/i.pc] are undetected for GN 42437 (thick light-blue arrow), but these non-detections still place GN 42437 within the traditional low-redshift star-forming locus (gray lines; Kewley et al. 2001; Kauffmann et al. 2003). Local low-metallicity galaxies with very-high-ionization emission lines are included in gold (circles and squares are Berg et al. 2019; Izotov et al. 2021, respectively), suggesting that GN 42437 is consistent with these local galaxies. High-redshift star-forming galaxies are included as dark-blue circles (Cameron et al. 2023; Sanders et al. 2023) to show that GN 42437 has similar strong-line non-detections as many previous observations within the epoch of reionization. High-redshift broad-line (Type /i.pc) AGN from Maiolino et al. (2023) are shown as light-blue diamonds. Broad-line AGN reside in a similar portion of this diagram as GN 42437 and other low-metallicity star-forming galaxies. A dark blue line in both panels shows the same intermediate mass black hole (IMBH) model as in Figure 6 with a metallicity of 10% solar, log 𝑈 = -1 . 5 and a varying contribution of the IMBH to the total number of hydrogen ionizing photons as shown by the colored points (Richardson et al. 2022). The gold lines show a mixture of shock and photoionized gas using an SMC metallicity (20% Z ⊙ ) and varying shock velocities (Allen et al. 2008). The upper limits on both the [S /i.pc/i.pc] and [N /i.pc/i.pc] lines are consistent with a model where an IMBH contributes about ∼ 30% of the ionizing photons (the light purple point, the forth point from the right in the [S /i.pc/i.pc] plot), as suggested by the [Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] ratio in Figure 6.  \n<!-- image -->  \n( > 20Myr; Garofali et al. 2024; Lecroq et al. 2024). Thus, the strength of the [Ne /v.pc] relative to other emission lines strongly suggests that HMXBs are unlikely to generate the observed very-high-ionization emission.  \nShocks have been invoked to explain the strong He /i.pc/i.pc and [Ne /v.pc] emission in lower redshift galaxies (Thuan & Izotov 2005; Izotov et al. 2012, 2021). Supernovae or stellar winds drive gas faster than the sound speed into adjacent gas. The fast moving gas creates a shock wave that develops thin and dense layers of gas that emits high-energy ionizing photons as a precursor to create high-ionization gas within the galaxy. The kinetic energy of the shock is radiated away in both the dense low-ionization gas in the shock front and the low-density highly-ionized precursor gas (Dopita et al. 2002).  \nWe use the shock plus precursor models of Allen et al. (2008) to test whether shocks could create the observed ionization structure. We use the lowest available metallicity (the SMC; 20% Z ⊙ ) and a wide range of shock velocities. Shocks do produce significant veryhigh-ionization emission, with He /i.pc/i.pc/H 𝛽 and [Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] emission line ratios approaching the observed values if the shock has a velocity near 600 km s -1 . This velocity is broadly consistent with previous studies of low-redshift low-metallicity galaxies (Izotov et al. 2012, 2021).  \nWhile fast radiative shocks could explain the He /i.pc/i.pc/H 𝛽 and [Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] ratios, the dense shocked regions produce significant lower ionization emission. The shock models predict that we should have strongly detected both [N /i.pc/i.pc]6585 Å and [S /i.pc/i.pc]6718+6733 Å emission. The shock models with a velocity of 600 km s -1 that  \nreproduce the strong [Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] predict log([S /i.pc/i.pc]/H 𝛼 ) = -0 . 36. This is 0.8 dex above the upper-limits of the observations. A similar over-prediction of the low-ionization emission lines occurs to drive down 𝑂 32 in Figure 9. We observe an 𝑂 32 ratio of 33 ± 14, but the 600 km s -1 shock model predicts much lower ionization of 𝑂 32 = 2.  \nTo reduce the strong low-ionization lines, some studies have suggested that only a fraction of the nebular emission lines come from a spatially distinct radiative shock and another fraction comes from gas ionized by massive stars. If we assume that 30% of the nebular emission comes from the 600 km s -1 radiative shock and the rest from the default photoionization model of subsection 7.1, we find that log([S /i.pc/i.pc]/H 𝛼 ) decreases to -0.9. We show this model in Figure 7 and it is still 0.2 dex above the [S /i.pc/i.pc]/H 𝛼 upper-limits. Similarly, the gold line on Figure 9 shows that this shock and stellar population mixture rises from O 32 = 2 to 14. While these mixture models still produce too much low-ionization emission (e.g. 𝑂 32 is below the observations), the mixture also produces too little [Ne /v.pc] compared to our observations. This is shown by the gold lines in Figure 6 and Figure 9. For instance, the [Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] ratio of the fastest shock plus massive star mixture model is 0.15, nearly a factor of two below the observations of GN 42437. The low-ionization lines and O 32 ratio requires even lower shock fractions that would reduce the very-highionization emission features even further below the observations. It could be possible that the shocks are driven into already ionized gas (Izotov et al. 2021) or that lower metallicity shock models could reproduce the line ratios, but the O 32 ratio will still be too large. There  \nFigure 8. Photoionization model predictions with massive stars only for the [Ne /i.pc/i.pc/i.pc] 3870 Å/[O /i.pc/i.pc] 3277+3730 Å versus [O /i.pc/i.pc/i.pc] 5008 Å/H 𝛽 flux ratio. The H 𝛽 flux is extrapolated from the H 𝛼 flux because the observations do not cover it (see section 5). Each line is a photoionization model with a single metallicity (0.05, 0.1, and 0.25 Z ⊙ ) and ionization parameter (log 𝑈 ) varying from-3(dark purple point shown on the 0.1 Z ⊙ curve) to -1 (yellow point). The GN 42437 observations lie close to the 0.1 Z ⊙ and log 𝑈 = -2 point. These parameters are roughly similar to the strong-line metallicity estimate of 57% Z ⊙ . These default parameters accurately reproduce the nebular structure of GN 42437 when including the IMBH. We use the 0.1 Z ⊙ and log 𝑈 = -2 model as our default throughout the paper.  \n<!-- image -->  \nis not a combination of current shock and massive star models that match the full nebular structure.  \nFinally, we are left with the contributions from accretion onto a black hole. GN 42437 has a low stellar (10 7 . 9 M ⊙ ) and dynamical mass (10 8 . 5 M ⊙ ). Any black hole in GN 42437 is unlikely to be > 10 7 M ⊙ , or a supermassive black hole (see section 8). Rather the black hole is likely on the order of 10 4 -7 M ⊙ . These black holes are called intermediate mass black holes (IMBH; Greene et al. 2020). We use the photoionization models discussed in subsection 7.1 from Richardson et al. (2022) to test whether IMBHs create the nebular structure. These models include the ''qsosed'' model for the IMBH SED (Kubota & Done 2018) that uses a self-consistent prescription for a 10 5 M ⊙ IMBH that accretes at sub-Eddington (0.1 M edd ) with a dimensionless spin parameter of 𝛼 ∗ = 0. IMBHs and massive stars can both contribute to the total number of ionizing photons of the models. We use models with IMBHs contributing 0, 4, 8, 16, 32, 50, 75, and 100% of the hydrogen ionizing photons (different point colors in Figure 6).  \nThe IMBH models successfully reproduce the entire ionization structure of GN 42437 using an IMBH fraction near 30%, log 𝑈 near -2, and 0.1 Z ⊙ . The blue line in Figure 6 shows the rise of the very-high-ionization lines ([Ne /v.pc] and He /i.pc/i.pc) along log-log tracks that depend on the fraction of hydrogen ionizing photons from the IMBH.Thisemphasizesthatthe[Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] ratio largely tracks the fraction of ionizing photons arising from the IMBH (a similar trend is also found from super massive black holes in Abel & Satyapal 2008). These tracks in the [Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] versus He /i.pc/i.pc/H 𝛽 plane are largely  \nFigure 9. The [Ne /v.pc]3427/[Ne /i.pc/i.pc/i.pc]3870 versus [O /i.pc/i.pc/i.pc]5008/[O /i.pc/i.pc]3727+3730 ( 𝑂 32 ) flux ratio demonstrates that IMBH photoionization models best match the observations of GN 42437 (blue square). Models from Richardson et al. (2022) with a fraction of the ionizing photons coming from intermediate mass black holes (IMBH) are shown as a solid blue line (for a 10 5 M ⊙ IMBH) and a dashed line (for a 10 3 M ⊙ IMBH). The two models are both consistent with the observed GN 42437 ionization structure. The gold line shows a model that combines 30% of the emission from shocks (with varying shock velocities) with 70% from massive stars. The shock mixture models produce too low of both [O /i.pc/i.pc/i.pc]/[O /i.pc/i.pc] and [Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] to match the observations of GN 42437.  \n<!-- image -->  \ninsensitive to metallicity: populations with 0.05 Z ⊙ and 0.25 Z ⊙ are indistinguishable from the default 0.1 Z ⊙ case (compare the blue shaded region). The gold shaded region illustrates the impact of log 𝑈 on the ratios: log 𝑈 = -1 generates larger [Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] at lower He /i.pc/i.pc/H 𝛽 , while lower log 𝑈 produces less [Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] at higher He /i.pc/i.pc/H 𝛽 . This is because He /i.pc/i.pc has a lower ionization potential that can be produced at lower log 𝑈 than [Ne /v.pc]. Unlike the BPT diagnostics discussed in subsection 7.1, the [Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] versus He/i.pc/i.pc/H 𝛽 diagram strongly diagnoses the contribution of AGNs even at low-metallicity.  \nThe IMBH plus massive star photoionization models matches more than just the very-high-ionization lines. The IMBH models match the upper-limits of the low-ionization emission lines, such as [N/i.pc/i.pc]6585/H 𝛼 and [S /i.pc/i.pc] 6718+6733 Å/H 𝛼 (see the points in Figure 7). Even the observed 𝑂 32 ratio, shown in Figure 9, is consistent with a model where IMBHs contribute approximately 30% of the hydrogen ionizing photons. The IMBH plus massive star model is the only model tested that matches the full ionization structure of GN 42437.  \nTheRichardson et al. (2022) models include multiple prescriptions for the IMBH ionizing spectrum, all of which are consistent with the observations of GN 42437. For instance, the dashed lines in Figure 6 and Figure 9 show a 10 3 M ⊙ black hole mass (M BH ) model instead of a 10 5 M ⊙ model. The 10 3 M ⊙ model with the default log 𝑈 and metallicity is ≈ 0.4 dex larger in [Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] than the observations of GN 42437. The Richardson et al. (2022) models also include a disk plus power law (''disk-plaw'') model that takes a disk model from Mitsuda et al. (1984) with a temperature set by the mass of  \nthe black hole (Peterson 1997), and combines it with a power-law with an exponent of 2.1. These models also match the full ionization structure of GN 42437 with the same default parameters, but an IMBH fraction of 16%.  \nHere, we have tested the ionizing source of the entire ionization structure of GN 42437. We have ruled out HMXBs as the dominate source of very-high-ionization photons because of the large H 𝛽 equivalent widths, the large He /i.pc/i.pc/H 𝛽 , and the elevated [Ne /v.pc]/H 𝛼 ratios. Shock models can produce very-high-ionization emission lines, but shocks over-predict the low-ionization emission lines by 0.8 dex. Models that require approximately 30% of the ionizing photons from an intermediate mass black hole and 70% from a massive star population reproduce the full nebular structure of GN 42437 without any more fine-tuning. Thus, we conclude that the combination of a burst of young stars and accretion onto an IMBH is the most likely origin of the strong very-high-ionization emission lines in GN 42437.", '8 BLACK HOLE MASS CONSTRAINTS': "In the previous section we discussed how a combination of ionizing photons from a massive star population and a black hole matches both the very-high-ionization emission and low-ionization emission of GN 42437. Now that we have revealed that GN 42437 hosts a black hole, we can ask: what are the properties of the black hole in GN 42437? The most fundamental parameter of a black hole is its mass (M BH ). Stellar and gaseous kinematics trace M BH because black holes have large gravitational potentials that cause significant high-velocity motions. In Saldana-Lopez et al. (in preparation) we explore the broad emission lines in the entire NIRSpec sample. We find that all of the emission lines in GN 42437 except [Ne /v.pc] 3427 Å are narrow. Even the strong [O /i.pc/i.pc/i.pc] 5008 Å (or H 𝛼 ) line has no detectable secondary broad component (see Figure 3). The Akaike Information Criterion (AIC; Akaike 1974) changes by less than 10 when including a second component. This suggests that there is less than a 2.5 𝜎 change in the residuals when adding a broader component. We do not detect a broad H 𝛼 component that would typically indicate the broad-line region of a Type /i.pc AGN.  \nThe non-detection of the broad-line region means that GN 42437 can be considered a narrow-line AGN, or Type /i.pc/i.pc AGN, within the Epoch of Reionization. At lower redshift Type /i.pc/i.pc AGN are usually separated using the traditional BPT diagrams. However, Figure 7 strongly shows that using the BPT diagrams will not be possible for all AGN at high-redshift. This is emphasized by the fact that both GN42437andmanyoftheType/i.pc AGN from the literature occupy the star-forming locus of the diagram. Rather, very-high-ionization line ratios, such as [Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] and He /i.pc/i.pc/H 𝛽 , are required to cleanly separate the AGN population at high-redshift (Figure 6). As with local Type /i.pc/i.pc AGN, it is challenging to robustly estimate M BH of GN 42437 because we do not have a direct probe of the sphere of influence of the black hole in GN 42437. Here, we explore a variety of different methods to approximate M BH in GN 42437 and Table 6 summarizes these different approximations.  \nThe dynamical mass (M dyn ) puts a fairly strict upper-limit on the total mass of the galaxy and the black hole itself. Using the intrinsic H 𝛼 line width (Table 5), the F444W size of the galaxy (the filter with H 𝛼 ; Table 2), and equation 2 from Maiolino et al. (2023) that assumes a Sérsic profile (Cappellari et al. 2006), we estimate log ( 𝑀 dyn / 𝑀 ⊙ ) = 8 . 5 ± 0 . 3. This is the maximum amount of mass possible within the MSA slit. We have already used the NIRSpec and NIRCam data to estimate that GN 42437 has a stellar mass of log(M ∗ / 𝑀 ⊙ ) = 7 . 9, or about 25% of the dynamical mass. Nearly  \nall of the observed stellar mass has formed within the last 3 Myr, indicating that the gas mass likely significantly contributes to the total mass. If we make the assumption that the local KennicuttSchmidt law holds for this galaxy (Kennicutt 1998; Gao & Solomon 2004; Kennicutt & Evans 2012) and use the observed star formation rate surface density ( Σ SFR = 10 M ⊙ yr -1 kpc -2 ), we obtain an order of magnitude estimate of the gas mass in the 300 pc observed region of GN 42437 to be near 1 -3 × 10 8 M ⊙ . One-third to an entirety of the dynamical mass is likely in gas. This means that the combination of the approximated gas mass and stellar mass in GN 42437 is at least two-thirds and up to 100% of the estimated dynamical mass of the galaxy. Therefore, the estimated M dyn , M ∗ and gas mass limits the M BH to be no more than a few times 10 7 M ⊙ .  \nThe black hole masses of narrow-line AGN at low-redshifts are usually estimated using the correlation between M BH and the stellar velocity dispersion of the central bulge ( 𝜎 ∗ ; Gebhardt et al. 2000; Ferrarese & Merritt 2000; Kormendy & Ho 2013; Greene et al. 2020). This M BH -𝜎 ∗ relation has been used to estimate M BH from thousands of Type /i.pc/i.pc AGN in the low-redshift universe (Heckman et al. 2004; Kauffmann & Heckman 2009; Greene et al. 2020). The high spectral resolution of the G395H grating allows us to spectroscopically resolve the H 𝛼 line. The intrinsic H 𝛼 velocity dispersion of 36 km s -1 estimates the velocity dispersion of the gas. However, there is a statistical offset between the velocity dispersion of the gas and the stars. We follow Übler et al. (2023) to correct for the stellar velocity dispersion by multiplying it by a factor of 1.4 (Bezanson et al. 2018) to estimate 𝜎 ∗ = 50 km s -1 . The left panel of Figure 10 shows local observations of the M BH -𝜎 ∗ in gold points and the gold line shows the best-fit relation from Greene et al. (2020) that uses a large compilation of local black holes and upper-limits from local galaxies (the 'All, limits' relation). We use this relation and the relation from Kormendy & Ho (2013) to estimate log ( M BH / M ⊙ ) = 5 . 7 and 5.9 respectively. We include the Greene et al. (2020) estimate in the left panel of Figure 10 as a light-blue circle.  \nThe intrinsic [Ne /v.pc] 3427 Å line width is 73 ± 15 km s -1 and 2.5 𝜎 broader than the H 𝛼 emission line (Table 5). The left panel of Figure 5 shows two Gaussian fits to the [Ne /v.pc] line: (1) fixing the [Ne /v.pc] velocity width to that of the H 𝛼 velocity width that has been corrected for the wavelength-dependent spectral resolution changes (gold region; de Graaff et al. 2023), and (2) an unconstrained velocity-width (light blue region). The light blue line fits the observed wings of the [Ne /v.pc] line better at the 2.5 𝜎 significance. This suggests that the [Ne /v.pc] emission line is broader than the other observed emission lines.  \nOnly the hard ionizing photons from the black hole accretion disk can produce [Ne /v.pc]. Observations of local AGN show that the line widths of very-high-ionization (e.g. [Ne /v.pc]) emission lines increase with increasing ionization potential (Dasyra et al. 2008, 2011). Additionally, [Ne /v.pc] 3427 Å has a larger critical density than other optical emission lines (10 7 cm -3 for [Ne /v.pc] and 10 5 cm -3 for [O /i.pc/i.pc/i.pc] 5008 Å, for instance). This suggests that [Ne /v.pc] may originate closer to the black hole than other rest-frame forbidden optical emission lines (Dasyra et al. 2008, 2011; Netzer 2013; Negus et al. 2021, 2023). Dasyra et al. (2008) and Dasyra et al. (2011) find that the width of the mid-infrared [Ne /v.pc] line scales significantly, but with appreciable scatter, with M BH for local AGN. We use the [Ne /v.pc] intrinsic 𝜎 = 73 km s -1 and equation 1 in Dasyra et al. (2008) to approximate the black hole mass to be log(M BH / M ⊙ ) = 5 . 9 ± 0 . 8 (Table 6). This is consistent with the estimate from the M BH -𝜎 ∗ above.  \nThe final method to estimate M BH is to assume that accretion onto the black hole creates a hot accretion disk that produces the observed emission in GN 42437. If we can estimate the total luminosity, or bolometric luminosity ( 𝐿 bol ), of the accretion disk, we can  \n<!-- image -->  \nFigure 10. Left Panel: The black hole mass (M BH ) versus the velocity dispersion of the galaxy ( 𝜎 ). Right Panel: The M BH versus the stellar mass (M ∗ ) of the galaxy. We estimate M BH for GN 42437 in a variety of ways (Table 6) using: the M BH -𝜎 ∗ relation (light-blue circle), the [O /i.pc/i.pc/i.pc] 5008 Å attenuation-corrected luminosity (light-blue square), and the attenuation-corrected H 𝛽 luminosity inferred from the H 𝛼 emission line (light-blue diamond). The error bar on the blue diamond shows how M BH changes if the black hole accretes at 10% of the Eddington Luminosity (upper error bar) or at 3 times the Eddington Luminosity (lower error bar). Other 𝑧 > 5 black holes from the literature are shown as light-blue diamonds and lower-limits (Furtak et al. 2023; Kokorev et al. 2023; Larson et al. 2023; Maiolino et al. 2023; Übler et al. 2023; Lambrides et al. 2024). For the Maiolino et al. (2023) sample we do not include the merging black hole candidates. Gold points show local AGN (Reines & Volonteri 2015; Greene et al. 2020) with the solid gold line showing the fit to the low-redshift sample that includes the upper limits (Greene et al. 2020). The dashed gold line in the right panel is the local M BH -𝑀 ∗ relation from Reines & Volonteri (2015). The gray dashed line in the right panel is the relation from Pacucci et al. (2023) using the stellar mass function and observed high-redshift black holes to predict the 𝑧 = 4 -7 M BH -𝑀 ∗ relation. In the left panel, GN 42437 and high-redshift galaxies have had their H 𝛼 velocity dispersions corrected for the observed difference between H 𝛼 and stellar velocity dispersions (an increase of a factor of 1.4; Bezanson et al. 2018). GN 42437 has a large range of estimated M BH , but it is at the low-end of the observed distribution.  \n<!-- image -->  \nTable 6. Various ways to infer the black hole mass (M BH ) for GN 42437. The first column lists the method used to determine the mass, either using a line width ( 𝜎 ) or assuming the bolometric luminosity of the black hole ( 𝐿 bol is equal to the Eddington Luminosity, 𝐿 Edd ). The details are discussed in section 8. The second column gives the observable used, the third column gives the value of the observable, and the forth column shows how each quantity is transformed into the M BH constraint. The fifth column highlights whether a low or high ionization parameter is assumed (log 𝑈 of -3.5 or -2.5). The final column gives the approximated black hole mass.  \nassume that 𝐿 bol is near the Eddington Luminosity ( 𝐿 Edd ). 𝐿 Edd is the luminosity at which the outward force of radiation pressure due to electron scattering is equal to the inward pull of gravity. This is defined as  \n𝐿 Edd = 1 . 26 × 10 38 𝑀 BH erg s -1 , (1)  \nsuch that if we assume 𝐿 bol = 𝐿 Edd we can estimate M BH . There are two challenges here: (1) GN 42437 has significant light from the burst of star formation and (2) the high ionization parameter (log 𝑈 ) of GN 42437 makes it challenging to compare to lower redshift AGN that have lower log 𝑈 . For the first challenge, we use the finding in subsection 7.2 that the black hole contributes approximately one-  \nionizing photons to correct for the contribution from star-formation.  \nWe first estimate L bol using one-third of the reddening corrected [O/i.pc/i.pc/i.pc] 5008 Å luminosity (L [ O/i.pc/i.pc/i.pc ] ). Low-redshift L bol / 𝐿 [ O/i.pc/i.pc/i.pc ] conversion factors are typically ∼ 600 (Kauffmann & Heckman 2009; Liu et al. 2009). Using this estimate leads to a log(M BH / 𝑀 ⊙ ) = 6 . 9, which is allowed by the M dyn of the system but an order of magnitude above the local M BH -𝜎 ∗ relation.  \nThese L [ O/i.pc/i.pc/i.pc ] corrections are made for local AGN that have significantly lower ionization parameters than GN 42437. Netzer (2009) shows that the local bolometric corrections can be obtained with photoionization models and relatively low log 𝑈 parameters of -3 . 5. When we estimate 𝐿 bol using the Netzer (2009) low log 𝑈  \ncorrections of 891 we obtain log(M BH / 𝑀 ⊙ ) = 7 . 0 (see Table 6). The Netzer (2009) photoionization models also include a higher log 𝑈 = -2 . 5 model, which is more similar to the log 𝑈 inferred for GN 42437 above. At higher log 𝑈 a larger fraction of the total light is emitted as nebular emission lines, such that the L bol / 𝐿 [ O/i.pc/i.pc/i.pc ] is lower. With the extreme [O /i.pc/i.pc/i.pc] equivalent widths, this high log 𝑈 is more appropriate for GN 42437. The higher log 𝑈 models estimate log(M BH / 𝑀 ⊙ ) = 5 . 9 (light-blue square in Figure 10). The higher log 𝑈 bolometric correction is consistent with the estimate using the local M BH -𝜎 ∗ relation.  \nThe photoionization models of Netzer (2009) also suggest that the H 𝛽 luminosity is less sensitive to the log 𝑈 of the AGN. Dividing the H 𝛼 luminosity by a factor of three to account for the impact of the starburst, we estimate log(M BH / 𝑀 ⊙ ) = 6 . 2 and 6.3, for high and low log 𝑈 nebulae. The high-ionization H 𝛽 estimate is shown on Figure 10 as a light-blue diamond. This highlights the consistency of using H 𝛽 as a bolometric calibrator for high-ionization AGN. However, these photoionization models assume a solar metallicity and large dust obscuration. Both of these properties appear untrue for GN 42437. High-redshift Type /i.pc/i.pc AGN likely require new bolometric corrections to accurately estimate 𝐿 bol . This is beyond the scope of this paper.  \nFinally, we can use the upper-limit of the X-ray Luminosity L 𝑋 to put an upper-limit of the M BH to be less than log(M BH / 𝑀 ⊙ ) < 5 . 6 using the L bol / 𝐿 𝑋 relations from (Duras et al. 2020). The estimate from L 𝑋 is roughly consistent with the previous estimates. Other estimates from the rest-frame optical continuum emission or midinfrared colors could also be made, but these are challenging to separate the star-forming from black hole components in Type /i.pc/i.pc AGN.  \nWe also include an error bar on the H 𝛽 estimate in Figure 10 with the upper-bound corresponding to a black hole accreting at 10% L Edd and the lower bound for a black hole accreting at 3L Edd . This range on L Edd roughly covers the expected range of accretion rates for low-mass black holes, and suggest that the black hole mass range is likely log(M BH / 𝑀 ⊙ ) ∼ 5 -7. These approximated M BH satisfy the constraints from the observed M dyn , M ∗ , and estimated gas mass.  \nTheblack hole in GN 42437 is possibly an intermediate mass black hole and among the lowest masses observed at either high or lowredshift. Figure 10 shows some of the previously measured 𝑧 > 5 black holes as light-blue diamonds (Furtak et al. 2023; Goulding et al. 2023; Kokorev et al. 2023; Übler et al. 2023; Maiolino et al. 2023; Lambrides et al. 2024). The left panel of Figure 10 shows that the high-redshift black holes generally follow the local M BH -𝜎 ∗ relation (Maiolino et al. 2023). GN 42437 is no exception. The range of M BH estimates that do not use the M BH -𝜎 ∗ relation are generally within the scatter of the local M BH -𝜎 ∗ relation. This suggests that the fundamental local relation between M BH and the gravitational potential of the galaxy may be in place as soon as 𝑧 = 5 . 58.  \nHowever, the M BH for GN 42437 is nearly two orders of magnitude larger than the extrapolated local M BH -𝑀 ∗ relation (right panel of Figure 10; Reines & Volonteri 2015; Greene et al. 2020). This extreme offset is similar, but slightly less, than other recently observed 𝑧 > 5 black holes (Furtak et al. 2023; Goulding et al. 2023; Kokorev et al. 2023; Maiolino et al. 2023; Übler et al. 2023). It is possible that the first black holes discovered at 𝑧 > 5 with JWST may be an extremely biased sample. However, a recent reanalysis of the M BH -𝑀 ∗ relation using the high-redshift stellar mass function agrees with the inferred 𝑀 ∗ of GN 42437 and possibly indicates that this offset is physical (gray line in Figure 10; Pacucci et al. 2023). This hints that the build-up of stellar mass lags behind the growth of the black hole, possibly because black hole growth occurs before the stellar mass growth. This early black hole growth may be naturally  \nexplained by the formation of a heavy black hole seed in the early universe (see below; Bhowmick et al. 2024). Future consumption of cold gas, either from gas that is already in the galaxy (Maiolino et al. 2023) or that will be accreted onto the galaxy, will likely preferentially grow the stellar mass of GN 42437 rather than black hole. This may naturally evolve GN 42437 back towards the local M BH -𝑀 ∗ relation at later times.", '9 NARROW-LINE IMBHS IN THE EARLY UNIVERSE': "Above we have shown that the strong very high-ionization emission lines are consistent with ∼ 30% of the ionizing photons coming from an IMBH. GN 42437 is reminiscent of lower redshift narrow-line, or Type /i.pc/i.pc, AGN, where the black hole is revealed using emission line ratios rather than from broad emission lines. Importantly, GN 42437 does not appear to be an AGN from the classic low-redshift diagnostics (Figure 7). While previous studies have found significant populations of Type /i.pc super massive black holes at 𝑧 > 5 (Labbé et al. 2023; Kocevski et al. 2023; Matthee et al. 2023; Harikane et al. 2023; Larson et al. 2023; Kokorev et al. 2024), there may be an even larger population of undiagnosed Type /i.pc/i.pc AGN with very-high-ionization emission lines waiting to be discovered in the epoch of reionization (Gilli et al. 2022; Scholtz et al. 2023).  \nUnlike low-redshift Type /i.pc/i.pc AGN which are separated from starforming galaxies using metal emission lines that are 10s of percent the strength of H 𝛼 (e.g. [N /i.pc/i.pc]/H 𝛼 ), the very-high-ionization emission lines that reveal the Type /i.pc/i.pc AGN in GN 42437 are much fainter and muchnarrower.Thevery-high-ionization emission lines in GN 42437 are extreme compared to local star-forming galaxies, but [Ne /v.pc] only has a rest-frame equivalent width of 11 Å. The 15 hour integration times in the high-resolution G235H configuration was required to detect this faint line. Additionally, with line widths of only 73 km s -1 , this narrow line would not likely be detected with similar integration times at lower spectral resolution. Observations covering very-highionization forbidden (e.g., [Ne /v.pc]) and permitted FUV lines (He /i.pc/i.pc, N/v.pc 1238+1242 Å, C /i.pc/v.pc 1548+1550, for example) must be used to describe the full IMBH population at high-redshift (e.g., Feltre et al. 2016). Establishing these samples will reveal the population demographics of black holes in the first billion years of cosmic history and describe the formation and growth of the first black holes.  \nWeobserve GN 42437 only 1 Gyr after the Big Bang. What formed the IMBHsoearlyincosmictime?GN42437isaunique 𝑧 > 5galaxy and some of the observations presented here may provide clues for the origin of the IMBH. The H 𝛼 equivalent width suggests that in the last 3 Myr GN 42437 has undergone a significant burst of star-formation. This means that the galaxy has formed 9 × 10 7 M ⊙ , or 25% of its dynamical mass, in the past 3 Myr. As we discussed in section 8, the observed M ∗ plus a rough estimate of the gas mass is consistent with the estimated M dyn . For there to be an unobserved, older stellar population (Papovich et al. 2023), it must be within a factor of a few from the observed M ∗ , or we are significantly under-estimating the M dyn .  \nWhat caused such a burst of star formation? The NIRCam imaging of GN 42437 suggests that there is a potential companion a few kpc from GN 42437 (Figure 1). If this apparent companion is confirmed to be at 𝑧 = 5 . 58, it is possible that the interaction has triggered a significant amount of cold gas to collapse to form one of the first major generations of stars in GN 42437. It is possible that the interaction either formed or has strongly fed the IMBH. The strong feeding of the IMBH is likely why we detect the [Ne /v.pc] and are able  \nto characterize GN 42437 as an AGN. Without this strong burst we would likely characterize GN 42437 as a star-forming galaxy.  \nThere are broadly three main formation mechanisms of IMBHs: steady growth by accretion of a stellar mass black hole, a large gas mass directly collapsing into a black hole without forming a star, and runaway stellar mergers within a dense cluster leading to a super massive star (Greene et al. 2020, and references within). Stellar mass black holes (10-100 M ⊙ ) are unlikely to grow sufficiently fast with Eddington-like accretion to match the observations of GN 42437 (Haiman & Loeb 2001; Bromm & Larson 2004; Madau et al. 2014). The black hole would need 300 Myr of constant Eddington accretion to grow from 100 M ⊙ to 10 5 M ⊙ (Madau et al. 2014). While this is possible within the age of the Universe, it would require the black hole to accrete at Eddington for 30% of the age of the Universe without forming an observable mass of stars before the past 3 Myr. If there is sufficient gas to constantly feed the black hole at Eddington for 300 Myr it is likely that some of the gas would collapse into a population of stars that we would observe. Thus, the star formation history of GN 42437 suggests that it is unlikely that the IMBH grew from a stellar mass black hole seed.  \nThe ongoing burst of star-formation could have also directly collapsed a large amount of gas into the IMBH. Metal-free gas clouds do not fragment into smaller clouds, like higher metallicity clouds. This creates a monolithic super massive (10 5 M ⊙ ) gas cloud without any force to counteract gravity. The entire 10 5 M ⊙ then directly collapses into a single IMBH without producing significant amounts of other stars or metals (Loeb & Rasio 1994; Bromm & Loeb 2003; Begelman et al. 2006; Mayer et al. 2010). This scenario could explain the IMBHinGN42437,but we detect significant oxygen and neon emission lines. For the IMBH in GN 42437 to be formed through direct collapse, metals must have been synthesized by a now invisible and earlier massive star population. It is possible that this first generation of stars is no longer observable because it has already exploded as supernovae and injected the metals into the ISM. Observations with MIRI searching for an older stellar population could shed light on this scenario (Papovich et al. 2023). The direct collapse scenario is a possibility and does not contradict with any GN 42437 observations.  \nFinally, if the current burst of star-formation in GN 42437 formed a very dense star cluster, it is possible that many of these massive stars could collide to form a single super massive star ( > 10 3 M ⊙ ) that collapses to form an IMBH. This 'runaway merger' process may naturally explain the presence of the IMBH, the sharp burst of star formation, and the higher observed metallicity of GN 42437 (Portegies Zwart & McMillan 2002; Gieles et al. 2018). N-body simulations suggest that if there are a large number of initial stars (10 6 -7 ) then a single super massive star can grow to 10 5 M ⊙ in just 2 Myr (Gieles et al. 2018).  \nOnce a super massive star forms, the only support against collapse into an IMBH is hydrogen fusion in its core. Once the helium abundance in the core raises significantly, fusion ceases to counteract gravity and the super massive star collapses into an IMBH (Denissenkov & Hartwick 2014). The super massive star is fully convective and efficiently mixes the helium fused in its core with the photospheric hydrogen. The luminosity of super massive stars exceeds the Eddington luminosity by Thomson scattering, which drives a helium and neutron capture element enriched stellar wind (Denissenkov & Hartwick 2014). This enrichment process is actually quicker than enrichment by supernovae because it occurs during the life time of the massive star. These super massive stars produce significant amounts of helium and lighter metals (e.g., N and C), including oxygen (Nagele &Umeda2023). Just one 10 4 M ⊙ super massive star injects 133 M ⊙ of oxygen into the ISM in 10 5 years. This modest oxygen production  \nby a single super massive star could account for significant portion of the observed metal enrichment within GN 42437. Thus, the formation of the super massive star, the metal enrichment of the ISM, and the presence of an IMBH can occur in the ongoing ∼ 3 Myr burst of star formation.  \nThis runaway merger and eventual super massive star scenario has been proposed to explain the abundance patterns of globular clusters (Carretta et al. 2009; Piotto et al. 2015) and the nitrogen enhancement of GNz-11 (Charbonnel et al. 2023; Marques-Chaves et al. 2024). If the IMBH in GN 42437 did form through gravitational runaway and the formation of a super massive star, this may suggest that gravitational runaway and super massive stars crucially drive the chemical, stellar, and mass evolution of the early universe.  \nThe three formation mechanisms would likely lead to different relationships between M BH -𝜎 ∗ and M BH -𝑀 ∗ (Greene et al. 2020). Steady accretion near Eddington requires a nearly constant source of cold gas that would likely fragment into stars. If the black hole in GN 42437 grew due to smooth accretion, the M ∗ would likely grow in tandem with the M BH . This is inconsistent with the current 𝑧 > 5 observations where black holes preferentially grow before their stellar mass (right panel of Figure 10). Both gravitational runaway or direct collapse would likely lead to preferential black hole growth, and are therefore both consistent with the observed M BH -𝜎 ∗ and M BH -𝑀 ∗ relations (Bhowmick et al. 2024; Natarajan et al. 2024). Further observations of black hole masses in the early universe are required to reveal a full, unbiased population of early black hole to trace the origin of these trends.  \nUltimately we cannot definitively distinguish the exact formation mechanism of the IMBH in GN 42437. Future evidence is required to solidify our understanding of the formation and evolution of the IMBH in GN 42437. First, for a direct collapse scenario, we would likely expect MIRI observations to reveal the first generation of stars that likely enriched the ISM to the observed level (Papovich et al. 2023). For the runaway merger scenario, super massive star models predict elevated helium, carbon, and nitrogen abundances (Denissenkov & Hartwick 2014; Nagele & Umeda 2023). GN 42437 does present abnormally strong He /i.pc (3187 Å and 5875 Å) and He /i.pc/i.pc emission lines (Stephenson et al. in preparation). While we do not detect [N/i.pc/i.pc] 6585 Å, GN 42437 is very highly ionized such that we might expect the FUV N and C lines to be very strong (Bunker et al. 2023; Isobe et al. 2023; Marques-Chaves et al. 2024; Topping et al. 2024). Rest-frame FUV observations of the C /i.pc/i.pc/i.pc], C /i.pc/v.pc, N /i.pc/i.pc/i.pc], and N/i.pc/v.pc] could shine light on the chemical abundances of GN 42437. These abundances would test whether gravitational runaway or a direct collapse created the observed IMBH in GN 42437.", '10 CONCLUSIONS': 'Here we presented JWST NIRSpec high-resolution G235H and G395Hobservations of an extreme starburst galaxy at 𝑧 = 5 . 59 called GN42437.Therest-frame optical spectra have high rest-frame equivalent width emission lines (H 𝛼 and [O /i.pc/i.pc/i.pc] 5008 equivalent width of 901 ± 145 and 1644 ± 260 Å) that indicate that the galaxy is currently undergoing a significant starburst. We measure 13 nebular emission lines at > 3 𝜎 significance within the spectra (Table 3, Figure 3) and place strong constraints on the emission line ratios within the galaxy (Table 4). The ratios of strong emission lines in GN 42437 (such as [O/i.pc/i.pc/i.pc] 5008 Å/H 𝛽 and [N /i.pc/i.pc] 6585 Å/H 𝛼 ) suggest that GN 42437 has similar low- and high-ionization emission lines as both star-forming galaxies and Type /i.pc (broad-line) AGN at 𝑧 > 5 (Figure 7).  \nGN 42437 has remarkable very-high-ionization emission lines  \n(ionization potentials above 54 eV). [Ne /v.pc] 3427 Å (7 𝜎 significance; Figure 4 and Figure 5) and He /i.pc/i.pc 4687 Å (3 𝜎 significance; Figure 5) exceed the strength of typical low-redshift star-forming galaxies and are reminiscent of local AGN (Figure 6). The [Ne /v.pc] rest-frame equivalent width is 11 ± 2Å,the[Ne/v.pc] flux is 0 . 26 ± 0 . 04the [Ne /i.pc/i.pc/i.pc] 3870 Å flux, and the [Ne /v.pc]/H 𝛼 ratio is 0 . 044 ± 0 . 007.  \n[Ne /v.pc] and He/i.pc/i.pc require energies above 97 and 54 eV, respectively, and cannot be produced by massive star populations alone (Figure 6). Neither shocks nor X-ray binaries match the weak lowionization emission lines and very-high-ionization ratios (Figure 7 and Figure 9). Using the [Ne /v.pc]/[Ne /i.pc/i.pc/i.pc] versus He /i.pc/i.pc/H 𝛽 diagram in Figure 6 we suggest that a combination of accretion onto an intermediate mass black hole (IMBH) and massive stars can reproduce the entire ionization structure of GN 42437. The models with IMBHs and massive stars can reproduce all of the emission lines observed in GN 42437. These models suggest that a significant fraction of the ionizing photons must come from accretion onto a black hole.  \nThe line widths of GN 42437 are spectroscopically resolved with intrinsic line widths near 36 km s -1 for the strong H 𝛼 and [O/i.pc/i.pc/i.pc] 5008 Å (Table 5). We do not detect broad components from either the H 𝛼 or [O /i.pc/i.pc/i.pc] 5008 Å emission lines (Figure 3), nor do we detect low-ionization emission line ratios that traditionally indicate an AGN at lower redshifts and higher metallicities. The common AGN diagnostics fail to reveal the presence of the IMBH. GN 42437 resembles a narrow-line AGN selected by the presence of its very-highionization emission lines, unlike low-redshift narrow-line galaxies that are selected based upon their low-ionization emission lines (e.g. [N/i.pc/i.pc]/H 𝛼 ). High-ionization emission lines may reveal the full population of AGN at high-redshift.  \nWe estimate the black hole mass (M BH ) using a variety of different methods (Table 6) with the inferred log(M BH / 𝑀 ⊙ ) ∼ 5 -7. This range of M BH is broadly consistent with the local M BH -𝜎 ∗ using the spectroscopically resolved velocity dispersion (left panel of Figure 10), but the black hole is ∼ 2 orders of magnitude larger than expected from the stellar mass of the galaxy (right panel of Figure 10). This may suggest that high-redshift black holes preferentially grow before the stars in the galaxy.  \nWe discuss various formation mechanisms of the IMBH (section 9). While we cannot definitively determine the formation mechanism, we speculate that runaway mergers of massive stars or a direct collapse black hole could have produced the IMBH within GN 42437.  \nThe [Ne /v.pc] 3427 Å emission line has the power to diagnose the population demographics of IMBHs in the early universe. Deep high spectral resolution observations are required to determine the occurrence of IMBHs in the early universe and how they shaped the evolution of the first galaxies. Deep surveys using the G235H grating may enable a powerful census of IMBHs in the early universe that will describe how these early black holes grew into the super massive black holes we observe in the local universe.', 'ACKNOWLEDGEMENTS': 'This work is based on observations made with the NASA/ESA/CSA James Webb Space Telescope. The data were obtained from the Mikulski Archive for Space Telescopes at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127 for JWST . These observations are associated with program #01871. Support for program #01871 was provided by NASA through a grant from the Space Telescope Science Institute, which is operated by the  \nAssociation of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127.  \nWe thank Karl Gebhardt for insightful conversations on the evolution of IMBHs. CTR acknowledges support from the Elon University Japheth E. Rawls Professorship. SRF acknowledges support from NASA/FINESST grant number 80NSSC23K1433. NGG and YII acknowledge support from the National Academy of Sciences of Ukraine by its project no. 0123U102248 and from the Simons Foundation. ASL acknowledges support from Knut and Alice Wallenberg Foundation', 'DATA AVAILABILITY': 'All data presented here are publicly available. Reasonable requests for the data to the authors will be accommodated.', 'REFERENCES': 'Carretta E., et al., 2009, A&A, 505, 117  \nThis paper has been typeset from a T E X/L A T E X file prepared by the author.'}