A redshifted excess in the broad emission lines after the flare of the $\gamma$-ray narrow-line Seyfert 1 PKS 2004-447

PKS 2004-447 is a narrow-line Seyfert 1 (NLS1) harbouring a relativistic jet with gamma-ray emission. On 2019-10-25, the Fermi-Large Area Telescope captured a $\gamma$-ray flare from this source, offering a chance to study the broad-line region (BLR) and jet during such violent events. This can provide insights to the BLR structure and jet interactions, which are important for active galactic nuclei and host galaxy coevolution. We report X-Shooter observations of enhancements in the broad line components of Balmer, Paschen and He I lines seen only during the post-flare and vanishing 1.5 years after. These features are biased redward up to $\sim$250 km s$^{-1}$ and are narrower than the pre-existing broad line profiles. This indicates a connection between the relativistic jet and the BLR of a young AGN, and how $\gamma$-ray production can lead to localised addition of broad emission lines.


Introduction
Since the first observations (Seyfert 1943) of broad features in active galactic nuclei (AGN) spectra, the nature of the emitting region in terms of its motions and structure has been the subject of much discussion and speculation (e.g., Mathews & Capriotti 1985;Leighly & Casebeer 2007;Bon et al. 2009;Shapovalova et al. 2009;Yong et al. 2017).The general view of the broad-line region (BLR) is that it is a region of photoionised gas that is subjected to high velocities -with a full width at half maximum (FWHM) > 2000 km s −1 -through random dispersion and/or virialised orbital rotational motions, as well as having a potential for bulk infall or outflow (e.g., Wanders et al. 1995;Done & Krolik 1996;Denney et al. 2009;Williams et al. 2018).
Broad-line-region outflows are of particular interest to the topic of AGN and host galaxy co-evolution (Zubovas & Nayakshin 2014;King & Pounds 2015;Harrison et al. 2018).In addition to cases of blue-shifted emission lines (e.g., Ge et al. 2019), broad absorption lines are observed in ∼20% of quasars (Hamann et al. 1993;Trump et al. 2006), and can reach velocities up to v ∼ 0.3c (Rogerson et al. 2016;Choi et al. 2020).This is sufficient energy (Scannapieco & Oh 2004;Hopkins & Elvis 2010;Harrison et al. 2018) to effectively cause feedback to the host galaxy as it eventually reaches and interacts with the interstellar medium.
In this work, we present a detailed analysis of the broad emission lines of X-shooter data of PKS 2004−447 surrounding the γ-ray flare.There are two 'low-state' spectra taken on 2017 May 28 (SP0) and 2017 June 19 (SP1, both in programme 099.B-0785).After the flare, a 'post-flare' spectrum was obtained on 2019 November 25 (SP2), and then a 'reverted-state' spectrum on 2021 May 05 (SP3, both in Director Discretionary Time programme 104.20UC).We deliberately observed the post-flare spectrum, SP2, when the γ-ray flare was predicted to reach the BLR (see Berton et al. 2021, for size estimations).We show, for the very first time, a temporal relation between a γ-ray event, jet flaring, and BLR changes localised to the redshifted part of the emitting region.This work adopts a standard ΛCDM cosmology with H 0 = 70 km s −1 Mpc −1 and Ω λ = 0.7.All magnitudes are provided in the AB system.

