REJ1034+396 is one of the few active galactic nuclei with a significant quasi-periodic oscillation (QPO). The QPO has been observed in over 1 Ms of XMM-Newton observations spanning over a decade. We investigate the power spectral density function (PSD) of 7 long ( 90 ks) XMM-Newton observations of the active galactic nucleus REJ1034+396 in two energy bands. The soft (0.3-0.5 keV) band targets emission from the disk, while the hard (2-7 keV) band isolates the primary X-ray continuum emission from the corona. The QPO is significantly detected in the hard band of 5 of the 7 observations. The best fitting models indicate that the QPO detection in both bands is entirely attributable to the coronal emission with no additional contribution from the disk. This explains the strong coherence between the hard and soft bands at the QPO frequency. The covariance spectrum is consistent with this picture as the variability at QPO frequencies is attributed solely to fluctuations in the hot corona. The time lag as a function of energy is well described by a 2000 s intrinsic soft lag, resulting from the disk responding to emission from the corona, that undergoes phase wrapping at approximately the QPO frequency. By demonstrating that in this system the QPO arises in the corona, we provide new insights into the mechanisms generating QPOs.

Ever since the discovery of the first active galactic nuclei (AGN), substantial observational and theoretical effort has been invested into understanding how massive black holes have evolved across cosmic time. Circum-nuclear obscuration is now established as a crucial component, with almost every AGN observed known to display signatures of some level of obscuration in their X-ray spectra. However, despite more than six decades of effort, substantial open questions remain: how does the accretion power impact the structure of the circum-nuclear obscurer? What are the dynamical properties of the obscurer? Can dense circum-nuclear obscuration exist around intrinsically weak AGN? How many intermediate mass black holes occupy the centers of dwarf galaxies? In this paper, we showcase a number of next-generation prospects attainable with the High-Energy X-ray Probe (HEX-P1) to contribute toward solving these questions in the 2030s. The uniquely broad (0.2–80 keV) and strictly simultaneous X-ray passband of HEX-P makes it ideally suited for studying the temporal co-evolution between the central engine and circum-nuclear obscurer. Improved sensitivities and reduced background will enable the development of spectroscopic models complemented by current and future multi-wavelength observations. We show that the angular resolution of HEX-P both below and above 10 keV will enable the discovery and confirmation of accreting massive black holes at both low accretion power and low black hole masses even when concealed by thick obscuration. In combination with other next-generation observations of the dusty hearts of nearby galaxies, HEX-P will be pivotal in paving the way toward a complete picture of black hole growth and galaxy co-evolution.

Future strategic X-ray astronomy missions will require unprecedentedly sensitive wide-field imagers providing high frame rates, low readout noise and excellent soft energy response. To meet these needs, our team is employing a multi-pronged approach to advance several key areas of technology. Our first focus is on advanced readout electronics, specifically integrated electronics, where we are collaborating on the VERITAS readout chip for the Athena Wide Field Imager, and have developed the Multi-Channel Readout Chip (MCRC), which enables fast readout and high frame rates for MIT-LL JFET (junction field effect transistor) CCDs. Second, we are contributing to novel detector development, specifically the SiSeRO (Single electron Sensitive Read Out) devices fabricated at MIT Lincoln Laboratory, and their advanced readout, to achieve sub-electron noise performance. Hardware components set the stage for performance, but their efficient utilization relies on software and algorithms for signal and event processing. Our group is developing digital waveform filtering and AI methods to augment detector performance, including enhanced particle background screening and improved event characterization. All of these efforts make use of an efficient, new X-ray beamline facility at Stanford, where components and concepts can be tested and characterized.

Bringing artificial intelligence (AI) alongside next-generation X-ray imaging detectors, including CCDs and DEPFET sensors, enhances their sensitivity to achieve many of the flagship science cases targeted by future X-ray observatories, based upon low surface brightness and high redshift sources. Machine learning algorithms operating on the raw frame-level data provide enhanced identification of background vs. astrophysical X-ray events, by considering all of the signals in the context within which they appear within each frame. We have developed prototype machine learning algorithms to identify valid X-ray and cosmic-ray induced background events, trained and tested upon a suite of realistic end-to-end simulations that trace the interaction of cosmic ray particles and their secondaries through the spacecraft and detector. These algorithms demonstrate that AI can reduce the unrejected instrumental background by up to 41.5 per cent compared with traditional filtering methods. Alongside AI algorithms to reduce the instrumental background, next-generation event reconstruction methods, based upon fitting physically-motivated Gaussian models of the charge clouds produced by events within the detector, promise increased accuracy and spectral resolution of the lowest energy photon events.

