HTTP_Request2_Exception Unable to connect to tcp://think-ws.ca.edps.org:85. Error: php_network_getaddresses: getaddrinfo failed: Name or service not known XMM-Newton – NuSTAR monitoring campaign of the Seyfert 1 galaxy IC 4329A | Astronomy & Astrophysics (A&A)
Open Access
Issue
A&A
Volume 687, July 2024
Article Number A51
Number of page(s) 15
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/202348760
Published online 26 June 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

Active galactic nuclei (AGN) are extremely luminous compact objects located at the center of massive galaxies. They are powered by the accretion of gas onto the central supermassive black hole (SMBH; MBH > 105M; Salpeter 1964; Ho 2008). The X-ray emission of AGN is due to high-energy processes that occur in the so-called hot corona, which contains high-energy electrons and is located near the SMBH (Haardt & Maraschi 1991; Haardt et al. 1994, 1997; Merloni 2003; Fabian et al. 2009; Zoghbi et al. 2012; De Marco et al. 2013). The hot electrons interact with thermal UV/optical photons emitted by the accretion disk through inverse-Compton scattering, which boosts the energy of the photons into the X-ray band. This process produces a broad power-law continuum in the X-ray spectrum with a cutoff at high energy (e.g., Sunyaev & Titarchuk 1980; Haardt & Maraschi 1993). Cold circumnuclear materials, such as the optically thick accretion disk and/or the molecular torus, absorb and reprocess the X-rays, which results in a feature in the X-ray spectrum that is known as Compton reflection, which peaks at around ∼30 keV (Pounds et al. 1990; George & Fabian 1991). Additionally, the reflection process produces absorption and emission lines, such as the prominent FeKα emission line at 6.4 keV (Nandra & Pounds 1994), which are the result of photoelectric absorption and fluorescence (Matt et al. 1997).

IC 4329A is a bright nearby (z = 0.01598, Koss et al. 2022a) AGN and it is classified in the optical as a Seyfert 1.5 galaxy (Oh et al. 2022). It is located at the center of an edge-on host galaxy, with a dust lane passing through the nucleus. It also shows a dusty lukewarm absorber in the UV band (i.e., T ∼ 3 × 104 K, Crenshaw & Kraemer 2001), characterized by saturated UV absorption lines (C IV, N V) near the systemic velocity of the host galaxy, which likely causes the reddening of the continuum and emission lines. The black hole mass of IC 4329A is log ( M BH ) [ M ] = 7 . 15 0.26 + 0.38 $ \log(M_{\mathrm{BH}})[M_{\odot}]=7.15^{+0.38}_{-0.26} $ (GRAVITY Collaboration 2024), and its Eddington ratio is λEdd = Lbol/LEdd = 0.17. This was computed using the bolometric luminosities calculated from the intrinsic luminosities in the 14–150 keV range, as shown in Ricci et al. (2017), with a bolometric correction of 8 (Koss et al. 2022b). Assuming Eq. (4) from Du et al. (2018), the dimensionless accretion rate of IC 4329A is M ˙ = 1.4 $ \dot{\mathscr{M}}=1.4 $ (GRAVITY Collaboration 2024).

IC 4329A has been extensively studied in the X-ray band in the past decades by all major X-rays satellites, including Einstein (Petre et al. 1984; Holt et al. 1989), EXOSAT (Singh et al. 1991), Ginga (Piro et al. 1990), ASCA/RXTE (Done et al. 2000), BeppoSAX (Perola et al. 2002; Dadina 2007), Chandra (McKernan & Yaqoob 2004; Shu et al. 2010), XMM-Newton (Steenbrugge et al. 2005; Nandra et al. 2007; Mehdipour & Costantini 2018), INTEGRAL (Beckmann et al. 2006), Swift/BAT (Winter et al. 2009; Ricci et al. 2017), Suzaku (Mantovani et al. 2014) NuSTAR (Brenneman et al. 2014a,b), and AstroSat (Tripathi et al. 2021; Dewangan et al. 2021). Its 2–10 keV flux ranges between 8.1 × 10−11 erg cm−2 s−1 (Winter et al. 2009) and 1.8 × 10−9 erg cm−2 s−1 (Shu et al. 2010). Previous BeppoSAX, ASCA+RXTE INTEGRAL and Swift/BAT observations of IC 4329A placed rough constraints on the high-energy cutoff at Ec ≥ 180 keV (Perola et al. 2002), Ec = 150 − 390 keV (Done et al. 2000), Ec = 60 − 300 keV (Molina et al. 2013) and Ec > 190 keV (Ricci et al. 2017), respectively. The NuSTAR hard X-ray spectrum of this AGN is characterized by a power law with a photon index of Γ = 1.73 ± 0.01 and a high-energy cutoff at Ec = 186 ± 14 keV (Brenneman et al. 2014b). Previous analysis also showed a neutral reflection component with a reflection fraction of Rrefl = 0.42 ± 0.02 (Mantovani et al. 2014).

This source is well known for its strong Fe Kα line (Piro et al. 1990). The Fe Kα line in IC 4329A shows a narrow core at 6.4 keV that is consistent with being produced by low-ionization material in the outer accretion disk, the broadline region (BLR), or in the torus. The narrow core component of the Fe Kα is also variable on timescales of weeks to months (Fukazawa et al. 2016; Andonie et al. 2022). From the analysis of the ASCA, RXTE, BeppoSAX, Suzaku, and NuSTAR observations (Done et al. 2000; Dadina 2007; Brenneman et al. 2014b; Mantovani et al. 2014), a broad component was reported that is most likely produced in the inner part of the accretion disk and is blurred by general relativistic effects. In other works (e.g., McKernan & Yaqoob 2004; Markowitz et al. 2006; Nandra et al. 2007; Tripathi et al. 2021), the observations suggested a modest or weak broad iron line, indicating that X-ray reflection from the inner disk is weak in this source. Analyzing the Suzaku and NuSTAR observations of IC 4329A with a relativistic reflection model, Ogawa et al. (2019) inferred a very low reflection fraction (R ∼ 3.2 × 10−3).

IC 4329A also shows an ionized warm absorber (WA) component with log(ξ/erg s−1 cm) in the range of −1.37 and 2.7 (Steenbrugge et al. 2005) and a highly ionized, that is, log(ξ/erg s−1 cm) = 5.34 ± 0.94, ultrafast outflowing component (UFO; Tombesi et al. 2015) with an outflow velocity v/c = 0.098 ± 0.004. This component was first detected by Markowitz et al. (2006) and was then confirmed by Tombesi et al. (2011).

Recently, IC 4329A was observed by the Imaging X-ray Polarimetry Explorer (IXPE) for ∼500 ks. From these observations, it appears that the source shows a 1σ confidence limit on a polarization degree of 3.3%±1.1% and a polarization angle of 78° ±10°, which is consistent with it being aligned with the radio jet (Ingram et al. 2023; Pal et al. 2023). These finding favor coronal geometries that are more asymmetric and might be outflowing. The coronal geometry is unconstrained within the 3σ confidence level, however.

We present the results of the analysis of the X-ray broadband spectra from the simultaneous XMM-Newton and NuSTAR observations of IC 4329A. These observations are part of a joint millimeter/X-ray study of this nearby source, which also includes observations from Swift (10 observations), NICER (20 observations), and from the Atacama Large Millimeter/submillimeter Array (ALMA; 10 observations). The aim of the campaign is to study the relation between the 100 GHz and X-ray continuum in AGN and to test the idea that the millimeter continuum is associated with self-absorbed synchrotron emission from the X-ray corona (e.g., Laor & Behar 2008; Inoue & Doi 2014), as suggested by the recent discovery of a tight correlation between the X-ray luminosity and the 100 GHz (Ricci et al. 2023) and 200 GHz (Kawamuro et al. 2022) luminosities of nearby AGN. We focus on the high-resolution X-ray spectroscopy of IC 4329A and therefore do not include Swift and NICER data, which will be presented together with the results of ALMA campaign in a dedicated forthcoming paper (Shablovinskaya et al., in prep.).

The paper is organized as follows: In Sect. 2, we present the dataset we analyzed in this work. In Sects. 3 and 4, we describe the timing and spectral data analysis processes, respectively. In Sect. 5, we discuss the results of our analysis, which are summarized in Sect. 6. Standard cosmological parameters (H = 70 km s−1 Mpc−1, ΩΛ = 0.73 and Ωm = 0.27) are adopted.

2. Data and data reduction

The dataset analyzed in this work consists of nine X-ray Multi-Mirror Mission (XMM-Newton; Jansen et al. 2001) observations, five of which were performed simultaneously with the Nuclear Spectroscopic Telescope Array (NuSTAR; Harrison et al. 2013). Details of the duration and exposure of the observations are reported in Table 1.

Table 1.

Summary of the 2021 XMM-Newton and NuSTAR observations of IC 4329A.

2.1. XMM-Newton data reduction

IC 4329A was observed once per day for ten consecutive days (from 10 August 2021 to 19 August 2021) by XMM-Newton (P.I. C. Ricci) during the XMM-Newton AO 19. The XMM-Newton observations were performed with the European Photon Imaging Camera (EPIC hereafter) detectors, and with the Reflection Grating Spectrometer (RGS; den Herder et al. 2001). The EPIC cameras were operated in small-window and thin-filter mode. Observation 0862090401 is not included in the analysis because due to a problem in the ground segment, the EPIC-pn exposure was lost.

The EPIC camera (Strüder et al. 2001) event lists were extracted with the EPPROC and EMPROC tools of the standard system analysis software (SAS v.18.0.0; Gabriel et al. 2004). The MOS detectors (Turner et al. 2001) and the RGS were not considered because due to the lower statistics of their spectra, they do not add information to the analysis.

The choice of the optimal time cuts for the flaring particle background and of the source and background extraction radii was performed via an iterative process that maximized the signal-to-noise ratio (S/N) as in Piconcelli et al. (2004). The resulting optimum extraction radius was 30″, and the background spectra were extracted from source-free circular regions with radii of ∼50″ for EPIC-pn for each observation. Response matrices and auxiliary response files were generated using the SAS tools RMFGEN and ARFGEN, respectively. EPIC-pn spectra were binned in order to oversample the instrumental energy resolution by a factor larger than three and to have no fewer than 20 counts in each background-subtracted spectral channel. No significant pile-up affected the EPIC data, as indicated by the SAS task EPATPLOT. EPIC-pn light curves were extracted by using the same circular regions for the source and the background as the spectra.

