Disc wind or disc line? The extraordinary Fe-K feature of Mrk 1513

We discuss the origin of a very unusual spectral structure observed in the Fe-K band of the Seyfert galaxy Mrk 1513, a local ($z$=0.063) active galactic nucleus (AGN) that is efficiently accreting matter onto its central supermassive black hole ($L_{\rm bol}/L_{\rm Edd}\sim$0.5). We consider the highest quality X-ray observation of this source available to date, performed in 2003 by XMM-Newton. The hard X-ray spectrum is characterised by a remarkable spectral drop at $\sim$7 keV, which can be interpreted as either the onset of a broad absorption trough or the blue wing of a relativistic emission line. Overall, this complex feature is significant at>5$\sigma$, and it is qualitatively reminiscent of a P-Cygni profile. A serendipitous spectrum of lower quality taken by XMM-Newton in 2015 qualitatively confirms the presence of similar Fe-K structures. Although it is not possible to distinguish between the two physical scenarios on sheer statistical grounds with the current data, several considerations lend weight to the possibility that Mrk 1513 is actually hosting a persistent outflow at accretion-disc scales, thus adding to the handful of known AGN in which a wide-angle X-ray wind has been identified so far.


Introduction
It is generally accepted that accreting supermassive black holes (SMBHs) in active galactic nuclei (AGN) play a key role in shaping the evolution of their host galaxies (Kormendy & Ho 2013).Feedback in the form of powerful mass outflows is thought to energetically connect these two spatially independent regions (Laha et al. 2021), as the propagation through the interstellar medium of the so-called ultra-fast outflows generates forward and backward shocks (e.g., Faucher-Giguère & Quataert 2012) that can enhance or quench the star formation in the host galaxy (e.g., Sturm et al. 2011;Zubovas et al. 2013;Zubovas & King 2014;Cresci et al. 2015).In AGN accreting with high efficiency, powerful winds with high column density (N H > 10 23 cm −2 ), ionisation (log ξ > 4) 1 , and velocity (−0.4c < ∼ v out < ∼ −0.1c) are expected from both theoretical arguments (e.g., King & Pounds 2003) and simulations (e.g., Proga & Kallman 2004).From an observational perspective, these outflows imprint X-ray absorption features at blueshifted energies, usually in the Fe-K band (i.e., Fe xxv and Fe xxvi).
Despite plenty of studies reporting on ultra-fast X-ray winds in nearby Seyferts and quasars, several open questions remain, such as the launching and driving mechanisms, and the geometry and duty cycle of these outflows.In particular, the detection ⋆ riccardo.middei@ssdc.asi.it 1 The ionisation parameter, defined as ξ = L ion /nr 2 (erg cm s −1 ), depends on the ratio between the ionising flux received by the gas and its density.For simplicity, its units are omitted throughout the text.rate of Fe-K absorption in relatively large samples of local AGN has been often used as an indirect measurement of the covering factor of disc winds, and these ensemble-based estimates in the range 0.4-0.6 (Tombesi et al. 2010;Gofford et al. 2013, see also Igo et al. 2020) suggest the wind to cover about half of the sky as seen by the X-ray continuum source.This is a key issue, as the degree of collimation of the ejected matter, combined with the sheer kinetic power, is a critical parameter to achieve an effective feedback on large scales.In the case of outflows with such a wide opening angle, a typical P-Cygni spectral feature would be expected (e.g., Castor & Lamers 1979).However, this distinctive emission-absorption profile is exceptionally rare in the X-ray spectra of AGN, so that a clear detection of a P-Cygni feature, or claims thereof, only exist for a few objects, including PG 1211+143 (Pounds & Reeves 2009), PDS 456 (Nardini et al. 2015), and I Zw 1 (Reeves & Braito 2019).Tentative identifications of this profile were also reported in some high-redshift objects (e.g., Vignali et al. 2015;Dadina et al. 2016).
The picture is further complicated by the possible spectral degeneracy between outflow and reflection signatures.Indeed, even in the presence of blueshifted Fe-K absorption, any emission counterpart can be also associated with the reflection of the primary X-ray radiation off the surface of the accretion disc.This makes it very difficult to disentangle the relative disc and wind contributions and so extract any sensible physical information from the observed line profile (Parker et al. 2022).It is worth noting that such a degeneracy works in both directions, as the extreme physical parameters implied by the strongest disc lines, among which a largely supersolar Fe abundance, can be significantly softened by allowing for the presence of a disc wind (e.g., Hagino et al. 2016).
A high Eddington rate is believed to imply favourable conditions for a fast X-ray wind to develop (e.g., Reynolds 2012), as the optical depth to Thomson scattering is high enough for continuum radiation pressure to sustain an outflow without requiring any substantial force multiplier associated to line absorption, which is negligible at high ionisation (but see Hagino et al. 2015).The importance of accretion rate and radiation driving is also revealed by the correlation between outflow velocity and Xray luminosity found by Matzeu et al. (2017) in PDS 456, while in NGC 2992 the wind signatures only emerge when the Eddington rate exceeds ∼2% (Marinucci et al. 2018).We refer to Giustini & Proga (2019) for a more general, holistic picture of the various wind manifestations as a function of the accretion rate.We note, however, that large uncertainties affect the estimate of the latter parameter, and it is therefore desirable to rely on BH masses derived via reverberation mapping techniques.
The ideal targets to gain deeper insights into the physics of accretion-disc winds are therefore AGN for which a sensible determination of BH mass and accretion rate is available, and possibly showing evidence of outflows also in other gas phases of lower ionisation.In this context, we revisit two archival XMM-Newton observations of the Seyfert 1 galaxy Mrk 1513 (RA = 21h 32m 27.8s, Dec = +10 • 08 ′ 19 ′′ ; also known as PG 2130+099 or UGC 11763), at redshift z = 0.063 (Springob et al. 2005).For this AGN, Grier et al. (2017) inferred a BH mass of log M BH = 6.92 +0.24  −0.23 M ⊙ by modelling the reverberation mapping data of the broad Hβ component (see also Hu et al. 2020).By applying the bolometric correction of Duras et al. (2020) to the hard X-ray luminosity of Mrk 1513, L 2−10 keV = 3.5 × 10 43 erg s −1 , we derive an accretion rate of L bol /L Edd ∼ 0.5.A nearly identical value is obtained by considering the 5100-Å luminosity with the standard bolometric correction factor of 9 (after Richards et al. 2006).Mrk 1513 can then be safely assumed to belong to the high-efficiency accretion regime.
Contrary to the extensive optical monitoring, only a few studies focussed on this AGN in the X-rays.Cardaci et al. (2009) reported on the presence of two warm-absorption components with moderate ionisation and column densities of a few ×10 21 cm −2 , mostly responsible for absorption by the Fe M-shell unresolved transition array (e.g., Sako et al. 2001).Both components appear to be inflowing at 500±300 km s −1 , yet the same authors caution against over-interpreting this result as the low signal-tonoise degraded the nominal resolution of the Reflection Grating Spectrometers (RGS, den Herder et al. 2001) onboard XMM-Newton.Interestingly, the Hubble Space Telescope (HST) Cosmic Origins Spectrograph data of Mrk 1513 show an intrinsic narrow absorption line (NAL) system with an outflow velocity v out = −1521 ± 20 km s −1 , seen in Lyα and high-ionisation species like N v, Si iv, and C iv (Tofany et al. 2014, and references therein), implying that a stream of hot gas is ejected at high velocity from the nuclear regions of this AGN.Here we discuss the intriguing possibility that Mrk 1513 actually hosts even an ultra-fast disc wind.This letter is organised as follows: Section 2 concerns the XMM-Newton observations and data reduction, Section 3 describes the results of the spectral analysis, while our conclusions are drawn in Section 4.

