Tracking the spectral properties of ESO 511-G030 across different epochs

The Type 1 active galactic nucleus (AGN) ESO 511-G030, a formerly bright and soft excess dominated source, was observed in 2019 in the context of a multi-wavelength monitoring campaign. In the new exposures, the source was found to be in a ﬂux state approximately ten times lower than archival exposures and without any trace of the soft excess. Interestingly, the X-ray weakening observed in the 2019 data corresponds to a comparable fading of the UV ﬂux, suggesting a strong link between these two components. The UV-X-ray spectral energy distribution (SED) of ESO 511-G030 shows remarkable variability. We tested both phenomenological and physically motivated models on the data, ﬁnding that the overall emission spectrum of ESO 511-G030 in this extremely low ﬂux state is due to the superposition of a power-law-like continuum ( Γ ∼ 1 . 7) and two reﬂection components emerging from hot and cold matter. Both the primary X-ray continuum and relativistic reﬂection are produced in the inner regions close to the supermassive black hole. The prominent variability of ESO 511-G030 and the lack of a soft excess can be explained by the dramatic change in the observed accretion rate, which dropped from an L / L Edd of 2% in 2007 to one of 0.2% in 2019. The X-ray photon index also became harder during the low ﬂux observations from 2019, perhaps as a result of a photon starved X-ray corona.


Introduction
The broadband emission in active galactic nuclei (AGNs) can be interpreted as an interplay between thermal and non-thermal processes taking place in the close surroundings of a central supermassive black hole (SMBH; see Padovani et al. 2017, for an overview on AGNs).Accretion onto the SMBH is responsible for the optical-UV emission, and a fraction of the thermal photons from the emission is intercepted and Comptonised up to the X-rays by the so-called hot corona (e.g.Galeev et al. 1979;Haardt & Maraschi 1991, 1993).Our understanding of this plasma is still lacking, though timing and microlensing arguments (e.g.Chartas et al. 2009;Morgan et al. 2012;De Marco et al. 2013;Kara et al. 2016) support the hypothesis of this component being compact and located in the inner regions of the accretion flow.
The presence of this hot plasma well explains the cutoff power-law-like continuum of AGNs, with the high-energy rollover being interpreted as a further signature of the inverse-Compton of the seed photons from the thermal plasma of electrons with E c , depending on the temperature of the relativistic electrons (kT; Rybicki & Lightman 1979).The characterisation of the primary X-ray continuum has been the focus of several studies (e.g.Perola et al. 2000;Dadina 2007;Molina et al. 2009Molina et al. , 2013;;Malizia et al. 2014;Ricci et al. 2018), and research in this area has been further boosted since the launch of NuSTAR (Harrison et al. 2013), as demonstrated by the increasing number of high-energy cut-off measurements (e.g.Baloković et al. 2020;Reeves et al. 2021;Kamraj et al. 2022).
Reflection off Compton-thin and Compton-thick matter imprints additional features onto the emerging X-ray spectrum.This is the case of the fluorescence Fe Kα emission line at 6.4 keV (e.g.Bianchi et al. 2009) and the Comptonhump (George & Fabian 1991;Matt et al. 1993;García et al. 2014).Noticeably, the analysis of the Fe Kα profile carries a wealth of information on the location of the reflecting materials.In fact, its intrinsically narrow profile can undergo distortions, resulting in a broader shape due to relativistic effects (e.g.Fabian et al. 1989Fabian et al. , 2000;;Tanaka et al. 1995;Nandra et al. 2007;de La Calle Pérez et al. 2010).
Finally, the soft X-ray band of AGNs ubiquitously shows a bump of counts below ∼2 keV that is not accounted by the high-energy power-law continuum (e.g Piconcelli et al. 2005;Bianchi et al. 2009;Gliozzi & Williams 2020).The origin of this component, the so-called soft excess, is still debated, and two possible scenarios have been tested using various datasets: blurred ionised reflection and warm Comptonisation (e.g.Crummy et al. 2006;Magdziarz et al. 1998;Jin et al. 2012;Done et al. 2012).The two models have been found to reproduce the data.Walton et al. (2013) tested the relativistic reflection origin of the soft excess on a broad number of Seyfert galaxies using Suzaku data (see also Crummy et al. 2006), while the two-coronae model was found to be the best-fit for the broadband emission spectrum of an increasing number of AGNs (e.g.Petrucci et al. 2018;Porquet et al. 2018;Kubota & Done 2018;Ursini et al. 2018Ursini et al. , 2020;;Mahmoud & Done 2020;Matzeu et al. 2020;Middei et al. 2020).
Variability is another key feature of AGN emission (Bregman 1990;Mushotzky et al. 1993;Wagner & Witzel 1995;Ulrich et al. 1997).Spectral variations are characterised by a softer-when-brighter behaviour that is commonly observed in nearby Seyfert galaxies and unobscured quasars (e.g.Sobolewska & Papadakis 2009;Serafinelli et al. 2017, respectively).X-ray amplitude variations have been witnessed over different time intervals, from hourly (Ponti et al. 2012) up to yearly changes (Vagnetti et al. 2016;Paolillo et al. 2017;Middei et al. 2017;Falocco et al. 2017;Gallo et al. 2018;Timlin et al. 2020).X-ray variability has been found to anti-correlate with the source luminosity; however, this can be naturally explained if changes result from the superposition of N randomly emitting sub-units.This scenario, already considered in optical studies (e.g.Pica & Smith 1983;Aretxaga et al. 1997), predicts a variability amplitude ∝N −0.5 ∝ L −0.5 and accounts for the sub-units being identical and the flare being independent (e.g. Green et al. 1993;Almaini et al. 2000).Interestingly, this behaviour has been reported by many authors for both local and high-redshift AGNs also in the X-rays (e.g.Barr & Mushotzky 1986;Lawrence & Papadakis 1993;Papadakis et al. 2008;Vagnetti et al. 2011Vagnetti et al. , 2016)).
In this context, we report on the X-ray spectral properties of ESO 511-G030, a nearby (z = 0.02239±0.00001,Theureau et al. 1998) spiral galaxy (de Vaucouleurs et al. 1991) hosting an unobscured Type 1 AGN (Véron-Cetty & Véron 2010).One of the first X-ray spectra of ESO 511-G030 was taken with ASCA (Advanced Satellite for Cosmology and Astrophysics) in 1998 (Turner et al. 2001b), and the source spectrum was best-fitted using an absorbed power-law.A subsequent INTEGRAL-based analysis allowed for measurement of a high-energy cut-off E c = 100 +101 −37 keV (Malizia et al. 2014).The AGN ESO 511-G030 is reported in the 105-month BAT catalogue (Oh et al. 2018) with a flux of ∼4×10 −11 erg s −1 cm −2 (14-195 keV), and it is one of Notes. (†) Exposures from the joint XMM-Newton-NuSTAR monitoring campaign.
the 13 objects in the FERO sample, an XMM-Newton-based collection of AGNs with a >5σ detection of a relativistic iron line (details in de La Calle Pérez et al. 2010).In a recent paper by Ghosh & Laha (2021), who studied a 2007 XMM-Newton exposure, the source flux was consistent with F 2−10 keV ∼ 2 × 10 −11 erg s −1 cm −2 .In the same paper, two Suzaku observations five years apart from the XMM-Newton one are also discussed.In those observations, the source flux was compatible with the XMM-Newton 2007 observation, and the second observation was about 50% higher.Moreover, ESO 511-G030 does not show the presence of cold and/or warm absorption components (e.g.Laha et al. 2014), leaving the soft X-ray band free from complex absorption features and with only a modest attenuation of the UV emission.This paper is organised as follows: Sect. 2 reports on data reduction, and in Sect. 3 the timing properties of the observations are discussed.In Sect.4, we describe the spectral analyses of the XMM-Newton-NuSTAR 2019 monitoring campaign, and we test the same spectral model on the archival data in Sect. 5. Then in Sects.6 and 7, we describe the broadband Swift data analysis, including XRT (Swift's X-Ray Telescope) and UVOT (Swift's Ultraviolet/Optical Telescope) data.In Sect.8, we test a self-consistent model on the 2019 and the 2007 XMM-Newton exposures also considering optical monitor (OM) data.Our conclusions and comments are reported in Sect.9.