2.2. NuSTAR data reduction

The five NuSTAR observations were performed simultaneously with five XMM-Newton observations during NuSTAR AO 7 (P.I. C. Ricci). The NuSTAR telescope observed the source with its two coaligned X-ray telescopes Focal Plane Modules A and B (FPMA and FPMB, respectively).

The NuSTAR high-level products were obtained using the NuSTAR data analysis software (NuSTARDAS) package (v2.1.1) within the HEASOFT package (version 6.29). Cleaned event files (level 2 data products) were produced and calibrated using standard filtering criteria with the NUPIPELINE task and the latest calibration files available in the NuSTAR calibration database (CALDB 20220802). For both FPMA and FPMB, the radii of the circular region used to extract source and background spectra were 75″; no other bright X-ray source is present within 75″ from IC 4329A, and no source was present in the background region. The low-energy (0.2–5 keV) effective area issue for FPMA (Madsen et al. 2020) does not affect our observations because no low-energy excess is found in the spectrum of this detector. The spectra were binned in order to oversample the instrumental resolution by a factor of 2.5 at least and to have an S/N higher than 3 in each spectral channel. Light curves were extracted using the NUPRODUCTS task, adopting the same circular regions as for the spectra.

3. Timing analysis

During the monitoring, IC 4329A varied within some observations (see Figs. 1 and 2). However, these variations did not affect the ratio of the flux in the 2.5–10 keV and the 0.2–2.5 keV band, (< 15%; see the bottom panels of Figs. 1 and 2).

thumbnail Fig. 1.

XMM-Newton EPIC-pn light curves (background subtracted) in the 0.2–10 keV, 0.2–2.5 keV and in the 2.5–10 keV are shown in the top, top middle, and bottom middle panels, respectively. In the bottom panels, the ratios of the XMM-Newton EPIC-pn background-subtracted light curves in the 2.5–10 keV and 0.2–2.5 keV energy bands are reported. All the light curves are extracted using a binning time of 500 s. The solid and dashed black lines indicate the average and the standard error of the mean, respectively. The color code used in this figure for different ObsID is applied throughout the paper.

thumbnail Fig. 2.

NuSTAR FPMA+B background-subtracted light curve in the 3–70 keV. 3–10 keV and 10–70 keV energy ranges are shown in the top, top middle, and bottom middle panels, respectively. In the bottom panels, the ratios of the NuSTAR NuSTAR background-subtracted light curves in the 10–70 keV and 3–10 keV energy bands are reported. All the light curves are extracted using a binning time of 500 s. The color code used in this figure for different ObsID is applied throughout the paper. The solid and dashed black lines indicate the average and the standard error of the mean, respectively.

We studied the variability spectra of the XMM-Newton EPIC-pn observations of IC 4329A using one of the common estimators for X-ray variability: the fractional variability Fvar (Vaughan et al. 2004; Papadakis 2004; McHardy et al. 2006; Ponti et al. 2012; Matzeu et al. 2016, 2017; De Marco et al. 2020). The Fvar, that is, the square root of the normalized excess variance ( σ NXS 2 $ \sigma^2_{\rm NXS} $ Nandra et al. 1997a; Edelson et al. 2002; Vaughan et al. 2003), is the difference between the total variance of the light curve and the mean squared error that is normalized for the average of the number of flux measurements squared. N is the number of good time intervals in a light curve, and xi and σi are the flux and the error in each interval. Fvar is defined (Vaughan et al. 2003) as

F var = σ NXS 2 = S 2 σ 2 ¯ x i ¯ 2 , $$ \begin{aligned} F_{\rm var}=\sqrt{\sigma ^2_{\rm NXS}}=\sqrt{\frac{S^2-\overline{\sigma ^2}}{\overline{x_i}^2}}, \end{aligned} $$(1)

where S 2 = 1 N 1 i = 1 N [ ( x i x i ¯ ) 2 ] $ S^2=\frac{1}{N-1}\sum_{i=1}^{N}[(x_i-\overline{x_i})^2] $ is the sample variance, that is, the integral of the power spectral density (PDS) between two frequencies, and σ 2 ¯ = 1 N i = 1 N [ σ i 2 ] $ \overline{\sigma^2}=\frac{1}{N}\sum_{i=1}^{N}[\sigma_i^2] $ is the mean square error.

We computed Fvar for each observation using the background-subtracted EPIC-pn light curves in different energy bands using a temporal bin of 1000 s on a timescale of 20 ks. The errors were computed using Eq. (B2) of Vaughan et al. (2003). The resulting Fvar spectra of each observation are shown in Fig. 3. To study the behavior of the variability spectra, we used the linmix code, a hierarchical Bayesian model for fitting a straight line to data with errors in both the x and y directions (Kelly 2007). To perform the fitting, we used a linear model of the data with the following fitting relation:

F var = α × E + β . $$ \begin{aligned} F_{\rm var}=\alpha \times E+\beta . \end{aligned} $$(2)

thumbnail Fig. 3.

Fractional variability spectra for the XMM-Newton EPIC-pn observations of IC 4329A as in Table 1. The spectra of this campaign are characterized by a mostly flat shape, with a fractional variability that never exceeds 5%. The dot-dashed black lines are the linear regressions, and the shaded regions represent the combined 3σ error on the slope and normalization.

We found the behavior of the variability spectra of this campaign to be nearly flat, with a correlation coefficient ρPearson spanning between −0.05 and 0.05.

The different spectral components observed in AGN can lead to spectral variability in different energy bands as each of them can dominate in a certain energy band. For example, the primary power-law component or the reflection component are dominant in the hard energy band (3–10 keV, Haardt & Maraschi 1991, 1993; Haardt & Matt 1993), while soft-excess and WA can dominate the soft (0.2–1 keV, Bianchi et al. 2009) and medium (1–3 keV, Blustin et al. 2005; Tombesi et al. 2013) energy bands. We therefore calculated the σ NXS 2 $ \sigma^2_{\rm NXS} $ from the 0.2–1 keV (soft), 1–3 keV (medium), and 3–10 keV (hard) light curves to obtain a complete picture of the X-ray variability of IC 4329A, and we compare these values in Fig. 4. The σ NXS 2 $ \sigma^2_{\rm NXS} $ in the soft, medium, and hard energy band for the XMM-Newton observations of IC 4329A in this campaign are very low in general. They are below 0.1% in all the observations, except for the σ NXS 2 $ \sigma^2_{\rm NXS} $ in the soft energy band of Obs.ID 0862090101 (0.15%). The σ NXS,1-3keV 2 $ \sigma^2_{\rm NXS, 1{-}3\,keV} $ vs. σ NXS,3-10keV 2 $ \sigma^2_{\rm NXS, 3{-}10\,keV} $ values are located in the vicinity or slightly below the one-to-one relation. The same holds for the σ NXS,0.2-1keV 2 $ \sigma^2_{\rm NXS, 0.2{-}1\,keV} $ vs. σ NXS,3-10keV 2 $ \sigma^2_{\rm NXS, 3{-}10\,keV} $ values, except again for Obs.ID 0862090101, which is located above the one-to-one relation. This could also be related to the small peak of the variability spectrum of this observation at around 0.5 keV. Thus, IC 4329A shows very weak variations overall, and the soft-excess and/or WA component variations are even weaker than those of the primary continuum and/or reflection component on a timescale of 20 ks (except for Obs.ID 0862090101, in which most likely the soft-excess or the absorbing components are more variable on timescales shorter than ∼20 ks with respect to the continuum), in agreement with what is expected for type 1 AGN (Ponti et al. 2012; Tortosa et al. 2023).

thumbnail Fig. 4.

Top panel: soft (0.2–1 keV) vs. hard (3–10 keV) σ NXS 2 $ \sigma^2_{\rm NXS} $. Bottom panel: medium (1–3 keV) vs. hard (3–10 keV) σ NXS 2 $ \sigma^2_{\rm NXS} $. The values are calculated for the XMM-Newton EPIC-pn observations of IC 4329A as in Table 1. The dashed black lines represent the one-to-one relation. The color code for the different observations is the same as in Fig. 1.

The NuSTAR light curves show no spectral variations(< 10%) in the ratio of the 3–10 and 10–70 keV count rates (see the bottom panels of Fig. 2).

Since the hardness ratios of the XMM-Newton and NuSTAR observations do not change strongly within each pointing (see the bottom panels of Figs. 1 and 2), we decided to use the time-averaged spectrum for our spectral analysis to improve the spectrum statistics.

4. Spectral analysis

We performed the spectral analysis using the XSPEC v.12.12.1 software package (Arnaud et al. 1996). The spectra obtained by the two NuSTAR modules (FPMA and FPMB) were fit simultaneously, with a typical cross-normalization constant lower than 5% (Madsen et al. 2015).

Errors and upper and/or lower limits were calculated using the Δχ2 = 2.71 criterion (corresponding to the 90% confidence level for one parameter), unless stated otherwise.

The Galactic column density at the position of the source (NH = 4.1 × 1020 cm−2; HI4PI Collaboration 2016) was always included in the fits, and it was modeled with the TBABS component (Wilms et al. 2000) in XSPEC, with NH kept frozen to the quoted value. We also assumed solar abundances unless stated otherwise.

Figure 5 shows the background-subtracted EPIC pn (0.3–10 keV), and FPMA-B spectra (3–75 keV), plotted with the corresponding X-ray background. In all figures, the spectra were corrected for the effective area of each detector.

thumbnail Fig. 5.

Background-subtracted XMM-Newton EPIC-pn and NuSTAR FPMA-B spectra and the corresponding background levels (shaded regions). The simultaneous XMM-Newton and NuSTAR observations are reported with the same color. The color code for different observations is the same as in Fig. 1. This color code is applied throughout the paper.

4.1. XMM-Newton data analysis