XMM-Newton observation and data reduction
Mrk 1513 was first observed with the EPIC pn and MOS cameras (Strüder et al. 2001;Turner et al. 2001, respectively) on-board XMM-Newton (Jansen et al. 2001) on 2003 May 17 for 37.5 ks (ObsID 0150470701).This remains the highest quality X-ray observation available for this AGN to date.Indeed, the source was also serendipitously observed on 2015 May 1 for 30 ks, yet in the latter data set (ObsID 0744370201) Mrk 1513 lies about 5 ′ off-axis and falls in the gap between two chips of the pn CCD.As this unfortunate position significantly affects the usefulness of the 2015 observation for the present analysis (Table 1), we will mainly focus on the 2003 spectrum.However, the main features discerned in the Fe-K band are remarkably similar in the two epochs, at least qualitatively (see Section 3.4).
The two XMM-Newton exposures were reprocessed using the standard pipeline within the Science Analysis System (sas, version 20.0.0) and the latest calibration files (March 2023).Except for the 2015 MOS data, all the other exposures are somewhat spoiled by some background flares.These periods were filtered via an iterative process aimed at maximising the S/N ratio in the 4-9 keV energy band (e.g., Piconcelli et al. 2004;Nardini et al. 2019).We extracted the source spectra from circular regions with radius of 30 ′′ (2003) or 35 ′′ (2015) for the pn, and 25 ′′ for both MOS detectors.The background regions were selected in order to avoid nearby sources and possible instrumental features (due to, e.g., Cu-K fluorescence), with a circular shape 2 and radii between 60 ′′ and 75 ′′ .The basic details of each observation are listed in Table 1.Response files were produced with the sas tools rmfgen and arfgen.Finally, the obtained spectra were binned in order to have at least 5 counts per bin.The sas routine epicspeccombine was used to sum the EPIC spectra together.The output files, however, are only employed for plotting purposes, unless stated otherwise.
In Fig. 1, we show the EPIC spectra of Mrk 1513 derived for the 2003 XMM-Newton observation, unfolded against a power law with Γ = 2 and normalisation equal to unity.A sharp spectral drop is clearly visible at ∼7 keV.A similar feature is present during the 2015 observation, despite the lower hard X-ray statistics and a different broadband spectral shape (see below).