X-shooter observations and data reductions
X-shooter observes simultaneously in the UVB, VIS, and NIR arms, spanning a total wavelength range of 3000-25 000 Å.All four spectra were taken with a 1 slit.Details of the four X-shooter spectra are listed in Table 1, and they are displayed over the full X-shooter range in the appendix (see Fig. A.2).
For the reduction of raw data, we used the default X-shooter pipeline provided by ESO 1 with bias, flat, lamp corrections, and response correction with a standard star.We also used ESO Molecfit2 to remove the telluric absorptions (Smette et al. 2015;Kausch et al. 2015).SP1 had sufficient signal in the NIR arm of X-shooter to accurately recover the features and the continuum that was absorbed.We de-reddened the spectra using a Fitzpatrick (1999) law with E(B − V) = 0.0846 (Schlafly & Finkbeiner 2011) and R v = 3.1.We measured the Na i doublet to obtain the redshift of the host galaxy, obtaining a value of z = 0.24035 ± 0.00009 for SP0, SP1, and SP3, and z = 0.24050 ± 0.00009 for SP2 due to the motion of the Earth.We assume the lower value redshift for calculations in this work, and an operating rest-frame wavelength of 2418-20 155 Å.
The airmasses and seeing vary for the four observations, suggesting a need to fine-tune flux calibration.This is usually done by assuming non-varying forbidden lines, but the different slit orientations sample different parts of the narrow-line region.This variation only affects the forbidden-line, as the continuum and permitted lines from the central region are well within the 1 slit.Thus, the flux per spectrum was fine-tuned based on the acquisition images (details in Appendix A).The calibration is accurate as the spectra match the spectral energy distribution (SED) from Gokus et al. (2021) in Fig. 1.It is not possible to compare the variation in the forbidden lines across these four spectra due to this systematic.
By averaging SP0, SP1, and SP3, we constructed SP013, a higher signal-to-noise ratio (S/N) stacked spectrum representing the long-term average of PKS 2004−447.This smooths out stochastic and systematic variation that is not characterised, and serves as the baseline for comparing with SP2.Statistically significant residuals in SP2−SP013 can be easily attributed to the γ-ray flare.