We started our spectral analysis with the 2022 EPIC-pn XMM-Newton spectra (in the 0.3–10 keV range). We fit the spectra separately and adopted a fitting model that was based on literature results. The model is composed of a neutral absorption component at the redshift of the source, a WA component (Steenbrugge et al. 2005), and an UFO component (Markowitz et al. 2006; Tombesi et al. 2011), XSPEC model: TBABS * ZTBABS * MTABLE(WACOMP) * MTABLE(UFOCOMP) * (BBODY + POWERLAW + ZGAUSS + ZGAUSS), HEREAFTER MODEL XMM. The two latter components were modeled with detailed grids computed using the photoionization code XSTAR (Kallman & Bautista 2001), with an input spectral energy distribution described by a power law with a photon index of Γ = 2. These tables consider standard solar abundances from Asplund et al. (2009), and take into account absorption lines and edges for all the metals characterized by an atomic number Z ≤ 30. Considering the typical values of turbulent velocity for WAs (Laha et al. 2014), we used a table computed considering a turbulent velocity of 100 km s−1 for the WA and a table computed considering a turbulent velocity of 1000 km s−1 for the UFO. We added a blackbody component to take into account any weak underlying soft X-ray spectral feature, such as the soft excess. The primary continuum was modeled by a simple power law. We focused on the XMM-Newton energy band range and therefore did not include reflection features at this stage. Two Gaussian lines were included to represent the iron Kα and Kβ lines. The line widths of both lines were free to vary, as well as their centroid energies. The fits obtained using this model are very good (see the left panels of Fig. 6). The values of the Fe Kα line parameters are shown in Table 2. The line width of the Kα line is consistent within the errors with the XMM-Newton resolution ∼200 eV. The line width of the Fe Kβ line is ∼80 − 90 eV because it fits a blend of lines that might also include Fe XXVI. The line width of the Fe Kα line in some observations is unconstrained, and we found only upper limits. When we tried to include a narrower core for the Fe Kα, by adding another ZGAUSS component with a line width fixed to 1 eV, we were not able to recover the line flux for both components. We either found upper limits for the narrower or for the broader component. This suggests that if a broad component of the Fe Kα line and a narrow core are present, they are unresolved during our X-ray monitoring campaign. To determine whether the broad unresolved Fe Kα component is due to relativistic reflection off the innermost regions of the accretion disk, we replaced the Gaussian component with the DISKLINE model. From this analysis, we found that a subdominant (EW ∼ 15 eV) relativistically broadened Fe Kα component it is required by the data in 50% of the observations of the monitoring, where we found an inner radius Rin ∼ 7Rg at which Rg = GMBH/c2.

thumbnail Fig. 6.

Data/model ratio for the EPIC-pn XMM-Newton observations (left panels), FPMA+B NuSTAR observations (central panels), and the X-ray broadband observations (right panels) of IC 4329A when Model XMM, Model Nu and Model A are applied, respectively (see Sect. 4).

Table 2.

Details of the iron Kα line parameters from the analysis of the XMM-Newton observations.

In the soft energy band (0.3–2 keV), we found a moderate soft-excess component with temperatures ranging from kTBB = 47 ± 15 eV to kTBB = 75 ± 6 eV and a neutral absorption component at the redshift of the source with NH = (2.67 ± 0.07)×1021 cm−2 over the 9 XMM-Newton observations. We also confirmed the presence of two X-ray ionized absorbing components, consistent with previous studies (Steenbrugge et al. 2005; Markowitz et al. 2006; Tombesi et al. 2011). These components are constant within their uncertainties during the monitoring. We found for the X-ray wind with turbulent velocity 100 km s−1, identified with a WA, an average column density NH = (1.49 ± 0.16)×1021 cm−2, an average ionization fraction log(ξ/erg s−1 cm) = 1.30 ± 0.07, and an average observed redshift of the absorption components z = −0.011 ± 0.005. The component with a turbulent velocity 1000 km s−1, identified with a UFO, shows an average column density NH = (1.97 ± 0.82)×1021 cm−2, an average ionization fraction log(ξ/erg s−1 cm) = 2.69 ± 0.03, and an average z = −0.177 ± 0.016.

4.2. NuSTAR data analysis

Before analyzing the broadband spectra, we first focused on the primary emission and on its reflected component using NuSTAR FPMA+B spectra (in the 3–75 keV range). During the fitting process, we left the NuSTAR FPMB cross-calibration constant free to vary. Again, we composed a model following the models presented in the literature (e.g., Brenneman et al. 2014a,b). We modeled the primary X-ray continuum with the CUTOFFPL model in XSPEC, which includes a power law with a high-energy exponential cutoff, and the reprocessed emission using the standard neutral reflection model XILLVER version [1.4.3] with log(ξ) fixed to zero (García et al. 2013) (XSPEC model: tbabs*(cutoffpl+ xillver), hereafter Model Nu). The cutoff energy, Ec, and the photon index, Γ, of the XILLVER component were linked those of the CUTOFFPL, and the reflection fraction, Rrefl, was forced to be negative. With these settings, XILLVER only reproduced the reflection component. The disk inclination angle was fixed to a value of 60° (Brenneman et al. 2014b). The data are very well reproduced by the model. The ratio of data and model are shown in the central panels of Fig. 6, in which we also report the χ r 2 = χ 2 $ \chi^2_r = \chi^2 $ = χ2/degrees of freedom (d.o.f.). In this analysis, we found a photon index of the primary power law ranging from Γ = 1.67 ± 0.04 to Γ = 1.74 ± 0.03 and a cutoff value Ec = 145 − 180 keV, in agreement with previous studies (e.g., Brenneman et al. 2014b). However, with respect to the results obtained by Brenneman et al. (2014b), we found a significantly lower value for the reflection fraction, Rrefl = 0.008 − 0.01, similar to the value found in the analysis of Suzaku and NuSTAR observations of IC 4329A by Ogawa et al. (2019). Together with the result of our analysis of the Fe Kα line with the XMM-Newton observations, which shows faint broad features (see Sect. 4.1), this suggests that the reflection from the inner disk is weak.

thumbnail Fig. 7.

Model A (see Sect. 4.3). The dotted black line represents the total model, which is composed of the primary X-ray emission (light blue line), the reprocessed component (dark blue line), and the soft excess component (green line). Absorption and emission features are included in the model.

We also tested a scenario in which the Fe Kα emission is slightly broadened by relativistic effects and comes from the accretion disk by substituting the XILLVER component with the RELXILL component (García et al. 2014; Dauser et al. 2014), which takes ionized relativistic reflection from an accretion disk illuminated by a hot corona into account. This model resulted in a significantly poorer fit ( χ r 2 $ \chi^2_r $ ∼ 1.4 − 1.7). Moreover, we were not able to constrain the black hole spin in any of the NuSTAR observations of this campaign, finding just a lower limit of a ≳ 0.1.

4.3. Broadband data analysis

Finally, we extended the analysis to include the whole 0.3–10 keV XMM-Newton EPIC-pn spectrum and the 3–75 keV NuSTAR FPMA/FPMB spectra for the simultaneous observations (see Table 1). For the spectral fitting, we used the following XSPEC model: TBABS * ZTBABS * MTABLE(WACOMP) * MTABLE(UFOCOMP) * (BBODY + CUTOFFPL + XILLVER) (hereafter Model A; see Fig. 7), which is the combination of the previous models. In Fig. 8 we show the unfolded data and ratio residuals for the X-ray broadband observations of IC 4329A when Model A was applied. The data/model ratio of this analysis and the resulting χ r 2 $ \chi^2_r $ are shown in the right panels of Fig. 6. We show the data, fitting model, and residuals ratio for all the simultaneous XMM-Newton and NuSTAR observations in Fig. 9.

thumbnail Fig. 8.

Unfolded data plus best-fitting model (upper panels) and ratio residuals (lower panels) for the X-ray broadband observations of IC 4329A when the best-fitting model (Model A) is applied. We show all the model components as in Fig. 7 (see Sect. 4 for more details).

thumbnail Fig. 9.

Left panels: data, fitting model, and residual ratio for XMM-Newton EPIC-pn (black) and NuSTAR FPMA (blue) and FPMB (magenta) spectra of the simultaneous observations of IC 4329A when Model A (see Sect. 4.3) is applied to the data. Right panels: Contour plots at 68%, 90%, and 99% resulting from the MCMC analysis of Model A applied to the broadband XMM-Newton and NuSTAR 0.3–75 keV spectra of IC 4329A. We show the outputs for the photon index (Γ), cutoff energy (Ecut [keV]), NuSTAR normalization (normnu [ph keV−1 cm−2 s−1]), iron abundance (AFe), and reflection fraction (Rrefl). We report from top to bottom the Obs.IDs 086209101+60702050002, 086209103+60702050004, 086209105+60702050006, 086209107+60702050008, and 086209109+60702050010.

Based on the analysis of the XMM-Newton data, we cannot distinguish different components of the Fe Kα emission. We therefore carried out additional tests to determine whether we might understand the origin of the Fe Kα line. The XILLVER model accounts for the narrow Fe Kα emission line, so that we included a Gaussian line for a broad component. This component appeared to be superfluous as its normalization was unconstrained, and we only found an upper limit of ∼10−7 ph cm−2 s−1 in all the observations, with an upper limit on the equivalent width in the range of 5–10 eV. We reported the best-fitting parameters for the X-ray primary continuum (Γ and Ec) and for the reflection component (AFe and Rrefl) in Table 3.

Table 3.

Best-fitting parameters for the X-ray broadband XMM-Newton plus NuSTAR spectra of IC 4329A.

The very low value of the reflection fraction and the faint broad component of the Fe Kα line in the fit of the broadband spectra support the hypothesis in which the disk reflection is weak and the narrow iron Kα signature comes from distant material. To further test this scenario, we replaced the XILLVER with the MYTORUS model (Murphy & Yaqoob 2009) to reproduce the reflection from distant neutral material (Model B1). The inclination angle of the torus was fixed at 30°. To reproduce the Compton-scattered continuum, we used a table1 in which the termination energy was fixed by default in the model and was 160 keV, the closest value to the cutoff energy found with Model A. To reproduce the fluorescent line spectra, we used another table2 with a line centroid offset Eoff = 0 keV and with the same termination energy as in the previous table. The column densities of the scattered and line components were linked and were free to vary. The normalization of the scattered and line components were tied to the normalization of the primary continuum, that is, the standard MYTORUS configuration (the so-called coupled reprocessor; see, e.g., Yaqoob 2012). The assumed geometry corresponds to a covering fraction of 0.5.

The statistical significance of the fits with Model B1 is similar to that of Model A. We found a value for the average column density of the torus of 2.6 × 1023 cm−2, consistent within the errors with the results obtained for all the observations of this campaign (see Table 3). The photon index values obtained for the primary continuum are consistent within the errors with those obtained with Model A, but we were not able to constrain the cutoff energies, most likely because of the degeneracies among the parameters. Moreover, MYTORUS does not include a cutoff energy, but a termination energy that was fixed to be 160 keV. We report the best-fitting parameters of Model B1 in Table 3.