Spectral analysis
The spectral analysis presented here was performed with the fitting software package xspec (Arnaud 1996).In all the fits the C-statistic (Cash 1979) was adopted and Galactic column density is always included via the tbabs model (N H = 3.7 × 10 20 cm −2 , HI4PI Collaboration et al. 2016).We also used a constant component free to vary to account for the cross-normalisation between the different EPIC detectors, which is within 5%.Cardaci et al. (2009) analysed the 2003 XMM-Newton observation of Mrk 1513 and found the RGS spectra to reveal the presence of warm-absorption features.Two mildly ionised components with a marginal redshift (hence inflowing) were statistically preferred over a single, blueshifted one (at v out ≃ −300 km s −1 ).As our aim is to investigate the nature of the spectral structures in the Fe-K band, and the velocity of these warm absorbers cannot be constrained at CCD resolution, in our analysis we assumed that both components are at rest in the local frame.Then, both the column density and the ionisation parameter ξ were calculated during the fitting procedure.In any case, the Fe-K band is insensitive to the best-fit properties of the warm absorbers.
The sharp drop observed at 7 keV in the XMM-Newton spectra of Mrk 1513 can be alternatively interpreted as due to either the blue horn/wing of a relativistic disc line, or an Fe-  K absorption edge (e.g., Boller et al. 2002).To better visualize the shape of this emission-absorption structure, we fitted the 0.3-10 keV spectrum with the following xspec model: tbabs × xabs × xabs × (bbody + powerlaw), where the two xabs components3 account for the warm absorber, and a phenomenological blackbody and a power law are used to reproduce, respectively, the soft excess and the hard nuclear continuum (compare with Fig. 1).All the parameters of the blackbody and power-law components were left free to vary.This simple model returns a statistics of C-stat = 3222 for 3131 degrees of freedom (d.o.f.) and, unsurprisingly, it leaves prominent residuals above 6 keV.We then carried out a blind line scan by adding to the underlying continuum model a Gaussian profile with fixed width of 150 eV (i.e., compatible with the CCD resolution).The line energy and normalisation were allowed to vary with 50 steps from 5.5 to 9.5 keV and from −1.5 to 1.5 × 10 −5 photons cm −2 s −1 , respectively.
The outcome of this procedure is shown in Fig. 2, which reveals two broad features, one in emission and the other one in absorption, whose combination is strongly reminiscent of a P-Cygni profile.However, if the continuum is more complex than assumed in our baseline model, or moderately bent at higher en-Fig.2. Results of the blind scan for the line normalisation and energy based on a phenomenological continuum model (see the text for more details).The contours correspond to the 68%, 90%, 99%, and 99.9% confidence levels.The proximity in energy of positive and negative residuals is strongly reminiscent of a P-Cygni profile.
ergies, the same residuals could also be explained by a single relativistic emission line.We thus explore both the wide-angle wind and the disc reflection scenarios in more depth below.