In this study, we examine 51 archival NICER observations and six archival NuSTAR observations of the neutron star ultracompact X-ray binary (UCXB) 4U 0614+091, which span over 5 yr. The source displays persistent reflection features, so we use a reflection model designed for UCXBs, with overabundant carbon and oxygen (xillverCO) to study how various components of the system vary over time. The flux of this source is known to vary quasiperiodically on a timescale of a few days, so we study how the various model components change as the overall flux varies. The flux of most components scales linearly with the overall flux, while the power law, representing coronal emission, is anticorrelated, as expected. This is consistent with previous studies of the source. We also find that during observations of the high-soft state, the disk emissivity profile as a function of radius becomes steeper. We interpret this as the corona receding to be closer to the compact object during these states, at which point the assumed power law illumination of xillverCO may be inadequate to describe the illumination of the disk.

Traditional cosmic ray filtering algorithms used in X-ray imaging detectors aboard space telescopes perform event reconstruction based on the properties of activated pixels above a certain energy threshold, within 3x3 or 5x5 pixel sliding windows. This approach can reject up to 98% of the cosmic ray background. However, the remaining unrejected background constitutes a significant impediment to studies of low surface brightness objects, which are especially prevalent in the high-redshift universe. The main limitation of the traditional filtering algorithms is their ignorance of the long-range contextual information present in image frames. This becomes particularly problematic when analyzing signals created by secondary particles produced during interactions of cosmic rays with body of the detector. Such signals may look identical to the energy deposition left by X-ray photons, when one considers only the properties within the small sliding window. Additional information is present, however, in the spatial and energy correlations between signals in different parts of the same frame, which can be accessed by modern machine learning (ML) techniques. In this work, we continue the development of an ML-based pipeline for cosmic ray background mitigation. Our latest method consist of two stages: first, a frame classification neural network is used to create class activation maps (CAM), localizing all events within the frame; second, after event reconstruction, a random forest classifier, using features obtained from CAMs, is used to separate X-ray and cosmic ray features. The method delivers > 40% relative improvement over traditional filtering in background rejection in standard 0.3-10 keV energy range, at the expense of only a small (< 2%) level of lost X-ray signal. Our method also provides a convenient way to tune the cosmic ray rejection threshold to adapt to a user's specific scientific needs.

Constraining the primary growth channel of supermassive black holes (SMBHs) remains one the most actively debated questions in the context of cosmological structure formation. Owing to the expected connection between SMBH spin parameter evolution and the accretion and merger history of individual black holes, population spin measurements offer a rare observational window into the cosmic growth of SMBHs. As of today, the most common method for estimating SMBH spin relies on modeling the relativistically broaden atomic profiles in the reflection spectrum observed in X-rays. In this paper, we study the observational requirements needed to confidently distinguish between the primary SMBH growth channels based on their distinct spin-mass distributions predicted by the Horizon-AGN cosmological simulation. Indoing so, we characterize outstanding limitations associated with the existing measurements and discuss the landscape of future observational campaigns which could be planned and executed with future X-ray observatories. We focus our attention on the High-Energy X-ray Probe (HEX-P), a proposed probe-class mission designed to serve the high-energy community in the 2030s.

The hard X-ray emission in active galactic nuclei (AGN) and black hole X-ray binaries is thought to be produced by a hot cloud of electrons referred to as the corona. This emission, commonly described by a power law with a high-energy cutoff, is suggestive of Comptonization by thermal electrons. While several hypotheses have been proposed to explain the origin, geometry, and composition of the corona, we still lack a clear understanding of this fundamental component. NuSTAR has been playing a key role improving our knowledge of X-ray coronae thanks to its unprecedented sensitivity above 10 keV. However, these constraints are limited to bright, nearby sources. The High Energy X-ray Probe (HEX-P) is a probe-class mission concept combining high spatial resolution X-ray imaging and broad spectral coverage (0.2-80 keV) with a sensitivity superior to current facilities. In this paper, we highlight the major role that HEX-P will play in further advancing our insights of X-ray coronae, notably in AGN. We demonstrate how HEX-P will measure key properties and track the temporal evolution of coronae in unobscured AGN. This will allow us to determine their electron distribution and test the dominant emission mechanisms. Furthermore, we show how HEX-P will accurately estimate the coronal properties of obscured AGN in the local Universe, helping address fundamental questions about AGN unification. In addition, HEX-P will characterize coronae in a large sample of luminous quasars at cosmological redshifts for the first time and track the evolution of coronae in transient systems in real time. We also demonstrate how HEX-P will enable estimating the coronal geometry using spectral-timing techniques. HEX-P will thus be essential to understand the evolution and growth of black holes over a broad range of mass, distance, and luminosity, and will help uncover the black holes' role in shaping the Universe.