Data reduction
In this section, we present the analysis of our multi-wavelength XMM-Newton/NuSTAR observational campaign from 2019, and we compare it with data taken using different facilities across a time interval of more than 20 years.In Table 1, the log of the observations is reported, while in Fig. 1 we show a snapshot of all spectra simply folded with a power-law (Γ = 2), from which it is possible to witness the variations of ESO 511-G030 in both the X-ray and optical-UV bands.
ASCA : We retrieved the already reduced data products of ESO 511-G030 from the Tartarus ASCA AGN database (Turner et al. 2001b).1. Variability is observed both in terms of spectral and amplitude changes.In the optical-UV band, the source faded by a factor of approximately ten.The underlying model folding the data is a power-law with Γ = 2 and a unitary normalisation.
Suzaku: ESO 511-G030 was observed with Suzaku (Mitsuda et al. 2007) on July 20 (OBSID: 707023010), July 22 (707023020), and August 6 (707023030), 2012, through the X-ray Imaging Spectrometer (XIS; Koyama et al. 2007) with net exposure times of 5.7 ks, 224 ks, and 51 ks, respectively.Following the processes described in the Suzaku data reduction guide1 , the XIS (X-ray Imaging Spectrometer) 0, 1, and 3 CCD spectra were extracted with HEASOFT (v6.29.1) and the latest version of the CALDB (November 2021).The cleaned event files were selected from the 3×3 and 5×5 edit modes and subsequently processed according to the suggested screening criteria.Both XIS source and background spectra were extracted from circular regions with radii of 3.0 arcmin.Care was also taken to avoid the chip corners containing Fe 55 calibration sources.The corresponding spectra and light curves were subsequently extracted using xselect for the second and third observations, as the first pointing was too short.For each detector, response matrices and ancillary response files were generated by running the xisrmfgen and xissimarfgen tasks.After verifying their consistency, we combined the front illuminated (XIS-FI) 0 and 3 spectra into a single XIS03-FI spectrum for observations 707023020 and 707023030.We used a cross-calibration constant to account for the inter-calibration between the XIS and Pin detectors.In the fits, the constant has its expected value of k ∼ 1.16 for the 707023030 data only, and it increases to a value of k ∼ 1.50 in the 707023020 observation.Such a particularly high value for this constant is explained by the non-simultaneity of the XIS-PIN exposures, which caused telemetry issues to occur and led to a shortening of the actual PIN exposure to about one-fifth of what was scheduled.Thus, the high value of the cross-correlation constant is straightforwardly explained by the intra-observation variability that the source had undergone.
XMM-Newton: We reduced and analysed both a 120 ks 2007 XMM-Newton orbit and five exposures of about 30 ksec, all of which were obtained simultaneously with NuSTAR in 2019.The exposures were made with the EPIC camera (European Photon Imaging Camera, Strüder et al. 2001;Turner et al. 2001a) operating in the Small Window mode.Data were processed using the XMM-Newton Science Analysis System (SAS; Version 19.0.0).Because of its larger effective area with respect to the MOS (Metal Oxide Semi-conductor) cameras, we only report the results for the pn instrument.Source spectra were derived using a circular region with a 40 arcsec radius centred on the source, while the background was extracted from a blank 50-A101, page 3 of 17 arcsec radius area near the source.The extraction regions were selected using an iterative process that maximises the S/N, similar to what is described in Piconcelli et al. (2004).The spectra were re-binned in order to have at least 30 counts for each bin and to not over sample the spectral resolution by a factor greater than three.Finally, from 'epatplot', pile up issues were found to not affect this dataset.
NuSTAR: We calibrated and cleaned raw NuSTAR (Harrison et al. 2013) data using the NuSTAR Data Analysis Software (NuSTARDAS; Perri et al., 2013 2 ) package (v.1.8.0).Level 2 cleaned products were obtained with the standard 'nupipeline' task, while third level science products were computed with the 'nuproducts' pipeline and the calibration database 20191219.A circular region with a radius of 50 arcsec was used to extract the source spectrum.The background was calculated using the same circular region but centred on a blank area near the source.To account for the inter-calibration of the two modules carried on the NuSTAR focal plane, we used a cross-normalisation constant in all the fits.The calibration constant was always found to be within 3%, indicating the FPMA/B spectra to be in good agreement.Spectra were binned so that each bin had at least 50 counts so as to not over sample the instrumental resolution by a factor greater than 2.5.

