A surprising finding of recent studies is the large number of Active Galactic Nuclei AGN associated with moderately massive black holes rm log Mbullet Modotsim 68 in the first billion years after the big bang z gt 5. In this context a relevant finding has been the large fraction of candidate dual AGN both at large separations several kpc and in close pairs less than a kpc likely in the process of merging. Frequent black hole merging may be a route for black hole growth in the early universe however previous findings are still tentative and indirect. We present JWSTNIRSpecIFU observations of a galaxy at z  7.15 in which we find evidence for a rm log Mbullet Modotsim 7.7 accreting black hole as traced by a broad component of H  emission associated with the Broad Line Region BLR around the black hole. This BLR is offset by 620 pc in projection from the centroid of strong restframe optical emission with a velocity offset of 40 km sSUP1SUP. The latter region is also characterized by narrow nebular emission features typical of AGN hence also likely hosting another accreting black hole although obscured Type 2 narrowline AGN. We exclude that the offset BLR is associated with Supernovae or massive stars and we interpret these results as two black holes in the process of merging. This finding may be relevant for estimates of the rate and properties of gravitationalwave signals from the early universe that will be detected by future observatories like LISA.
['Astrophysics - Astrophysics of Galaxies']
GANIFS JWST discovers an offset AGN 740 million years after the big bang
{'No Header': ', 1-11 (2023)', 'GA-NIFS: JWST discovers an offset AGN 740 million years after the Big Bang': "Hannah Übler, 1 , 2 ★ Roberto Maiolino, 1 , 2 Pablo G. Pérez-González, 3 Francesco D'Eugenio, 1 , 2 Michele Perna, 3 Mirko Curti, 4 Santiago Arribas, 3 Andrew Bunker, 5 Stefano Carniani, 6 Stéphane Charlot, 7 Bruno Rodríguez Del Pino, 3 William Baker, 1 , 2 Torsten Böker, 8 Giovanni Cresci, 9 James Dunlop, 10 Norman A. Grogin, 11 Gareth C. Jones, 5 Nimisha Kumari, 12 Isabella Lamperti, 3 Nicolas Laporte, 1 , 2 , 13 Madeline A. Marshall, 14 , 15 Giovanni Mazzolari, 16 , 17 , 1 Eleonora Parlanti, 7 Tim Rawle, 18 Jan Scholtz, 1 , 2 Giacomo Venturi, 7 Joris Witstok 1 , 2  \n- 1 Kavli Institute for Cosmology, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK\n- 2 Cavendish Laboratory, University of Cambridge, 19 JJ Thomson Avenue, Cambridge, CB3 0HA, UK\n- 3 Centro de Astrobiología (CAB), CSIC-INTA, Ctra. de Ajalvir km 4, Torrejón de Ardoz, E-28850, Madrid, Spain\n- 4 European Southern Observatory, Karl-Schwarzschild-Straße 2, 85748, Garching, Germany\n- 5 University of Oxford, Department of Physics, Denys Wilkinson Building, Keble Road, Oxford OX13RH, UK\n- 6 Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy\n- 7 Sorbonne Université, CNRS, UMR 7095, Institut d' Astrophysique de Paris, 98 bis bd Arago, 75014 Paris, France\n- 8 European Space Agency, c/o STScI, 3700 San Martin Drive, Baltimore, MD 21218, USA\n- 9 INAF - Osservatorio Astrofisco di Arcetri, largo E. Fermi 5, 50127 Firenze, Italy\n- 10 Institute for Astronomy, University of Edinburgh, Royal Observatory, Edinburgh EH9 3HJ, UK\n- 11 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA  \n12  \nAURA for the European Space Agency, Space Telescope Science Institute, Baltimore, MD, USA  \n- 13 Aix-Marseille Université, CNRS, CNES, LAM (Laboratoire d'Astrophysique de Marseille), UMR 7326, 13388 Marseille, France\n- 14 National Research Council of Canada, Herzberg Astronomy & Astrophysics Research Centre, 5071 West Saanich Road, Victoria, BC V9E 2E7, Canada\n- 15 ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), Australia\n- 16 Dipartimento di Fisica e Astronomia, Universit'a di Bologna, Via Gobetti 93/2, I-40129 Bologna, Italy\n- 17 INAF - Osservatorio di Astrofisica e Scienza dello Spazio di Bologna, Via Gobetti 93/3, I-40129 Bologna, Italy\n- 18 European Space Agency, ESAC, Villanueva de la Cañada, E-28692 Madrid, Spain  \nAccepted XXX. Received YYY; in original form ZZZ", 'ABSTRACT': 'A surprising finding of recent studies is the large number of Active Galactic Nuclei (AGN) associated with moderately massive black holes (log ( M · / M ⊙ ) ∼ 6 -8), in the first billion years after the Big Bang ( 𝑧 > 5). In this context, a relevant finding has been the large fraction of candidate dual AGN, both at large separations (several kpc) and in close pairs (less than a kpc), likely in the process of merging. Frequent black hole merging may be a route for black hole growth in the early Universe; however, previous findings are still tentative and indirect. We present JWST/NIRSpec-IFU observations of a galaxy at 𝑧 = 7 . 15 in which we find evidence for a log ( M · / M ⊙ ) ∼ 7 . 7 accreting black hole, as traced by a broad component of H β emission, associated with the Broad Line Region (BLR) around the black hole. This BLR is offset by 620 pc in projection from the centroid of strong rest-frame optical emission, with a velocity offset of ∼ 40 km/s. The latter region is also characterized by (narrow) nebular emission features typical of AGN, hence also likely hosting another accreting black hole, although obscured (type 2, narrow-line AGN). We exclude that the offset BLR is associated with Supernovae or massive stars, and we interpret these results as two black holes in the process of merging. This finding may be relevant for estimates of the rate and properties of gravitational wave signals from the early Universe that will be detected by future observatories like LISA.  \nKey words: galaxies: high-redshift - galaxies: active - galaxies: supermassive black holes', '1 INTRODUCTION': "★ hu215@cam.ac.uk  \nAfter one year of operations, one of the most surprising findings in the field of extragalactic astronomy produced by the James Webb Space Telescope (JWST; Gardner et al. 2023; Rigby et al. 2023) is the unexpectedly high abundance of active galactic nuclei (AGN)  \nduring the first billion years of cosmic history ( 𝑧 ≳ 5 . 5; Harikane et al. 2023; Maiolino et al. 2023b; Greene et al. 2023; Matthee et al. 2023; Kocevski et al. 2023; Yang et al. 2023). These are associated with accreting black holes and are identified through the properties of the spectral energy distribution (SED; e.g. Yang et al. 2023), through the presence of (narrow) nebular line emission typical of AGN (e.g. Übler et al. 2023; Maiolino et al. 2023a; Scholtz et al. 2023), and through the detection of a broad component of the Balmer lines, associated with the so-called Broad Line Region (BLR) of AGN (e.g. Kocevski et al. 2023; Übler et al. 2023; Harikane et al. 2023; Maiolino et al. 2023b; Greene et al. 2023; Matthee et al. 2023). With the strongest lines of the hydrogen Balmer series, H α and H β , being redshifted out of atmospheric windows for 𝑧 > 4 galaxies, many of these accreting black holes are now identified for the first time with JWST spectroscopy (primarily NIRSpec-MSA and NIRCam-slitless Jakobsen et al. 2022; Böker et al. 2023; Rieke et al. 2023a; Ferruit et al. 2022) - not least due to a remarkable lack of accompanying 𝑋 -ray emission, previously taken as the main indicator to constrain the cosmic black hole accretion rate density (Madau & Dickinson 2014; but see also Bogdán et al. 2023; Goulding et al. 2023).  \nThe JWST finding of the many AGN associated with moderately massive black holes has sparked the exploration of a large number of models, simulations and theories in order to explain their properties and abundances (e.g. Volonteri et al. 2023; Trinca et al. 2023; Schneider et al. 2023; Pacucci et al. 2023; Bennett et al. 2024; Zhang et al. 2023; Weller et al. 2023; Di Matteo et al. 2023a; Mayer et al. 2023). These models involve various potential early black hole seeding scenarios, including: heavy seeds resulting from the direct collapse of primordial clouds; remnants of early massive stars; rapid merging of stars and (stellar mass) black holes in dense nuclear clusters. Episodes of super-Eddington accretion have also been suggested (Schneider et al. 2023; Trinca et al. 2023; Bennett et al. 2024).  \nThe role of merging between intermediate mass black holes has also been extensively investigated by models and simulations (Barai et al. 2018; Barausse et al. 2020; Volonteri et al. 2020, 2022; Di Matteo et al. 2023a,b; Valentini et al. 2021; Vito et al. 2022; Di Mascia et al. 2021; Mannerkoski et al. 2022; Chen et al. 2023). This is not only a potential channel for the growth of black holes in the remote Universe, but also relevant for the detection of gravitational wave signals from early cosmic epochs with future gravitational waves experiments. Within this context, significant progress has been made observationally in recent years through the detection of dual quasars and AGN separated by a few/several kpc at high redshift (z ∼ 1-3) (Mannucci et al. 2022; Ciurlo et al. 2023; Perna et al. 2023a), in a number significantly larger than predicted by cosmological simulations (e.g. Rosas-Guevara et al. 2019; De Rosa et al. 2019; Volonteri et al. 2022; Chen et al. 2023). At even higher redshift (z > 4), JWST multi-object spectroscopic data have suggested the presence of a significant population of dual AGN, on scales smaller than 1 kpc, based on the complex profile of their BLR, suggesting the presence of two separate BLRs associated with two accreting black holes (Maiolino et al. 2023b). However, in the latter case no spatial information was available to confirm the presence of two AGN. Moreover, the complex profile of the broad lines may potentially also be associated with some peculiar properties of an individual BLR.  \nIn this paper, we report the detection of a broad H β line in the 𝑧 = 7 . 15 galaxy system ZS7, which identifies it as a (type 1) AGN. ZS7 is also known by its CANDELS ID (Grogin et al. 2011; Koekemoer et al. 2011; Nayyeri et al. 2017) COSMOS13679 (Pentericci et al. 2016, 2018), COSY-0237620370 (Roberts-Borsani et al. 2016; Laporte et al. 2017; Witten et al. 2023b), and COS-zs7-1 (Stark et al. 2017). This galaxy was already identified to host an AGN, based on the  \ndetection of N /v.pc] λ 1239andHe/i.pc/i.pc λ 1640in the rest-frame UV (Laporte et al. 2017). However, the peculiarity of this system is that the BLR is offset from the primary (narrow) nebular line emitter. Additionally, also the latter emitter displays nebular line signatures typical of an AGN. We interpret our results in the context of supermassive black hole mergers, which will be detectable with future gravitational waves observatories like LISA.  \nThroughout this work, we use the AB magnitude system and assume a flat Λ CDM cosmology with Ω 𝑚 = 0 . 315 and 𝐻 0 = 67 . 4 km/s/Mpc (Planck Collaboration et al. 2020). With this cosmology, 1 '' corresponds to a transverse distance of 5.27 proper kpc at 𝑧 = 7 . 15.", '2 NIRSPEC-IFU OBSERVATIONS AND DATA PROCESSING': "ZS7 was observed with the NIRSpec instrument (Jakobsen et al. 2022) in Integral Field Unit (IFU) mode (Böker et al. 2022) on board JWST as part of the NIRSpec-IFU GTO programme 'Galaxy Assembly with NIRSpec IFS' (GA-NIFS) under programme 1217 (PI: Nora Lützgendorf). The NIRSpec data were taken on May 9, 2023, with a medium cycling pattern of eight dithers and a total integration time of 4.1 h with the medium-resolution grating/filter pair G395M/F290LP, covering the wavelength range 2 . 87 -5 . 27 𝜇 m (spectral resolution 𝑅 ∼ 700 -1300; Jakobsen et al. 2022). It was also imaged with the NIRCam instrument (Rieke et al. 2023b) in eight filters though the PRIMER programme (PID: 1837, PI: James Dunlop).  \nRaw data files were downloaded from the Barbara A. Mikulski Archive for Space Telescopes (MAST) and subsequently processed with the JWST Science Calibration pipeline 1 version 1.8.2 under the Calibration Reference Data System (CRDS) context jwst\\_1068.pmap. We made several modifications to the default reduction steps to increase data quality, which are described in detail by Perna et al. (2023a) and which we briefly summarise here. Countrate frames were corrected for 1 / 𝑓 noise through a polynomial fit. During calibration in Stage 2, we removed regions affected by failed open MSA shutters. We also removed regions with strong cosmic ray residuals in several exposures. Remaining outliers were flagged in individual exposures using an algorithm similar to /l.pc/a.pc/c.pc/o.pc/s.pc/m.pc/i.pc/c.pc (van Dokkum 2001): we calculated the derivative of the count-rate maps along the dispersion direction, normalised it by the local flux (or by three times the rms noise, whichever was highest), and rejected the 96 th percentile of the resulting distribution (see D'Eugenio et al. 2023, for details). The final cube was combined using the 'drizzle' method, for which we used an official patch to correct a known bug. 2 The main analysis in this paper is based on the combined cube with a pixel scale of 0 . 05 '' . We used spaxels away from the central source and free of emission features to perform a background subtraction.", '3.1 Spectral fitting and detection of broad H β': "Fig.1 shows two spectra (zoomed around H β ) extracted from the IFU cube at the location of the [O /i.pc/i.pc/i.pc] λ 5007 centroid and at a location  \nFigure 1. Top left: Color-composite image of ZS7 and its immediate environment created using public NIRCam data from the PRIMER program (PI: James Dunlop). The approximate NIRCam PSF is indicated by the white filled (F200W, ∼ 0 . 06 '' ) and open (F444W, ∼ 0 . 14 '' ) circles. The effective NIRSpec-IFU field of view (FOV) is indicated by the dashed white box. Top right: Zoom-in on the environment of ZS7, as indicated by the rectangular solid white box in the top left panel. The color scale shows a line map integrated over 90 Å around the [O /i.pc/i.pc/i.pc] λ 5007 position, where we indicate the positions of two faint companions detected spectrally at the same redshift with dashed grey circles. The centroids of the [O /i.pc/i.pc/i.pc] λ 5007 emission and the BLR emission are indicated by a white plus and a golden star, respectively, where the small white circles on top of the symbols indicate five times the fitting errors on the centroid positions. The approximate NIRSpec PSF at 4 𝜇 m is indicated by the white circle. The centroid of the [O /i.pc/i.pc/i.pc] λ 5007 emission coincides with the NIRCam F200W emission, while the BLR location coincides with the F444W emission, reflected in the color gradient in the NIRCam false-color image. Bottom: Zoom-in on the spectra extracted from a single spaxel at the [O /i.pc/i.pc/i.pc] λ 5007 centroid position (left) and at the BLR location (right). We show our full fit in blue, with narrow emission lines in green and the BLR component in purple. At the BLR location, the broad H β emission is evident, while the [O /i.pc/i.pc/i.pc] λλ 4959,5007 emission is narrow at all locations. Some flux of the BLR point-source is spread through the point-spread function (PSF FWHM ∼ 0 . 13 -0 . 14 '' at 4 𝜇 m) and therefore is still visible at the [O /i.pc/i.pc/i.pc] λ 5007 centroid position.  \n<!-- image -->  \nabout 0.12 '' to the NW, where a broad component of H β is clearly visible. We fit the spectra with a combination of Gaussian functions for the line emission, and a power-law continuum. We include a narrow component for each of the following lines: [O /i.pc/i.pc] λ 3727, [O/i.pc/i.pc] λ 3729, [Ne /i.pc/i.pc/i.pc] λ 3869, He I 𝜆 3889+H8, [Ne III] 𝜆 3967, H ε , H δ , H γ , [O /i.pc/i.pc/i.pc] λ 4363, H β , [O /i.pc/i.pc/i.pc] λ 4959, [O /i.pc/i.pc/i.pc] λ 5007. The intensity ratio of the [O III] doublet is fixed to [O /i.pc/i.pc/i.pc] λ 5007/[O /i.pc/i.pc/i.pc] λ 4959 = 2.98, based on atomic physics (Osterbrock & Ferland 2006). We note that the [O II] doublet as well as [Ne III] 𝜆 3967 and H ε are blended. We show our full fits to spectra extracted from the BLR location and the centroid of the [O /i.pc/i.pc/i.pc] λ 5007 emission in Appendix A.  \nAs mentioned, for H β , a broad component is observed in some spaxels. We see no evidence for broad components in other lines. As it will be discussed later on, we interpret this broad H β component as a BLR, offset relative to the [O /i.pc/i.pc/i.pc] λ 5007 centroid.  \nFollowing Übler et al. (2023), we re-scale the formal noise based on the 'ERR' extension with a measurement of the standard deviation  \nin regions free of line emission to take into account correlations due to the non-negligible size of the PSF relative to the spaxel size. This increases the noise by a factor of up to two for individual spaxels.  \nWe first fit the integrated spectrum extracted from 3 × 3 spaxels centred on the BLR location to robustly constrain the shape of the BLR emission. Measuring the BLR FWHM on a larger aperture including all spaxel with 𝑆 / 𝑁 H 𝛽, BLR ≥ 5(see Fig. 2) gives consistent results within the uncertainties. Here and elsewhere, the 𝑆 / 𝑁 is measureddirectly from the best fit. Then, we fit the emission extracted from individual spaxels by fixing the BLR FWHM to the previous best-fit value. We list the resulting emission line fluxes extracted at the BLR location and [O /i.pc/i.pc/i.pc] λ 5007 centroid position in Table 1. Zoom-ins on the corresponding spectral fits extracted from individual spaxels are shown in Fig. 1.  \nIf we integrate over a larger aperture centred on the [O /i.pc/i.pc/i.pc] λ 5007 emission in ZS7, we can detect additional fainter lines. A circular aperture with radius 0 . 3 '' (encompassing both the [O /i.pc/i.pc/i.pc] λ 5007  \nemitter and the BLR, similar to the region shown in Fig. 2) reveals tentative line emission at the expected wavelengths of H11, H10, H9, and a detection of He I 𝜆 5876. While He /i.pc/i.pc λ 1640 has been detected in previous observations from the ground (Laporte et al. 2017), we do not detect He /i.pc/i.pc λ 4686 in ZS7. This is not necessarily surprising, because He /i.pc/i.pc λ 1640 is the brighter line. If we extract a spectrum collapsed from a pseudo-slit corresponding to the previous X-Shooter observations by Laporte et al. (2017), we find an upper limit on the He/i.pc/i.pc λ 4686 flux which is consistent within the uncertainties with their measurement for He /i.pc/i.pc λ 1640 within a factor of 8-13, consistent with theoretical expectations.", '3.2 Maps and velocity field': 'We use the fluxes resulting from full spectral fitting to obtain emission line maps, for the regions where the lines are detected at more than 5 𝜎 . In Fig.2 we show specifically the flux maps of [O /i.pc/i.pc/i.pc] λ 5007 (which is also shown in Fig.1), the broad H β component, and [O/i.pc/i.pc/i.pc] λ 4363.  \nThe map of the broad component of H β is consistent with the PSF, i.e. a spatially unresolved source. The maps confirm what was already visible from Fig.1, namely the H β broad component is clearly offset from the [O /i.pc/i.pc/i.pc] λ 5007 line centroid.  \nThe auroral line [O /i.pc/i.pc/i.pc] λ 4363 peaks at the location of the H β broad, not unexpectedly, as strong [O /i.pc/i.pc/i.pc] λ 4363 emission was already pointed out to be likely associated with AGN activity (Dors et al. 2020; Brinchmann 2023; Perna et al. 2017). However it is interesting to see that it also extend to the [O /i.pc/i.pc/i.pc] λ 5007 peak; we will discuss that the strength of the [O /i.pc/i.pc/i.pc] λ 4363 line also at this location suggest the presence of an additional (type 2) AGN.  \nWe also derive the [O /i.pc/i.pc/i.pc] λ 5007 velocity map, as shown in Fig. 2. We use the same fitting approach as described in Section 3.1,but fit only the [O /i.pc/i.pc/i.pc] λλ 4959,5007 emission which has the highest signalto-noise ratio. The velocity is derived from the line centroid after subtracting the systemic velocity measured at the [O /i.pc/i.pc/i.pc] λ 5007 centroid position. From our fits, we find an average uncertainty of 6 km/s on the velocities shown in Fig. 2. The region around the BLR location appears clearly separated in velocity space, as indicated by the cluster of positive relative velocities. The absolute difference in velocity to the centroid of the [O /i.pc/i.pc/i.pc] λ 5007 emitter is only about 20 km/s, but significant. Being at the the same redshift, the small projected distance and the small velocity difference suggest that we are witnessing the imminent merger of two systems, one of which hosts a BLR AGN.  \nFinally, we also extract and map the continuum in the bluest (2.93 𝜇 m) and reddest (4.9-5 𝜇 m) parts of the spectrum. As illustrated in Fig. 2, the red continuum peaks at the location of the broad H β , while the blue continuum peaks on the [O /i.pc/i.pc/i.pc] λ 5007 centroid. This essentially confirms what was already seen in the NIRCam colors shown in Fig. 1. This further indicates that we are likely observing two distinct objects in close proximity. The blue continuum centered on the [O /i.pc/i.pc/i.pc] λ 5007 centroid is likely tracing a young stellar population. Instead, it is not straight-forward to interpret the nature of the compact red source at the location of the broad H β : it could be a reddened or aged stellar population, but it may also be associated with the reddened continuum of the AGN accretion disc. It is not possible to easily distinguish between these two scenarios, because the continuum spectral information is limited to 3600 Å < 𝜆 rest < 5200 Å and has modest S/N.  \nTable 1. Observed narrow emission line fluxes extracted from one spaxel in units of 10 -20 erg/s/cm 2 from full spectral fits to the BLR location and [O/i.pc/i.pc/i.pc] λ 5007 centroid position and derived quantities.', '3.3 Stellar mass estimate for the [O /i.pc/i.pc/i.pc ] λ 5007 emitter': "Weestimate the stellar mass of the [O /i.pc/i.pc/i.pc] λ 5007 emitter through a full SED modelling using imaging data from the PRIMER programme on the pixel level, following Pérez-González et al. (2023) (see also García-Argumánez et al. 2023; D'Eugenio et al. 2023). SEDs were constructed for each one of the 262 30 mas (150 pc, with the PSF FWHM of the reddest image being around 5 pixels) pixels of the system counting with more then three NIRCam filters with 𝑆 / 𝑁 > 5. The NIRCam data covers the wavelength range 0.9-4.4 𝜇 mwith eight filters (one medium-band filter, the rest being broad-band). The SEDs were fitted to stellar population models by Bruzual & Charlot (2003), assuming a delayed exponential for the star formation history (SFH) with timescale 𝜏 taking values from 1 Myr (i.e., an instantaneous burst) to 10 Gyr (i.e., a constant SFH). The attenuation was treated with the Calzetti et al. (2000) law, with A(V) values from 0 mag to 5 mag. The metallicity was left to vary between 0 . 005 𝑍 ⊙ and 𝑍 ⊙ . The IMF was set to Chabrier (2003). With a total of four free parameters, the best fitting model (using a Monte Carlo method to estimate errors) is scaled to the photometry to get an estimation of the (surviving) stellar mass.  \nTo estimate the stellar mass of the [O /i.pc/i.pc/i.pc] λ 5007 emitter, we exclude the Northern region where the BLR is detected with 𝑆 / 𝑁 > 5 to mitigate contamination from the AGN (see Fig. 2). This leads to a total stellar mass of 𝑀 ★ ∼ 2 . 9 ± 0 . 9 × 10 9 M ⊙ . If we fill in the excluded region with values from the Southern part of the galaxy at similar distance from the [O /i.pc/i.pc/i.pc] λ 5007 centroid, we find 𝑀 ★ ∼ 3 . 8 ± 1 . 1 × 10 9 M ⊙ . This estimate is consistent with what has been found in previous work from SED fitting to the CANDELS and also IRAC data (Pentericci et al. 2016; Laporte et al. 2017). This shows that the [O /i.pc/i.pc/i.pc] λ 5007 emitter is a moderately massive galaxy for its redshift (see e.g. Curti et al. 2023; Nakajima et al. 2023), expected to host an intermediate-mass to massive black hole in its centre.  \nIn Section 7 we will argue that the narrow line emission associated with [O /i.pc/i.pc/i.pc] λ 4363 may also trace recent or ongoing AGN activity. Therefore, we make another estimate of the stellar mass of the [O/i.pc/i.pc/i.pc] λ 5007 emitter by excluding not only the region where the H 𝛽 BLR is detected with 𝑆 / 𝑁 > 5, but also the region where the narrow [O/i.pc/i.pc/i.pc] λ 4363 is detected with 𝑆 / 𝑁 > 5. This would lead to a stellar mass of 𝑀 ★ ∼ 1 . 7 ± 0 . 5 × 10 9 M ⊙ (without filling in the excluded regions with other values). We consider this a lower limit on the total stellar mass of the [O /i.pc/i.pc/i.pc] λ 5007 emitter, and note that a reduction of the stellar mass by a factor 2-3 does not impact the results and conclusion presented in the following sections.  \n<!-- image -->  \n<!-- image -->  \nFigure 2. Top left: continuum flux density map at 𝜆 = 2 . 9 -3 𝜇 m. Bottom left: continuum flux density map at 𝜆 = 4 . 9 -5 𝜇 m. Top centre, bottom centre, and bottom right: integrated line flux maps with 𝑆 / 𝑁 ≥ 5 from our spectral fitting for the narrow [O /i.pc/i.pc/i.pc] λ 5007 component, the H β BLR component, and the narrow [O/i.pc/i.pc/i.pc] λ 4363 component, respectively. Top right: velocity map of the narrow [O /i.pc/i.pc/i.pc] λλ 4959,5007 emission. In each panel, the black cross and star indicate the centroid of the [O /i.pc/i.pc/i.pc] λ 5007 emission and the BLR location as shown in Fig. 1, respectively. The black scale bar to the left indicates 0 . 2 '' ∼ 1 kpc, and the grey outline indicates the mask used for our full spectral fitting. There is a clear shift in the location of the continuum emission from bluer to redder wavelengths, aligning with the centroids of the [O /i.pc/i.pc/i.pc] λ 5007 emission and the BLR position, reflecting the color gradient seen in the NIRCam image (Fig. 1). The flux of the auroral [O /i.pc/i.pc/i.pc] λ 4363 line peaks at the BLR position, but covers also the [O /i.pc/i.pc/i.pc] λ 5007 flux peak, indicating a clear gradient of the [O /i.pc/i.pc/i.pc] λ 5007/[O /i.pc/i.pc/i.pc] λ 4363 line ratio from the position of the AGN to the [O /i.pc/i.pc/i.pc] λ 5007 emitter.  \n<!-- image -->", '4 EVIDENCE FOR AN OFFSET TYPE 1 AGN AND ENVIRONMENT OF ZS7': "At a projected distance of 620 pc from the centroid of the [O /i.pc/i.pc/i.pc] λ 5007 emission in ZS7, we detect a broad component at the wavelength of H β (Fig.1). As discussed, the emission is consistent with stemming from a point source, spread by the point-spread function (PSF; PSF FWHM ∼ 0 . 13 -0 . 14 '' at 4 𝜇 m; see Fig. 2). The broad component does not have a counterpart in the forbidden lines, in particular not in the strong [O /i.pc/i.pc/i.pc] λλ 4959,5007 doublet, hence it cannot be associated with emission of high velocity ISM or CGM, for instance due to galactic outflows. Therefore, the most plausible explanation is that the broad H β is tracing the Broad Line Region of an accreting black hole, i.e. a type 1 AGN (we discuss and discard other interpretations in Section 6).  \nThe finding of an AGN in ZS7 is not so surprising, as evidence for an AGN was already found in the past based on the detection of N/v.pc] λ 1239 and He /i.pc/i.pc λ 1640 in the rest-frame UV (Laporte et al. 2017). The novelty of our finding is that the type 1 AGN is offset relative to the [O /i.pc/i.pc/i.pc] λ 5007 emitter. Also based on the kinematic offset around the BLR (Fig.S2 and Sect.S3.2) and the markedly different colors, these results suggests that the type 1 AGN traced by the broad H β is hosted in a system separate from the galaxy traced by the bulk of the [O /i.pc/i.pc/i.pc] λ 5007 emission, and in the process of merging with it. The ongoing merger may be channelling gas to the 5 × 10 7 M ⊙ black hole as a consequence of tidal interactions (Di Matteo et al. 2005; Johansson et al. 2009; see also e.g. Dougherty et al. 2024).  \nWe consider the possibility that the BLR is actually at the centre of a single system but partially obscured by dust, while the [O/i.pc/i.pc/i.pc] λ 5007 brighter emission region is tracing the less obscured part of the galaxy (or a less obscured part of an extended Narrow Line Region). However, the dust extinction towards nebular lines as traced by the H β /H γ Balmer decrement does not show a significant gradient between the two locations and, if anything, slightly lower  \nreddening towards the narrow nebular lines at the BLR location (see Appendix B and Fig. B1). In addition, the fact the [O /i.pc/i.pc/i.pc] λ 4363 line, which is clearly resolved across the two locations, peaks at the BLR location, provides further evidence that the BLR location is not the dust-obscured centre of the [O /i.pc/i.pc/i.pc] λ 5007 emitter. In the case of differential reddening, one would expect that [O /i.pc/i.pc/i.pc] λ 4363 would peak at the same location as [O /i.pc/i.pc/i.pc] λ 5007, whereas it peaks at the BLR location.  \nIt is also implausible that the BLR is tracing the centre of a single system and that the [O /i.pc/i.pc/i.pc] λ 5007 centroid is tracing a burst of star formation in its disc, at ∼ 600 pc from the nucleus. Indeed, the [O/i.pc/i.pc/i.pc] λ 5007 emission is not an unresolved clump of star formation, it is a well resolved structure extending over ∼ 1 . 5 kpc. Additionally, if the BLR was the centre of a single system then the velocity field should be approximately symmetric relative to it; on the contrary, the velocity field around the BLR stands out for having a clear offset relative to the surrounding region, just as expected for a separate system. Furthermore, there is no evidence of significant extended line or continuum emission North-West of the BLR. Finally, we will show that also the [O /i.pc/i.pc/i.pc] λ 5007 centroid is likely hosting another (type 2) AGN, indicating that it is associated with a separate nucleus.  \nAt the same redshift as ZS7, we detect three fainter companions through emission lines of H β and [O /i.pc/i.pc/i.pc] λλ 4959,5007: two to the North at projected separations of 2.5 and 4.1 kpc, and one to the South-East at a projected distance of 3.9 kpc. Only the closest satellite is visible in the imaging data (see Figure 1), and has been previously photometrically identified by Witten et al. (2023a). The high number of companions, in addition to the offset BLR, suggests that we are witnessing a snapshot of the complex assembly of a larger system just 740 million years after the Big Bang.  \nTable 2. Measurements of the broad H β component extracted from an integrated spectrum with 𝑆 / 𝑁 H β , BLR > 5 (see bottom centre panel of Fig. 2), and corresponding black hole properties.", '5 BLACK HOLE MASS AND LUMINOSITY OF THE OFFSET TYPE 1 AGN': 'By assuming standard virial relations, we can estimate the black hole mass using our best-fit values of the H β BLR flux and FWHM, based on the 20-spaxel integrated spectrum around the BLR location, where 𝑆 / 𝑁 H β , BLR ≥ 5. The FWHM is corrected for the instrumental spectral resolution which is approximately 𝑅 ∼ 1000 at the observed wavelength of H β , corresponding to 𝜎 ∼ 127 km/s (Jakobsen et al. 2022). Using the single-epoch calibration by Greene & Ho (2005) derived from data of nearby AGN, we find a black hole mass of log ( M · / M ⊙ ) = 7 . 7 ± 0 . 4. The uncertainties are dominated by the scatter of local scaling relations, which we take into account by adding 0.4 dex in quadrature to our error budget (e.g. Ho & Kim 2014). We also estimate the AGN bolometric luminosity to be 𝐿 bol = 8 . 5 (± 0 . 8 ) × 10 44 erg/s based on the H β BLR flux following Dalla Bontà et al. (2020), indicating that the black hole is accreting at about 14 per cent of its Eddington limit. We report these measurements in Table 2.', '6 EXCLUDING ALTERNATIVE SCENARIOS FOR THE BROAD H β : SUPERNOVAE, OUTFLOW AND MASSIVE STARS': "We briefly discuss alternative interpretations of the broad H β component: a local outflow, a supernova (SN) in the galaxy disc, or a cluster of Wolf-Rayet or very massive stars.  \nWhile we recall that outflows as fast as a few 1000 km/s are only associated with AGN feedback (e.g. Fabian 2012; Heckman & Best 2014; Carniani et al. 2015; Förster Schreiber et al. 2019; Veilleux et al. 2020; see also discussion by Maiolino et al. 2023a), the outflow scenario should result in even more prominent broad emission in [O/i.pc/i.pc/i.pc] λλ 4959,5007, as observed in galactic outflows locally and at high redshift, and especially because outflows are typically metal enriched with high [O /i.pc/i.pc/i.pc] λ 5007/H β (Origlia et al. 2004). 3  \nCore collapse SNe do show broad Balmer line emission; at their peak luminosity their broad H β luminosity is of the order of ∼ 10 40 erg/s, although some extremely luminous SNe can reach broad H β luminosities 4 of a few times ∼ 10 41 erg/s (Pastorello et al. 2002; Tartaglia et al. 2016; Kokubo et al. 2019; Dickinson et al. 2023; Fransson et al. 2014), still lower than what is observed in ZS7 ( 𝐿 H 𝛽, broad ∼ 10 42 erg/s; Table 2). The additional important  \ndifference is that when SNe reach such high luminosities, the broad component of H α and H β is generally much larger than observed in ZS7, typically several 1000 km/s and even exceeding 10,000 km/s. Moreover, core-collapse SNe, especially the luminous ones, have an asymmetric blue-shifted profile (from a few 100 km/s to a few 1000 km/s) and/or P-Cygni profiles (e.g. Gal-Yam 2012; Taddia et al. 2013; Gutiérrez et al. 2017), while the broad H β component in ZS7 is symmetric (Fig. 1). Since we detect the continuum, we can also rule out strong absorption lines. We also note that the spectrum of corecollapse SNe (especially in the early phases) is further generally characterised by prominent [O /i.pc] λ 6300 emission and/or Fe bumps near H β , which are not seen in ZS7.  \nAdditional compelling evidence against the SN scenario comes from the fact that in this case the source at the location of the broad H β should display some flux variation. We can leverage the fact that the NIRCam and NIRSpec-IFU observations were obtained at two different epochs separated by 5 months, i.e. 20 days rest-frame. While this is not a long time lag, depending on the phase of the putative SN light curve, we would expect a flux variation ranging from a factor of several to at least 20%. We have consistently extracted the flux from the same aperture both from the NIRCam and NIRSpec data. For NIRCam we have used the filters overlapping with our NIRSpec spectrum, i.e. F444W, F410M and F356W, and in NIRSpec we have extracted the spectra matching the same wavelength range as the three filters. At least a 5% cross-calibration uncertainty is estimated due to the fact that the PSF of the two instruments are not identical (i.e. slight differential losses, despite the large adopted aperture), as well as slight background disuniformities in the IFU data. Nonetheless we obtain that the flux in the three bands are fully consistent within 3%, 2% and 0.07%, in the three band respectively. Such a good match, and the implied constancy of the flux (both of the continuum and of the broad H β , which is included in all three filters), makes the single SN scenario simply impossible.  \nHowever, one can explore whether the broad H β is not due to a single SN, but the overlap of multiple, co-located SNe associated with star formation within the small region associated with the BLR. However, even assuming the extreme case of no AGN contribution to the narrow H β , this would give an upper limit on the star formation rate of < 6 𝑀 ⊙ / yr, i.e. < 10 -1 𝑆𝑁 / yr; even in the extreme case that the peak of the broad line emission of each SN lasts for one year, this would imply an average broad H β luminosity of 10 39 -40 erg/s, much lower than observed in ZS7. Also, one would expect the peak of the core collapse SNe to pile up at the location of the peak of star formation, i.e. at the [O /i.pc/i.pc/i.pc] λ 5007 centroid.  \nAlthough it is not a common feature, broad Balmer lines have been associated with clusters of Wolf-Rayet stars in some individual cases (e.g. in NGC 4214, see Sargent & Filippenko 1991). More distinctive features of Wolf-Rayet stars are for instance the 'blue bump' around 4600-4720 Å, associated with broad He II 𝜆 4686, N /i.pc/i.pc/i.pc 𝜆 4640 and/or C /i.pc/i.pc/i.pc//i.pc/v.pc 𝜆 4650, and the 'red bump' around 5800 Å, associated with C /i.pc/v.pc 𝜆 5808 (e.g. Sargent & Filippenko 1991; Terlevich et al. 1996; Schaerer et al. 1997, 1999; Brinchmann et al. 2008; MonrealIbero et al. 2017). These spectral signatures are not seen in ZS7. We calculate an upper limit on the flux in the 'blue bump' region, and on He II 𝜆 4686, from our aperture where the H 𝛽 BLR is detected with 𝑆 / 𝑁 > 5. From this, we find 'blue bump'/H 𝛽 < 0 . 033 and He II 𝜆 4686/H 𝛽 < 0 . 006. For comparison, Sargent & Filippenko (1991) find 'blue bump'/H 𝛽 ∼ 0 . 91 in the region where they also detect broad H 𝛼 , and values of He II 𝜆 4686/H 𝛽 ∼ 0 . 01 -0 . 1 are detected in nearby Wolf-Rayet galaxies (e.g. Vacca & Conti 1992; Thuan et al. 1996; Fernandes et al. 2004). To further constrain from the H 𝛽 flux and the upper limit on He II 𝜆 4686 the possible contribu-  \n7, we follow Equation 15 by Vacca & Conti (1992) to calculate the ratio of the number of WNL stars, the dominant sub-type of Wolf-Rayet stars, and the 'equivalent' number of O stars, and find 0.019. Considering the detection of He II 𝜆 1640 in ZS7, the two independent observations discussed by Laporte et al. (2017) reveal line widths that are consistent with the spectral resolution of the respective instruments, while again generally broad emission would be expected for WolfRayet stars (e.g. Chisholm et al. 2019; Leitherer et al. 2019; Vanzella et al. 2020). In summary, while we cannot exclude the presence of Wolf-Rayet stars, we can conclude that their potential contribution to any spectral features in ZS7 is minor. This specifically holds for the broad H 𝛽 emission.  \nVery massive stars are also known to show broad Balmer lines, but in these cases the Equivalent Width of the broad component of H β is only a few Å (Martins et al. 2020), while in ZS7 the Equivalent Width of the broad component of H β is about 60 Å. Furthermore, in the case of very massive stars the broad H β is also accompanied by prominent and broad He /i.pc/i.pc λ 4686 (Martins & Palacios 2022), not seen in ZS7. Finally, we note that an outflow driven by a large numbers of massive stars has been invoked as a possible origin of broad offnuclear H α emission without a counterpart in [O /i.pc/i.pc/i.pc] λλ 4959,5007 in a nearby galaxy (Rodríguez Del Pino et al. 2019). However, in the latter case the width of the broad H α is only ∼ 800 -1000 km/s. Moreover, in the latter scenario, as well as all previous scenarios in which the broad H β is associated with clusters of young stars, one would expect the broad H β to peak at the location of the peak of star formation, i.e. at the location of the peak of (narrow) nebular emission, and not offset as in the case of ZS7.  \nGiven the several arguments outlined above, and also the fact that previous studies already found evidence for an AGN in ZS7 (Laporte et al. 2017), we conclude that the most plausible explanation is that the broad H β component is tracing the BLR of a massive, accreting black hole. Furthermore, the H β BLR and the [O /i.pc/i.pc/i.pc] λ 5007 distribution are likely associated with two different systems in the process of merging. This was already discussed in Section 4, and will be further discussed in the next two sections.", '7 INDICATION OF A SECOND, TYPE 2 AGN IN THE [O /i.pc/i.pc/i.pc ] λ 5007 EMITTER': 'As discussed by various authors, the detection of type 2 AGN, i.e. for which the BLR is obscured along our line of sight, and only the Narrow Line Region (NLR) is visible, is more challenging at high redshift. Indeed, the narrow line ratios of AGN at 𝑧 > 4 are displaced from the local AGN locus in classical optical diagnostic diagrams, such as the BPT (Baldwin et al. 1981), involving [O /i.pc/i.pc/i.pc] λ 5007/H β vs [NII]6584/H α , and the diagrams proposed by Veilleux & Osterbrock (1987), in which [NII]6584 is replaced by the [SII]6720 doublet (Übler et al. 2023; Kocevski et al. 2023; Harikane et al. 2023; Maiolino et al. 2023b; Scholtz et al. 2023). This is likely a consequence of the gas in the NLR becoming metal-poor, which results in an overlap of AGN with the sequence of star-forming galaxies on these same diagrams. Some of the alternative diagnostic diagrams to identify AGN based on UV lines proposed by, for instance, Scholtz et al. (2023), are not usable in our case, due the limited spectral coverage of our data.  \nHowever, an additional diagnostic can be provided by the strong [O/i.pc/i.pc/i.pc] λ 4363 auroral line, which is found to peak at the same location of the broad H β (Fig. 2), with an intensity comparable to H γ . Such high [O /i.pc/i.pc/i.pc] λ 4363 is extremely rare and, in combination with  \nFigure 3. Narrow line ratio diagnostic diagram of [O /i.pc/i.pc/i.pc] λ 4363/H γ vs. [O/i.pc/i.pc/i.pc] λ 5007/[O /i.pc/i.pc/i.pc] λ 4363. Different regions within ZS7 are indicated by golden symbols: the BLR location (star), the [O /i.pc/i.pc/i.pc] λ 5007 centroid (diamond), and the surrounding region (circle; see Sec. 7 for details). Local star-forming galaxies (grey shading and points) and AGN (red shading and points), dust-corrected following Cardelli et al. (1989), separate into two parallel sequences in this diagram, with AGN occupying a region of higher electron temperatures and having elevated ratios of [O /i.pc/i.pc/i.pc] λ 4363/H γ . The red diamond shows a galaxy from the JWST Early Release Observations (Pontoppidan et al. 2022) at similar redshift, which is identified as an AGN through a [Ne IV] 𝜆𝜆 2422 , 2424 detection (Brinchmann 2023; B23), and the red star indicates the narrow line ratios of a BLR AGN at 𝑧 = 5 . 55 (Übler et al. 2023; Ü23). The BLR centroid of ZS7 is located at the extreme end of the local AGN distribution, and also the [O /i.pc/i.pc/i.pc] λ 5007 centroid falls on top of the local AGN cloud, while the host emission is located in a region where local star-forming galaxies and AGN overlap. The narrow line ratios extracted from the [O /i.pc/i.pc/i.pc] λ 5007 centroid position could indicate recent AGN activity also in this part of ZS7.  \n<!-- image -->  \n[O/i.pc/i.pc/i.pc] λ 5007, indicates very high gas temperatures, of the order of 𝑇 𝑒, [O/i.pc/i.pc/i.pc] λ 5007 = 22700 ± 1900 K (following Dors et al. 2020). Such a high ratio of [O /i.pc/i.pc/i.pc] λ 4363/H γ is not seen in nearby star-forming galaxies, but it is seen in AGN (Perna et al. 2017; Brinchmann 2023). This is further illustrated in Figure 3, where we plot nearby AGN (red contours) and star-forming galaxies (grey contours) from SDSS (Abazajian et al. 2009) in the [O /i.pc/i.pc/i.pc] λ 4363/H γ vs. [O/i.pc/i.pc/i.pc] λ 5007/[O /i.pc/i.pc/i.pc] λ 4363 line ratio diagram. In this diagram, local star-forming galaxies and AGN separate into two parallel sequences, with AGN occupying a region of higher electron temperatures and having elevated ratios of [O/i.pc/i.pc/i.pc] λ 4363/H γ (Fig. 3; also Mazzolari et al., in prep.). Indeed, the narrow emission line ratios at the location of the broad H β in ZS7 lie at the extreme end of the nearby AGN distribution (golden star in Fig. 3). This is yet another confirmation that the broad H β is tracing a (type 1) AGN.  \nIntriguingly, as illustrated by the bottom right panel of Fig. 2, the auroral [O /i.pc/i.pc/i.pc] λ 4363 line is spatially extended across the locations of both the BLR and [O /i.pc/i.pc/i.pc] λ 5007 centroid. We can measure the [O /i.pc/i.pc/i.pc] λ 4363/H γ and [O /i.pc/i.pc/i.pc] λ 5007/[O /i.pc/i.pc/i.pc] λ 4363 line ratios at the [O/i.pc/i.pc/i.pc] λ 5007 centroid position, and we find that also the emission line  \nproperties at the location of the [O /i.pc/i.pc/i.pc] λ 5007 centroid are consistent with local AGN, as shown by golden diamond in Fig. 3. This finding strongly suggest that also the galaxy associated with the [O /i.pc/i.pc/i.pc] λ 5007 centroid is hosting a (type 2) AGN.  \nWefinally note that, if we extract these emission line ratios from an integrated spectrum excluding the region where the BLR is detected with 𝑆 / 𝑁 > 5, and excluding the 3 × 3 spaxels around the [O /i.pc/i.pc/i.pc] λ 5007 centroid, then we find values falling into the region which is locally occupied by both AGN and star-forming galaxies (golden circle in Fig. 3).', '8 AN IMMINENT GALAXY AND MASSIVE BLACK HOLE MERGER': "As discussed in the previous Sections, the NIRSpec data and the NIRCam data provide evidence that ZS7 is composed of two components in close vicinity. The projected separation of 620 pc and velocity difference of the narrow line emission between the BLR location and the centroid of the [O /i.pc/i.pc/i.pc] λ 5007 emitter ( ∼ 20 km/s), likely hosting another AGN, are small, and likely imply an imminent merger.  \nIn Section 5, we measured the mass of the black hole associated with the type 1 AGN, traced by the broad H β , by using virial relations. Estimating the mass of the putative second black hole, associated with the [O /i.pc/i.pc/i.pc] λ 5007 emitter, is more difficult. Using our stellar mass estimate of the [O /i.pc/i.pc/i.pc] λ 5007 emitter ( 𝑀 ★ ∼ 2 . 9 × 10 9 𝑀 ⊙ , i.e. excluding the region where the BLR is significantly detected, see Section 3.3) and assuming a 𝑀 ★ -𝑀 · relation calibrated on nearby AGN (Reines & Volonteri 2015), results in a secondary black hole mass associated with the [O /i.pc/i.pc/i.pc] λ 5007 emitter of log ( M · / M ⊙ ) ∼ 5 . 8 ± 0 . 6. If we assume a lower limit for the stellar mass of 𝑀 ★ ∼ 1 . 7 × 10 9 𝑀 ⊙ (see discussion in Section 3.3), we find log ( M · / M ⊙ ) ∼ 5 . 6 ± 0 . 6. Comparing this to the black hole mass estimate for the BLR AGN (log ( M · / M ⊙ ) = 7 . 7 ± 0 . 4), we find a ratio that is very similar to the dual AGN candidates identified by Maiolino et al. (2023b). If instead we assume a 𝑀 ★ -𝑀 · relation calibrated on local S/S0 and elliptical galaxies, we get a black holes mass higher by about one order of magnitude. Alternatively, if we assume a black hole-to-stellar mass ratio between 0.01 and 0.1, as found for many high-z BLR-AGN discovered with JWST, the corresponding mass of the putative second black hole would be comparable to or even larger than the mass of the BLR-identified AGN in ZS7. We note that simulations find that in 'offset' AGN systems (i.e. interacting galaxies in which only one of the two black holes is accreting) the active black hole is typically more massive (Steinborn et al. 2016; Chen et al. 2023), although these results are only available at 𝑧 ≲ 3.  \nContinuing with the interpretation of an imminent massive black hole merger, at a projected separation of 620 pc the black hole pair is likely in the dynamical friction phase when the interaction with stars and gas assists the black holes in reducing their angular momenta to sink to the common centre of mass (Chandrasekhar 1943; Ostriker 1999), before three-body interactions with individual stars (Begelman et al. 1980; Quinlan 1996) and eventually the emission of gravitational waves lead to coalescence of the black hole pair (Peters &Mathews 1963; Peters 1964). Using our black hole mass estimate, the stellar mass estimate of the [O /i.pc/i.pc/i.pc] λ 5007 emitter summed within roughly 3 kpc, and assuming an absolute separation of the BLR and the [O /i.pc/i.pc/i.pc] λ 5007 centroid of 1 kpc, we find an approximative inspiral time of 100-200 Myr (Binney & Tremaine 2008).", '9 IMPLICATIONS': "It is not possible to provide any reliable statistics on the number and fraction of merging black holes at high redshift, based on this single case. We note, however, that this galaxy was observed as part of a sample of 11 galaxies at 𝑧 > 6 that were observed with the NIRSpecIFU mode within the GA-NIFS survey. Of these, the ZS7 system and a second galaxy (Cresci et al., in prep.) are found to clearly host an AGN based on the detection of broad Balmer lines. The other AGN at 𝑧 > 6 in GA-NIFS is also found in an interacting system, on a few kpc separation. Also the two 𝑧 ∼ 6 . 8 QSOs analysed so far within GA-NIFS undergo galaxy-galaxy mergers (Marshall et al. 2023; see also e.g. Yue et al. 2021). Although these are low-number statistics, these results suggest that AGN (and more broadly black hole activity) in the early Universe are associated with galaxy interactions. In the case of ZS7, the separation of the individual systems is so small that the galaxies and their supermassive black holes will most likely merge within the next few hundred million years. At slightly lower redshifts (4 < 𝑧 < 6), a recent study by Maiolino et al. (2023b) identified a significant number (3/11) of candidate dual AGN based on the complex profile of broad Balmer lines in the NIRSpec multiobject JWST Deep Extragalactic Survey (JADES; Eisenstein et al. 2023, without 2d spatial information). At 𝑧 ∼ 3, Perna et al. (2023b) identify 20-30% dual AGN on separation scales of 3-30 kpc within GA-NIFS. Summarising, the findings of our and other studies using different approaches to assess black hole activity in the early Universe are consistent with a potentially large number of black holes in the process of merging.  \nThere are no simulations available that provide the expected number of binary black holes, with separations less than 1 kpc, in this mass range and at 𝑧 ≳ 7. However, the mentioned finding of a large fraction of dual AGN at lower redshift (z ∼ 3) on scales of a few/several kpc already highlighted tension with the expectation of simulations, which predict a lower fraction of dual AGN at 𝑧 < 5 (Rosas-Guevara et al. 2019; De Rosa et al. 2019; Volonteri et al. 2022; Chen et al. 2023) than observed (Koss et al. 2012; Spingola et al. 2019; Perna et al. 2023b). However, while some simulations find an increasing fraction of massive black hole pairs with decreasing redshift due to long coalescence timescales of massive black hole mergers in these models (Chen et al. 2023), most simulations predict the fraction of massive black hole pairs to increase with increasing redshift due to higher galaxy-galaxy merger rates (Volonteri et al. 2016; RosasGuevara et al. 2019; De Rosa et al. 2019; Volonteri et al. 2022; see also discussion by Sesana et al. 2009).  \nWealsomentionanother possible scenario explaining our observations of ZS7, namely that the displaced accreting black hole actually results from the recoil of a recent supermassive black hole merger. An-isotropic emission of gravitational waves during a black hole merger event may give the resulting black hole a velocity kick (Peres 1962; Bekenstein 1973). Some simulations predict these velocities could be as high as several hundred to a few thousand km/s for premerger black holes with comparable masses, which could even lead to the ejection of the merged black hole from the host galaxy (e.g. Madau & Quataert 2004; Merritt et al. 2004; Blecha et al. 2016). Recoiling black holes may be observable through both a spatial and a velocity offset to the host galaxy (see e.g. Castignani et al. 2022), and this could be the case in ZS7: a spatial offset of 620 pc between the host galaxy emission and the BLR is evident in this case. In addition, we register a small velocity offset between the narrow line emission of ZS7 and the BLR emission of ∼ 40 km/s. While recoil velocities depend on the masses and relative spins of the initial black holes, the fact that we do not observe a large velocity offset of the  \nBLR makes the interpretation of a recoiling supermassive black hole less appealing, although we caution that we cannot measure the absolute velocity difference. Furthermore, a recoiled black hole would probably not have a NLR around it and, if it did, it would likely be the interstellar medium of the host galaxy, hence it should share the same velocity field as the host. In the case of ZS7 we clearly see that the narrow component of the nebular emission around the accreting black hole, i.e. its NLR, has a different velocity (Fig. 2). This suggests a more likely scenario in which the accreting black hole has yet to merge and is still embedded in the interstellar medium of its own galaxy. Finally, the indication that the [O /i.pc/i.pc/i.pc] λ 5007 centroid is also likely to host another (type 2) AGN, implies that the merging has yet to happen, unless the narrow lines signatures are actually the fossil record of past accretion at that location. However, even the recoil scenario would be pointing at a black hole merging event; the difference between the two scenarios would be that in one case the merging has yet to happen, while in the other case it has just happened.  \nFinally, we note that another possible explanation for the offset BLR could be that it is tracing one of the so-called 'wandering intermediate-mass black holes (IMBH)'. Recent simulations have shown that (at least on scales of several kpc) the dynamical friction during galaxy interactions/mergers may not be efficient at sinking black holes at high redshift (Di Matteo et al. 2023b). As a consequence, this may result into a potentially large population of IMBHs in the haloes of galaxies. However, such population of wandering IMBHs is expected to be on scales much larger (tens of kpc) than the separation observed in ZS7 (where dynamical friction is much more efficient). Additionally, a 'wandering' black hole, even if accreting, should be essentially 'naked', i.e. should not have its own NLR. This is in contrast with what is observed for the offset BLR in ZS7, which is clearly exciting at least [O /i.pc/i.pc/i.pc] λ 4363 tracing low density gas outside the BLR.", '10 CONCLUSIONS': 'The results presented in this work highlight the power of NIRSpecIFU observations in the study of high -𝑧 galaxies, mergers, and the detection of moderately massive black holes through imagingspectroscopy. Our observations provide clear and robust evidence for a massive black hole involved in a merger with another galaxy, likely hosting another accreting black hole, at 𝑧 = 7 . 15, only 740 Myr after the Big Bang. Overall, our results seem to support a scenario of an imminent massive black hole merger in the early Universe, highlighting this as an additional important channel for the early growth of black holes. Together with other recent findings in the literature, this suggests that massive black hole merging in the distant Universe is common. Our observations may be used as guidance for the modelling of gravitational wave events originating from massive black hole mergers that will be detectable with future observatories like LISA (e.g. Haehnelt 1994; Jaffe & Backer 2003; Sesana et al. 2005; Valiante et al. 2021; Amaro-Seoane et al. 2023).', 'ACKNOWLEDGEMENTS': "We are grateful to the anonymous referee for a constructive report that helped to improve the quality of this manuscript. We thank Martin Bourne, Massimo Dotti, Sophie Koudmani, Alberto Sesana, Debora Sĳacki, Sandro Tacchella and Naoki Yoshida for useful discussions. AJB and GCJ acknowledge funding from the 'FirstGalaxies' Advanced Grant from the European Research Council (ERC)  \nunder the European Union's Horizon 2020 research and innovation program (Grant agreement No. 789056). BRdP, MP and SA acknowledge grant PID2021-127718NB-I00 funded by the Spanish Ministry of Science and Innovation/State Agency of Research (MICIN/AEI/ 10.13039/501100011033). FDE, JS and JW acknowledge support by the Science and Technology Facilities Council (STFC), by the ERC through Advanced Grant 695671 'QUENCH', and by the UKRI Frontier Research grant RISEandFALL. HÜ gratefully acknowledges support by the Isaac Newton Trust and by the Kavli Foundation through a Newton-Kavli Junior Fellowship. IL acknowledges support from PID2022-140483NB-C22 funded by AEI 10.13039/501100011033 and BDC 20221289 funded by MCIN by the Recovery, Transformation and Resilience Plan from the Spanish State, and by NextGenerationEU from the European Union through the Recovery and Resilience Facility. JSD acknowledges the support of the Royal Society through a Royal Society Research Professorship. MAM acknowledges the support of a National Research Council of Canada Plaskett Fellowship, and the Australian Research Council Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), through project number CE170100013. PGP-G acknowledges support from grant PID2022139567NB-I00 funded by Spanish Ministerio de Ciencia e Innovación MCIN/AEI/10.13039/501100011033, FEDER, UE. RM acknowledges support by the Science and Technology Facilities Council (STFC), from the ERC Advanced Grant 695671 'QUENCH', and funding from a research professorship from the Royal Society. SCa and GV acknowledge support from the European Union (ERC, 'WINGS', 101040227).", 'DATA AVAILABILITY': 'The NIRSpec data used in this research has been obtained within the NIRSpec-IFU GTO programme GA-NIFS and will become publicly available in May 2024. The NIRCam observations of ZS7 are public through the PRIMER programme (PI: James Dunlop). Data presented in this paper will be shared upon reasonable request to the corresponding author.', 'REFERENCES': '```\nAbazajian K. N., et al., 2009, ApJS, 182, 543 Amaro-Seoane P., et al., 2023, Living Reviews in Relativity, 26, 2 Baldwin J. A., Phillips M. M., Terlevich R., 1981, PASP, 93, 5 Barai P., Gallerani S., Pallottini A., Ferrara A., Marconi A., Cicone C., Maiolino R., Carniani S., 2018, MNRAS, 473, 4003 Barausse E., Dvorkin I., Tremmel M., Volonteri M., Bonetti M., 2020, ApJ, 904, 16 Begelman M. C., Blandford R. D., Rees M. J., 1980, Nature, 287, 307 Bekenstein J. D., 1973, ApJ, 183, 657 Bennett J. S., Sĳacki D., Costa T., Laporte N., Witten C., 2024, MNRAS, 527, 1033 Binney J., Tremaine S., 2008, Galactic Dynamics: Second Edition. Princeton University Press Blecha L., et al., 2016, MNRAS, 456, 961 Bogdán Á., et al., 2023, Nature Astronomy, Böker T., et al., 2022, A&A, 661, A82 Böker T., et al., 2023, PASP, 135, 038001 Brinchmann J., 2023, MNRAS, 525, 2087 Brinchmann J., Pettini M., Charlot S., 2008, MNRAS, 385, 769 Bruzual G., Charlot S., 2003, MNRAS, 344, 1000 Calzetti D., Armus L., Bohlin R. C., Kinney A. L., Koornneef J., StorchiBergmann T., 2000, ApJ, 533, 682\n```  \nVolonteri M., et al., 2020, MNRAS, 498, 2219  \nVolonteri M., Pfister H., Beckmann R., Dotti M., Dubois Y., Massonneau W.,  \nMusoke G., Tremmel M., 2022, MNRAS, 514, 640  \nVolonteri M., Habouzit M., Colpi M., 2023, MNRAS, 521, 241  \nWeller E. J., Pacucci F., Ni Y., Chen N., Di Matteo T., Siwek M., Hernquist  \nL., 2023, MNRAS, 520, 3955  \nWitten C., et al., 2023a, arXiv e-prints, p. arXiv:2303.16225  \nWitten C. E. C., Laporte N., Katz H., 2023b, ApJ, 944, 61  \nYang G., et al., 2023, ApJ, 950, L5  \nYue M., Fan X., Yang J., Wang F., 2021, ApJ, 921, L27  \nZhang H., et al., 2023, arXiv e-prints, p. arXiv:2309.07210  \nvan Dokkum P. G., 2001, PASP, 113, 1420', 'APPENDIX A: FULL SPECTRAL FITS': "In Fig. A1 we show full spectral fits to spectra extracted from individual spaxels at the BLR location, and at the centroid location of the [O /i.pc/i.pc/i.pc] λ 5007 emitter. At the BLR location, the broad H β emission is evident. Some flux of the BLR point-source is spread through the point-spread function (PSF FWHM ∼ 0 . 13 -0 . 14 '' at 4 𝜇 m) and therefore is still visible at the [O /i.pc/i.pc/i.pc] λ 5007 centroid position.", 'APPENDIX B: DUST EXTINCTION': "In Fig. B1 we show a map of the PSF-matched H γ to H β narrow emission line flux, derived using the empirical fitting functions for the NIRSpec-IFU PSF provided by D'Eugenio et al. (2023). The centroid locations of the [O /i.pc/i.pc/i.pc] λ 5007 and BLR are indicated by a cross and star, respectively. The black line in the colorbar shows the theoretical line ratio of H γ /H β =0.468, whereas lower values indicate dust extinction. The data show that there is no enhanced dust extinction at the location of the BLR. This rules out a scenario where the BLR centroid location is the true, dust-obscured centre of the ZS7 system.  \nThis paper has been typeset from a T E X/L A T E X file prepared by the author.  \nFigure A1. Full spectral fits to spectra extracted from individual spaxels at the BLR location (top), and at the centroid location of the [O /i.pc/i.pc/i.pc] λ 5007 emitter (bottom). Residuals show (data - model)/uncertainties. We show our full fit in blue, with narrow emission lines in green and the BLR component in purple. At the BLR location, the broad H β emission is evident, while the [O /i.pc/i.pc/i.pc] λλ 4959,5007 emission is narrow at all locations. We note the strong auroral [O /i.pc/i.pc/i.pc] λ 4363 emission line in both spectra.  \n<!-- image -->  \nFigure B1. Map of the PSF-matched H γ to H β narrow emission line flux where both lines have 𝑆 / 𝑁 > 5, with the centroid location of the [O /i.pc/i.pc/i.pc] λ 5007 and the location of the BLR indicated by a cross and star, respectively.  \n<!-- image -->"}