Galactic cosmic ray (GCR) particles have a significant impact on the particle-induced background of X-ray observatories, and their flux exhibits substantial temporal variability, potentially influencing background levels. In this study, we present one-day binned high-energy reject rates derived from the Chandra-ACIS and XMM-Newton EPIC-pn instruments, serving as proxies for GCR particle flux. We systematically analyze the ACIS and EPIC-pn reject rates and compare them with the AMS proton flux. Our analysis initially reveals robust correlations between the AMS proton flux and the ACIS/EPIC-pn reject rates when binned over 27-day intervals. However, a closer examination reveals substantial fluctuations within each 27-day bin, indicating shorter-term variability. Upon daily binning, we observe finer. temporal structures in the datasets, demonstrating the presence of recurrent variations with periods of ∼ 25 days and 23 days in ACIS and EPIC-pn reject rates, respectively, spanning the years 2014 to 2018. Notably, during the 2016–2017 period, we additionally detect periodicities of ∼13.5 days and 9 days in the ACIS and EPIC-pn reject rates, respectively. Intriguingly, we observe a time lag of ∼ 6 days between the AMS proton flux and the ACIS/EPIC-pn reject rates during the second half of 2016. This time lag is not visible before 2016 and aftern2017. The underlying physical mechanisms responsible for this time lag remain a subject of ongoing investigation.

Aims. We present an analysis of the Hubble Space Telescope Cosmic Origins Spectrograph spectrum of I Zw 1 aiming to probe the absorbing medium associated with the active galactic nucleus (AGN). Methods. We fitted the emission spectrum and performed spectral analysis of the identified absorption features to derive the corresponding ionic column densities and covering fractions of the associated outflows. We employed photoionisation modelling to constrain the total column density and the ionisation parameter of four detected kinematic components. By investigating the implications of the results together with the observed kinematic properties of both emission and absorption features, we derived constraints on the structure and geometry of the absorbing medium in the AGN environment. Results. We find and characterise absorption line systems from outflowing ionised gas in four distinct kinematic components, located at −60, −280, −1950, and −2900 km s −1 with respect to the source rest frame. While the two slower outflows are consistent with a full covering of the underlying radiation source, the well-constrained doublet line ratios of the faster two, higher column density, outflows suggest partial covering, with a covering fraction of C f ∼ 0.4. The faster outflows show also line-locking in the N V doublet, a signature of acceleration via line absorption. This makes I Zw 1 possibly the closest object that shows evidence for hosting line-driven winds. The observed −1950 km s −1 absorption is likely due to the same gas as an X-ray warm absorber. Furthermore, the behaviour in UV and X-ray bands implies that this outflow has a clumpy structure. We find that the highly asymmetric broad emission lines in I Zw 1, indicative of a collimated, outflowing broad line region, are covered by the absorbing gas. Finally, the strongest UV–X-ray absorber may be connected to some of the blueshifted line emission, indicative of a more spatially extended structure of this ionised medium.