Swift:
The Swift satellite observed ES0511-G030 from 2018 to 2021, and we reduced the data acquired with XRT and UVOT.The XRT observed the source in photon counting mode, and we derived the corresponding scientific products using the facilities provided by the Space Science Data Center (SSDC3 ) of the Italian Space Agency.In particular, spectra were extracted adopting a circular region of ∼60 arcsec centred on the source, and a concentric annulus was used for the background.The spectra were binned in order to have at least five counts in each bin.The UVOT aperture photometry was used to obtain the monochromatic fluxes for all the available filters.A source extraction region of 5 arcsec radius was adopted, and an appropriate blank annular region concentric with the source was adopted for the background.
All the errors quoted in our tables and text account for 90% uncertainties, while 68% errors are shown in the figures.Fits were performed using XSPEC (Arnaud 1996) and assuming the standard cosmological framework given by H 0 = 70 Km s −1 Mpc −1 , Ω Λ = 0.73, and Ω m = 0.27.

Timing properties
We extracted the background-subtracted light curves for all the observations.In the top panels of Fig. 2, we show the soft and hard X-ray light curves of the 2019 monitoring campaign time series, while those of archival observations (XMM-Newton and Suzaku) are shown in the bottom panels.From Fig. 2, both short-and long-term variability characterises the ESO 511-G030 light curves.During 2019, ESO 511-G030 had a quite stable behaviour, with both the soft and hard X-rays being fairly constant within each exposure and amongst the different pointings.This constancy, apart from a small fraction in the fifth exposure, can be observed through the ratios of the panels.This behaviour suggests that the balance between soft and hard X-rays did not change during the campaign.The flat shape of the 2019 hardness ratios can be quantitatively (qualitatively) compared with those computed from XMM-Newton (Suzaku) archival exposures.The ratios between the 0.5-2 and 2-10 keV bands were more variable in the 2007 XMM-Newton exposure, with changes of about 25%.Hardness ratios from Suzaku show moderate variations within the same exposure, fluctuating by approximately 10% on daily rather than monthly timescales.
Short-term X-ray variations are related to the intrinsic properties of the AGN, such as its SMBH mass or its luminosity (e.g.Vaughan et al. 2003;Papadakis 2004;McHardy et al. 2006).Ponti et al. (2012) computed the normalised excess variance of ESO 511-G030 using the 2007 XMM-Newton light curves, and the authors derived the black hole mass of ESO 511-G030 to be log M BH = 7.89 +0.30  −0.20 M .Another commonly adopted estimator suitable for X-ray variability characterisation is the fractional root mean square variability amplitude (F var ; e.g.Edelson et al. 2002;Vaughan et al. 2003;Ponti et al. 2004).The F var tool calculates the square root of the normalised excess variance, and it has been widely used to characterise the variability properties of AGN in X-rays (e.g.Vaughan et al. 2004;Ponti et al. 2006;Matzeu et al. 2016Matzeu et al. , 2017;;Alston et al. 2019;Parker et al. 2020;De Marco et al. 2020;Igo et al. 2020;Middei et al. 2020).We studied the variability properties of ESO 511-G030 by computing the F var spectra for each of the 2019 XMM-Newton observations.This spectral sample contained variability on timescales ranging between ∼1−40 ks.We used the background-subtracted light curves calculated in different energy intervals and adopted a temporal bin of 1000 sec.Following the same procedure, we also computed the F var spectrum of the 2007 XMM-Newton observation.These samples had longer timescales ranging between ∼ 1 − 120 ks and enabled us to identify variable spectral components contributing to the time-averaged spectrum.The resulting F var spectra of each observation are shown in Fig. 3.The errors were computed using Eq.(B2) of Vaughan et al. (2003) and account only for the uncertainty caused by Poisson noise.Aside from some excess towards the high-energy region of the spectra (e.g., due to residual background variability), all the observations from the 2019 monitoring campaign show a rather flat F var spectrum, therefore implying a similar variability power across different energy bands.In contrast, the 2007 F var spectrum clearly shows a divergence from the 2019 observations in the soft X-rays, which is consistent with the presence of an additional variability component.To better highlight this difference, we overlayed the 2007 F var spectrum with the grouped F var spectra from the 2019 campaign (see Fig. 4).We used publicly available table models4 to describe the 2007 F var spectrum in terms of combined contribution from flux variability of a power-law-like continuum and a soft excess.While the model well reproduces the 2007 F var spectrum and the high-energy part of the 2019 data, it clearly overestimates the soft band part of the 2019 F var spectrum.One possible explanation for this is the presence of an additional soft variability component in the 2007 data that is not present in the 2019 data.Fig. 3. F var spectra for the XMM-Newton observations in Table 1.Spectra from the 2019 monitoring campaign are all characterised by a fairly flat shape, while data from 2007 clearly show a variability bump in the soft X-rays.

The Fe Kα complex
We begun our investigation by focusing on the 2019 EPIC-pn spectra (between 4 and 8 keV) and the properties of the Fe Kα emission line.We adopted two components, a power-law and a Gaussian line that were assumed to be narrow (σ = 0 eV) and centred at 6.4 keV.We simultaneously fitted all the spectra, computing the photon index (we assumed its value to be the same amongst the pointings) and the power-law normalisation.
For the Gaussian component, we calculated its normalisation in all the exposures.This simple fit returned a Fe Kα flux consistent with being constant, Norm FeKα = (6.8± 0.8) × 10 −6 photons cm −2 s −1 .We then considered the 2007 XMM-Newton data on the 4-8 keV energy range.We tested the data with the same model and further assumed the Fe Kα to have the same normalisation from 2007 to 2019.In other words, we only fitted the photon index and the normalisation of the continuum for the newly added data (see top panel of Fig. 5).The narrow emission line only reproduced the 2019 spectra, and an additional broader component was required for the 2007 data.We thus froze the narrow Gaussian component to its best-fit value and added a new broad Gaussian emission line.In this fit, the line energy centroid and width were computed, and we assumed these values to be the same among the all observations.We only allowed the line's normalisation to vary between the 2019 dataset and the 2007 exposure.This new model resulted in the fit shown in the bottom panel of Fig. 5 with χ 2 /d.o.f.= 400/310.We subsequently re-fit the data, also allowing the narrow Gaussian normalisation to vary between the datasets.This test resulted in a slight benefit in terms of statistic, with ∆χ 2 /∆d.o.f.= −10/−1.In this case, however, only an upper limit was returned for the flux of the narrow Gaussian components.We further tested the origin of the Fe Kα emission line by computing the fit once again but only including a broad Gaussian component.Again, we assumed the line's energy centroid and width to be the same across the years, and only its normalisation was calculated separately for the 2019 and 2007 data.Interestingly, this step led to a statistically equivalent fit χ 2 /d.o.f.=400/310, which suggests the Fe Kα in ESO 511-G030 to be consistent either with a superposition of a narrow and constant core plus a broad and variable component or with a single and moderately broad Gaussian that varies in time and becomes stronger at a higher continuum flux.This result is illustrated in Fig. 6, whereby the line normalisation for a single broadened Gaussian significantly decreases between the high flux (2007) and low flux (2019) observations.This behaviour, which is quite at odds with what is commonly observed in AGNs, was also seen in NGC 2992 (Marinucci et al. 2020).