Additionally, we tested the reflected emission from a dusty torus with a variable covering fraction (Model B2) using the RXTORUSD3 model (Paltani & Ricci 2017; Ricci & Paltani 2023), the first torus model to consider dusty gas. In these tables, the cutoff energy was fixed to 200 keV, which is consistent with the value of Ec we found in IC 4329A (see Table 3). The level of the statistical significance observed in the fits using Model B2 is comparable to that observed using Model B1. The values of the photon index and the equatorial column density of the torus are also comparable within the errors with those of Model B (see Table 3). With Model B2, we were able to measure the ratio of the inner to outer radius of the torus, that is, r/R = 0.69 ± 0.02, and the viewing angle i ∼ 45°. These values were constant in the observations. All the parameters values are reported in Table 3.

The results obtained with Models A, B1, and B2 are consistent with the iron Kα line originating from neutral Compton-thin (NH = 1022 − 1023 cm−2) material that does not produce a prominent Compton reflection.

We searched for possible degeneracies between the fitting parameters in Model A, performing a Monte Carlo Markov chain (MCMC) model using the XSPEC-EMCEE tool4. This is an implementation of the emcee code (Foreman-Mackey et al. 2013) to analyze X-ray spectra within XSPEC. We used 50 walkers with 10 000 iterations each, burning the first 1000. The walkers started at the best-fit values found in XSPEC, following a Gaussian distribution in each parameter, with the standard deviation set to the delta value of that parameter. In the right panels of Fig. 9 we show the contour plots resulting from the MCMC analysis of Model A applied to the broadband 0.3–75 keV spectra of IC 4329A.

To directly measure the coronal temperature parameter, we used the XILLVERCP versions of the XILLVER tables (Model C). This model assumed that the primary emission is due to the Comptonization of thermal disk photons in a hot corona, and the reflection spectrum was calculated by using a more physically motivated primary continuum, implemented with the analytical Comptonization model NTHCOMP (Zdziarski et al. 1996; Życki et al. 1999) instead of a simple cutoff power law. The seed photon temperature was fixed at 50 eV by default in this model. This value is consistent with the maximum disk temperature expected for a source with M BH = 6 . 8 1.1 + 1.2 × 10 7 M $ M_{\mathrm{BH}}=6.8^{+1.2}_{-1.1}\times10^7\,M_{\odot} $ (Bentz et al. 2023) (i.e., TWien = 32 ± 10 eV). As reported in Table 3, we were able to constrain the coronal temperature for only two observations. The electron temperature of the corona is thought to be related to the spectral cutoff energy, Ec ∼ 2 − 3 × kTe (Petrucci et al. 2000, 2001). Assuming this relation, the values we obtained with Model C for the electron temperature of the corona are consistent with the cutoff values that we obtained with Model A.

4.4. Reanalysis of previous observations

To better understand the origin of the Fe Kα line, the evolution of the spectral parameters of the primary continuum and of the soft energy components (e.g., soft-excess or outflowing components), we compared the results of our campaign with those obtained by previous observations of IC 4329A (see Table 4). This was done by reanalyzing the archival XMM-Newton (2001, 2003, 2018), NuSTAR (2012), and Suzaku (2007, 2012) data. In Fig. 11 we show the 2–10 keV flux of IC 4329A during the years, including the literature values from observations carried out in the past four decades. During this period, the 2–10 keV flux of IC 4329A has changed by a factor ∼2.1. The highest state was in 2001, when the flux was F2 − 10 = 1.75 × 10−10 erg cm−2 s−1, while in 2021, the source reached the lowest state with F2 − 10 = 8.15 × 10−11 erg cm−2 s−1. This is lower by ∼70% than the value of 2001.

Table 4.

Summary of previous XMM-Newton, NuSTAR, and Suzaku observations of IC 4329A.

thumbnail Fig. 11.

Historical 2–10 keV light curve of IC 4329A as observed with various X-ray satellites. In chronological order: Einstein, 1979 (Holt et al. 1989); EXOSAT, 1985 (Singh et al. 1991); Ginga, 1989 (Piro et al. 1990); BeppoSAX, 1998 (Dadina 2007); Chandra, 2001 (Shu et al. 2010); XMM-Newton 2001, 2003 (Nandra et al. 2007); Suzaku, 2007 (Mantovani et al. 2014); Swift/XRT, 2007 (Winter et al. 2009); AstroSat, 2017 (Tripathi et al. 2021); XMM-Newton 2022; medium flux of the observations analyzed in this work.

We applied the same fitting approach as outlined in Sect. 4.1 to the archival XMM-Newton observations and to the 2007 Suzaku observations, and we used Model A (see Sect. 4.3) to fit the simultaneous NuSTAR and Suzaku observations of 2012.

The continuum emission of all previous XMM-Newton observations and of the Suzaku observations of 2007 is well fit by a power law and a blackbody component. The latter was included to take the soft excess into account. The simultaneous NuSTAR and Suzaku observations of 2012 are well fit by Model A (see Sect. 4.3). From the analysis of these simultaneous observations, we found a primary power law characterized by a photon index Γ = 1.74 ± 0.01 and an exponential cutoff at Ec = 173 ± 13 keV, in agreement with what was found by Brenneman et al. (2014b). We also found an iron abundance AFe = 1.08 ± 0.22 and a low value of the reflection fraction: Rrefl = 0.009 ± 0.004. This agrees with the results obtained from our monitoring campaign (see Sect. 4.3 and Table 3). We report the values of the best-fitting parameters of the primary continuum for the XMM-Newton and Suzaku observation in Table 5. The data, fitting models, and residuals ratio are shown in Fig. 9. Moreover, all these observations, as well as the Suzaku observation from 2012, show a WA and a UFO component, except for the XMM-Newton observation from 2001, for which ionized absorbers could not be confirmed due to the low S/N. The values of the parameters of these components are reported in Table 6. During the 2007 Suzaku monitoring (five observations, each ∼5 ks long) the spectral parameter values are constant within the errors. We therefore report the average values. For a consistent analysis, we also applied the same model to the XMM-Newton observation of 2003 as in Tombesi et al. (2011) (i.e., a power-law component absorbed by Galactic column density, and an XSTAR table with a turbulent velocity of 1000 km s−1 to the 4–10 keV band). We found values similar to those reported in Tombesi et al. (2011): NH = (1.69 ± 0.23)×1022 cm−2, log(ξ/erg s−1 cm) = 4.23 ± 0.82, vw/c = 0.095 ± 0.004.

Table 5.

Best-fitting parameters of the primary continuum and details of the iron Kα line parameters for the analysis of the spectra of IC 4329A from previous XMM-Newton and Suzaku observations as in Table 4.

Table 6.

Spectral parameters, velocity, upper and lower limits on the radial location, and upper limits on the energetic parameters for each X-ray absorber during the years.

Additionally, we found the values of the equivalent width (EW) of the FeKα line obtained from previous observations and from the observations of the 2021 campaign to be inversely correlated with the X-ray luminosity in the 2–10 keV band. Fitting these two parameters with a linear regression, we found a marginally significant (2.8σ) anticorrelation with a Pearson correlation coefficient ρ = −0.664 (see Fig. 10). This trend is known as the X-ray Baldwin or Iwasawa–Taniguchi effect (hereafter, IT effect) (Iwasawa & Taniguchi 1993). This inverse correlation is similar to the classical Baldwin effect, in which the equivalent width of the [C IV] λ1550 emission line is inversely correlated with the UV continuum luminosity (Baldwin 1977). The IT effect has been found in a large sample of different objects studied with different instruments (Page et al. 2004; Jiang et al. 2006; Bianchi et al. 2007; Shu et al. 2010; Ricci et al. 2014). Several explanations have been proposed for the IT effect: a decrease in the covering factor of the material forming the FeKα line (Page et al. 2004; Ricci et al. 2013), the dependence of the luminosity on the ionization state of the material that produces the line (Nandra et al. 1997b; Nayakshin 2000), and the variability related to the nonsimultaneous reaction to flux changes in the continuum of the reprocessing material (Jiang et al. 2006). One notable distinction we found in the reanalysis of previous observations of IC 4329A is the lack of a distinct relativistically broadened FeKα component. Excluding the XMM-Newton observation of 2001 in which it was not possible to detect the FeKα line because the S/N ratio of the observation was too low, the other observations show an unresolved symmetric FeKα with a line width that is mostly consistent with the XMM-Newton resolution at 6.5 keV and a low value of the line width and equivalent width, similar to what was found in the 2021 XMM-NewtonNuSTAR monitoring (see Table 5).

thumbnail Fig. 10.

Equivalent width (EW) of the neutral iron line plotted against the 2–10 keV X-ray luminosity (L2 − 10). The black line shows the linear regressions, and the shaded black regions represent the combined 3σ error on the slope and normalization.

5. Discussion

5.1. Spectral variability

The continuum 2–10 keV flux of IC 4329A showed fluctuations over the past 40 years (Fig. 11), and we also observed significant variability during our campaign. The highest value of the observed 2–10 keV flux in our monitoring is higher by ∼30% than the lowest observed value (see Table 3). Consequently, we aimed to investigate whether the spectral properties of the sources are connected to its flux level. The photon index of the power law shows some evidence for variability between the different observations. Including the values from archival observations, we found that Γ shows a moderate correlation with the 2–10 keV flux with a Pearson correlation coefficient of 0.63, corresponding to a 1 − Pvalue of 98% (see the top left panel of Fig. 12). This agrees with the softer-when-brighter behavior that is typically observed in AGN (Sobolewska & Papadakis 2009; Trakhtenbrot et al. 2017). Specifically, it was found that as the 2–10 keV flux increases, the photon index tends to become steeper (Shemmer et al. 2006). This could be related to changes in the physical conditions of the accretion disk as the flux increases (Haardt & Maraschi 1991), or to the effect of pair production in the X-ray corona (Ricci et al. 2018).

thumbnail Fig. 12.

Top panels: Γ vs. F2 − 10 keV (left) and Γ vs. λEdd (right) relation for IC 4329A. The black line shows the linear regression, and the shaded black regions represent the combined 1σ error on the slope and normalization. The dotted, dotted dashed, and dashed lines in the right panel represent the relations from Shemmer et al. (2008), Brightman et al. (2013) and Trakhtenbrot et al. (2017). Bottom left panel: Ec vs. λEdd. Bottom right panel: Rrefl vs. λEdd. The plots include the data from our campaign and the data obtained from archival XMM-Newton, NuSTAR and Suzaku observations. Ec and Rrefl cannot be measured when ony XMM-Newton data are available. The Eddington ratio was computed using the bolometric luminosity estimated using the bolometric correction to the 2–10 keV X-ray luminosity (κ2 − 10) from Lusso et al. (2010).