Results
The continuum of PKS 2004−447 is blended.There is the power law from the AGN accretion disk; Fe ii multiplet emission; a host galaxy that is resolved in the aperture (Kreikenbohm et al. 2016); and synchrotron emission from the jet.Its UV-NIR frequencies are dominated by synchrotron emission, resulting in a negative slope instead of the usual positive power law (Gokus et al. 2021).During the flare, the flux is boosted across all X-shooter wavelengths, with a greater increase for lower frequencies.
To characterise the differences between the four datasets, we constructed root mean square (rms) spectra using Eq.(3) from Peterson et al. (2004).Figure 2 shows the rms spectra with and without SP2.Without SP2, no features appear above the noise around Hδ, Hγ, Hβ, and He iλ5875.With SP2, those lines show noticeable enhancements.We investigated each spectrum's variability relative to SP013 and we find that the enhancement only appears in SP2, suggesting that the broad emission lines in the Balmer and helium series changed temporarily in SP2, likely related to the γ-ray flare.
We can also see the effects of mismatched slit orientations here.While SP1−SP013 has somewhat clean subtractions in the forbidden lines, SP0−SP013 and SP3−SP013 do not.This makes interpreting the changes in Hα difficult.While this emission line is also affected, it is hard to disentangle observed and systematic variation on the [N ii] line.
A closer look at the enhancement in Hβ as well as a comparison of all epochs can be found in Fig. 3 (see Fig. B.2 for the full version).Outside of SP2, the Hβ line remains constant with the subtracted spectrum, consistent with a ∼0 mean value within the noise level.All of the Hβ enhancement is significant3 at >6σ.We also discover that the peak of the Hβ excess is redshifted by4 240±50 km s −1 (4.8σ), which is consistent across all of the subtracted spectrum.We refer to this enhancement with a redward bias as the 'red excess'.All of the Balmer lines and He iλ10830 exhibit the red excess, but it is less obvious in He iλ5875 (see Line fitting is not trivial for PKS 2004−447 as the continuum is not a power law.Approximations can work for individual epochs, but they are inconsistent across all four spectra.This is because both the emission lines and continuum can vary, so it is   not possible to disentangle the source of any measured variation.Therefore, we only fitted SP013 to provide average BLR properties and SP2−SP013 to estimate red excess properties (details in Appendix B).Enhancements and red excesses are present independently and they are not due to continuum modelling.
The width of all lines are around FWHM = 1600 km s −1 .All Balmer and Paschen lines have a <50 km s −1 velocity offset, while all He i lines are blueshifted to 85 km s −1 .In addition, O iλ8446 is redshifted to 180 km s −1 , suggesting an inflow.
For the red excess, we measured a FWHM ∼ 1000 km s −1 for all Balmer and Paschen lines, but the offset differs for different line series.Furthermore, Hδ, Hγ, and Hβ are redshifted by ∼250 km s −1 , and Hα is redshifted by 180 km s −1 .Also He i red excesses are comparatively weaker for this analysis, but we observe roughly the same width and up to 340 km s −1 of redshift.As mentioned, Paschen lines do not have a noticeable redshift, suggesting that the red excess is more prominent towards the blue.

The impact of observational effects on our findings
We explore the possibility that the red excesses are caused by observational effects.This is a possibility since the observing conditions between all spectra are different.Namely, SP2, which is the spectrum closest to the flare and the spectrum that is showing the red excess, was observed under the worse conditions relative to the rest (see Table 1).
The immediate cause of concern with SP2 is flux dimming caused by high ending airmass at ∼2.1 and seeing 1.2 .The drop in flux caused by poor seeing or field differential dispersion is uniform, and this was corrected when we calibrated the flux photometrically (see Appendix A).Also, as seen in Table 1, SP0 and SP1 have a similar S/N given similar exposure times, followed by SP3 and SP2 with the best S/N since PKS 2004−447 is at its brightest.This demonstrates that the quality of spectra is more dependant on exposure times and source brightness.
Moreover, PKS 2004−447 has a scale of 3.8 kpc/ at z = 0.24.For a BLR radius of 15 light days (Berton et al. 2021), this corresponds to a diameter of 2.5×10 −5 kpc, which covers 1 .With a slit of 1 , this BLR scale is too small for airmass and seeing to have any relevant effect on the emission lines.Even if they did, that should cause the BLR to be dimmer, not brighter, and that would also affect the NLR and Mg ii line, which is not observed.
Finally, VLT/X-shooter is equipped with atmospheric differential correctors (ADCs), which prevent spectral distortions for airmass lower than 2.5.Therefore, it is unlikely for the red excess to be an observational artefact.
L14, page 3 of 19 Grey regions represents sections of the X-shooter spectra with bad pixels.SP013 in both panels are smoothed with a median filter with a 11 pixel window, while the other lines are smoothed with a 31 pixel window.

BLR geometry
The difference in peak offset between H (close to systematic redshift) and He (almost 100 km s −1 blueshifted) indicates a stratified BLR, while the larger FWHM in He i indicates that He is at a smaller radial distance.This is consistent with the theory (Korista & Goad 2004) and observations from velocity-resolved reverberation mapping (e.g., Kollatschny 2003;Bentz et al. 2010).However, the large uncertainty associated with He i (see Table B.1) makes it impossible to meaningfully conclude the BLR geometry based on this information.
The kinematics of the Balmer lines were well measured and so was the velocity offset of the redshifted O i line.The latter implies the presence of an inflow.We know that this jetted source is observed face on, with a viewing angle of <50 deg, as constrained by Schulz et al. (2016).Therefore, under a diskwind model, for Balmer lines to not have a significant peak offset implies the emitting region is moving perpendicular to the line of sight.Building on this assumption, two scenarios exist for placing the other emitting regions.The first is a single stream of material that follows the locally optimally emitting clouds (LOC, Baldwin et al. 1995) model, with He i, Paschen, and Balmer lines, ending with O i.This implies a greater push at the start of the stream, thus the blueshifted peaks of He i.The second is a dual-stream disk model (e.g., Yong et al. 2017) with He wind radially closer and at a higher inclination than the H wind, not placing Oi inflow at any required position.Reverberation mapping will be needed to fully investigate the structure, and a higher resolution spectrum with better S/N focusing on the Paschen, He i, and He ii lines is required.
Our results reveal perplexing line fluxes and ratios.The Balmer decrement is very steep at Hα:Hβ:Hγ:Hδ = 4.67:1.00:0.31:0.15.Also, Paα/Paβ = 1.08, but Paα/Hβ = 0.50 is very low (all values were derived from the fluxes presented in Table B.1).A reddening effect does not explain this as a consistent range of case B values for Hγ/Hβ would worsen L14, page 4 of 19 Fig. 3. Variation of the Hβ region in the low-state, post-flare and reverted-state.All lines are smoothed with a median filter with a 5 pixel window, which is much much narrower than any features that could be present.In each panel, a relevant spectrum, shown in orange, is subtracted against SP0, which is acting as the baseline in these plots and shown in blue.Their subtracted result is shown in green.The observed flux is the y-axis and velocity offset is x-axis.We also indicate the 0 value in the y-and x-axes for every panel with vertical and horizontal dashed lines.These visual aids emphasise the non-zero redshift in the emission line excess in Hβ during the post-flare.The statistical significance of the excess is also given in the relevant panels, see footnote 3 for details.We see that the green line in the low-state and reverted-state are consistent with 0, indicating no broad-line variation outside of the flare.This figure is a truncated version of the Paα/Hβ ratio and result in Hα/Hβ being too low.A differential in optical depth could explain the Balmer ratios, yet this would require an unusually high optical depth of 1000 (see Fig. 9 of Davidson & Netzer 1979).Balmer line ratios are known to be correlated with a velocity offset, where the Hα/Hβ increases towards the peak of a line profile (Stirpe 1991).This could be a feasible explanation as the line profile observed in PKS 2004−447 is highly dominated by the peak.Alternatively, the emission lines may not be completely recombination driven, and a partial local thermodynamical equilibrium exists and the Saha-Boltzmann relations hold for part of the BLR (Popović 2003).This is supported by the Balmer ratios in PKS 2004−447 being similar to 3C 390.3 and 3C 382 from that work.

Understanding of red excess
The red excess and the flare are correlated events.This is the first observed event where a γ-ray flare is associated with a local change in the BLR with distinct kinematic properties.If multiple spectra were taken as soon as the flare started, we could assess if the red excess appeared before or after the flare.If they were taken before, we would observe the red excess fading, implying a link to the cause of flare.If they were taken after, we would see the excess grow, suggesting a direct contribution of the jet to the red excess.
The most curious aspect of the red excess is that it is redshifted.Here we speculate about four possible scenarios on how it is formed.All of the speculated scenarios would be extremely interesting to study, and we encourage frequent observation immediately after a flare for all γ-NLS1s.
(1) The simplest scenario is a bipolar outflow, with receding gas visible and approaching gas partially obscured.This is motivated by the red excess being more prominent towards the blue and absent for the IR Paschen lines.The O i inflow also supports this model.The flare would reverberate throughout both flows equally, reaching most of the material within a month.However, shorter wavelength emissions from the approaching ions would be obscured, biasing a redward enhancement.
(2) The γ-ray flare directly creates the red excess.If the event occurred above the H and He emitting region, a shock wave from the event would push the gas down, introducing more material into the emitting region.In this scenario, we should also observe red excesses in Mg ii and O i, but that is not seen.This is a point of contention for this scenario.Using the BLR radius calculated by Berton et al. (2021) of ∼15 light days, a FWHM ∼ 1500 km s −1 , and assuming t dyn ≈ R BLR /FWHM gives a dynamical timescale of 8 years.This is the time needed for the introduced material to re-equilibrate, which is longer than the 1.5 years between SP2 and SP3.Reconciling this requires that the event occurred at a distance <3 light days, which is consistent with the predicted location of γ-ray production at ∼10 3 r g ∼ 1 light day (Ghisellini et al. 2010).
(3) The flaring event led to an increased illumination onto a region close enough to the black hole to experience strong gravitational redshift.A shift of 250 km s −1 yields an emitting radius of ∼1 light day (with M bh = 1.5 × 10 7 M from Berton et al. 2021), which is consistent with an increased illumination from the γ-ray flare itself.The emergence of the red excess should be instant due to its proximity to the black hole and γ-ray source.However, its visibility a month after the flare poses an issue, even accounting for gravitational time dilation.
(4) The accretion disk or BLR undergo a structural change that leads to the flare and increased illumination in a specific part of the BLR.This region could be close to or far from the black hole.If it is close, the duration of the illumination is not tied to the short-lived flare unlike in scenario (3).If it is far, since O i is an established inflow, that might imply the outermost parts of the BLR are always infalling.This region might contribute to part of the red wing observed in the permitted line profiles, which can be verified with velocity-resolved reverberation mapping.The relevant timescale for this interaction is not well understood, as the infalling timescale based on the viscosity parameter L14, page 5 of 19 has been shown to be flawed as demonstrated by changing-look AGN (Lawrence 2018).

Appendix B: Spectral analysis and line fitting
We used the python script PyQSOFit by Guo et al. (2018) that utilises the least squares fitting package kmpfit.The residuals were weighted based on the variance spectra.PyQSOFit has the functionality to fit a continuum and lines.For continuum fitting, while optimising the fitting range to feasibility, we approximated the rest-frame range from 3300-7500Å.This encompasses half of the UVB and the entire VIS arm.We fitted, with a free power law, a third-order polynomial, and the Fe ii template from Boroson & Green (1992).The NIR arm was fitted separately and can be approximated by a power law in its entirely.
For line fitting, we used a modified version of PyQSOFit called PyQSOFit_SBL5 .This version allows for skewed Gaussian and Voigt profile fitting.It also has the functionality to link the parameters of an unlimited number of fitted components.This reduced the degree of arbitrariness as we de-blended the Hα line.It also reduced the number of individual parameters to vary per fitting, as well as increased the sampling of the parameter of each related component for a more accurate result.While the script is able to perform bootstrapping for error estimation, the more realistic errors would be from the variance of each measured value throughout the three spectra for the SP013.
One key finding through trial and error is that the forbidden lines need to be fitted with three components.These are a broad skewed Voigt for the outflows, a narrow skewed Gaussian intermediate component (NLR A), and a very narrow Gaussian component with very little offset relative to the host galaxy (NLR B).The NLR A and NLR B components are defined from [O iii]λ5007; however, the broad component is different for every line.Either coincidentally, or physically, the narrow components of the permitted lines can be associated with NLR A and NLR B, so these were used to simplify the fitting procedure.The line profiles are subsequently described here and are also shown Helium lines -Of the observed He lines, we could only reliably measure He iλ10830.The only He ii line is at λ4686, but it sits in the region with heavy Fe ii emission, which we were unable to constrain.We note that He i λ3889 and λ7065 are very weak; the former suffers from continuum definition issues and the latter suffers from profile blending from multiple lines.Furthermore, He i λ5876 has continuum definition issues, but we were able to acquire a reliable fit by assuming that it has the exact same profile as He iλ10830.Both He i λ5876 and λ10830 have a significant NLR A component.
Red excess All of the red excesses are fitted with a single non-skewed Gaussian component, with no width or offset constraints.The measurement for He iλ5876 was extremely problematic as the continuum is not well defined in that region, resulting in contamination from continuum variation.Also He iλ10830 suffers from contamination from telluric correction variation.As such, their widths and peak offsets are not as reliable as the Balmer lines.We also observed a red excess on He iλ6678, but this is heavily contaminated by [S ii]λ6725 and Hα variations.While not listed because of the inaccurate measurement, it is a significant red excess with a FWHM= 1765 ± 425 km s −1 , a redshift of 620±600 km s −1 , and a flux of 21.52±8.52× 10 −17 erg/s/cm/cm.L14, page 8 of 19    3, showing all of the lines with enhancements or red excesses.This is a grid of plots showing an emission line region from the four X-SHOOTER spectra.All lines have been smoothed with a median filter with a 5 pixel window, which is significantly narrower than any features that could be present.In each panel, two spectra (blue and orange line) and their subtracted result (green line) are plotted with the observed flux in the y-axis and the velocity offset in the x-axis.The two spectra in each panel are plotted depending on the row and column as indicated in large font.For example, the bottom right corner plots SP2 and SP3.In the top left corner, the panel with a black box plots SP013 and SP2.We also indicate the 0 value in the y-and x-axes for every panel with vertical and horizontal dashed lines.These visual aids emphasise the non-zero redshift in the emission line excess in panels involving SP2.The statistical significance of the excess is also given in the relevant panels (see footnote 3 for details).We see that in all other panels without SP2, the green line is consistent with 0, indicating almost no variation outside of the flare.

Fig
Fig. B.2). Weaker enhancements are seen in Paβ and Paα; however, they are not significantly redshifted.Also O i and Mg ii have no enhancements.Line fitting is not trivial for PKS 2004−447 as the continuum is not a power law.Approximations can work for individual epochs, but they are inconsistent across all four spectra.This is because both the emission lines and continuum can vary, so it is Notes.Columns: (1) spectrum name; (2) observation date; (3) observation date in modified Julian date (MJD); (4) airmass at the beginning of the observation; (5) airmass at the end of the observation; (6) DIMM seeing during the observation; (7) exposure time in seconds; (8) slit rotation as measured with respect to the first spectrum, SP0; (9) signal-to-noise ratio measured from the featureless region between 5300-5400 Å; (10) SDSS magnitudes obtained from the flux calibrated acquisition images.