Spectral modelling
We investigated the broadband spectral properties of ESO 511-G030 by testing a purely phenomenological model.We A101, page 6 of 17 Fig. 5. Zoom of the Fe Kα complex.Top panels: EPIC-pn data in the 4-8 keV energy band fitted using a power-law and a narrow Gaussian emission line assumed to be constant among the exposures.Residuals in the Fe K energy range are present.Bottom panels: Same as above but with the addition of a broad Fe Kα emission line for which the line normalisation was assumed to be constant among the 2019 pointings but free to vary in 2007.modelled the XMM-Newton/NuSTAR data with a cut-off powerlaw absorbed for the Galaxy (N H = 4.33 × 10 20 cm −2 , HI4PI Collaboration 2016), a moderately broad Gaussian component for the Fe Kα, and a thermal component to account for curvature in the soft band (Model A).We fitted each observation separately.We report the inferred best-fit values and the statistic associated to each fit in Table 2.This procedure revealed that no significant spectral variations occurred during the campaign.The Fe Kα was constant in terms of normalisation, and in all observations except two and four, the line profile was consistent with being broad.The high-energy cut-off was constrained in observations one and five, while only lower limits were obtained in the remaining exposures.At lower energies, a weak and constant black-body-like component did not vary amongst the different observations.Given the little variability amongst the parameters, we fitted all the observations simultaneously, tying the photon index, the cut-off energy for the primary continuum, the energy centroid and width of the Gaussian component, and the temperature and normalisation of the black body.This resulted in a fit with χ 2 = 1470 for 1416 d.o.f and an associated null probability of 0.1.Moreover we found a Γ = 1.62 ± 0.02, while the high- Fig. 6.Confidence regions at 99%, 90%, and 68% for the line's width and normalisation.The contours were obtained assuming the model power-law+zGauss broad and using a single broad Gaussian component.
We subsequently tested a more reliable physical framework for the ESO 511-G030 2019 spectra.We started by considering two scenarios: (i) one accounting for a narrow Fe Kα, signature of distant reflecting material and (ii) another in which this emission feature is a blend of a relativistically broadened component plus a narrow one.The model Borus (e.g.Baloković et al. 2018) was used to account for the distant reflection and the relxill model (e.g.García et al. 2013García et al. , 2014) ) for the relativistic component.Within Borus, the toroidal X-ray reprocessor was assumed to have a spherical shape with conical cutouts at both poles, and the X-ray source was assumed to be at its centre.We used the table borus01_v161215a.ftz.Relxill (e.g.Dauser et al. 2016) is part of a model suite that accounts for ionised reflection from an accretion disc illuminated by a hot corona.In XSPEC notation, we thus tested the following models: tbabs G ×(cutoffpl + Borus), referred to as Model B, and tbabs G × (cutoffpl + relxill + Borus), referred to as Model C, for cases (i) and (ii), respectively.These models were applied to each XMM-Newton and NuSTAR dataset, and we fitted the Γ, the high-energy cut-off, and the normalisation of the primary continuum, tying these values with those of the Borus table.The column density and the normalisation of the Borus table were also computed in each exposure.We proceeded similarly when testing Model C, but in this case, we assumed the iron abundance to be solar (A Fe = 1) and computed the ionisation parameter ξ and the inner radius r in .Model C better reproduces the data, and Fig. 7 reports the corresponding best-fit values.
The ESO 511-G030 spectra are well described by a primary continuum with Γ = 1.73±0.02.Lower limits for the high-energy cut-off were inferred in all observations but observation five, for which E c = 75 ± 20 keV was obtained.Then, the narrow core of the Fe Kα emerges from a Compton-thin medium with an averaged column density N H ∼ 1.8 × 10 23 cm −2 , which is also responsible for the moderate high-energy curvature of the spectra.A relativistic reflection component likely originates from the mildly ionised matter, as the ionisation parameter is consistent amongst the exposures.However, the inner radius of this A101, page 7 of 17 Notes.The corresponding statistics are also shown.The cross-correlation constants used for XMM-Newton and the FPMA&B data are always within a 10%.The L39 quantity refers to the luminosity in units of 10 39 erg s −1 , and D 10 is the distance in units of 10 kpc, as defined in the XSPEC manual (see https://heasarc.gsfc.nasa.gov/xanadu/xspec/manual/node137.html).
component is poorly constrained in this model.The best-fit quantities inferred from the fits are shown in Table 2.
As a final test, we added a black body component in each observation to account for any weak underlying soft X-ray spectral features.The addition of this soft component did not provide any significant improvement to the fit.This test accords with the scenario of an absent or weak soft X-ray excess in the source.

Archival X-ray observations of ESO 511-G030
Thus far, ESO 511-G030 AGN has been observed by different facilities.Detailed studies on ASCA, Suzaku as well as the 2007 XMM-Newton exposure have already been published (e.g.Turner et al. 2001b;Ghosh & Laha 2021).Thus, we performed straightforward fits on these data to extract information on the spectral properties of ESO 511-G030 in previous years.
We tested our Model C by adopting the same fitting procedure described in Sect.4.2 but also slightly modifying the model when necessary.As clearly shown in Fig. 1, all archival exposures have larger fluxes than the data in 2019 and different spectral shapes, especially in the soft band.For this reason, we fitted separately the ASCA data, all the Suzaku observations, and XMM-Newton data.Concerning Suzaku, the additional crosscalibration constant (k) was used for the XIS and PIN data.
Interestingly, our baseline Model C failed to reproduce all the archival spectra, and thus additional components to this model were needed.In particular, a neutral absorption was required by ASCA data for which no relativistic reflection was necessary.For the Suzaku data, we added a single black body component to account for the soft excess observed below ∼1 keV, while two black body components were added to fit the 2007 XMM-Newton exposure (see also Ghosh & Laha 2021).These steps yielded the fits in Fig. 8, and in Table 3 we report the corresponding best-fit values.
Apart from the 1999 absorption event observed in the ASCA data, due to matter with constant N H = 1.2±0.3×10 21cm −2 , ESO 511-G030 archival spectra are consistent with a variable powerlaw that is, on average, softer than the power-law in the 2019 data (∆Γ ∼ 0.2).The continuum flux in the 2-10 keV energy range is up to a factor of six larger than in 2019, while the soft X-ray flux in the 0.5-2 keV range is even larger, by up to a factor of ten.The large changes in the soft flux can be ascribed to the presence(absence) of the soft excess component.The soft excess clearly plays a major role in shaping the soft band in the 2007 XMM-Newton and 2011 Suzaku observations.The relativistic reflection observed in the archival observations is about ten times larger than the one inferred in 2019, which is in agreement with the hot reflection responding to the primary changes.On the other hand, the reflected spectrum due to cold matter and that also contributes to the Fe Kα has a more constant behaviour across the years.Values derived for the normalisation of the Borus table are, in fact, rather constant across the epochs.Thus, a lower (higher) primary flux would correspond to a larger (smaller) reflection fraction.Finally, the Compton-thin nature of the reprocessor in ESO 511-G030 is further confirmed by this data.

The Swift Monitoring campaign
We present in this section the Swift X-ray and UV data that were taken in the context of two different observational campaigns, one of which is still ongoing, performed between 2020 and 2022.The objective of the campaigns is to keep track of the broadband properties of ESO 511-G030 and, possibly, observe a revenant soft excess.Up to fall 2021, however, the source flux A101, page 8 of 17 state had been consistent with the state observed in 2019 with XMM-Newton and NuSTAR.In Fig. 9, we show the light curves derived for different X-ray bands (0.3-0.5, 0.5-0.8,0.8-1.2,1.2-2, 2-5, and 5-10 keV) and optical-UV filters (which are possibly dominated by the host galaxy).The Swift-XRT data show a min-to-max variability of about a factor of five across the light curve.Once converted into fluxes, both X-rays and UVs were about a factor of approximately 10 fainter than in the XMM-Newton observation of 2007.Thus ESO 511-G030 as observed with Swift appears to be an extension of the quiescent state of ESO 511-G030 observed during 2019.By analysing the Swift spectra, our main objective was to establish whether a simple power-law component (plus Galactic absorption) was enough to explain the data.Our test confirmed that a simple power-law well represents the 0.3-10 keV energy range of our source and no additional components are required.In Table A.1, the inferred best-fit quantities are listed, and the corresponding best-fit values are shown as a function of the observing time in Fig. 10.The Swift data are consistent with a fairly flat spectral shape We show these spectra in black, magenta and blue, respectively.Bottom panels report on the data-model ratios.
of Γ = 1.62 ± 0.09, and only small variability was observed in the power-law normalisation.This may suggest the source still remains in a quiescent state with no soft excess awakening.
To test whether the soft excess was or was not present in XRT spectra, we used all the observations cited in Table A.1 to produce a stacked spectrum of ESO 511-G030.The obtained spectrum has no signature of this component.A simple powerlaw absorbed by the Galaxy (Γ = 1.62 ± 0.02) is in fact enough to account for the data (χ 2 = 337 for 339 d.o.f.), thus ruling out any additional component.The Swift fluxes in the 0.5-2 and 2-10 keV bands are compatible with what was observed during the 2019 XMM-Newton and NuSTAR monitoring campaigns, with F 0.5−2 keV = (2.18 ± 0.02) × 10 −12 erg cm −2 s −1 and F 2−10 keV = (5.1 ± 0.1) × 10 −12 erg cm −2 s −1 .Moreover, the stacked spectrum does not show any evidence of a Fe Kα emission line.The lack of this feature can be likely explained by the coupled effects of the source low flux and the small effective area of the Swift-XRT telescope in the hard X-rays.Finally, we note that the last A101, page 9 of 17 360/403 p.null 0.9 Flux 0.5−2 keV 6.0±0.1 ×10 −12 erg cm −2 s −1 Flux 2−10 keV 14.0±0.5 Notes.The L39 stands for the luminosity in units of 10 39 ergs s −1 , and D 10 is the distance in units of 10 kpc.
four Swift exposures of ESO 511-G030 were consistent with a flux increase of the source in both the X-ray and UV energy bands.

Relation between X-rays and UVs
A viable way to quantify the actual relation between X-rays and UVs band is provided by the α OX parameter.The non-linear relation between the UV and X-ray luminosity in AGNs was discovered in late 1970s by Tananbaum et al. (1979), and several other authors have investigated the physical meaning of such a relation (e.g.Tananbaum et al. 1986;Zamorani et al. 1981;Vagnetti et al. 2010;Martocchia et al. 2017;Chiaraluce et al. 2018) respectively), providing rich information on the UV continuum slope.We therefore derived the X-ray luminosity at 2 keV in order to compute the α OX values.We relied on our best-fit to the EPIC-pn spectra discussed in Sects.4 and 5, while the XRT spectra were modelled with a simple power-law with Galactic absorption (see Appendix A).
A101, page 10 of 17 (MJD < 59580), after which a moderate increase of ultraviolet emission can be observed, the optical-UV curves could be highly affected by the host emission.We noticed that to quantify the actual amplitude of the optical-UV emission, an estimate of the contribution of the host galaxy is required (see Sect. 6).
The total AGN plus Galaxy luminosity at log ν * = 15.08 corresponding to λ = 2500 Å was estimated at each epoch as an interpolation between the data in the closest filters, which for these observations are UVW1 and UVM2.We computed the optical-UV spectral energy distributions (SEDs) to get an estimate of the host galaxy fraction at ν * , which may be significant, and of the AGN luminosity, L UV , following the prescriptions by Vagnetti et al. (2013).We assumed each optical-UV SED to be the sum of an AGN spectrum proportional to the average SED by Richards et al. (2006) with a typical slope of α AGN = −0.57(with F ν ∼ ν α AGN ) and a host galaxy contribution whose spectrum was modelled with a typical slope α galaxy = −3 (see e.g., Lusso et al. 2010;Vagnetti et al. 2013).The host galaxy fraction f g at ν * was then derived at each epoch as a function of the sole spectral index.Indeed, expressing the total luminosity as The fraction f g was thus estimated by inverting the previous relation as f g = −(0.57+ α)/2.43.Therefore, when f g tends to zero, then the slope α tends to −0.57, as per the Richards et al. (2006) AGN SED.However, when f g = 1, α tends to −3, which is consistent with a pure host galaxy spectrum.
The monochromatic fluxes for the optical and UV filters were computed by converting the observed rates using the appropriate conversion factors.The spectral index was derived by applying least squares on these data.Figure 11 shows the optical-UV SEDs for the XMM-Newton and Swift observations.From 2007 to 2022 pointings, the ESO511-G030 SEDs underwent dramatic spectral changes.The optical-UV ESO 511-G030 SEDs appear quite steep, with the exception of the 2007 exposure and those from Swift in 2022.The corresponding spectral index lies in the range between −2.7 and −3.These values are far steeper than the more typical slope of −0.57derived from the average spectral energy distribution of a statistically significant number of AGNs by Richards et al. (2006) around ν * .The variation in the SEDs suggests that the nuclear emission changes, while a substantial constant contribution from the host galaxy is also present.In particular, the steep slopes derived for the observations taken after 2019 but before 2022 can be ascribed to the dominant shape of the host galaxy.
The host galaxy luminosity at 2500 Å was defined as the average value of its estimates at different epochs, which was found to be log L galaxy = 27.68 (erg s −1 Hz −1 ) with a small dispersion, σ = 0.05.This average value was subtracted from the total luminosity at 2500 Å to obtain the AGN luminosity, L UV , at each epoch.This luminosity was often smaller than the host galaxy luminosity, as shown in the top panel of Fig. 12.For most observations, the ratio between the AGN and the total monochromatic luminosity (AGN+host) is in the range of 0-50%.In many cases, the SED slope is − 3, and the AGN fraction is negligible and thus not plotted in Figs. 12 and 13.Regarding the SED of the 2007 archival observation, it is flatter than the one reported by Richards et al. (2006), with α = −0.44.Thus, for this 2007 observation, we assumed the total monochromatic luminosity at 2500 Å to be due to the AGN only.
The α OX derived for our source can be also compared with the well-known L UV − α OX anti-correlation (Vignali et al. 2003;Just et al. 2007;Vagnetti et al. 2010).Figure 13 shows the track of ESO 511-G030 in the log L UV − α OX plane.