Several studies in the past decades have shown that Γ is strongly correlated with λEdd (Lu & Yu 1999; Wang et al. 2004; Shemmer et al. 2006, 2008; Brightman et al. 2013; Trakhtenbrot et al. 2017). We searched for this correlation between the values of the photon index of IC 4392A found in our monitoring and from archival observations, with the Eddington ratio computed using the bolometric luminosity estimated using the bolometric correction to the 2–10 keV X-ray luminosity (κ2 − 10) from Lusso et al. (2010). We found a moderate Γ vs. λEdd correlation (Pearson correlation coefficient ρcorr = 0.62, corresponding to a 1 − Pvalue of 97%; see the top right panel of Fig. 12). The linear regression has a slope value of 1.69 ± 0.80, which is consistent within the errors with the values expected from the literature (e.g., Shemmer et al. 2008; Brightman et al. 2013; Trakhtenbrot et al. 2017), but with a lower normalization value of 1.52 ± 0.11.

This result implies that IC 4329A fits the picture in which the mass accretion rate drives the physical conditions of the AGN corona well, which causes the shape of the X-ray spectrum. Therefore, the cutoff energy is expected to show a correlation with the 2–10 keV flux/luminosity and with the Eddington ratio. In Tortosa et al. (2018), we did not find clear evidence of a significant correlation between Ec and λEdd in the sample of bright AGN observed at the time with NuSTAR. However, the sample was still small and did not allow us to exclude a relation like this. Ricci et al. (2018), using the values from the BAT AGN Spectroscopic Survey (BASS) X-ray catalog (Ricci et al. 2017), found a negative correlation in the Ec versus λEdd relation, suggesting that the Eddington ratio is the main parameter driving the cutoff energy in AGN. The high-energy cutoff values of IC 4329A are marginally different among the observations, but they are still consistent within the errors with previous literature values (Brenneman et al. 2014b) and with the median values of bright nearby AGN (Ricci et al. 2018). We observe no clear correlation in this monitoring between this parameter and λEdd (see the lower panels of Fig. 12).

In the different observations of our campaign, the reflection fraction is consistent with being constant within the errors (see the lower right panel of Fig. 12). This parameter shows a very low value, consistent with what has been observed in the analysis of Suzaku and NuSTAR observations of IC 4329A by Ogawa et al. (2019). This result, in addition to our analysis of the Fe Kα line based on XMM-Newton observations, which shows faint broad features (see Sect. 4.1), suggests that the reflection from the inner disk is weak. This interpretation is further supported by the fact that the narrow features of the iron line seen in the XMM-Newton observations are consistent with being associated with emission from distant neutral material.

IC 4329A also clearly shows a moderate soft-excess component. However, this component does not vary in the observations. Its temperature ranges from kTBB = 68 ± 5 eV to kTBB = 75 ± 6 eV, consistent within the error in the observations of this campaign and with the temperature found from analyzing the archival observations (see Table 4). The flux of the soft-excess component is constant in the observations of our campaign. It is (2.31 − 2.89)±0.75 × 10−4 ph cm−2 s−1.

5.2. Fe Kα line

The central energy of the Gaussian line in our analysis is always consistent within the errors with low-ionization or neutral material in all the observations of this campaign (see Table 2). IC 4329A is well known to exhibit a strong Fe Kα line (Piro et al. 1990) with a narrow core at 6.4 keV. A broad component, most likely produced in the inner part of the accretion disk and blurred by general relativistic effects, was found in some X-ray observations (Done et al. 2000; Dadina 2007; Brenneman et al. 2014b; Mantovani et al. 2014). From our modeling, just 50% of the observations of this campaign seems to require a relativistic broad component of the Fe Kα line.

Using the mean value of the measured σFe Kα without considering the upper limits (see Table 2), we computed a full width at half-maximum (FWHM) of the iron Kα line: FWHMFe Kα = 9078 ± 250 km s−1. With this value, we estimated the radius of the iron Kα line emission, assuming that the line width represents the Keplerian velocity emitting material, using the following equation:

R Fe K α = G M BH v 2 , $$ \begin{aligned} R_{\rm Fe~K\alpha }=\frac{G\,M_{\rm BH}}{v^2}, \end{aligned} $$(3)

where G is the gravitational constant, and MBH is the mass of the central SMBH (Netzer et al. 1990; Peterson et al. 2004), and assuming a spherical geometry and an isotropic velocity distribution. This results in the assumption of v = 3 / 2 v FWHM $ v=\sqrt{3}/2v_{FWHM} $. We obtained RFe Kα = (4.68 ± 0.14)×10−5 pc, which corresponds to ∼7Rg, consistent with the inner radius recovered in 50% of the observations when fitting the line with the DISKLINE model. However, with the mean value of the FWHM of the iron Kα line obtained by the reanalysis of previous observations, we obtain R Fe K α past = ( 0.015 ± 0.005 ) $ R_{\mathrm{Fe~K\alpha}}^{\mathrm{past}}=(0.015\pm 0.005) $ pc. As in Nandra (2006), Gandhi et al. (2015), we can compare the radial position of the iron Kα line emitting region with the optical BLR radius RHβ, inferred from Hβ reverberation studies and with the value of the dust sublimation radius Rsub (i.e., the inner radius of the dusty torus; Nenkova et al. 2008). Rsub is defined as

R sub = 0.4 ( L bol 10 45 erg s 1 ) 0.5 ( 1500 K T sub ) 2.6 [ pc ] , $$ \begin{aligned} R_{\rm sub}=0.4\left(\frac{L_{\rm bol}}{10^{45} \mathrm{erg\,s^{-1}}}\right)^{0.5}\left(\frac{1500\,K}{T_{\rm sub}}\right)^{2.6}\,[\mathrm{pc}], \end{aligned} $$(4)

where Lbol is the bolometric luminosity, and Tsub is the dust sublimation temperature, assumed to be the sublimation temperature of graphite grains (T ∼ 1500 K, Kishimoto et al. 2007). Considering a bolometric luminosity of log(Lbol/erg s−1) = 45.04 (Koss et al. 2017), we obtain Rsub = 0.42 pc while RHβ = 0.012 pc from both Bentz et al. 2009, 2023. This value is consistent with the BLR radius resolved by GRAVITY Collaboration (2024) (RBLR = 0.013 pc). GRAVITY Collaboration (2024) inferred also an inner BLR radius of RBLR, min = 0.0037 pc. Based on these values, the bulk of the Fe Kα emission line of IC 4329A appears to likely originate in the inner accretion disk or in the BLR.

5.3. X-ray outflows

During the fitting process, we took ionized absorption into account by including two XSTAR tables, as outlined in Sect. 4.1. The resulting parameters for the absorption components obtained by our analysis of the 2021 campaign, together with the results obtained for the archival observations, are listed in Table 6.

Following the same approach outlined in different studies (e.g., Tombesi et al. 2013; Serafinelli et al. 2019; Tortosa et al. 2022), we computed upper and lower bounds on the radial position of the absorbers. The lower limit on the radial location of the outflows can be placed by estimating the radius at which the observed velocity equals the escape velocity: rmin = 2GMBHv−2. The upper limit is placed using the definition of the ionization parameter ξ = Lion/nr2 (Tarter et al. 1969), where Lion is the unabsorbed ionizing luminosity emitted by the source between 1 Ryd and 1000 Ryd (1 Ryd = 13.6 eV), n is the number density of the absorbing material, and r is the distance from the central source. To apply this definition, we need to assume that the thickness of the absorber does not exceed the distance from the SMBH, that is, that the absorbers are somewhat compact (Crenshaw & Kraemer 2012; Tombesi et al. 2013), r max = L ion N H 1 ξ 1 , $ r_{\mathrm{max}} = L_{\mathrm{ion}}N_{\mathrm{H}}^{-1}\xi^{-1}, $ with NH and ξ the column density and ionization fraction of the absorber, respectively. We computed the ionizing luminosity using the LUMINOSITY task in XSPEC on the unabsorbed best-fit spectral model (see Sects. 4 and 4.4) between 13.6 eV and 13.6 keV for all the observations.

The values of the estimated radial location of the absorber components in IC 4329A are reported in Table 6. By comparing these values with the literature, it is possible to see that the outflow component with turbulent velocity 100 km s−1 (hereafter Wind 1) shows the typical upper and lower limits of the distance of the WAs for the type 1 Seyfert galaxies (Tombesi et al. 2013), while the outflow component with a turbulent velocity 1000 km s−1 (hereafter Wind 2) is within the range of the average locations of UFOs (∼3 × 10−4 − 3 × 10−2 pc; see Tombesi et al. 2012).

Regarding the outflow energetics, we computed the mass outflow rate = 4πrNHμmpCgvw (Crenshaw & Kraemer 2012), where r is the radial location of the absorber, NH is the equivalent hydrogen column density, μ is the mean atomic mass per proton (=1.4 for solar abundances), mp is the mass of the proton, Cg is the global covering factor (≃0.5, Tombesi et al. 2010), and vw is the radial velocity centroid. It is interesting to compare this quantity with the mass accretion rate of the source (see Table 6). The typical value of the mass outflow rate for sources accreting below or close to the Eddington limit is out ≳ 5 − 10%acc both for UFOs and slower outflows (Tombesi et al. 2012). Both X-ray absorbers agree with this scenario, in which the mechanical power is enough to exercise a significant feedback on the surrounding environment.

We also computed the value of the momentum rate of the outflow, that is, the rate at which the outflow transports momentum into the environment of the host galaxy, = out/vw. This quantity is compared in Table 6 with the momentum of the radiation of the source, which is defined as the ratio of the observed luminosity and the speed of light. For IC 4329A, this is log(rad/erg cm−1) = 35.58. Outflows accelerated through the continuum radiation pressure are expected to have a out/rad ≳ 1 (King & Pounds 2015). The median value of this ratio for UFOs is ∼0.96 after the relativistic correction and ∼0.64 without the relativistic corrections (Luminari et al. 2020). The values of out/rad we found for IC 4329A are in the range between 0.80 and 5.35. This suggests that radiation pressure accelerates the material to the escape velocity.