Fig. 1 .
Fig. 1.Comparision of the spectral energy distribution of PKS 2004−447 to the photometrically calibrated X shooter spectra.Left and middle panels are adopted from Fig. 6 of Gokus et al. (2021), refer to original paper for detailed legend.Right panel shows SP0, the low state in black, and SP2, the post-flare in red.All three plots have the same axes units.The red and blue horizontal and vertical lines in the left and middle panel indicates the y-and x-axis limit plotted in the right panel.The match between the SED and spectrum for both states are very close.

Fig. 2 .
Fig. 2. Broad line enhancements in SP2, the post-flare spectrum.All plotted lines are indicated by the respective legends, while vertical dashdotted line indicates the key Balmer and He iλ5875 emission line.Blue region represents the UVB arm of X-shooter while the green region is VIS.
Fig.3.Variation of the Hβ region in the low-state, post-flare and reverted-state.All lines are smoothed with a median filter with a 5 pixel window, which is much much narrower than any features that could be present.In each panel, a relevant spectrum, shown in orange, is subtracted against SP0, which is acting as the baseline in these plots and shown in blue.Their subtracted result is shown in green.The observed flux is the y-axis and velocity offset is x-axis.We also indicate the 0 value in the y-and x-axes for every panel with vertical and horizontal dashed lines.These visual aids emphasise the non-zero redshift in the emission line excess in Hβ during the post-flare.The statistical significance of the excess is also given in the relevant panels, see footnote 3 for details.We see that the green line in the low-state and reverted-state are consistent with 0, indicating no broad-line variation outside of the flare.This figure is a truncated version of Fig. B.2.