Galactic cosmic-ray (GCR) particles have a significant impact on the particle-induced background of X-ray observatories, and their flux exhibits substantial temporal variability, potentially influencing background levels. In this study, we present 1 day binned high-energy reject rates derived from the Chandra-ACIS and XMM-Newton EPIC-pn instruments, serving as proxies for the GCR particle flux. We systematically analyze the ACIS and EPIC-pn reject rates and compare them with the AMS proton flux. Our analysis initially reveals robust correlations between the AMS proton flux and the ACIS/EPIC-pn reject rates when binned over 27 day intervals. However, a closer examination reveals substantial fluctuations within each 27 day bin, indicating shorter-term variability. Upon daily binning, we observe finer temporal structures in the data sets, demonstrating the presence of recurrent variations with periods of ∼25 days and 23 days in the ACIS and EPIC-pn reject rates, respectively, spanning the years 2014–2018. Notably, during the 2016–2017 period, we additionally detect periodicities of ∼13.5 days and 9 days in the ACIS and EPIC-pn reject rates, respectively. Intriguingly, we observe a time lag of ∼6 days between the AMS proton flux and the ACIS/EPIC-pn reject rates during the second half of 2016. This time lag is not visible before 2016 and after 2017. The underlying physical mechanisms responsible for this time lag remain a subject of ongoing investigation.

Traditional image segmentation methods employed with X-ray imaging detectors aboard X-ray space telescopes consist of two stages: first, a low energy threshold is applied; groups of activated pixels are then classified according to their shapes and identified as valid X-ray events or rejected as being possibly induced by cosmic rays. This method is fast and removes up to 98% of the cosmic ray-induced background. However, these traditional methods fail to address two important problems: first, they struggle to recover the true energies of, and sometimes fail to detect entirely, low-energy photons (photon energies less than 0.5keV); second, they consider only the shape of the active pixel regions, ignoring the longer-range context within the image frames. This limits their sensitivity to a specific type of cosmic ray signal: "islands" created by secondary particles produced by cosmic rays hitting the body of the telescope (the shapes of which are often indistinguishable from X-ray photon signals). Together, these limitations hinder investigations of faint, diffuse targets, such as the outskirts of galaxies and galaxy clusters, and of "low energy" targets such as individual stars, galaxies and high redshift systems. Both limitations can, however, be addressed with machine learning (ML) models. This work is part of our effort to develop fast and efficient background reduction methods for future astronomical X-ray missions using ML methods. We highlight several significant improvements in the classification and semantic segmentation of our background filtering pipeline. Our more realistic training and test data now incorporate the effects of readout noise and charge diffusion. In the presence of charge diffusion, our model is able to obtain an 80% relative improvement in lost signal recovery compared to the traditional background reduction techniques. We identify several directions for further development of the model.

We present the first results from a 100-day Swift, NICER and ground-based X-ray/UV/optical reverberation mapping campaign of the Narrow-Line Seyfert 1 Mrk 335, when it was in an unprecedented low X-ray flux state. Despite dramatic suppression of the X-ray variability, we still observe UV/optical lags as expected from disk reverberation. Moreover, the UV/optical lags are consistent with archival observations when the X-ray luminosity was >10 times higher. Interestingly, both low- and high-flux states reveal UV/optical lags that are 6-11 times longer than expected from a thin disk. These long lags are often interpreted as due to contamination from the broad line region, however the u band excess lag (containing the Balmer jump from the diffuse continuum) is less prevalent than in other AGN. The Swift campaign showed a low X-ray-to-optical correlation (similar to previous campaigns), but NICER and ground-based monitoring continued for another two weeks, during which the optical rose to the highest level of the campaign, followed 10 days later by a sharp rise in X-rays. While the low X-ray countrate and relatively large systematic uncertainties in the NICER background make this measurement challenging, if the optical does lead X-rays in this flare, this indicates a departure from the zeroth-order reprocessing picture. If the optical flare is due to an increase in mass accretion rate, this occurs on much shorter than the viscous timescale. Alternatively, the optical could be responding to an intrinsic rise in X-rays that is initially hidden from our line-of-sight.

UV and optical continuum reverberation mapping is a powerful tool for probing the accretion disk and inner broad-line region. However, recent reverberation mapping campaigns in the X-ray, UV, and optical have found lags consistently longer than those expected from the standard disk reprocessing picture. The largest discrepancy to date was recently reported in Mrk 335, where UV/optical lags are up to 12 times longer than expected. Here, we perform a frequency-resolved time lag analysis of Mrk 335, using Gaussian processes to account for irregular sampling. For the first time, we compare the Fourier frequency-resolved lags directly to those computed using the popular interpolated cross-correlation function method applied to both the original and detrended light curves. We show that the anticipated disk reverberation lags are recovered by the Fourier lags when zeroing in on the short-timescale variability. This suggests that a separate variability component is present on long timescales. If this separate component is modeled as reverberation from another region beyond the accretion disk, we constrain a size scale of roughly 15 lt-days from the central black hole. This is consistent with the size of the broad-line region inferred from H β reverberation lags. We also find tentative evidence for a soft X-ray lag, which we propose may be due to light travel time delays between the hard X-ray corona and distant photoionized gas that dominates the soft X-ray spectrum below 2 keV.