Ultraviolet to X-ray modelling
In accordance with Sect.6, the multi-wavelength properties of ESO 511-G030 varied dramatically with time.In 2007, in fact, the source UV-to-X-rays were consistent with that of a bare Type 1 AGN.Later, once re-observed in 2019, neither the UVs nor the X-rays where compatible with their historical fluxes.Moreover, the UV emission, which in 2007 was fully ascribable to the accretion process, turned out to be Galaxy dominated, with only a few percent of the flux being due to the AGN.
To better understand the interplay amongst the different emission components in ESO 511-G030 across the years, we modelled the XMM-Newton spectrum taken in 2007 and those obtained in 2019.At this stage, we included the corresponding OM data and tested AGNSED (Kubota & Done 2018).This model allowed us to self-consistently reproduce the UV-to-Xray spectra of ESO 511-G030.In accordance with Done et al. ( 2012), AGNSED accounts for three distinct emitting regions: an outer standard disc region, a warm Comptonising corona, and an inner hot Comptonising plasma.The flow is radially stratified and emits as a standard black-body-like disc from R out to R warm , as warm Comptonisation from R warm to R hot (adopting the passive disc scenario by Petrucci et al. 2018), and as the typical hot Comptonisation component below R hot down to R ISCO .We thus tested the model5 : tbabs G × redden × (galaxy + AGNSED + relxillcp + Borus).(1) The tbabs G was fixed to the value N H = 4.33 × 10 20 cm −2 (HI4PI Collaboration 2016), while the 'redden' component was fixed to a value of E(B − V) = 0.056, in agreement with Schlafly & Finkbeiner (2011).We tested this model by fitting each observation separately, and within AGNSED, we allowed the photon indices of the warm and hot coronae to vary as well as the R hot and R warm parameters.The hot coronal temperature was fixed to 100 keV to fit the 2007 data because there is no constraint above 10 keV, while it was left free to vary in the 2019 data.The warm coronal temperature was computed in both datasets.We used a co-moving distance of 96 Mpc (derived from the redshift) 6 and adopted the mass for the SMBH in ESO 511-G030 by Ponti et al. (2012).The models Relxillcp and Borus (borus11_v190815a.fits), both accounting for a Comptonised continuum, were set similarly to what we described in Sect.4.2.These different flavours of relxill and Borus were used in order to tie the model parameters accounting for a Comptonised continuum (e.g.kT e instead of E cut ) with the corresponding parameters in AGNSED, which assumes an underlying Comptonised continuum and not simply a cut-off power-law.In the fit procedure, we assumed the spin of the central SMBH to be maximally rotating in both relxillcp and AGNSED, and the same disc inclination of 30 • was also assumed.Within AGNSED, the upper limit of the scale height for the hot Comptonisation component was set to be 10 R g , which mimics a spherical Comptonisation region of similar radius.In relxillcp, we set the emissivity profile to its default value of three.The inner radius for the relxillcp reflection was fixed to R hot .This assumption is discussed later.Then the outer disc radius was set to be the same (R out = 400R g ) between the model component AGNSED and relxillcp.After preliminary tests, we fixed the relxillcp parameter R in to 100 R g , as it could not be constrained in the 2019 data.Finally, we added a galaxy template accounting for the host galaxy (matching the morphological type of ESO 511-G030; i.e.Sc; Lauberts 1982) contributing to the UV flux.The table was computed following Ezhikode et al. (2017) and included within XSPEC as a template named 'hostpol' (Polletta et al. 2007).We note that our galaxy model has a spectral shape between 3-10 eV of ∼ Γ = 4.This shape is fully consistent with our assumption in Sect.7 for the 6 https://www.astro.ucla.edu/~wright/CosmoCalc.html.