Fig. A. 2 .
Fig. A.2. Full rest-frame wavelength range of SP0, SP1, SP2, and SP3.The colour sections indicate the three arms, while the two grey sections in the NIR arm mark the most problematic telluric regions that Molecfit was not able to recover.
in Figure B.1, including [O iii], and the results of the fitting is summarised in Table B.1.Balmer lines -In the SP013 spectrum, the broad component of Hδ, Hγ, and Hβ are similar.In particular, Hγ and Hβ have a detectable NLR A component.Furthermore, Hγ fitting is considered problematic as it is heavily blended with [O iii]λ4363 and the continuum approximation over-subtracts the blue wing.We note that Hα is complicated as it is directly blended with [N ii], and potentially blended with He iλ6678 and [S ii] on the red wing.The final fitting to measure Hα broad emission lines used a total of seven components; four for the [N ii] doublet with two NLR A and two broad components set to the theoretical ratio; three for Hα with NLR A, NLR B, and a broad component that is associated with Hβ (not fixed, but has to be similar).The blending from He iλ6678 and [S ii] was not significant, and the final component of [N ii] was comparatively too weak to alter the Hα broad line.Paschen lines -Paschen lines observed ranges from Pa to Paα; however, we have only been able to acquire reliable linefitting measurements for Paβ and Paα.We note that Pa is heavily blended with [S iii]λ9531, and Paδ and Paγ are both dominated by noise that was introduced by telluric absorption.We find a broad component and a NLR A component for both Paβ and Paα.Broad Paα is problematic due to telluric noise, causing an unreliable continuum definition.