The last parameter we computed was the instantaneous kinetic power of the outflow, K ˙ = 1 2 M ˙ out v w 2 $ \dot{K}=\frac{1}{2}\dot{M}_{\mathrm{out}}v_{\mathrm{w}}^2 $, and we compared it with the outflowing observed bolometric luminosity of the source. Tombesi et al. (2012) showed that for a significant feedback in the environment of an AGN, a minimum ratio of the mechanical power of the outflow and the bolometric luminosity of ∼0.3% for UFOs and ∼0.8% for WAs is required. Based on the values found for IC 4329A, the source fits this scenario well, in which the outflowing winds can impress a strong feedback.

6. Conclusions

We have presented the detailed broadband analysis of the simultaneous XMM-NewtonNuSTAR monitoring of the Seyfert 1 galaxy IC 4329A that was carried out in 2021. The results of our analysis are listed below.

  • Our spectral analysis shows that the X-ray broadband spectra of IC 4329A from the 2021 X-ray monitoring campaign are well fit by a model that includes a soft-excess component, a WA, and an UFO, primary emission modeled by a power law with a cutoff at high energy, and a distant neutral reflection component accounting for the Fe fluorescence (see Fig. 7 and Sect. 4).

  • We found no significant variations in the XMM-Newton and NuSTAR hardness ratios of the source (lower than 10%; see Figs. 1, 2 and Sect. 3) within the single observations (∼20 ks). The variability spectra of IC 4329A during the campaign are nearly flat (see Fig. 3). The analysis of the excess variance values of the observations of this monitoring also suggests that IC 4329A shows very weak spectral variations overall (see Fig. 4).

  • We found a photon index of the primary power law ranging from Γ = 1.65 ± 0.02 to Γ = 1.77 ± 0.02 and a cutoff value Ec = 140 − 250 keV, which agrees with the results of Brenneman et al. (2014b). However, with respect to the results from Brenneman et al. (2014b), we found a lower value for the reflection fraction, Rrefl = 0.006 − 0.009, as found by the analysis of Suzaku and NuSTAR observations of IC 4329A by Ogawa et al. (2019). Together with the very faint broad component of the Fe Kα (see Sect. 4.1), this suggests that the reflection from the inner disk is weak.

  • Our analysis of the broadband spectra of IC 4329A suggests that if it is present, the narrow core of the iron Kα line likely originates from neutral Compton-thin (NH = 1022 − 1023 cm−2) material that does not produce a prominent Compton reflection. This is consistent with the low reflection fraction found in the analysis. Moreover, in this 2021 XMM-NewtonNuSTAR monitoring, we found an indication of a relativistically broadened iron Kα line in IC 4329A in 50% of the observations of the monitoring. However, the data from previous IC 4329A observations did not require an asymmetrical Fe Kα component. From an estimate of the position of the iron K emitting region (see Sect. 5.2), the Fe Kα emission line of IC 4329A appears to likely originate in the inner accretion disk or in the BLR. At the XMM-Newton resolution, the different components of the Fe Kα line are unfortunately unresolved. Future observations with a higher-resolution spectrometer (e.g., XRISM-Resolve) will be pivotal for distinguishing the different components of the low-ionization Fe Kα emission that is indicated in the current data.

  • Small spectral changes are observed that follow the softer-when-brighter behavior typically observed in unobscured AGN. This is most likely related to changes in the physical conditions of the accretion flow (see Sect. 5.1). We found a statistically significant Γ versus λEdd correlation (see the upper right panel of Fig. 12), indicating that the Eddington ratio may drive the physical conditions of the X-ray emitting corona.

  • We estimated the distance from the SMBH of the iron Kα line emitting region from the mean FWHM of the line. We find that the Fe Kα emission line of IC 4329A is consistent with originating in the inner disk (see Sect. 5.2).

  • From the analysis of the XMM-Newton and NuSTAR observations of IC 4329A from this campaign and from archival data, we found X-ray outflows composed of two phases: one UFO, and one WA with a lower velocity (see Sect. 5.3). From the energetics of these X-ray winds, we conclude that both components are powered by radiation pressure and that they could exert a significant feedback on the AGN environment.


Acknowledgments

C.R. acknowledges support from the Fondecyt Regular grant 1230345 and ANID BASAL project FB210003. E.S. acknowledges support from ANID BASAL project FB210003 and Gemini ANID ASTRO21-0003. T.K. is grateful for support from RIKEN Special Postdoctoral Researcher Program and is supported by JSPS KAKENHI grant number JP23K13153. This work is based on observations obtained with the ESA science mission XMM-Newton, with instruments and contributions directly funded by ESA Member States and the USA (NASA), the NuSTAR mission, a project led by the California Institute of Technology, managed by the Jet Propulsion Laboratory and funded by NASA. This research has made use of the NuSTAR Data Analysis Software (NuSTARDAS) jointly developed by the ASI Space Science Data Center (SSDC, Italy) and the California Institute of Technology (Caltech, USA). We thank the anonymous referee for his/her valuable suggestions.