Spectral timing analyses based upon wavelet transforms provide a new means to study the variability of the X-ray emission from accreting systems, including AGN, stellar mass black holes and neutron stars, and can be used to trace the time variability of X-ray reverberation from the inner accretion disc. The previously-missing iron K reverberation time lags in the AGN IRAS 13224-3809 and MCG-6-30-15 are detected and found to be transitory in nature. Reverberation can be hidden during periods in which variability in the iron K band becomes dominated by ultrafast outflows (UFO). Following the time evolution of the reverberation lag between the corona and inner accretion disc, we may observe the short-timescale increase in scale height of the corona as it is accelerated away from the accretion disc during bright X-ray flares in the AGN I Zw 1. Measuring the variation of the reverberation lag that corresponds to the continuous, stochastic variations of the X-ray luminosity sheds new light on the disc-corona connection around accreting black holes. Hysteresis is observed between the X-ray count rate and the scale height of the corona, and a time lag of 10~40ks is observed between the rise in luminosity and the increase in reverberation lag. This correlation and lag are consistent with viscous propagation through the inner accretion disc, leading first to an increase in the flux of seed photons that are Comptonised by the corona, before mass accretion rate fluctuations reach the inner disc and are able to modulate the structure of the corona.

The standard alpha-disc model predicts an anticorrelation between the density of the inner accretion disc and the black hole mass times square of the accretion rate, as seen in higher mass (M-BH > 10(6) M-circle dot) active galactic nuclei (AGNs). In this work, we test the predictions of the alpha-disc model and study the properties of the inner accretion flow for the low-mass end (M-BH approximate to 10(5-6) M-circle dot) of AGNs. We utilize a new high-density disc reflection model where the density parameter varies from n(e) = 10(15) to 10(20) cm(-3) and apply it to the broad-band X-ray (0.3-10 keV) spectra of the low-mass AGN sample. The sources span a wide range of Eddington fractions and are consistent with being sub-Eddington or near-Eddington. The X-ray spectra reveal a soft X-ray excess below similar to 1.5 keV which is well modelled by high-density reflection from an ionized accretion disc of density n(e) similar to 10(18) cm(-3) on average. The results suggest a radiation pressure-dominated disc with an average of 70 per cent fraction of the disc power transferred to the corona, consistent with that observed in higher mass AGNs. We show that the disc density higher than 10(15) cm(-3) can result from the radiation pressure compression when the disc surface does not hold a strong magnetic pressure gradient. We find tentative evidence for a drop in black hole spin at low-mass regimes.

ABSTRACT The narrow-line Seyfert 1 galaxy $\rm I\,\,Zwicky\, 1$ shows a unique and complex system of ionized gas in outflow, which consists of an ultra-fast wind and a two-component warm absorber. In the last two decades, XMM–Newton monitored the source multiple times enabling the study of the long-term variability of the various outflows. Plasma in photoionization equilibrium with the ionizing source responds and varies accordingly to any change of the ionizing luminosity. However, detailed modelling of the past Reflection Grating Spectrometer (RGS) data has shown no correlation between the plasma ionization state and the ionizing continuum, revealing a complex long-term variability of the multiphase warm absorber. Here, we present a new observation of $\rm I\,\,Zwicky\, 1$ by XMM–Newton taken in early 2020 characterized by a lower X-ray flux state. The soft X-ray spectrum from the RGS reveals the two components of the warm absorber with log ξ ∼ −1.0 and log ξ ∼ 1.7. Comparing our results with the previous observations, the ionization state of the two absorbing gas components is continuously changing, following the same unpredictable behaviour. The new results strengthen the scenario in which the ionization state of the warm absorber is driven by the density of the gas rather than the ionizing luminosity. In particular, the presence of a radiation driven, inhomogeneous clumpy outflow may explain both the variability in ionization throughout the years and the line-locked N v system observed in the ultraviolet band. Finally, the EPIC-pn spectrum reveals an ultra-fast wind with an outflow velocity of ∼0.26c and ionization parameter of log ξ ∼ 3.8.