XMM-Newton 2007
Swift XMM-Newton 2019 Richards et al. 2006 Fig. 11.ESO 511-G030 optical-UV SEDs derived using XMM-Newton and Swift data (dark green).The corresponding least squares fits are shown as straight dashed lines.XMM-Newton and Swift data are in blue and dark green, respectively.The solid, thick black line is the average SED by Richards et al. (2006) for Type 1 objects in the SDSS.The vertical dashed line is at log ν * = 15.08,corresponding to 2500 Å.We noticed a remarkable change in spectral shape from the 2007 observation compared to the later ones.
host galaxy.We tied the normalisation of hostpol across the five spectra.
A101, page 12 of 17 Fig.12.Light curves for the AGN luminosity at 2500 Å, at 2 keV, and the inferred α OX .Blue dots identify values derived using XMM-Newton observations while in green the same quantities were estimated using Swift.The dashed horizontal line accounts for the adopted host galaxy luminosity at 2500 Å .Finally, both L UV and L 2 keV are in units of erg s −1 Hz −1 .
This procedure led us to the best-fit and model for data ratios, shown in Fig. 14.The fit information as well as the corresponding yielded quantities are shown in Table 4.
According to these fits, the SED of ESO 511-G030 varied dramatically from 2007 to 2019.As for the Eddington ratio, it varied from a value of L/L Edd ∼ 2% (in the 2007 spectrum soft excess+power-law) to a rather low and radiatively inefficient value of ∼0.2% in 2019.The standard configuration assumed within AGNSED, that is, the presence of a hot plasma (for R between R isco and R hot ∼ 27R g ), a warm plasma (from R hot to R warm ∼ 150R g ), and an accretion disc radially segregated, agrees with the 2007 XMM-Newton observation.As stated before, the reflection in our model is produced beyond R hot .However, in AGNSED a warm corona is present between R hot and R warm , and the standard disc starts only beyond R warm .Though we believe that the impact on our best-fit results should be limited, we note that our modelling does not take into account the presence of the warm corona in the reflection computation and thus does not provide a fully self-consistent physical picture of the emission emerging from ESO511-G030.To our knowledge, the model REXCOR (Xiang et al. 2022) is the sole publicly available model that self-consistently computes the reflection spectrum from the combination of an outer standard disc and inner warm corona.However, this model cannot be extended to the UV energy range and cannot be used to simultaneously fit the OM data.
The lack of a substantial soft excess during the 2019 campaign leads towards a different layout where the hot corona is now more extended than before and no clear indication of a warm Comptonising region is found.This warm region has shrunk, and the disc already extends from 50R g , close to the hot component.The relativistic reflection component is less prominent, and no constraints on the inner radius were obtained.The cold reflection component is compatible amongst the 2007 and 2019 data, in agreement with its distant origin from the central engine.Finally, the change of the hot corona is also accompanied by the spectra evolving from a softer to a harder state, from Γ = 1.91 ± 0.03 to an average value of Γ = 1.75 ± 0.02.