An ultra-fast disc wind?
According to the line scan, the 7-keV spectral drop is apparently due to the presence of a broad absorption trough immediately blueward of a prominent emission feature, whose peak is compatible with fluorescence from iron in a low-ionisation state.The trough, which is tempting to associate with an accretion-disc wind, possibly has a complex structure with two dips astride 8 keV, although the higher energy one is only marginally significant.Based on Fig. 2, we added to the baseline continuum two Gaussian components with, respectively, positive and negative amplitude, which yield a statistical improvement of ∆Cstat/∆d.o.f.≈ −45/−6.The probability of chance improvement is < 10 −7 , and the aggregate feature is significant at > 5σ.Yet, the centroids and widths of both lines are hard to constrain simultaneously in a reasonable way, as they are all highly degenerate with each other; in fact, the two profiles substantially overlap in order to reproduce the steepness of the flux drop as well as the mild asymmetry of the emission and absorption tails.
Scattering by outflowing matter that covers a large fraction of the solid angle can give rise to intense Fe-K emission at lower energies than the absorption trough (e.g., Nardini et al. 2015).As the Gaussian-line fit is poorly informative, we thus tested a self-consistent P-Cygni profile adopting the pcygx model, customized for xspec by Done et al. (2007) after (Lamers et al. 1987).This model mimics the simplified scenario of an expanding, spherically symmetric wind with a velocity field of the form w = w 0 + (1 − w 0 )(1 − 1/x) γ , where w = v/v inf is the velocity of the wind normalised to the terminal velocity, w 0 is the dimen-sionless velocity of the wind photosphere, x represents the ratio between the radial distance r and the radius of the photosphere, and γ sets the scaling relation between distance and velocity of the wind.Both w 0 and γ have negligible impact on the fits, hence we fixed their values to 0.01 and 1, respectively.Under these assumptions, the line optical depth can be written as τ(w) ∝ τ tot w α 1 (1− w) α 2 , where α 1 and α 2 regulate the P-Cygni profile smoothness (see Ramírez et al. 2005).We also assumed for simplicity that α 1 = α 2 = 0, and considered as free parameters the rest-frame energy of the line E 0 , the terminal velocity of the wind v inf , and the total optical depth of the gas τ tot .
This model provides a satisfactory representation of the Fe-K spectrum of Mrk 1513, with a statistics comparable to the double-Gaussian model but three more degrees of freedom, i.e., C-stat/d.o.f.= 3176/3128.The line rest energy is fully consistent with Fe xxv Heα at E 0 = 6.72 +0.11  −0.14 keV, the gas optical depth is τ tot ∼ 0.2, and the terminal velocity is v inf ∼ 0.34c.All the bestfit parameters are reported in Table 4, while best-fit model and residuals are shown in Fig. 3 (top and middle panels).