The Wide Field Imager (WFI) flying on Athena will usher in the next era of studying the hot and energetic Universe. WFI observations of faint, diffuse sources will be limited by uncertainty in the background produced by high-energy particles. These particles produce easily identified "cosmic-ray tracks" along with signals from secondary photons and electrons generated by particle interactions with the instrument. The signal from these secondaries is identical to the X-rays focused by the optics, and cannot be filtered without also eliminating these precious photons. As part of a larger effort to understand the WFI background, we here present results from a study of background-reduction techniques that exploit the spatial correlation between cosmic-ray particle tracks and secondary events. We use Geant4 simulations to generate a realistic particle background, sort this into simulated WFI frames, and process those frames in a similar way to the expected flight and ground software to produce a WFI observation containing only particle background. The technique under study, Self Anti-Coincidence or SAC, then selectively filters regions of the detector around particle tracks, turning the WFI into its own anti-coincidence detector. We show that SAC is effective at improving the systematic uncertainty for observations of faint, diffuse sources, but at the cost of statistical uncertainty due to a reduction in signal. If sufficient pixel pulse-height information is telemetered to the ground for each frame, then this technique can be applied selectively based on the science goals, providing flexibility without affecting the data quality for other science. The results presented here are relevant for any future silicon-based pixelated X-ray imaging detector, and could allow the WFI and similar instruments to probe to truly faint X-ray surface brightness.

ABSTRACTWe report on X-ray flares that were observed from the active galactic nucleus I Zwicky 1 (I Zw 1) in 2020 January by the NuSTAR and XMM-Newton observatories. The X-ray spectrum is well-described by a model comprised of the continuum emission from the corona and its reflection from the accretion disc around a rapidly spinning (a > 0.94) black hole. In order to model the broadband spectrum, it is necessary to account for the variation in ionization across the disc. Analysis of the X-ray spectrum in time periods before, during and after the flares reveals the underlying changes to the corona associated with the flaring. During the flares, the reflection fraction drops significantly, consistent with the acceleration of the corona away from the accretion disc. We find the first evidence that during the X-ray flares, the temperature drops from $140_{-20}^{+100}$ keV before to $45_{-9}^{+40}$ keV during the flares. The profile of the iron K line reveals the emissivity profile of the accretion disc, showing it to be illuminated by a compact corona extending no more than $7_{-2}^{+4}$ rg over the disc before the flares, but with tentative evidence that the corona expands as it is accelerated during the flares. Once the flares subsided, the corona had collapsed to a radius of $6_{-2}^{+2}$ rg. The rapid time-scale of the flares suggests that they arise within the black hole magnetosphere rather than in the accretion disc, and the variation of the corona is consistent with the continuum arising from the Comptonization of seed photons from the disc.

Actively accreting supermassive black holes significantly impact the evolution of their host galaxies, truncating further star formation by expelling large fractions of gas with wide-angle outflows. The X-ray band is key to understanding how these black hole winds affect their environment, as the outflows have high temperatures ($\sim$10$^{5-8}$K). We have developed a Bayesian framework for characterizing Active Galactic Nuclei (AGN) outflows with an improved ability to explore parameter space and perform robust model selection. We applied this framework to a new 700 ks and an archival 315 ks Chandra High Energy Transmission Gratings observation of the Seyfert galaxy NGC 4051. We have detected six absorbers intrinsic to NGC 4051. These wind components span velocities from 400 km s$^{-1}$ to 30,000 km s$^{-1}$. We have determined that the most statistically significant wind component is purely collisionally ionized, which is the first detection of such an absorber. This wind has $T\approx10^7$ K and $v\approx880$ km s$^{-1}$ and remains remarkably stable between the two epochs. Other slow components also remain stable across time. Fast outflow components change their properties between 2008 and 2016, suggesting either physical changes or clouds moving in and out of the line of sight. For one of the fast components we obtain one of the tightest wind density measurements to date, log $n/$[cm$^{-3}$]=13.0$^{+0.01}_{-0.02}$, and determine that it is located at $\sim$240 gravitational radii. The estimated total outflow power surpasses 5% of the bolometric luminosity (albeit with large uncertainties) making it important in the context of galaxy-black hole interactions.