Discussion and conclusions
We reported on the spectral and temporal properties of ESO 511-G030 that showed significant variability in both the optical-UV and X-ray bands.Our analysis revealed the ESO 511-G030 spectrum to be consistent with a primary power-law Γ = 1.73 ± 0.02 and accompanied by a poorly constrained high-energy rollover.The reflected flux we observed in the X-rays of ESO 511-G030 emerges from regions of different densities (Compton-thin and Compton-thick).In Fig. 7, we find the high-energy spectrum of the source to be dominated by reflection off a Compton-thin medium N H ∼ 1.8 × 10 23 cm −2 that does not produce a relevant Compton-hump but accounts for a narrow Fe Kα line.The relativistic reflection component reproduces the moderate broad shape of the same emission line and contributes to the overall spectral curvature.
Testing our Model C on archival data (see Sect. 5 and the work by Ghosh & Laha 2021) revealed the primary component to be harder in 2019 than in past observations.We found these harder states to correspond to lower flux levels, suggesting the commonly observed softer-when-brighter trend.The reflected flux also varied.In particular, the relativistic reflection was found to follow the variations of the primary emission (in agreement with this flux being released in the close surroundings of the central engine), while a less variable behaviour was observed for the cold reflection.
One of the main features of the 2019 observational campaign is a lack of substantial soft excess.A simple power-law, in fact, dominates the soft-to-hard X-ray spectrum of ESO 511-G030, at least from 2019.The lack of a soft excess, or its negligible contribution to the overall emission spectrum, is further supported by the analysis of the 2007 and 2019 excess variance spectra.Two different components are in fact needed to account for the 2007 F var spectrum: one responsible for the changes in the X-ray A101, page 13 of 17 continuum and a second accounting for the soft excess.In contrast, the F var computed for the 2019 XMM-Newton exposures only requires a single component accounting for variance due to the nuclear continuum.It is worth noting that while the F var spectra of the Fig. 4 sample has slightly different timescales (see Sect. 3), we verified that when cutting the 2007 observation into shorter segments (so as to sample similar timescales as for the 2019 F var spectra), a soft excess component still appears in the F var of the lowest flux segment.
The absence of a strong soft excess is quite unusual since it is ubiquitously observed in AGNs (e.g.Piconcelli et al. 2005;Bianchi et al. 2009;Gliozzi & Williams 2020).We note that ESO 511-G030 data do not require any absorbing component, as also discussed in Laha et al. (2014), thus the lack of the soft excess in this source cannot be explained with an absorption process.
The case of ESO 511-G030 is peculiar, and we can only compare its behaviour with that of Mrk 1018.This AGN had been studied in depth by Noda & Done (2018), who observed very different spectral shapes corresponding to different Eddington ratios.From a typical Type 1 spectrum with a strong soft excess, the source dimmed down, became harder, and showed a weaker soft excess.This spectral transition corresponded to a change in the Eddington ratio from L/L Edd ∼2% to L/L Edd ∼0.4%, which is very similar to what is observed here for ESO 511-G030.As the soft excess is responsible for most of the ionising photons, the dramatic drop in the X-rays also led to the disappearance of the Broad Line Region, producing the 'changing-look' phenomenon.In other words, the presence (lack) of the soft excess corresponded to a general softening (hardening) of the X-ray continuum emission with an accompanying dramatic change in the disc emission and a disappearance of the optical broad lines.
Similar to Mrk 1018, ESO 511-G030 had a dramatic change in its accretion rate, passing from L/L Edd ∼ 2% in 2007 down to L/L Edd ∼ 0.2% in 2019.This dimming was also accompanied by a dramatic change in the UV SEDs (see Figs. 11 and 14), though it did not lead to a 'complete' changing-look process.An optical FLOYDS spectrum was in fact taken quasi-simultaneously with the XMM-Newton-NuSTAR campaign to check whether broad lines were present or not.From a quick comparison between the FLOYDS spectrum and an 6dF archival one taken in 2000, the Hβ line does not disappear in 2019, and its velocity width was similar between the spectra (FWHM Hβ ∼ 4500 km s −1 , private communications with Keith Horne and Juan V. Hernández Santisteban).
The strong decrease of the accretion rate between 2007 and 2019 seems to be the crucial element for explaining the observed spectral UV-X-ray behaviour.Indeed this decrease naturally explains the strong decrease of the UV emission (see Figs. 12 and 13).To the 0th order, the decrease of the UV flux would also mean a reduction of the soft photons flux entering and cooling the hot corona.So we would expect an increase of the hot plasma temperature and a hardening of the X-ray spectrum with respect to the observation from 2007.Spectral hardening was observed, and the absence of a stringent high-energy cut-off signature in 2019 agrees with a high corona temperature, much higher than the usual values observed in Seyfert galaxies (e.g.Fabian et al. 2015Fabian et al. , 2017;;Tamborra et al. 2018;Middei et al. 2019).The absence of high signal-high energy observations in the archives prevents any comparison with past observations that would help to support this scenario.But the strong decrease of the accretion rate could also explain the absence of the soft X-ray excess in 2019, at least in the case of the warm corona model.The observations agree with the warm corona being the upper layers of the accretion disc (Petrucci et al. 2018, and references therein).More importantly, to reproduce the soft X-ray spectral shape, simulations show that a large enough accretion power has to be released inside this warm corona and not in the accretion disc underneath (Różańska et al. 2015;Petrucci et al. 2020;Ballantyne 2020).So if the accretion power becomes too low, the warm corona cannot be energetically sustained.It is less obvious to understand why the soft X-ray excess would disappear if it is due to relativistically blurred ionised reflection.It is possible, however, that at a low accretion rate, the disc becomes more optically thin (or even recedes) and produces less reflection, as is indeed observed in 2019.However, other possible explanations for the lack of the soft excess in ESO 511-G030 may be viable, and new exposures, possibly performed during the awakening of this component, are needed in order to shed light onto the engine of this Seyfert galaxy.The increasing UV and X-ray fluxes observed by Swift in the first quarter of 2022 encourage us A101, page 14 of 17 Notes.The symbol is used for those values that were kept fixed in the fitting procedure.
to ask for more observing time to, possibly, observe the revenant soft excess of ESO 511-G030.