A reflection-dominated spectrum?
Fluorescence, scattering, and recombination are known to occur also from iron inside the disc, hence the emission feature at 6-7 keV is not an unambiguous signature of a wide-angle wind (see Parker et al. 2022).We therefore replaced the P-Cygni component with an emission line from a relativistic accretion disc, modelled through the diskline profile in xspec (e.g., Fabian et al. 1989).The free parameters of this component are the line energy and normalisation, and the disc inner radius and inclination.We eventually assumed a disc emissivity scaling with radius as (1 − √ 6/R)/R 3 , where R is in units of gravitational radii (r g = GM BH /c 2 ), as this is preferred by the fit.
As shown in Fig. 3 (top and bottom panels), this model delivers an equally good fit to the XMM-Newton spectrum of Mrk 1513, with C-stat/d.o.f.= 3173/3127.The best-fit parameters are listed in Table 4.The line rest energy of E 0 = 6.48 +0.06 −0.03 keV suggests that the gas in the disc is only moderately ionised.The disc is seen at intermediate inclination, i ∼ 40 • , and extends down to the innermost stable circular orbit of a non-rotating black hole.Even without taking the line and disc parameters at face value, it is clear that also a relativistic line can effectively describe the complex 7-keV spectral structures observed in Mrk 1513.We note, however, that any disc line is required to be very strong, with an equivalent width of EW = 330 ± 120 eV, whereas the average value for broad Fe-K lines is ∼100 eV (e.g., Nandra et al. 2007;Patrick et al. 2012).This entails strong light bending to efficiently illuminate the disc and suppress the fraction of direct continuum that reaches the observer at infinity (Miniutti & Fabian 2004).Hence, we would expect the broadband X-ray spectrum of Mrk 1513 to be reflection-dominated, and we should be able, in principle, to account for both the Fe-K feature and the excess emission in the soft band within a relativistic reflection framework.To this aim, we performed several tests by taking advantage of the relxill suite (e.g., García et al. 2014;Dauser et al. 2016), which allows one to reproduce the primary and reflected emission for various physical and geometrical configurations of the disc-corona system.Rest-frame energy (keV) Fig. 3. Top panel: pcygx (magenta) and diskline (cyan) best-fit models, superimposed for reference to the combined EPIC spectrum (rebinned for clarity).Note that the effects of the warm absorber on the Fe-K band are negligible in either case.Middle panel: fit residuals with respect to pcygx model, computed as (data−model)/error for the individual EPIC spectra (colour-coded as in Fig. 1).Bottom panel: same as above for the diskline model.

Broadband relativistic reflection models
The relxill package provides a set of spectral models in which the relativistic reflection from an accretion disc illuminated by the primary X-ray is self-consistently calculated during the fitting procedure.Given the strength of the disc line required to explain the 7-keV drop, we tested whether any model in the relxill family allows us to account for the entire 0.3-10 keV spectrum of Mrk 1513.We started by adopting the lamppost configuration for the primary X-ray source, as rendered through the following model expression: tbabs × xabs × xabs × relxilllp, where the latter component is meant to incorporate both the direct and the reprocessed flux under the assumed geometry.We performed the fit computing the photon index Γ of the primary continuum, the height of the corona above the disc (h, which is also set to control the reflection fraction), and the ionisation parameter ξ d , inclination i, and iron abundance A Fe of the disc.The SMBH spin a * was also left as a free parameter under the hypothesis that the disc extends down to the innermost stable circular orbit.The continuum high-energy cutoff was frozen to the value of 300 keV, well beyond the XMM-Newton operating energy range, and the two warm-absorption components were treated as described above.The relxilllp model yields a statistics of C-stat/d.o.f.= 3245/3128, and it is moderately worse than the phenomenological continuum model.In particular, P-Cygni-like residuals remain in the Fe-K band (Fig. 4, top panel).The inferred best-fit quantities are reported in Table 3.
This poor result can be ascribed to the lack of physical and geometrical flexibility of the relxilllp model.We then replaced it with the relxillD flavour, where the density of the disc is no more fixed at log (n/cm −3 ) = 15 but becomes a free parameter.Moreover, the lamppost configuration is dismissed, allowing the disc emissivity to follow a broken-power-law dependence on radius (for simplicity, we assumed a trend ∝ r −q with no break) and the reflection fraction R (see Dauser et al. 2016) to vary.With the emissivity exponent and reflection fraction as new free parameters in place of the coronal height, and the addition of the disc density, the fit improves down to Cstat/d.o.f.= 3194/3126 (Table 3).The residuals in the Fe K band are somewhat attenuated but not completely offset (Fig. 4, bottom panel).This suggests that, for the 7-keV feature to have a reflection origin, a different physical process must be responsible for the bulk of the soft X-ray emission in Mrk 1513.Indeed, the fit is invariably driven by the soft excess, which calls for a degree of relativistic blurring that is far too large to simultaneously produce a sharp drop as the one seen around 7 keV in Mrk 1513.We then argue that none of the broadband reflection-based models, while largely acceptable in terms of overall statistics, are able to account for the Fe-K spectral structure, leaving positive and negative residuals with a P-Cygni-like shape.