References

  1. Andonie, C., Bauer, F. E., Carraro, R., et al. 2022, A&A, 664, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Arnaud, K. A. 1996, in Astronomical Data Analysis Software and Systems V, eds. G. H. Jacoby, & J. Barnes, ASP Conf. Ser., 101, 17 [NASA ADS] [Google Scholar]
  3. Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
  4. Baldwin, J. A. 1977, ApJ, 214, 679 [NASA ADS] [CrossRef] [Google Scholar]
  5. Beckmann, V., Gehrels, N., Shrader, C. R., & Soldi, S. 2006, ApJ, 638, 642 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bentz, M. C., Peterson, B. M., Netzer, H., Pogge, R. W., & Vestergaard, M. 2009, ApJ, 697, 160 [Google Scholar]
  7. Bentz, M. C., Onken, C. A., Street, R., & Valluri, M. 2023, ApJ, 944, 29 [NASA ADS] [CrossRef] [Google Scholar]
  8. Bianchi, S., Guainazzi, M., Matt, G., & Fonseca Bonilla, N. 2007, A&A, 467, L19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Bianchi, S., Guainazzi, M., Matt, G., Fonseca Bonilla, N., & Ponti, G. 2009, A&A, 495, 421 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Blustin, A. J., Page, M. J., Fuerst, S. V., Branduardi-Raymont, G., & Ashton, C. E. 2005, A&A, 431, 111 [CrossRef] [EDP Sciences] [Google Scholar]
  11. Brenneman, L. W., Madejski, G., Fuerst, F., et al. 2014a, ApJ, 781, 83 [CrossRef] [Google Scholar]
  12. Brenneman, L. W., Madejski, G., Fuerst, F., et al. 2014b, ApJ, 788, 61 [NASA ADS] [CrossRef] [Google Scholar]
  13. Brightman, M., Silverman, J. D., Mainieri, V., et al. 2013, MNRAS, 433, 2485 [Google Scholar]
  14. Crenshaw, D. M., & Kraemer, S. B. 2001, ApJ, 562, L29 [NASA ADS] [CrossRef] [Google Scholar]
  15. Crenshaw, D. M., & Kraemer, S. B. 2012, ApJ, 753, 75 [NASA ADS] [CrossRef] [Google Scholar]
  16. Dadina, M. 2007, A&A, 461, 1209 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Dauser, T., Garcia, J., Parker, M. L., Fabian, A. C., & Wilms, J. 2014, MNRAS, 444, L100 [Google Scholar]
  18. De Marco, B., Ponti, G., Cappi, M., et al. 2013, MNRAS, 431, 2441 [Google Scholar]
  19. De Marco, B., Adhikari, T. P., Ponti, G., et al. 2020, A&A, 634, A65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. den Herder, J. W., Brinkman, A. C., Kahn, S. M., et al. 2001, A&A, 365, L7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Dewangan, G. C., Tripathi, P., Papadakis, I. E., & Singh, K. P. 2021, MNRAS, 504, 4015 [NASA ADS] [CrossRef] [Google Scholar]
  22. Done, C., Madejski, G. M., & Życki, P. T. 2000, ApJ, 536, 213 [Google Scholar]
  23. Du, P., Zhang, Z.-X., Wang, K., et al. 2018, ApJ, 856, 6 [NASA ADS] [CrossRef] [Google Scholar]
  24. Edelson, R., Turner, T. J., Pounds, K., et al. 2002, ApJ, 568, 610 [CrossRef] [Google Scholar]
  25. Fabian, A. C., Zoghbi, A., Ross, R. R., et al. 2009, Nature, 459, 540 [Google Scholar]
  26. Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [Google Scholar]
  27. Fukazawa, Y., Furui, S., Hayashi, K., et al. 2016, ApJ, 821, 15 [NASA ADS] [CrossRef] [Google Scholar]
  28. Gabriel, C., Denby, M., Fyfe, D. J., et al. 2004, in Astronomical Data Analysis Software and Systems (ADASS) XIII, eds. F. Ochsenbein, M. G. Allen, & D. Egret, ASP Conf. Ser., 314, 759 [NASA ADS] [Google Scholar]
  29. Gandhi, P., Hönig, S. F., & Kishimoto, M. 2015, ApJ, 812, 113 [NASA ADS] [CrossRef] [Google Scholar]
  30. García, J., Dauser, T., Reynolds, C. S., et al. 2013, ApJ, 768, 146 [Google Scholar]
  31. García, J., Dauser, T., Lohfink, A., et al. 2014, ApJ, 782, 76 [Google Scholar]
  32. George, I. M., & Fabian, A. C. 1991, MNRAS, 249, 352 [Google Scholar]
  33. GRAVITY Collaboration (Amorim, A.et al.,) 2024, A&A, 684, A167 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Haardt, F., & Maraschi, L. 1991, ApJ, 380, L51 [Google Scholar]
  35. Haardt, F., & Maraschi, L. 1993, ApJ, 413, 507 [Google Scholar]
  36. Haardt, F., & Matt, G. 1993, MNRAS, 261, 346 [NASA ADS] [CrossRef] [Google Scholar]
  37. Haardt, F., Maraschi, L., & Ghisellini, G. 1994, ApJ, 432, L95 [NASA ADS] [CrossRef] [Google Scholar]
  38. Haardt, F., Maraschi, L., & Ghisellini, G. 1997, ApJ, 476, 620 [CrossRef] [Google Scholar]
  39. Harrison, F. A., Craig, W. W., Christensen, F. E., et al. 2013, ApJ, 770, 103 [Google Scholar]
  40. HI4PI Collaboration (Ben Bekhti, N., et al.) 2016, A&A, 594, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Ho, L. C. 2008, ARA&A, 46, 475 [Google Scholar]
  42. Holt, S. S., Turner, T. J., Mushotzky, R. F., & Weaver, K. 1989, in Two Topics in X-Ray Astronomy, Volume 1: X Ray Binaries. Volume 2: AGN and the X Ray Background, eds. J. Hunt, & B. Battrick, ESA Special Publ., 296, 1105 [NASA ADS] [Google Scholar]
  43. Ingram, A., Ewing, M., Marinucci, A., et al. 2023, MNRAS, 525, 5437 [CrossRef] [Google Scholar]
  44. Inoue, Y., & Doi, A. 2014, PASJ, 66, L8 [NASA ADS] [CrossRef] [Google Scholar]
  45. Iwasawa, K., & Taniguchi, Y. 1993, ApJ, 413, L15 [NASA ADS] [CrossRef] [Google Scholar]
  46. Jansen, F., Lumb, D., Altieri, B., et al. 2001, A&A, 365, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Jiang, P., Wang, J. X., & Wang, T. G. 2006, ApJ, 644, 725 [NASA ADS] [CrossRef] [Google Scholar]
  48. Kallman, T., & Bautista, M. 2001, ApJS, 133, 221 [Google Scholar]
  49. Kawamuro, T., Ricci, C., Imanishi, M., et al. 2022, ApJ, 938, 87 [NASA ADS] [CrossRef] [Google Scholar]
  50. Kelly, B. C. 2007, ApJ, 665, 1489 [Google Scholar]
  51. King, A., & Pounds, K. 2015, ARA&A, 53, 115 [NASA ADS] [CrossRef] [Google Scholar]
  52. Kishimoto, M., Hönig, S. F., Beckert, T., & Weigelt, G. 2007, A&A, 476, 713 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Koss, M., Trakhtenbrot, B., Ricci, C., et al. 2017, ApJ, 850, 74 [Google Scholar]
  54. Koss, M. J., Ricci, C., Trakhtenbrot, B., et al. 2022a, ApJS, 261, 2 [NASA ADS] [CrossRef] [Google Scholar]
  55. Koss, M. J., Trakhtenbrot, B., Ricci, C., et al. 2022b, ApJS, 261, 1 [NASA ADS] [CrossRef] [Google Scholar]
  56. Laha, S., Guainazzi, M., Dewangan, G. C., Chakravorty, S., & Kembhavi, A. K. 2014, MNRAS, 441, 2613 [NASA ADS] [CrossRef] [Google Scholar]
  57. Laor, A., & Behar, E. 2008, MNRAS, 390, 847 [NASA ADS] [CrossRef] [Google Scholar]
  58. Lu, Y., & Yu, Q. 1999, ApJ, 526, L5 [NASA ADS] [CrossRef] [Google Scholar]
  59. Luminari, A., Tombesi, F., Piconcelli, E., et al. 2020, A&A, 633, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Lusso, E., Comastri, A., Vignali, C., et al. 2010, A&A, 512, A34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Madsen, K. K., Harrison, F. A., Markwardt, C. B., et al. 2015, ApJS, 220, 8 [Google Scholar]
  62. Madsen, K. K., Grefenstette, B. W., Pike, S., et al. 2020, ArXiv e-prints [arXiv:2005.00569] [Google Scholar]
  63. Mantovani, G., Nandra, K., & Ponti, G. 2014, MNRAS, 442, L95 [Google Scholar]
  64. Markowitz, A., Reeves, J. N., & Braito, V. 2006, ApJ, 646, 783 [NASA ADS] [CrossRef] [Google Scholar]
  65. Matt, G., Fabian, A. C., & Reynolds, C. S. 1997, MNRAS, 289, 175 [NASA ADS] [CrossRef] [Google Scholar]
  66. Matzeu, G. A., Reeves, J. N., Nardini, E., et al. 2016, MNRAS, 458, 1311 [Google Scholar]
  67. Matzeu, G. A., Reeves, J. N., Nardini, E., et al. 2017, MNRAS, 465, 2804 [NASA ADS] [CrossRef] [Google Scholar]
  68. McHardy, I. M., Koerding, E., Knigge, C., Uttley, P., & Fender, R. P. 2006, Nature, 444, 730 [NASA ADS] [CrossRef] [Google Scholar]
  69. McKernan, B., & Yaqoob, T. 2004, ApJ, 608, 157 [NASA ADS] [CrossRef] [Google Scholar]
  70. Mehdipour, M., & Costantini, E. 2018, A&A, 619, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  71. Merloni, A. 2003, MNRAS, 341, 1051 [NASA ADS] [CrossRef] [Google Scholar]
  72. Molina, M., Bassani, L., Malizia, A., et al. 2013, MNRAS, 433, 1687 [NASA ADS] [CrossRef] [Google Scholar]
  73. Murphy, K. D., & Yaqoob, T. 2009, MNRAS, 397, 1549 [Google Scholar]
  74. Nandra, K. 2006, MNRAS, 368, L62 [NASA ADS] [CrossRef] [Google Scholar]
  75. Nandra, K., & Pounds, K. A. 1994, MNRAS, 268, 405 [Google Scholar]
  76. Nandra, K., George, I. M., Mushotzky, R. F., Turner, T. J., & Yaqoob, T. 1997a, ApJ, 476, 70 [Google Scholar]
  77. Nandra, K., George, I. M., Mushotzky, R. F., Turner, T. J., & Yaqoob, T. 1997b, ApJ, 488, L91 [NASA ADS] [CrossRef] [Google Scholar]
  78. Nandra, K., O’Neill, P. M., George, I. M., & Reeves, J. N. 2007, MNRAS, 382, 194 [NASA ADS] [CrossRef] [Google Scholar]
  79. Nayakshin, S. 2000, ApJ, 534, 718 [NASA ADS] [CrossRef] [Google Scholar]
  80. Nenkova, M., Sirocky, M. M., Ivezić, Ž., & Elitzur, M. 2008, ApJ, 685, 147 [Google Scholar]
  81. Netzer, H. 1990, in Active Galactic Nuclei, eds. R. D. Blandford, H. Netzer, L. Woltjer, T. J. L. Courvoisier, & M. Mayor, 57 [Google Scholar]
  82. Ogawa, S., Ueda, Y., Yamada, S., Tanimoto, A., & Kawaguchi, T. 2019, ApJ, 875, 115 [NASA ADS] [CrossRef] [Google Scholar]
  83. Oh, K., Koss, M. J., Ueda, Y., et al. 2022, ApJS, 261, 4 [NASA ADS] [CrossRef] [Google Scholar]
  84. Page, K. L., O’Brien, P. T., Reeves, J. N., & Turner, M. J. L. 2004, MNRAS, 347, 316 [NASA ADS] [CrossRef] [Google Scholar]
  85. Pal, I., Stalin, C. S., Chatterjee, R., & Agrawal, V. K. 2023, J. Astrophys. Astron., 44, 87 [NASA ADS] [CrossRef] [Google Scholar]
  86. Paltani, S., & Ricci, C. 2017, A&A, 607, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  87. Papadakis, I. E. 2004, MNRAS, 348, 207 [NASA ADS] [CrossRef] [Google Scholar]
  88. Perola, G. C., Matt, G., Cappi, M., et al. 2002, A&A, 389, 802 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  89. Peterson, B. M., Ferrarese, L., Gilbert, K. M., et al. 2004, ApJ, 613, 682 [Google Scholar]
  90. Petre, R., Mushotzky, R. F., Krolik, J. H., & Holt, S. S. 1984, ApJ, 280, 499 [NASA ADS] [CrossRef] [Google Scholar]
  91. Petrucci, P. O., Haardt, F., Maraschi, L., et al. 2000, ApJ, 540, 131 [Google Scholar]
  92. Petrucci, P. O., Haardt, F., Maraschi, L., et al. 2001, ApJ, 556, 716 [Google Scholar]
  93. Piconcelli, E., Jimenez-Bailón, E., Guainazzi, M., et al. 2004, MNRAS, 351, 161 [Google Scholar]
  94. Piro, L., Yamauchi, M., & Matsuoka, M. 1990, ApJ, 360, L35 [NASA ADS] [CrossRef] [Google Scholar]
  95. Ponti, G., Papadakis, I., Bianchi, S., et al. 2012, A&A, 542, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  96. Pounds, K. A., Nandra, K., Stewart, G. C., George, I. M., & Fabian, A. C. 1990, Nature, 344, 132 [NASA ADS] [CrossRef] [Google Scholar]
  97. Ricci, C., & Paltani, S. 2023, ApJ, 945, 55 [NASA ADS] [CrossRef] [Google Scholar]
  98. Ricci, C., Paltani, S., Awaki, H., et al. 2013, A&A, 553, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  99. Ricci, C., Ueda, Y., Paltani, S., et al. 2014, MNRAS, 441, 3622 [NASA ADS] [CrossRef] [Google Scholar]
  100. Ricci, C., Trakhtenbrot, B., Koss, M. J., et al. 2017, ApJS, 233, 17 [Google Scholar]
  101. Ricci, C., Ho, L. C., Fabian, A. C., et al. 2018, MNRAS, 480, 1819 [NASA ADS] [CrossRef] [Google Scholar]
  102. Ricci, C., Chang, C.-S., Kawamuro, T., et al. 2023, ApJ, 952, L28 [NASA ADS] [CrossRef] [Google Scholar]
  103. Salpeter, E. E. 1964, ApJ, 140, 796 [NASA ADS] [CrossRef] [Google Scholar]
  104. Serafinelli, R., Tombesi, F., Vagnetti, F., et al. 2019, A&A, 627, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  105. Shemmer, O., Brandt, W. N., Netzer, H., Maiolino, R., & Kaspi, S. 2006, ApJ, 646, L29 [Google Scholar]
  106. Shemmer, O., Brandt, W. N., Netzer, H., Maiolino, R., & Kaspi, S. 2008, ApJ, 682, 81 [Google Scholar]
  107. Shu, X. W., Yaqoob, T., & Wang, J. X. 2010, ApJS, 187, 581 [NASA ADS] [CrossRef] [Google Scholar]
  108. Singh, K. P., Rao, A. R., & Vahia, M. N. 1991, ApJ, 377, 417 [NASA ADS] [CrossRef] [Google Scholar]
  109. Sobolewska, M. A., & Papadakis, I. E. 2009, MNRAS, 399, 1597 [NASA ADS] [CrossRef] [Google Scholar]
  110. Steenbrugge, K. C., Kaastra, J. S., Sako, M., et al. 2005, A&A, 432, 453 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  111. Strüder, L., Briel, U., Dennerl, K., et al. 2001, A&A, 365, L18 [Google Scholar]
  112. Sunyaev, R. A., & Titarchuk, L. G. 1980, A&A, 86, 121 [NASA ADS] [Google Scholar]
  113. Tarter, C. B., Tucker, W. H., & Salpeter, E. E. 1969, ApJ, 156, 943 [Google Scholar]
  114. Tombesi, F., Cappi, M., Reeves, J. N., et al. 2010, A&A, 521, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  115. Tombesi, F., Cappi, M., Reeves, J. N., et al. 2011, ApJ, 742, 44 [Google Scholar]
  116. Tombesi, F., Cappi, M., Reeves, J. N., & Braito, V. 2012, MNRAS, 422, L1 [Google Scholar]
  117. Tombesi, F., Cappi, M., Reeves, J. N., et al. 2013, MNRAS, 430, 1102 [Google Scholar]
  118. Tombesi, F., Meléndez, M., Veilleux, S., et al. 2015, Nature, 519, 436 [NASA ADS] [CrossRef] [Google Scholar]
  119. Tortosa, A., Bianchi, S., Marinucci, A., Matt, G., & Petrucci, P. O. 2018, A&A, 614, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  120. Tortosa, A., Ricci, C., Tombesi, F., et al. 2022, MNRAS, 509, 3599 [Google Scholar]
  121. Tortosa, A., Ricci, C., Arévalo, P., et al. 2023, MNRAS, 526, 1687 [NASA ADS] [CrossRef] [Google Scholar]
  122. Trakhtenbrot, B., Ricci, C., Koss, M. J., et al. 2017, MNRAS, 470, 800 [NASA ADS] [CrossRef] [Google Scholar]
  123. Tripathi, P., Dewangan, G. C., Papadakis, I. E., & Singh, K. P. 2021, ApJ, 915, 25 [NASA ADS] [CrossRef] [Google Scholar]
  124. Turner, M. J. L., Abbey, A., Arnaud, M., et al. 2001, A&A, 365, L27 [CrossRef] [EDP Sciences] [Google Scholar]
  125. Vaughan, S., Edelson, R., Warwick, R. S., & Uttley, P. 2003, MNRAS, 345, 1271 [Google Scholar]
  126. Vaughan, S., Fabian, A. C., Ballantyne, D. R., et al. 2004, MNRAS, 351, 193 [NASA ADS] [CrossRef] [Google Scholar]
  127. Wang, J.-M., Watarai, K.-Y., & Mineshige, S. 2004, ApJ, 607, L107 [NASA ADS] [CrossRef] [Google Scholar]
  128. Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914 [Google Scholar]
  129. Winter, L. M., Mushotzky, R. F., Reynolds, C. S., & Tueller, J. 2009, ApJ, 690, 1322 [NASA ADS] [CrossRef] [Google Scholar]
  130. Yaqoob, T. 2012, MNRAS, 423, 3360 [Google Scholar]
  131. Zdziarski, A. A., Johnson, W. N., & Magdziarz, P. 1996, MNRAS, 283, 193 [NASA ADS] [CrossRef] [Google Scholar]
  132. Zoghbi, A., Fabian, A. C., Reynolds, C. S., & Cackett, E. M. 2012, MNRAS, 422, 129 [NASA ADS] [CrossRef] [Google Scholar]
  133. Życki, P. T., Done, C., & Smith, D. A. 1999, MNRAS, 309, 561 [Google Scholar]