Fig. 1 .
Fig.1.Unfolded spectra of ESO 511-G030 corresponding to the observations listed in Table1.Variability is observed both in terms of spectral and amplitude changes.In the optical-UV band, the source faded by a factor of approximately ten.The underlying model folding the data is a power-law with Γ = 2 and a unitary normalisation.

Fig. 2 .
Fig.2.ESO511-G030 light curves.Top panels: Background-subtracted light curves from the 2019 multi-wavelength campaign.The first and second row refer to the XMM-Newton 0.5-2 and 2-10 keV bands, respectively, while the last row shows the ratios between the two bands.The NuSTAR light curves extracted in the 10-79 keV band are also shown.Bottom panels: Background-subtracted time series for ASCA (black), the 2007 XMM-Newton orbit (red) and the Suzaku visits two and three (in blue).

Fig. 4 .
Fig. 4. Comparison between the F var spectrum of the 2007 and 2019 data.The 2019 observations (magenta) were obtained by grouping together the five exposures via the standard command setplotgroup within XSPEC.The spectrum derived from 2007 data is in black.The magenta curve represents a theoretical F var spectrum, which depicts the combined contribution from a variable power-law-like continuum and a variable soft excess.

Fig. 7 .
Fig. 7. Model C fitting to the 2019 spectra.Borus and relxill account for the cold and hot reflection components, respectively.

Fig. 8 .
Fig.8.Fit to the ASCA, 2007 XMM-Newton, and Suzaku data.We show these spectra in black, magenta and blue, respectively.Bottom panels report on the data-model ratios.

Fig. 9 .
Fig.9.Multi-filter light curves derived from 2018 up to 2022 for ESO 511-G030 using Swift-XRT and -UVOT.Until the end of 2021 (MJD < 59580), after which a moderate increase of ultraviolet emission can be observed, the optical-UV curves could be highly affected by the host emission.We noticed that to quantify the actual amplitude of the optical-UV emission, an estimate of the contribution of the host galaxy is required (see Sect. 6).

Fig. 10 .
Fig. 10.Best-fit parameters and fluxes in the 0.3-2 (Soft) and 2-10 keV (hard) bands derived from Swift.Straight solid lines show the mean value of the parameter, while dotted lines represent the corresponding standard deviation.Fluxes are in units of 10 −12 erg cm −2 s −1 .

Fig. 13 .
Fig. 13.ESO 511-G030 α OX estimates compared with the sample by Chiaraluce et al. (2018).The solid black line represents logL UV − α ox linear relation as derived in the same work.The α OX computed of ESO 511-G030 shows significant variability in agreement with the drop in the OM data shown in Fig. 1.

Table 1 .
Log of the observations.

Table 2 .
Best-fit parameters obtained using models A and B (see Sect. 4.2 for details).

Table 3 .
Parameters derived by adapting the baseline Model C to the archival data.

Table 4 .
Best-fit parameters derived adopting AGNSED on the XMM-Newton 2007 exposure and data from the 2019 monitoring campaign.