The 2015 XMM-Newton observation
In the 2015 observation, the source lies in the gap between two chips of the EPIC-pn CCD.This translates into an extremely poor data quality above 5 keV, where the pn camera collected nearly four times less counts than in the 2003 exposure (see Table 1).In practice, the effective area of the pn in the Fe-K band  in 2015 was equivalent to that a single MOS detector.Still, also this second XMM-Newton observation provides some qualitative information about the possible Fe-K structure in this source.
In Fig. 5 we show the 2003 and 2015 EPIC spectra obtained by combining the individual pn and MOS spectra using the sas routine epicspeccombine.The source is found in a strikingly different state, denoted by a much more prominent soft excess.
On the other hand, the high-energy part of the 2015 spectrum shows milder changes with respect to the 2003 one.Remarkably, a deficit of net counts can be tentatively observed around ∼8 keV also in the second pointing of Mrk 1513.There are again hints of a P-Cygni profile, although the emission and absorption component are suggested to be more detached by the blind line search, performed on the combined EPIC spectrum (Fig. 6).Owing to the limited data quality, the statistical significance of both features is marginal, approaching the ≈95% or 2σ level.Despite the low statistics in the Fe-K band, we attempted to apply the pcygx and diskline models to the 2015 EPIC spectrum.The enhanced emission in the soft band was accounted for by adding a second phenomenological blackbody to the (red) observations.The soft X-ray spectrum shows a significant change, while the source is found in a similar spectral and flux state in the Fe-K band.The edge-like feature at ∼7-8 keV seems to be persistent, with an energy that possibly evolves with time.4.

Discussion and Conclusions
We have reported on the X-ray spectral properties of the highly accreting AGN Mrk 1513, with specific focus on the peculiar edge-like feature discernible at ∼7 keV (Fig. 1).Despite the short duration (37.5 ks) of the 2003 XMM-Newton observation, such a complex Fe-K emission and/or absorption structure is seen in the data from all the three EPIC detectors, and is tentatively present also in a lower quality spectrum that was serendipitously ac- We have shown that, on sheer statistical grounds, it is not possible to discriminate between a wide-angle, ultra-fast wind and a relativistic accretion-disc line as the origin of the 7-keV feature.This notwithstanding, we can rely on some physical considerations to lean towards the most likely scenario.Indeed, in our analysis we have adopted a simplified phenomenological continuum, consisting of a soft X-ray blackbody and a hard X-ray power law, only modified by the known warm absorbers.However, the involved strength of any disc line, irrespective of its exact parameters, would suggest that X-ray reflection heavily contributes to the broadband spectral shape.When we switch to a physically motivated model for the disc-line interpretation, we find that the sharpness of the 7-keV feature is at odds with the overall spectral smoothness.This argument is not yet sufficient in itself to completely rule out the disc-line scenario, though, since a different physical component, such as warm Comptonisation, can partly or entirely account for the soft excess (e.g., Middei et al. 2019;Matzeu et al. 2020).
The nature of the broadband continuum is not as severe an issue in the disc-wind interpretation instead, as any emission/absorption wind signatures would be imprinted on the spectrum from the background X-ray source regardless of its shape and/or physical origin.In other words, the source of the continuum and that of the Fe-K feature(s) do not need to be one and the same.Hence, the P-Cygni solution remains fully viable at this stage, even in the absence of a self-consistent characterization of the entire spectrum.Interestingly, there is the hint of a change in the power-law photon index between the pcygx and diskline models, which is slightly higher (1.70 against 1.66; Table 2) in the latter case.The steeper slope partly compensates for the lack of an absorption component, and suggests that the continuum trend above 10 keV could help discriminating between the two models.Unfortunately, with the current data quality, the number of significant spectral bins beyond the 7-keV edge is scarce, and it is not possible to reliably anchor the continuum level blueward of the apparent absorption trough.A similar situation was encountered after the first XMM-Newton observation of the quasar PDS 456, where a broad trough was detected extending across the upper end of the EPIC operating bandpass (Reeves et al. 2003).Only with the advent of NuSTAR, which eventually probed the high-energy continuum, PDS 456 has been conclu-sively established as the prototype of wide-angle disc winds in AGN (Nardini et al. 2015).
The disc-wind explanation is also appealing because of the unmistakable indications of a P-Cygni profile.Although the detection rate of blueshifted Fe-K features in the X-ray spectra of AGN is fairly high, leading to the notion that such outflows must cover a sizeable fraction of the solid angle, only a handful of P-Cygni profiles have been actually discovered so far.This, on the one hand, casts some doubts on the commonly accepted covering factor of ultra-fast X-ray winds among AGN and, on the other hand, prevents a better understanding of how the matter and energy flux from the central regions couples with the largescale environment.At the time of the 2003 XMM-Newton observation analysed here, a persistent outflow of highly ionised gas at v out ∼ −1500 km s −1 was present in Mrk 1513, as revealed by the 1995, 1999, and 2010 HST observations (Tofany et al. 2014).The 2015 XMM-Newton observation, despite its poor quality, provisionally suggests that also the X-ray wind might be persistent on similar timescales.Confirming the physical connection between the UV and X-ray outflows (e.g., Mehdipour et al. 2023) would shed new light on the relation between efficient accretion, disc winds, and large-scale feedback in AGN.Our analysis has demonstrated that Mrk 1513 is potentially a very promising source to investigate the AGN outflow phenomenon and, as such, it also stands out as a suitable target for the newly launched XRISM mission.Besides probing any possible fine structure in the Fe-K emission/absorption profiles, the high-resolution view would pinpoint the properties of the warm absorbers and reveal any residual continuum curvature, helping us to disentangle the contribution from disc (or even distant) reflection.Longer X-ray exposures and, ideally, broader spectral coverage of this source are also desirable to achieve an accurate characterization of the disc wind properties and so determine its possible impact on the host galaxy.