All Tables

Table 1.

Summary of the 2021 XMM-Newton and NuSTAR observations of IC 4329A.

Table 2.

Details of the iron Kα line parameters from the analysis of the XMM-Newton observations.

Table 3.

Best-fitting parameters for the X-ray broadband XMM-Newton plus NuSTAR spectra of IC 4329A.

Table 4.

Summary of previous XMM-Newton, NuSTAR, and Suzaku observations of IC 4329A.

Table 5.

Best-fitting parameters of the primary continuum and details of the iron Kα line parameters for the analysis of the spectra of IC 4329A from previous XMM-Newton and Suzaku observations as in Table 4.

Table 6.

Spectral parameters, velocity, upper and lower limits on the radial location, and upper limits on the energetic parameters for each X-ray absorber during the years.

All Figures

thumbnail Fig. 1.

XMM-Newton EPIC-pn light curves (background subtracted) in the 0.2–10 keV, 0.2–2.5 keV and in the 2.5–10 keV are shown in the top, top middle, and bottom middle panels, respectively. In the bottom panels, the ratios of the XMM-Newton EPIC-pn background-subtracted light curves in the 2.5–10 keV and 0.2–2.5 keV energy bands are reported. All the light curves are extracted using a binning time of 500 s. The solid and dashed black lines indicate the average and the standard error of the mean, respectively. The color code used in this figure for different ObsID is applied throughout the paper.

In the text
thumbnail Fig. 2.

NuSTAR FPMA+B background-subtracted light curve in the 3–70 keV. 3–10 keV and 10–70 keV energy ranges are shown in the top, top middle, and bottom middle panels, respectively. In the bottom panels, the ratios of the NuSTAR NuSTAR background-subtracted light curves in the 10–70 keV and 3–10 keV energy bands are reported. All the light curves are extracted using a binning time of 500 s. The color code used in this figure for different ObsID is applied throughout the paper. The solid and dashed black lines indicate the average and the standard error of the mean, respectively.

In the text
thumbnail Fig. 3.

Fractional variability spectra for the XMM-Newton EPIC-pn observations of IC 4329A as in Table 1. The spectra of this campaign are characterized by a mostly flat shape, with a fractional variability that never exceeds 5%. The dot-dashed black lines are the linear regressions, and the shaded regions represent the combined 3σ error on the slope and normalization.

In the text
thumbnail Fig. 4.

Top panel: soft (0.2–1 keV) vs. hard (3–10 keV) σ NXS 2 $ \sigma^2_{\rm NXS} $. Bottom panel: medium (1–3 keV) vs. hard (3–10 keV) σ NXS 2 $ \sigma^2_{\rm NXS} $. The values are calculated for the XMM-Newton EPIC-pn observations of IC 4329A as in Table 1. The dashed black lines represent the one-to-one relation. The color code for the different observations is the same as in Fig. 1.

In the text
thumbnail Fig. 5.

Background-subtracted XMM-Newton EPIC-pn and NuSTAR FPMA-B spectra and the corresponding background levels (shaded regions). The simultaneous XMM-Newton and NuSTAR observations are reported with the same color. The color code for different observations is the same as in Fig. 1. This color code is applied throughout the paper.

In the text
thumbnail Fig. 6.

Data/model ratio for the EPIC-pn XMM-Newton observations (left panels), FPMA+B NuSTAR observations (central panels), and the X-ray broadband observations (right panels) of IC 4329A when Model XMM, Model Nu and Model A are applied, respectively (see Sect. 4).

In the text
thumbnail Fig. 7.

Model A (see Sect. 4.3). The dotted black line represents the total model, which is composed of the primary X-ray emission (light blue line), the reprocessed component (dark blue line), and the soft excess component (green line). Absorption and emission features are included in the model.

In the text
thumbnail Fig. 8.

Unfolded data plus best-fitting model (upper panels) and ratio residuals (lower panels) for the X-ray broadband observations of IC 4329A when the best-fitting model (Model A) is applied. We show all the model components as in Fig. 7 (see Sect. 4 for more details).

In the text
thumbnail Fig. 9.

Left panels: data, fitting model, and residual ratio for XMM-Newton EPIC-pn (black) and NuSTAR FPMA (blue) and FPMB (magenta) spectra of the simultaneous observations of IC 4329A when Model A (see Sect. 4.3) is applied to the data. Right panels: Contour plots at 68%, 90%, and 99% resulting from the MCMC analysis of Model A applied to the broadband XMM-Newton and NuSTAR 0.3–75 keV spectra of IC 4329A. We show the outputs for the photon index (Γ), cutoff energy (Ecut [keV]), NuSTAR normalization (normnu [ph keV−1 cm−2 s−1]), iron abundance (AFe), and reflection fraction (Rrefl). We report from top to bottom the Obs.IDs 086209101+60702050002, 086209103+60702050004, 086209105+60702050006, 086209107+60702050008, and 086209109+60702050010.

In the text
thumbnail Fig. 11.

Historical 2–10 keV light curve of IC 4329A as observed with various X-ray satellites. In chronological order: Einstein, 1979 (Holt et al. 1989); EXOSAT, 1985 (Singh et al. 1991); Ginga, 1989 (Piro et al. 1990); BeppoSAX, 1998 (Dadina 2007); Chandra, 2001 (Shu et al. 2010); XMM-Newton 2001, 2003 (Nandra et al. 2007); Suzaku, 2007 (Mantovani et al. 2014); Swift/XRT, 2007 (Winter et al. 2009); AstroSat, 2017 (Tripathi et al. 2021); XMM-Newton 2022; medium flux of the observations analyzed in this work.

In the text
thumbnail Fig. 10.

Equivalent width (EW) of the neutral iron line plotted against the 2–10 keV X-ray luminosity (L2 − 10). The black line shows the linear regressions, and the shaded black regions represent the combined 3σ error on the slope and normalization.

In the text
thumbnail Fig. 12.

Top panels: Γ vs. F2 − 10 keV (left) and Γ vs. λEdd (right) relation for IC 4329A. The black line shows the linear regression, and the shaded black regions represent the combined 1σ error on the slope and normalization. The dotted, dotted dashed, and dashed lines in the right panel represent the relations from Shemmer et al. (2008), Brightman et al. (2013) and Trakhtenbrot et al. (2017). Bottom left panel: Ec vs. λEdd. Bottom right panel: Rrefl vs. λEdd. The plots include the data from our campaign and the data obtained from archival XMM-Newton, NuSTAR and Suzaku observations. Ec and Rrefl cannot be measured when ony XMM-Newton data are available. The Eddington ratio was computed using the bolometric luminosity estimated using the bolometric correction to the 2–10 keV X-ray luminosity (κ2 − 10) from Lusso et al. (2010).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.