Fig. B. 1 .
Fig. B.1.Line fitting on continuum subtracted SP013.Each panel has a top plot to show the components used to fit them, and a bottom plot to show the residuals.All lines are colour coded based on the legend of [O iii].The last plot shows the red-excess fitting for Hβ, which has the best and isolated signal.It is not conclusive if the feature can be approximated with a Gaussian.

Fig. B. 2 .
Fig. B.2. Full version of Figure3, showing all of the lines with enhancements or red excesses.This is a grid of plots showing an emission line region from the four X-SHOOTER spectra.All lines have been smoothed with a median filter with a 5 pixel window, which is significantly narrower than any features that could be present.In each panel, two spectra (blue and orange line) and their subtracted result (green line) are plotted with the observed flux in the y-axis and the velocity offset in the x-axis.The two spectra in each panel are plotted depending on the row and column as indicated in large font.For example, the bottom right corner plots SP2 and SP3.In the top left corner, the panel with a black box plots SP013 and SP2.We also indicate the 0 value in the y-and x-axes for every panel with vertical and horizontal dashed lines.These visual aids emphasise the non-zero redshift in the emission line excess in panels involving SP2.The statistical significance of the excess is also given in the relevant panels (see footnote 3 for details).We see that in all other panels without SP2, the green line is consistent with 0, indicating almost no variation outside of the flare.

Fig. B. 2 .
Fig. B.2. Continued.For the grid plots for Hγ, the additional vertical red dashed line indicates the wavelength position of [O iii]λ4364, and they appear to be properly cancelled in the subtracted line for all.

Fig. B. 2 .
Fig. B.2. Continued.Statistical significance measured using RMS value taken from SP2 featureless continuum between 18000-18500Å.This is different to Paβ as the quality from X-Shooter degrades towards the red end of the IR arm.

Table B .
1. Summary of permitted lines observed and fitted, and the red excesses observed and fitted.Ctm refers to lines for which we are unable to accurately define the continuum for a fit.For telluric, problematic, and blend, refer to text in Appendix B; (4) FWHM of the line profile from SP013 in km/s, with empty entries indicating that the feature was not fitted; (5) velocity offset of the line profile from SP013 in km/s, with negative values indicating a blueshift; (6) kurtosis or skew of the line profile from SP013; (7) flux of the emission line from SP013 in 10 −17 erg/s/cm/cm; (8) FWHM of the red excess in km/s; (9) velocity offset of the red excess in km/s; and (10) flux of the red excess in 10 −17 erg/s/cm/cm.