Fig. 1 .
Fig. 1.EPIC-pn (black), MOS1 (red), and MOS2 (green) spectra from the 2003 XMM-Newton observation of Mrk 1513, unfolded against a power law with Γ = 2.Besides the soft excess component below 1 keV, all the spectra are characterized by a sharp drop around 7 keV.The emission deficit up to ∼9 keV is clearly noticeable in the pn data, but also confirmed by the two MOS detectors despite their rapidly declining effective area at these energies.

Fig. 4 .
Fig. 4. Fit residuals for the relxilllp (top) and relxillD (bottom) relativistic reflection models applied to the 2003 XMM-Newton spectrum of Mrk 1513.The colour code is the same as in Fig. 1.

Fig. 5 .
Fig. 5. Unfolded EPIC spectra of Mrk 1513 in the 2003 (black) and 2015(red) observations.The soft X-ray spectrum shows a significant change, while the source is found in a similar spectral and flux state in the Fe-K band.The edge-like feature at ∼7-8 keV seems to be persistent, with an energy that possibly evolves with time.

Fig. 6 .
Fig. 6.Same as Fig. 2 for the 2015 XMM-Newton observation of Mrk 1513.The emission component and the secondary dip in the absorption profile appear to be rather stable in energy compared to 2003.The contours correspond to the 68% and 90% confidence levels, and the significance of both features is ≈2σ.

Table 2 .
Best-fit parameters derived for the alternative models adopted to reproduce the 7-keV spectral feature(s) of Mrk 1513.Errors are given at the 90% confidence level for the single parameter of interest.

Table 3 .
Best-fit parameters derived for the different relativistic reflection models applied to the 2003 EPIC spectra of Mrk 1513.

Table 4 .
Best-fit parameters derived for the pcygx and diskline models adopted to reproduce the hard X-ray spectral feature(s) of Mrk 1513 as observed in 2015.Errors are given at the 90% confidence level for the single parameter of interest.Appendix 3.4).To date, these are the only observations of Mrk 1513 carried out with the latest generation of X-ray facilities.