Search

The authors wish to thank J.M. Torrejón for the organization of the Athena-multi-messenger Workshops, held on November 27 - 29, 2018, in Alicante, Spain and on 5 May 2020 - 6 May 2020 in videoconference., In this paper we explore the scientific synergies between Athena and some of the key multi-messenger facilities that should be operative concurrently with Athena. These facilities include LIGO A+, Advanced Virgo+ and future detectors for ground-based observation of gravitational waves (GW), LISA for space-based observations of GW, IceCube and KM3NeT for neutrino observations, and CTA for very high energy observations. These science themes encompass pressing issues in astrophysics, cosmology and fundamental physics such as: the central engine and jet physics in compact binary mergers, accretion processes and jet physics in Super-Massive Binary Black Holes (SMBBHs) and in compact stellar binaries, the equation of state of neutron stars, cosmic accelerators and the origin of Cosmic Rays (CRs), the origin of intermediate and high-Z elements in the Universe, the Cosmic distance scale and tests of General Relativity and the Standard Model. Observational strategies for implementing the identified science topics are also discussed. A significant part of the sources targeted by multi-messenger facilities is of transient nature. We have thus also discussed the synergy of Athena with wide-field high-energy facilities, taking THESEUS as a case study for transient discovery. This discussion covers all the Athena science goals that rely on follow-up observations of high-energy transients identified by external observatories, and includes also topics that are not based on multi-messenger observations, such as the search for missing baryons or the observation of early star populations and metal enrichment at the cosmic dawn with Gamma-Ray Bursts (GRBs)., Open Access funding provided thanks to the CRUE-CSIC agreement with Springer Nature., This research has been supported by the European Union’s Horizon 2020 programme under the AHEAD (grant agreement number 654215) and AHEAD2020 project (grant agreement n. 871158). FJC and SMN acknowledges funding under project RTI2018-096686-B-C21 funded by MCIN/AEI/10.13039/501100011033 and by “ERDF A way of making Europe”, and by the Unidad de Excelencia María de Maeztu, ref. MDM-2017-0765.

Previously, we modelled the X-ray spectra of the narrow-line Seyfert 1 galaxy IRAS 13224$-$3809 using a disc reflection model with a fixed electron density of $10^{15}$ cm$^{-3}$. An additional blackbody component was required to fit the soft X-ray excess below 2 keV. In this work, we analyse simultaneously five flux-resolved XMM-Newton spectra of this source comprising data collected over 2 Ms. A disc reflection model with an electron density of $n_{\rm e}\approx10^{20}$ cm$^{-3}$ and an iron abundance of $Z_{\rm Fe}=3.2\pm0.5Z_{\odot}$ is used to fit the broad-band spectra of this source. No additional component is required to fit the soft excess. Our best-fit model provides consistent measurements of black hole spin and disc inclination angle as in previous models where a low disc density was assumed. In the end, we calculate the average illumination distance between the corona and the reflection region in the disc of IRAS 13224$-$3809 based on best-fit density and ionisation parameters, which changes from 0.43$\sqrt{f_{\rm AD}/f_{\rm INF}}$ $r_{\rm g}$ in the lowest flux state to 1.71$\sqrt{f_{\rm AD}/f_{\rm INF}}$ $r_{\rm g}$ in the highest flux state assuming a black hole mass of $2\times10^{6}M_{\odot}$. $f_{\rm AD}/f_{\rm INF}$ is the ratio between the flux of the coronal emission that reaches the accretion disc and infinity. This ratio depends on the geometry of the coronal region in IRAS 13224$-$3809. So we only discuss its value based on the simple `lamp-post' model, although detailed modelling of the disc emissivity profile of IRAS 13224$-$3809 is required in future to reveal the exact geometry of the corona., Comment: 17 pages, 15 figures including appendices; accepted for publication in MNRAS

The standard alpha-disc model predicts an anti-correlation 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_{\rm BH}>10^{6} M_{\odot}$) 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_{\rm BH}\approx 10^{5-6}M_{\odot}$) of AGNs. We utilize a new high-density disc reflection model where the density parameter varies from $n_{\rm e}=10^{15}$ to $10^{20}$ cm$^{-3}$ and apply it to the broadband 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 $\sim 1.5$ keV which is well modeled by high-density reflection from an ionized accretion disc of density $n_{\rm e}\sim 10^{18}$ cm$^{-3}$ on average. The results suggest a radiation pressure-dominated disc with an average of 70% 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., 20 pages, 10 figures, 6 tables. Accepted for publication in MNRAS

Ultra-fast outflows (UFOs) have been detected in the high-quality X-ray spectra of a number of active galactic nuclei (AGN) with fairly high accretion rates and are thought to significantly contribute to the AGN feedback. After a decade of dedicated study, their launching mechanisms and structure are still not well understood, but variability techniques may provide useful constraints. In this work, therefore, we perform a flux-resolved X-ray spectroscopy on a highly accreting and variable NLS1 AGN, 1H 0707-495, using all archival XMM-Newton observations to study the structure of the UFO. We find that the wind spectral lines weaken at higher luminosities, most likely due to an increasing ionization parameter as previously found in a few similar sources. Instead, the velocity is anticorrelated with the luminosity, which is opposite to the trend observed in the NLS1 IRAS 13224-3809. Furthermore, the detection of the emission lines, which are not observed in IRAS 13224-3809, indicates a wind with a larger opening angle in 1H 0707-495, presumably due to a higher accretion rate. The emitting gas is found to remain broadly constant with the luminosity. We describe the variability of the wind with a scenario where the strong radiation extends the launch radius outwards and shields the outer emitting gas, similarly to super-Eddington compact objects, although other possible explanations are discussed. Our work provides several hints for a multi-phase outflow in 1H 0707-495., Comment: 20 pages, 9 figures, accepted for publication in MNRAS

We present a detailed analysis of the velocity structure of the Virgo cluster using {\it XMM-Newton} observations. Using a novel technique which uses uses the Cu K$\alpha$ instrumental line to calibrate the EPIC-pn energy scale, we are able to obtain velocity measurements with uncertainties down to $\Delta v \sim 100$ km/s. We created 2D projected maps for the velocity, temperature, metallicity, density, pressure and entropy with an spatial resolution of 0.25$'$. We have found that in the innermost gas there is a high velocity structure, most likely indicating the presence of an outflow from the AGN while our analysis of the cluster cool core using RGS data indicates that the velocity of the gas agrees with the M87 optical redshift. An overall gradient in the velocity is seen, with larger values as we move away from the cluster core. The hot gas located within the western radio flow is redshifted, moving with a velocity $\sim 331$ km/s while the hot gas located within the eastern radio flow is blueshifted, with a velocity $\sim 258$ km/s, suggesting the presence of backflows. Our results reveal the effects of both AGN outflows and gas sloshing, in the complex velocity field of the Virgo cluster., Comment: 19 pages, 21 figures; This paper is part of a series on the ICM velocity structure using XMM-Newton observations. Related series papers: arXiv:2203.12635, arXiv:2303.17556

Soft Ultra-Luminous X-ray (ULXs) sources are a subclass of the ULXs that can switch from a supersoft spectral state, where most of the luminosity is emitted below 1 keV, to a soft spectral state with significant emission above 1 keV. In a few systems, dips have been observed. The mechanism behind this state transition and the dips nature are still debated. To investigate these issues, we obtained a long XMM-Newton monitoring campaign of a member of this class, NGC 247 ULX-1. We computed the hardness-intensity diagram for the whole dataset and identified two different branches: the normal branch and the dipping branch, which we study with four and three hardness-intensity resolved spectra, respectively. All seven spectra are well described by two thermal components: a colder ($kT_{\rm bb}$ $\sim$ 0.1-0.2 keV) black-body, interpreted as emission from the photo-sphere of a radiatively-driven wind, and a hotter ($kT_{\rm disk}$ $\sim$ 0.6 keV) multicolour disk black-body, likely due to reprocessing of radiation emitted from the innermost regions. In addition, a complex pattern of emission and absorption lines has been taken into account based on previous high-resolution spectroscopic results. We studied the evolution of spectral parameters and the flux of the two thermal components along the two branches and discuss two scenarios possibly connecting the state transition and the dipping phenomenon. One is based on geometrical occultation of the emitting regions, the other invokes the onset of a propeller effect., 13 pages, accepted by MNRAS on 2021/08/16

SS433 is the only Galactic binary system known to persistently accrete at highly super-critical (or hyper-critical) rates, similar to those in tidal disruption events, and likely needed to explain the rapid growth of those very high redshift quasars containing massive SMBHs. Probing the inner regions of SS433 in the X-rays is crucial to understanding this system, and super-critical accretion in general, but is highly challenging due to obscuration by the surrounding wind, driven from the accretion flow. NuSTAR observed SS433 in the hard X-ray band across multiple phases of its 162 day super-orbital precession period. Spectral-timing tools allow us to infer that the hard X-ray emission from the inner regions is likely being scattered towards us by the walls of the wind-cone. By comparing to numerical models, we determine an intrinsic X-ray luminosity of $\ge$ 2$\times$10$^{37}$ erg/s and that, if viewed face on, we would infer an apparent luminosity of $>$ 1$\times$10$^{39}$ erg/s, confirming SS433's long-suspected nature as an ultraluminous X-ray source (ULX). We present the discovery of a narrow, $\sim 100$s lag due to atomic processes occurring in outflowing material travelling at least 0.14-0.29c, which matches absorption lines seen in ULXs and -- in the future -- will allow us to map a super-critical outflow for the first time., 14 pages, 10 figures, accepted for publication in MNRAS

The intracluster medium (ICM) contains the vast majority of the baryonic matter in galaxy clusters and is heated to X-ray radiating temperatures. X-ray spectroscopy is therefore a key to understanding both the morphology and the dynamics of galaxy clusters. Here we recall crucial evolutionary problems of galaxy clusters unveiled by 19 years of high-resolution X-ray spectroscopy with the reflection grating spectrometer (RGS) on board XMM-Newton. Its exquisite combination of effective area, spectral resolution, and excellent performance over two decades has enabled transformational science and important discoveries such as the lack of strong cooling flows, the constraints on ICM turbulence, and cooling-heating balance. The ability of RGS to resolve individual ICM spectral lines reveals in great detail the chemical enrichment in clusters by supernovae and asymptotic giant branch stars. RGS spectra clearly show that the ICM plasma is overall in thermal equilibrium, which is unexpected given the wealth of energetic phenomena such as jets from supermassive black holes and mergers.

We present a detailed analysis of the spectral properties of the black hole transient GRS 1716-249, based on the archival Swift and NuSTAR observations taken during the outburst of this source in 2016-2017. The first six NuSTAR observations show that the source is in a canonical hard state, where the spectrum is dominated by a power-law continuum. The seventh NuSTAR observation is taken during the intermediate state where both a disc thermal component and a power-law continuum are shown. All of our observations show a broad emission line feature in the iron band and a Compton hump above 10 keV. We model the broad band spectra using a high density disc reflection model, where the soft X-ray emission in the hard state is interpreted as part of the disc reflection component. This model enables us to constrain the disc density parameter of GRS 1716-249 in the range of $10^{19}$-$10^{20}$ cm$^{-3}$. We only obtain an upper limit of the inner disc radius using high density disc reflection spectroscopy and the results indicate either a non-truncated disc or a slightly truncated disc with $R_{\rm in}, 10 pages, 7 figures, accepted by MNRAS

We present new X-ray observations of luminous heavily dust-reddened quasars (HRQs) selected from infrared sky surveys. HRQs appear to be a dominant population at high redshifts and the highest luminosities, and may be associated with a transitional "blowout" phase of black hole and galaxy co-evolution models. Despite this, their high-energy properties have been poorly known. We use the overall sample of $10$ objects with XMM-Newton coverage to study the high-energy properties of HRQs at $\left< L_{\rm bol} \right> = 10^{47.5}$ erg/s and $\left< z \right>= 2.5$. For the seven sources with strong X-ray detections, we perform spectral analyses. These find a median X-ray luminosity of $\left< L_{\rm 2-10\,keV} \right> = 10^{45.1}$ erg/s, comparable to the most powerful X-ray quasars known. The gas column densities are $N_{\rm H}=(1$-$8)\times 10^{22}$ cm$^{-2}$, in agreement with the amount of dust extinction observed. The dust to gas ratios are sub-Galactic, but are higher than found in local AGN. The intrinsic X-ray luminosities of HRQs are weak compared to the mid-infrared ($L_{\rm 6\mu m}$) and bolometric luminosities ($L_{\rm bol}$), in agreement with findings for other luminous quasar samples. For instance, the X-ray to bolometric corrections range from $\kappa_{\rm bol}\approx 50$-$3000$. The moderate absorption levels and accretion rates close to the Eddington limit ($\left< \lambda_{\rm Edd} \right>=1.06$) are in agreement with a quasar blowout phase. Indeed, we find that the HRQs lie in the forbidden region of the $N_{\rm H}$-$\lambda_{\rm Edd}$ plane, and therefore that radiation pressure feedback on the dusty interstellar medium may be driving a phase of blowout that has been ongoing for a few $10^{5}$ years. The wider properties, including [OIII] narrow-line region kinematics, broadly agree with this interpretation., Comment: In press. Accepted for publication in MNRAS

1ES 1927+654 is a paradigm-defying AGN and one of the most peculiar X-ray nuclear transients. In early 2018, this well-known AGN underwent a changing-look event, in which broad optical emission lines appeared and the optical flux increased. Yet, by July 2018, the X-ray flux had dropped by over two orders of magnitude, indicating a dramatic change to the inner accretion flow. With three years of observations with NICER, XMM-Newton, and NuSTAR, we present the X-ray evolution of 1ES 1927+654, which can be broken into three phases-(1) an early super-Eddington phase with rapid variability in X-ray luminosity and spectral parameters, (2) a stable super-Eddington phase at the peak X-ray luminosity, and (3) a steady decline back to the pre-outburst luminosity and spectral parameters. For the first time, we witnessed the formation of the X-ray corona, as the X-ray spectrum transitioned from thermally-dominated to primarily Comptonized. We also track the evolution of the prominent, broad 1 keV feature in the early X-ray spectra and show that this feature can be modeled with blueshifted reflection (z = -0.33) from a single-temperature blackbody irradiating spectrum using xillverTDE, a new flavor of the xillver models. Thus, we propose that the 1 keV feature could arise from reflected emission off the base of an optically thick outflow from a geometrically thick, super-Eddington inner accretion flow, connecting the inner accretion flow with outflows launched during extreme accretion events (e.g. tidal disruption events). Lastly, we compare 1ES 1927+654 to other nuclear transients and discuss applications of xillverTDE to super-Eddington accretors., Comment: 25 pages, 10 figures, accepted for publication in ApJ

We present a spectral analysis of NuSTAR and NICER observations of the luminous, persistently accreting neutron star (NS) low-mass X-ray binary Cygnus X-2. The data were divided into different branches that the source traces out on the Z-track of the X-ray color-color diagram; namely the horizontal branch, normal branch, and the vertex between the two. The X-ray continuum spectrum was modeled in two different ways that produced a comparable quality fit. The spectra showed clear evidence of a reflection component in the form of a broadened Fe K line, as well as a lower energy emission feature near 1 keV likely due to an ionized plasma located far from the innermost accretion disk. We account for the reflection spectrum with two independent models (relxillns and rdblur*rfxconv). The inferred inclination is in agreement with earlier estimates from optical observations of ellipsoidal light curve modeling (relxillns: $i=67^{\circ}\pm4^{\circ}$, rdblur*rfxconv: $i=60^{\circ}\pm10^{\circ}$). The inner disk radius remains close to the NS ($R_{\rm in}\leq1.15\ R_{\mathrm{ISCO}}$) regardless of the source position along the Z-track or how the 1 keV feature is modeled. Given the optically determined NS mass of $1.71\pm0.21\ M_{\odot}$, this corresponds to a conservative upper limit of $R_{\rm in}\leq19.5$ km for $M=1.92\ M_{\odot}$ or $R_{\rm in}\leq15.3$ km for $M=1.5\ M_{\odot}$. We compare these radius constraints to those obtained from NS gravitational wave merger events and recent NICER pulsar light curve modeling measurements., Comment: Accepted for publication in ApJ, 13 pages, 5 tables, 5 figures

Molecular emission was imaged with ALMA from numerous components near and within bright H2-emitting knots and absorbing dust globules in the Crab Nebula. These observations provide a critical test of how energetic photons and particles produced in a young supernova remnant interact with gas, cleanly differentiating between competing models. The four fields targeted show contrasting properties but within them, seventeen distinct molecular clouds are identified with CO emission; a few also show emission from HCO+, SiO and/or SO. These observations are compared with Cloudy models of these knots. It has been suggested that the Crab filaments present an exotic environment in which H2 emission comes from a mostly-neutral zone probably heated by cosmic rays produced in the supernova surrounding a cool core of molecular gas. Our model is consistent with the observed CO J=3-2 line strength. These molecular line emitting knots in the Crab present a novel phase of the ISM representative of many important astrophysical environments., Comment: ApJ in press, comments welcome

Accretion and ejection of matter in active galactic nuclei (AGN) are tightly connected phenomena and represent fundamental mechanisms regulating the growth of the central supermassive black hole and the evolution of the host galaxy. However, the exact physical processes involved are not yet fully understood. We present a high-resolution spectral analysis of a simultaneous \xmm\ and \nustar\ observation of the narrow line Seyfert 1 (NLS1) AGN 1H 1934-063, during which the X-ray flux dropped by a factor of $\sim6$ and subsequently recovered within 140 kiloseconds. By means of the time-resolved and flux-resolved X-ray spectroscopy, we discover a potentially variable warm absorber and a relatively stable ultra-fast outflow (UFO, $v_\mathrm{UFO}\sim-0.075\,c$) with a mild ionization state ($\log(\xi/\mathrm{erg\,cm\,s^{-1})}\sim1.6$). The detected emission lines (especially a strong and broad feature around 1\,keV) are of unknown origin and cannot be explained with emission from plasmas in photo- or collisional-ionization equilibrium. Such emission lines could be well described by a strongly blueshifted ($z\sim-0.3$) secondary reflection off the base of the equatorial outflows, which may reveal the link between the reprocessing of the inner accretion flow photons and the ejection. However, this scenario although being very promising is only tentative and will be tested with future observations., Comment: 15 pages, 17 figures, 1 table, accepted for publication in MNRAS

The brightest Fast Blue Optical Transients (FBOTs) are mysterious extragalactic explosions that may represent a new class of astrophysical phenomena. Their fast time to maximum brightness of less than a week and decline over several months and atypical optical spectra and evolution are difficult to explain within the context of core-collapse of massive stars which are powered by radioactive decay of Nickel-56 and evolve more slowly. AT2018cow (at redshift of 0.014) is an extreme FBOT in terms of rapid evolution and high luminosities. Here we present evidence for a high-amplitude quasi-periodic oscillation (QPO) of AT2018cow's soft X-rays with a frequency of 224 Hz (at 3.7$\sigma$ significance level or false alarm probability of 0.02%) and fractional root-mean-squared amplitude of >30%. This signal is found in the average power density spectrum taken over the entire 60-day outburst and suggests a highly persistent signal that lasts for a billion cycles. The high frequency (rapid timescale) of 224 Hz (4.4 ms) argues for a compact object in AT2018cow, which can be a neutron star or black hole with a mass less than 850 solar masses. If the QPO is the spin period of a neutron star, we can set limits on the star's magnetic field strength. Our work highlights a new way of using high time-resolution X-ray observations to study FBOTs., Comment: Published in Nature astronomy on 13th December 2021

We present an analysis of the NuSTAR and Swift spectra of the black hole candidate MAXI J1813-095 in a failed-transition outburst in 2018. The NuSTAR observations show evidence of reflected emission from the inner region of the accretion disc. By modelling the reflection component in the spectra, we find a disc inner radius of $R_{\rm in}, 9 pages, 12 figures, accepted for publication in MNRAS

We present a new 2 Ms Chandra observation of the cluster merger Abell 2146, which hosts two huge M~2 shock fronts each ~500 kpc across. For the first time, we resolve and measure the width of cluster merger shocks. The best-fit width for the bow shock is 17+/-1 kpc and for the upstream shock is 10.7+/-0.3 kpc. A narrow collisionless shock will appear broader in projection if its smooth shape is warped by local gas motions. We show that both shock widths are consistent with collisionless shocks blurred by local gas motions of 290+/-30 km/s. The upstream shock forms later on in the merger than the bow shock and is therefore expected to be significantly narrower. From the electron temperature profile behind the bow shock, we measure the timescale for the electrons and ions to come back into thermal equilibrium. We rule out rapid thermal equilibration of the electrons with the shock-heated ions at the 6 sigma level. The observed temperature profile instead favours collisional equilibration. For these cluster merger shocks, which have low sonic Mach numbers and propagate through a high $\beta$ plasma, we find no evidence for electron heating over that produced by adiabatic compression. Our findings are expected to be valid for collisionless shocks with similar parameters in other environments and support the existing picture from the solar wind and supernova remnants. The upstream shock is consistent with this result but has a more complex structure, including a ~2 keV increase in temperature ~50 kpc ahead of the shock., Comment: 18 pages, 17 figures, 1 table, accepted to MNRAS

Analyses of absorption from disk winds and atmospheres in accreting compact objects typically treat the central emitting regions in these systems as point sources relative to the absorber. This assumption breaks down if the absorbing gas is located within $few \times 1000\cdot GM/{c}^{2}$, in which case a small component of the absorber's Keplerian motion contributes to the velocity-width of absorption lines. Here, we demonstrate how this velocity-broadening effect can be used to constrain the sizes of central engines in accreting compact objects via a simple geometric relationship, and develop a method for modeling this effect. We apply this method on the Chandra/HETG spectra of three ultra-compact and short period neutron star X-ray binaries in which evidence of gravitationally redshifted absorption, owing to an inner-disk atmosphere, has recently been reported. The significance of the redshift is above $5\sigma$ for XTE J1710$-$281 (this work) and 4U 1916$-$053, and is inconsistent with various estimates of the relative radial velocity of each binary. For our most sensitive spectrum (XTE J1710$-$281), we obtain a 1$\sigma$ upper bound of 310 $\text{km}$ $\text{s}^{-1}$ on the magnitude of this geometric effect and a central engine of size ${R}_{CE} < 60 ~ GM/{c}^{2}$ (or, $< 90 ~ GM/{c}^{2}$ at the $3\sigma$ level). These initial constraints compare favorably to those obtained via microlensing in quasars and approach the sensitivity of constraints via relativistic reflection in neutron stars. This sensitivity will increase with further exposures, as well as the launch of future microcalorimeter and grating missions., Comment: Accepted to ApJ, 32 Pages (3 Appendices), 17 Figures

The radiative cooling time of the hot gas at the centres of cool cores in clusters of galaxies drops down to 10 million years and below. The observed mass cooling rate of such gas is very low, suggesting that AGN feedback is very tightly balanced or that the soft X-ray emission from cooling is somehow hidden from view. We use an intrinsic absorption model in which the cooling and coolest gas are closely interleaved to search for hidden cooling flows in the Centaurus, Perseus and A1835 clusters of galaxies. We find hidden mass cooling rates of between 10 to 500 Msunpyr as the cluster mass increases, with the absorbed emission emerging in the Far Infrared band. Good agreement is found between the hidden cooling rate and observed FIR luminosity in the Centaurus Cluster. The limits on the other two clusters allow for considerable hidden cooling. The implied total mass of cooled gas is much larger than the observed molecular masses. We discuss its fate including possible further cooling and collapse into undetected very cold clouds, low mass stars and substellar objects, Comment: 11 pages. 7 figures, 2 tables MNRAS in press

In this paper, we discuss atomic processes modifying the soft X-ray spectra from optical depth effects like photoelectric absorption and electron scattering suppressing the soft X-ray lines. We also show the enhancement in soft X-ray line intensities in a photoionized environment via continuum pumping. We quantify the suppression/enhancement by introducing a “line modification factor (f mod).” If 0 ≤ f mod ≤ 1, the line is suppressed, which could be the case in both collisionally ionized and photoionized systems. If f mod ≥ 1, the line is enhanced, which occurs in photoionized systems. Hybrid astrophysical sources are also very common, where the environment is partly photoionized and partly collisionally ionized. Such a system is V1223 Sgr, an Intermediate Polar binary. We show the application of our theory by fitting the first-order Chandra Medium Energy Grating (MEG) spectrum of V1223 Sgr with a combination of Cloudy-simulated additive cooling-flow and photoionized models. In particular, we account for the excess flux for O vii, O viii, Ne ix, Ne x, and Mg xi lines in the spectrum found in a recent study, which could not be explained with an absorbed cooling-flow model.

[Context] It is thought that ultraluminous X-ray sources (ULXs) are mainly powered by super-Eddington accreting neutron stars or black holes as shown by the recent discovery of X-ray pulsations and relativistic winds. [Aims] This work presents a follow-up study of the spectral evolution over two decades of the pulsing ULX NGC 1313 X-2 in order to understand the structure of the accretion disc. The primary objective is to determine the shape and nature of the dominant spectral components by investigating their variability with the changes in the source luminosity. [Methods[ We performed a spectral analysis over the canonical 0.3-10.0 keV energy band of all the high signal-to-noise XMM-Newton observations (96% of the available data), and we tested a number of different spectral models, which should approximate super-Eddington accretion discs. The baseline model consists of two thermal blackbody components with different temperatures plus an exponential cutoff powerlaw. [Results] The baseline model provides a good description of the X-ray spectra. In particular, the hotter and brighter (LXa., a., 6-9a ×a 1039 erg s-1) thermal component describes the emission from the super-Eddington inner disc and the cutoff powerlaw describes the contribution from the accretion column of the neutron star. Instead, the cooler component describes the emission from the outer region of the disc close to the spherisation radius and the wind. The luminosity-temperature relation for the cool component follows a negative trend, which is not consistent with La, T4, as is expected from a sub-Eddington thin disc of Shakura-Sunayev. This is not consistent with La, T2 either, as is expected for an advection-dominated disc. However, this would rather agree with a wind-dominated X-ray emitting region. Instead, the (Lx, Tdisk) relation for the hotter component is somewhere in between the first two theoretical scenarios. [Conclusions] Our findings agree with the super-Eddington scenario and provide further detail on the disc structure. The source spectral evolution is qualitatively similar to that seen in NGC 1313 X-1 and Holmberg IX X-1, indicating a common structure and evolution among archetypal ULXs.

We report on a spectroscopic analysis of the X-ray emission from IGR J17062-6143 in the aftermath of its June 2020 intermediate duration Type I X-ray burst. Using the Neutron Star Interior Composition Explorer, we started observing the source three hours after the burst was detected with MAXI/GSC, and monitored the source for the subsequent twelve days. We observed the tail end of the X-ray burst cooling phase, and find that the X-ray flux is severely depressed relative to its historic value for a three day period directly following the burst. We interpret this intensity dip as the inner accretion disk gradually restoring itself after being perturbed by the burst irradiation. Superimposed on this trend we observed a $1.5$ d interval during which the X-ray flux is sharply lower than the wider trend. This drop in flux could be isolated to the non-thermal components in the energy spectrum, suggesting that it may be caused by an evolving corona. Additionally, we detected a 3.4 keV absorption line at $6.3\sigma$ significance in a single $472$ s observation while the burst emission was still bright. We tentatively identify the line as a gravitationally redshifted absorption line from burning ashes on the stellar surface, possibly associated with ${}^{40}{\rm Ca}$ or ${}^{44}{\rm Ti}$., Comment: 14 pages, 7 figures, 3 tables. Accepted for publication in ApJ

Ultra-luminous X-ray pulsars (ULXPs) provide a unique opportunity to study super-Eddington accretion. We present the results of a monitoring campaign of ULXP NGC 7793 P13. Over our four-year monitoring campaign with Swift, XMM-Newton, and NuSTAR, we measured a continuous spin-up with $\dot P$ ~ -3.8e-11 s/s. The strength of the spin-up is independent of the observed X-ray flux, indicating that despite a drop in observed flux in 2019, accretion onto the source has continued at largely similar rates. The source entered an apparent off-state in early 2020, which might have resulted in a change in the accretion geometry as no pulsations were found in observations in July and August 2020. We used the long-term monitoring to update the orbital ephemeris and the periodicities seen in both the observed optical/UV and X-ray fluxes. We find that the optical/UV period is very stable over the years, with $P_\text{UV}$ = 63.75 (+0.17, -0.12) d. The best-fit orbital period determined from our X-ray timing results is 64.86 +/- 0.19 d, which is almost a day longer than previously implied, and the X-ray flux period is 65.21+/- 0.15 d, which is slightly shorter than previously measured. The physical origin of these different flux periods is currently unknown. We study the hardness ratio to search for indications of spectral changes. We find that the hardness ratios at high energies are very stable and not directly correlated with the observed flux. At lower energies we observe a small hardening with increased flux, which might indicate increased obscuration through outflows at higher luminosities. We find that the pulsed fraction is significantly higher at low fluxes. This seems to imply that the accretion geometry already changed before the source entered the deep off-state. We discuss possible scenarios to explain this behavior, which is likely driven by a precessing accretion disk., Comment: 10 pages, 6 figures, accepted for publication in A&A

Ultraluminous X-ray sources (ULXs) are mainly powered by accretion in neutron stars or stellar-mass black holes. Accreting at rates exceeding the Eddington limit by factors of a few up to hundreds, radiation pressure is expected to inflate the accretion disc, and drive fast winds that have in fact been observed at significant fractions of the speed of light. Given the super-Eddington luminosity, the accretion disc will be thicker than in sub-Eddington accretors such as common active galactic nuclei and X-ray binaries, leading to a different spectral energy distribution and, possibly, a different thermal status of the wind. Here we show the first attempt to calculate the photoionization balance of the winds driven by strong radiation pressure in thick discs with a focus on ULXs hosting black holes or non-magnetic neutron stars. We find that the winds are generally in thermally stable equilibrium, but long-term variations in the accretion rate and the inclination due to precession may have significant effects on the wind appearance and stability. Our model trends can explain the observed correlation between the spectral residuals around 1 keV and the ULX spectral state. We also find a possible correlation between the spectral hardness of the ULX, the wind velocity and the ionization parameter in support of the general scenario., Comment: 17 pages, 9 figures, 5 tables, accepted for publication in MNRAS

MAXI J1820+070 (optical counterpart ASASSN-18ey) is a black hole candidate discovered through its recent very bright outburst. The low extinction column and long duration at high flux allow detailed measurements of the accretion process to be made. In this work, we compare the evolution of X-ray spectral and timing properties through the initial hard state of the outburst. We show that the inner accretion disc, as measured by relativistic reflection, remains steady throughout this period of the outburst. Nevertheless, subtle spectral variability is observed, which is well explained by a change in coronal geometry. However, characteristic features of the temporal variability - low-frequency roll-over and QPO frequency - increase drastically in frequency, as the outburst proceeds. This suggests that the variability timescales are governed by coronal conditions rather than solely by the inner disc radius. We also find a strong correlation between X-ray luminosity and coronal temperature. This can be explained by electron pair production with a changing effective radius and a non-thermal electron fraction of ~20%., Re-submitted to MNRAS after reviewers comments. 14 pages, 14 figures

We present new Karl G. Jansky Very Large Array (VLA, 1.5 GHz) radio data for the giant elliptical galaxy IC 4296, supported by archival radio, X-ray (Chandra, XMM-Newton) and optical (SOAR, HST) observations. The galaxy hosts powerful radio jets piercing through the inner hot X-ray emitting atmosphere, depositing most of the energy into the ambient intra-cluster medium (ICM). Whereas the radio surface brightness of the A configuration image is consistent with a Fanaroff-Riley Class I (FR I) system, the D configuration image shows two bright, relative to the central region, large (~160 kpc diameter), well-defined lobes, previously reported by Killeen et al., at a projected distance r~>230 kpc. The XMM-Newton image reveals an X-ray cavity associated with one of the radio lobes. The total enthalpy of the radio lobes is ~7x10^59 erg and the mechanical power output of the jets is ~10^44 erg/s. The jets are mildly curved and possibly re-brightened by the relative motion of the galaxy and the ICM. The lobes display sharp edges, suggesting the presence of bow shocks, which would indicate that they are expanding supersonically. The central entropy and cooling time of the X-ray gas are unusually low and the nucleus hosts a warm H\alpha+[NII] nebula and a cold molecular CO disk. Because most of the energy of the jets is deposited far from the nucleus, the atmosphere of the galaxy continues to cool, apparently feeding the central supermassive black hole and powering the jet activity., Comment: Accepted for publication in MNRAS

We present a detailed X-ray timing analysis of the highly variable NLS1 galaxy, IRAS 13224-3809. The source was recently monitored for 1.5 Ms with XMM-Newton which, combined with 500 ks archival data, makes this the best studied NLS1 galaxy in X-rays to date. We apply standard time- and Fourier-domain in order to understand the underlying variability process. The source flux is not distributed lognormally, as would be expected for accreting sources. The first non-linear rms-flux relation for any accreting source in any waveband is found, with $\mathrm{rms} \propto \mathrm{flux}^{2/3}$. The light curves exhibit significant strong non-stationarity, in addition to that caused by the rms-flux relation, and are fractionally more variable at lower source flux. The power spectrum is estimated down to $\sim 10^{-7}$ Hz and consists of multiple peaked components: a low-frequency break at $\sim 10^{-5}$ Hz, with slope $\alpha < 1$ down to low frequencies; an additional component breaking at $\sim 10^{-3}$ Hz. Using the high-frequency break we estimate the black hole mass $M_\mathrm{BH} = [0.5-2] \times 10^{6} M_{\odot}$, and mass accretion rate in Eddington units, $\dot m_{\rm Edd} \gtrsim 1$. The non-stationarity is manifest in the PSD with the normalisation of the peaked components increasing with decreasing source flux, as well as the low-frequency peak moving to higher frequencies. We also detect a narrow coherent feature in the soft band PSD at $0.7$ mHz, modelled with a Lorentzian the feature has $Q \sim 8$ and an $\mathrm{rms} \sim 3$ %. We discuss the implication of these results for accretion of matter onto black holes., Comment: Accepted to MNRAS. 19 pages, 17 figures

We analyze 5 epochs of NICER data of the black hole X-ray binary MAXI J1820+070 during the bright hard-to-soft state transition in its 2018 outburst with both reflection spectroscopy and Fourier-resolved timing analysis. We confirm the previous discovery of reverberation lags in the hard state, and find that the frequency range where the (soft) reverberation lag dominates decreases with the reverberation lag amplitude increasing during the transition, suggesting an increasing X-ray emitting region, possibly due to an expanding corona. By jointly fitting the lag-energy spectra in a number of broad frequency ranges with the reverberation model reltrans, we find the increase in reverberation lag is best described by an increase in the X-ray coronal height. This result, along with the finding that the corona contracts in the hard state, suggests a close relationship between spatial extent of the X-ray corona and the radio jet. We find the corona expansion (as probed by reverberation) precedes a radio flare by ~5 days, which may suggest that the hard-to-soft transition is marked by the corona expanding vertically and launching a jet knot that propagates along the jet stream at relativistic velocities., 16 pages, 10 figures, accepted to be published in ApJL

KUG 1141+371 is a Seyfert 1 galaxy that shows a simultaneous flux increase in the optical and UV bands in the past decade. For instance, the latest Swift observation in 2019 shows that the UVW2 flux of the AGN in KUG 1141+371 has increased by over one order of magnitude since 2009. Meanwhile, the soft X-ray flux of KUG 1141+371 also shows a steady increase by one order of magnitude since 2007. The significant multi-wavelength luminosity change is likely due to a boost of mass accretion rate from approximately 0.6% of the Eddington limit to 3.2%, assuming a black hole mass of $10^{8}M_{\odot}$. In this work, we conduct detailed multi-epoch X-ray spectral analysis focusing on the variability of the X-ray continuum emission and the puzzling soft excess emission. In addition, our SED models also suggest a simultaneous increase of disc temperature and a decreasing inner disc radius along with the increasing accretion rate. Finally, we discuss possible connection between KUG 1141+371 and black hole transients in outburst., 18 pages, 14 figures, accepted by MNRAS

We present new deep, high-resolution, 1.5 GHz observations of the prototypical nearby Perseus galaxy cluster from the Karl G. Jansky Very Large Array. We isolate for the first time the complete tail of radio emission of the bent-jet radio galaxy NGC 1272, which had been previously mistaken to be part of the radio mini-halo. The possibility that diffuse radio galaxy emission contributes to mini-halo emission may be a general phenomenon in relaxed cool-core clusters, and should be explored. The collimated jets of NGC 1272 initially bend to the west, and then transition eastward into faint, 60 kpc-long extensions with eddy-like structures and filaments. We suggest interpretations for these structures that involve bulk motions of intracluster gas, the galaxy's orbit in the cluster including projection effects, and the passage of the galaxy through a sloshing cold front. Instabilities and turbulence created at the surface of this cold front and in the turbulent wake of the infalling host galaxy most likely play a role in the formation of the observed structures. We also discover a series of faint rings, south-east of NGC 1272, which are a type of structure that has never been seen before in galaxy clusters., 16 pages, 9 figures, Submitted to ApJ

We present deep (265 ks) Chandra X-ray observations of PSO J352.4034$-$15.3373, a quasar at z=5.831 that, with a radio-to-optical flux ratio of R>1000, is one of the radio-loudest quasars in the early universe and is the only quasar with observed extended radio jets of kpc-scale at $z \gtrsim 6$. Modeling the X-ray spectrum of the quasar with a power law, we find a best fit of $\Gamma = 1.99^{+0.29}_{-0.28}$, leading to an X-ray luminosity of $L_{2-10} = 1.26^{+0.45}_{-0.33} \times 10^{45}\ {\rm erg}\ {\rm s}^{-1}$ and an X-ray to UV brightness ratio of $\alpha_{\rm OX} = -1.36 \pm 0.11$. We identify a diffuse structure 50 kpc (${\sim}8^{\prime\prime}$) to the NW of the quasar along the jet axis that corresponds to a $3\sigma$ enhancement in the angular density of emission and can be ruled out as a background fluctuation with a probability of P=0.9985. While with few detected photons the spectral fit of the structure is uncertain, we find that it has a luminosity of $L_{2-10}\sim10^{44}\ {\rm erg}\ {\rm s}^{-1}$. These observations therefore potentially represent the most distant quasar jet yet seen in X-rays. We find no evidence for excess X-ray emission where the previously-reported radio jets are seen (which have an overall linear extent of $0.^{\prime\prime}28$), and a bright X-ray point source located along the jet axis to the SE is revealed by optical and NIR imaging to not be associated with the quasar., Comment: 16 pages, 7 Figures. Accepted for publication the Astrophysical Journal

Most Ultraluminous X-ray sources (ULXs) are thought to be powered by super-Eddington accretion onto stellar-mass compact objects. Accretors in this extreme regime are naturally expected to ionise copious amounts of plasma in their vicinity and launch powerful radiation-driven outflows from their discs. High spectral resolution X-ray observations (with RGS gratings onboard XMM-Newton) of a few ULXs with the best datasets indeed found complex line spectra and confirmed such extreme (0.1-0.3c) winds. However, a search for plasma signatures in a large ULX sample with a rigorous technique has never been performed, thereby preventing us from understanding their statistical properties such as the rate of occurrence, to constrain the outflow geometry and its duty cycle. We developed a fast method for automated line detection in X-ray spectra and applied it to the full RGS ULX archive, rigorously quantifying the statistical significance of any candidate lines. Collecting the 135 most significant features detected in 89 observations of 19 objects, we created the first catalogue of spectral lines detected in soft X-ray ULX spectra. We found that the detected emission lines are concentrated around known rest-frame elemental transitions and thus originate from low-velocity material. The absorption lines instead avoid these transitions, suggesting they were imprinted by blueshifted outflows. Such winds therefore appear common among the ULX population. Additionally, we found that spectrally hard ULXs show fewer line detections than soft ULXs, indicating some difference in their accretion geometry and orientation, possibly causing over-ionisation of plasma by the harder spectral energy distributions of harder ULXs., Accepted for publication in MNRAS. 22 pages, 6 figures and 3 tables

NGC1275 is the Brightest Cluster Galaxy (BCG) in the Perseus cluster and hosts the active galactic nucleus (AGN) that is heating the central 100\,kpc of the intracluster medium (ICM) atmosphere via a regulated feedback loop. Here we use a deep 490ks Cycle-19 Chandra High-Energy Transmission Grating (HETG) observation of NGC1275 to study the anatomy of this AGN. The X-ray continuum is adequately described by an unabsorbed power-law with photon index $\Gamma\approx 1.9$, creating strong tension with the detected column of molecular gas seen via HCN and HCO$^+$ line absorption against the parsec-scale core/jet. This tension is resolved if we permit a composite X-ray source; allowing a column of $N_H\sim 8\times 10^{22}\,{\rm cm}^{-2}$ to cover $\sim 15$% of the X-ray emitter does produce a significant improvement in the statistical quality of the spectral fit. We suggest that the dominant unabsorbed component corresponds to the accretion disk corona, and the sub-dominant X-ray component is the jet working surface and/or jet cocoon that is expanding into clumpy molecular gas. We suggest that this may be a common occurence in BCG-AGN. We conduct a search for photoionized absorbers/winds and fail to detect such a component, ruling out columns and ionization parameters often seen in many other Seyfert galaxies. We detect the 6.4keV iron-K$\alpha$ fluorescence line seen previously by XMM-Newton and Hitomi. We describe an analysis methodology which combines dispersive HETG spectra, non-dispersive microcalorimeter spectra, and sensitive XMM-Newton/EPIC spectra in order to constrain (sub)arcsec-scale extensions of the iron-K$\alpha$ emission region., Comment: 13 pages, accepted for publication in MNRAS

In order to connect galaxy clusters to their progenitor protoclusters, we must constrain the star formation histories within their member galaxies and the timescale of virial collapse. In this paper we characterize the complex star-forming properties of a $z=2.5$ protocluster in the COSMOS field using ALMA dust continuum and new VLA CO(1-0) observations of two filaments associated with the structure, sometimes referred to as the "Hyperion" protocluster. We focus in particular on the protocluster "core" which has previously been suggested as the highest redshift bona fide galaxy cluster traced by extended X-ray emission in a stacked Chandra/XMM image. We re-analyze this data and refute these claims, finding that at least 40 $\pm$ 17% of extended X-ray sources of similar luminosity and size at this redshift arise instead from Inverse Compton scattering off recently extinguished radio galaxies rather than intracluster medium. Using ancillary COSMOS data, we also constrain the SEDs of the two filaments' eight constituent galaxies from the rest-frame UV to radio. We do not find evidence for enhanced star formation efficiency in the core and conclude that the constituent galaxies are already massive (M$_{\star} \approx 10^{11} M_{\odot}$), with molecular gas reservoirs $>10^{10} M_{\odot}$ that will be depleted within 200-400 Myr. Finally, we calculate the halo mass of the nested core at $z=2.5$ and conclude that it will collapse into a cluster of 2-9 $\times 10^{14} M_{\odot}$, comparable to the size of the Coma cluster at $z=0$ and accounting for at least 50% of the total estimated halo mass of the extended "Hyperion" structure., 30 pages, 12 figures, 1 appendix. Accepted for publication in ApJ

We present a recent ALMA observation of the CO(1-0) line emission in the central galaxy of the Zw 3146 galaxy cluster ($z=0.2906$). We also present updated X-ray cavity measurements from archival Chandra observations. The $5\times 10^{10}\,M_{\odot}$ supply of molecular gas, which is confined to the central 4 kpc, is marginally resolved into three extensions that are reminiscent of the filaments observed in similar systems. No velocity structure that would be indicative of ordered motion is observed. The three molecular extensions all trail X-ray cavities, and are potentially formed from the condensation of intracluster gas lifted in the wakes of the rising bubbles. Many cycles of feedback would be require to account for the entire molecular gas reservoir. The molecular gas and continuum source are mutually offset by 2.6 kpc, with no detected line emission coincident with the continuum source. It is the molecular gas, not the continuum source, that lies at the gravitational center of the brightest cluster galaxy. As the brightest cluster galaxy contains possible tidal features, the displaced continuum source may correspond to the nucleus of a merging galaxy. We also discuss the possibility that a gravitational wave recoil following a black hole merger may account for the displacement., Accepted to ApJ. 12 pages, 9 figures

We present the analysis of X-ray observations of the black hole binary 4U~1630$-$47 using relativistic reflection spectroscopy. We use archival data from the RXTE, Swift, and NuSTAR observatories, taken during different outbursts of the source between $1998$ and $2015$. Our modeling includes two relatively new advances in modern reflection codes: high-density disks, and returning thermal disk radiation. Accretion disks around stellar-mass black holes are expected to have densities well above the standard value assumed in traditional reflection models (i.e., $n_{\rm e}\sim10^{15}~{\rm cm^{-3}}$). New high-density reflection models have important implications in the determination of disk truncation (i.e., the disk inner radius). This is because one must retain self-consistency in the irradiating flux and corresponding disk ionization state, which is a function of disk density and system geometry. We find the disk density is $n_{\rm e}\ge10^{20}~{\rm cm^{-3}}$ across all spectral states. This density, combined with our constraints on the ionization state of the material, implies an irradiating flux impinging on the disk that is consistent with the expected theoretical estimates. Returning thermal disk radiation -- the fraction of disk photons which bend back to the disk producing additional reflection components -- is expected predominantly in the soft state. We show that returning radiation models indeed provide a better fit to the soft state data, reinforcing previous results which show that in the soft state the irradiating continuum may be blackbody emission from the disk itself., Comment: 20 pages, 6 figures, 6 tables, Accepted by ApJ

Most ultraluminous X-ray sources (ULXs) are believed to be stellar mass black holes or neutron stars accreting beyond the Eddington limit. Determining the nature of the compact object and the accretion mode from broadband spectroscopy is currently a challenge, but the observed timing properties provide insight into the compact object and details of the geometry and accretion processes. Here we report a timing analysis for an 800 ks XMM-Newton campaign on the supersoft ultraluminous X-ray source, NGC 247 ULX-1. Deep and frequent dips occur in the X-ray light curve, with the amplitude increasing with increasing energy band. Power spectra and coherence analysis reveals the dipping preferentially occurs on $\sim 5$ ks and $\sim 10$ ks timescales. The dips can be caused by either the occultation of the central X-ray source by an optically thick structure, such as warping of the accretion disc, or from obscuration by a wind launched from the accretion disc, or both. This behaviour supports the idea that supersoft ULXs are viewed close to edge-on to the accretion disc., Comment: Accepted MNRAS version following minor comments. 8 pages, 6 figures

We present results for the first observed outburst from the transient X-ray binary source MAXI J0637-430. This study is based on eight observations from the Nuclear Spectroscopic Telescope Array (NuSTAR) and six observations from the Neil Gehrels Swift Observatory X-Ray Telescope (Swift/XRT) collected from 2019 November 19 to 2020 April 26 as the 3-79 keV source flux declined from 8.2e-10 to 1.4e-12 erg/cm^2/s. We see the source transition from a soft state with a strong disk-blackbody component to a hard state dominated by a power-law or thermal Comptonization component. NuSTAR provides the first reported coverage of MAXI J0637-430 above 10 keV, and these broadband spectra show that a two-component model does not provide an adequate description of the soft state spectrum. As such, we test whether blackbody emission from the plunging region could explain the excess emission. As an alternative, we test a reflection model that includes a physical Comptonization continuum. Finally, we also test a spectral component based on reflection of a blackbody illumination spectrum, which can be interpreted as a simple approximation to the reflection produced by returning disk radiation due to the bending of light by the strong gravity of the black hole. We discuss the physical implications of each scenario and demonstrate the value of constraining the source distance., Comment: 17 pages, 10 figures, accepted to The Astrophysical Journal (ApJ)

The actual mechanism(s) powering galactic outflows in active galactic nuclei (AGN) is still a matter of debate. At least two physical models have been considered in the literature: wind shocks and radiation pressure on dust. Here we provide a first quantitative comparison of the AGN radiative feedback scenario with observations of galactic outflows. We directly compare our radiation pressure-driven shell models with the observational data from the most recent compilation of molecular outflows on galactic scales. We show that the observed dynamics and energetics of galactic outflows can be reproduced by AGN radiative feedback, with the inclusion of radiation trapping and/or luminosity evolution. The predicted scalings of the outflow energetics with AGN luminosity can also quantitatively account for the observational scaling relations. Furthermore, sources with both ultra-fast and molecular outflow detections are found to be located in the `forbidden' region of the $N_\mathrm{H} - \lambda$ plane. Overall, an encouraging agreement is obtained over a wide range of AGN and host galaxy parameters. We discuss our results in the context of recent observational findings and numerical simulations. In conclusion, AGN radiative feedback is a promising mechanism for driving galactic outflows that should be considered, alongside wind feedback, in the interpretation of future observational data., Comment: accepted for publication in MNRAS

We present new XMM-Newton observations extending the mosaic of the Perseus cluster out to the virial radius to the west. Previous studies with ROSAT have reported a large excess in surface brightness to the west, possibly the result of large scale gas sloshing, but lacked the spatial resolution and depth to determine if this excess lay behind a cold front. In our new XMM observations we have found that there is a sharp edge in X-ray surface brightness near the cluster virial radius (1.7Mpc) to the west, with a width comparable to the mean free path. The temperature measurements obtained with Suzaku data across this edge show that the temperature increases sharply outside this edge. All of these properties are consistent with this edge being the largest cold front ever observed in a galaxy cluster. We also find a second edge to the west, 1.2Mpc from the core, which also appears to be a cold front. Our results indicate that magnetic fields are able to stabilize the cold fronts against turbulence all the way out to the cluster virial radius. By comparing with numerical simulations, we find that these large cold fronts require large impact parameter, low mass ratio mergers which can produce fast gas motions without destroying the cluster core., 10 pages, 7 figures, submitted to ApJ

We report on the discovery of a new, transient ultraluminous X-ray source (ULX) in the galaxy NGC 7090. This new ULX, which we refer to as NGC 7090 ULX3, was discovered via monitoring with $Swift$ during 2019-20, and to date has exhibited a peak luminosity of $L_{\rm{X}} \sim 6 \times 10^{39}$ erg s$^{-1}$. Archival searches show that, prior to its recent transition into the ULX regime, ULX3 appeared to exhibit a fairly stable luminosity of $L_{\rm{X}} \sim 10^{38}$ erg s$^{-1}$. Such strong long-timescale variability may be reminiscent of the small population of known ULX pulsars, although deep follow-up observations with $XMM$-$Newton$ and $NuSTAR$ do not reveal any robust X-ray pulsation signals. Pulsations similar to those seen from known ULX pulsars cannot be completely excluded, however, as the limit on the pulsed fraction of any signal that remains undetected in these data is $\lesssim$20\%. The broadband spectrum from these observations is well modelled with a simple thin disc model, consistent with sub-Eddington accretion, which may instead imply a moderately large black hole accretor ($M_{\rm{BH}} \sim 40 ~ M_{\odot}$). Similarly, though, more complex models consistent with the super-Eddington spectra seen in other ULXs (and the known ULX pulsars) cannot be excluded given the limited signal-to-noise of the available broadband data. The nature of the accretor powering this new ULX therefore remains uncertain., 11 pages, 7 figures; accepted for publication in MNRAS

We present the first high-resolution 230-470 MHz map of the Perseus cluster obtained with the Karl G. Jansky Very Large Array. The high dynamic range and resolution achieved has allowed the identification of previously-unknown structures in this nearby galaxy cluster. New hints of sub-structures appear in the inner radio lobes of the brightest cluster galaxy NGC 1275. The spurs of radio emission extending into the outer X-ray cavities, inflated by past nuclear outbursts, are seen for the first time at these frequencies, consistent with spectral aging. Beyond NGC 1275, we also analyze complex radio sources harbored in the cluster. Two new distinct, narrowly-collimated jets are visible in IC 310, consistent with a highly-projected narrow-angle tail radio galaxy infalling into the cluster. We show how this is in agreement with its blazar-like behavior, implying that blazars and bent-jet radio galaxies are not mutually exclusive. We report the presence of filamentary structures across the entire tail of NGC 1265, including two new pairs of long filaments in the faintest bent extension of the tail. Such filaments have been seen in other cluster radio sources such as relics and radio lobes, indicating that there may be a fundamental connection between all these radio structures. We resolve the very narrow and straight tail of CR 15 without indication of double jets, so that the interpretation of such head-tail sources is yet unclear. Finally, we note that only the brightest western parts of the mini-halo remain, near NGC 1272 and its bent double jets., Comment: 17 pages, 12 figures, Accepted for publication in MNRAS

We continue the analysis of NuSTAR data from the recent discovery outburst of MAXI J1820+070 (optical counterpart ASASSN-18ey), focussing on an observation including unusual flaring behaviour during the hard to soft state transition, which is a short phase of outbursts and so comparatively rarely observed. Two plateaus in flux are separated by a variable interval lasting ~10 ks, which shows dipping then flaring stages. The variability is strongest (with fractional variability up to $F_{\rm Var}\sim10\%$) at high energies and reduces as the contribution from disc emission becomes stronger. Flux-resolved spectra show that the variability is primarily due to the power law flux changing. We also find a long soft lag of the thermal behind the power law emission, which is $20_{-1.2}^{+1.6}$ s during the flaring phase. The lag during the dipping stage has a different lag-energy spectrum, which may be due to a wave passing outwards through the disc. Time resolved spectral fitting suggests that the lag during the flaring stage may be due to the disc re-filling after being disrupted to produce the power law flare, perhaps related to the system settling after the jet ejection which occurred around 1 day before. The timescales of these phenomena imply a low viscosity parameter, $\alpha\sim10^{-3}$, for the inner region of the disc., Comment: MNRAS, resubmitted after minor revision. 12 pages, 10 figures, 1 table

The very small accretion disks in ultra-compact X-ray binaries (UCXBs) are special laboratories in which to study disk accretion and outflows. We report on three sets of new (250 ks total) and archival (50 ks) Chandra/HETG observations of the "dipping" neutron-star X-ray binary 4U 1916$-$053, which has an orbital period of $P\simeq 50$~minutes. We find that the bulk of the absorption in all three spectra originates in a disk atmosphere that is redshifted by $v\simeq 220-290$ $\text{km}$ $\text{s}^{-1}$, corresponding to the gravitational redshift at radius of $R \sim 1200$ $GM/{c}^{2}$. This shift is present in the strongest, most highly ionized lines (Si XIV and Fe XXVI), with a significance of 5$\sigma$. Absorption lines observed during dipping events (typically associated with the outermost disk) instead display no velocity shifts and serve as a local standard of rest, suggesting that the redshift is intrinsic to an inner disk atmosphere and not due to radial motion in the galaxy or a kick. In two spectra, there is also evidence of a more strongly redshifted component that would correspond to a disk atmosphere at $R \sim 70$ $GM/{c}^{2}$; this component is significant at the 3$\sigma$ level. Finally, in one spectrum, we find evidence of disk wind with a blue shift of $v = {-1700}^{+1700}_{-1200}$ $\text{km}$ $\text{s}^{-1}$. If real, this wind would require magnetic driving., Comment: 15 Pages, 6 Figures, Accepted to ApJ Letters

The black hole candidate and X-ray binary MAXI J1535-571 was discovered in September 2017. During the decay of its discovery outburst, and before returning to quiescence, the source underwent at least four reflaring events, with peak luminosities of $\sim$10$^{35-36}$ erg s$^{-1}$ (d/4.1 kpc)$^2$. To investigate the nature of these flares, we analysed a sample of NICER observations taken with almost daily cadence. In this work we present the detailed spectral and timing analysis of the evolution of the four reflares. The higher sensitivity of NICER at lower energies, in comparison with other X-ray detectors, allowed us to constrain the disc component of the spectrum at $\sim$0.5 keV. We found that during each reflare the source appears to trace out a q-shaped track in the hardness-intensity diagram similar to those observed in black hole binaries during full outbursts. MAXI J1535-571 transits between the hard state (valleys) and softer states (peaks) during these flares. Moreover, the Comptonised component is undetected at the peak of the first reflare, while the disc component is undetected during the valleys. Assuming the most likely distance of 4.1 kpc, we find that the hard-to-soft transitions take place at the lowest luminosities ever observed in a black hole transient, while the soft-to-hard transitions occur at some of the lowest luminosities ever reported for such systems., Accepted for publication in MNRAS main journal

arXiv:2007.15660v1, Cosmological simulations, as well as mounting evidence from observations, have shown that supermassive black holes play a fundamental role in regulating the formation of stars throughout cosmic time. This has been clearly demonstrated in the case of galaxy clusters in which powerful feedback from the central black hole is preventing the hot intracluster gas from cooling catastrophically, thus reducing the expected star formation rates by orders of magnitude. These conclusions, however, have been almost entirely based on nearby clusters. Based on new Chandra X-ray observations, we present the first observational evidence for massive, runaway cooling occurring in the absence of supermassive black hole feedback in the high-redshift galaxy cluster SpARCS104922.6 + 564032.5 (z = 1.709). The hot intracluster gas appears to be fueling a massive burst of star formation (≈900 M⊙ yr−1) that is offset by dozens of kpc from the central galaxy. The burst is co-spatial with the coolest intracluster gas but not associated with any galaxy in the cluster. In less than 100 million years, such runaway cooling can form the same amount of stars as in the Milky Way. Therefore, intracluster stars are not only produced by tidal stripping and the disruption of cluster galaxies, but can also be produced by runaway cooling of hot intracluster gas at early times. Overall, these observations show the dramatic impact when supermassive black hole feedback fails to operate in clusters. They indicate that in the highest overdensities, such as clusters and protoclusters, runaway cooling may be a new and important mechanism for fueling massive bursts of star formation in the early universe., J.H.-L. acknowledges support from NSERC via the Discovery grant program, as well as the Canada Research Chair program. C.R. acknowledges financial support from the physics department of the Université de Montréal. G.W. acknowledges support from the National Science Foundation through grant AST-1517863, by HST program number GO15294, and by grant number 80NSSC17K0019 issued through the NASA Astrophysics Data Analysis Program (ADAP). Support for program number GO-15294 was provided by NASA through a grant from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Incorporated, under NASA contract NAS5-26555. M.J.J. acknowledges support for the current research from the National Research Foundation of Korea under the programs 2017R1A2B2004644 and 2020R1A4A2002885. H.R.R. acknowledges support from an STFC Ernest Rutherford Fellowship and an Anne McLaren Fellowship.

Galactic outflows are known to consist of several gas phases; however, the connection between these phases has been investigated little and only in a few objects. In this paper, we analyse Multi Unit Spectroscopic Explorer (MUSE)/Very Large Telescope (VLT) data of 26 local (U)LIRGs and study their ionized and neutral atomic phases. We also include objects from the literature to obtain a sample of 31 galaxies with spatially resolved multiphase outflow information. We find that the ionized phase of the outflows has on average an electron density three times higher than the disc (ne,disc ∼145 cm-3 versus ne,outflow ∼500 cm-3), suggesting that cloud compression in the outflow is more important than cloud dissipation. We find that the difference in extinction between outflow and disc correlates with the outflow gas mass. Together with the analysis of the outflow velocities, this suggests that at least some of the outflows are associated with the ejection of dusty clouds from the disc. This may support models where radiation pressure on dust contributes to driving galactic outflows. The presence of dust in outflows is relevant for potential formation of molecules inside them. We combine our data with millimetre data to investigate the molecular phase. We find that the molecular phase accounts for more than 60 per cent of the total mass outflow rate in most objects and this fraction is higher in active galactic nuclei (AGN)-dominated systems. The neutral atomic phase contributes of the order of 10 per cent, while the ionized phase is negligible. The ionized-to-molecular mass outflow rate declines slightly with AGN luminosity, although with a large scatter. © 2021 The Author(s) Published by Oxford University Press on behalf of Royal Astronomical Society., RM and AF acknowledge ERC Advanced Grant 695671 ‘QUENCH’ and support by the Science and Technology Facilities Council (STFC). SA acknowledges support from the Spanish Ministerio de Ciencia e Innovación through grant ESP2017-83197-P. S Cazzoli acknowledges financial support from the State Agency for Research of the Spanish MCIU through the ‘Center of Excellence Severo Ochoa’ award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709). AM, GC, and FM acknowledge the grant PRIN 2017PH3WAT_003 from the Italian Ministry for University and Research. MP is supported by the Programa Atracción de Talento de la Comunidad de Madrid via grant 2018-T2/TIC-11715. GV acknowledges support from ANID programs FONDECYT Postdoctorado 3200802 and Basal-CATA AFB-170002. EB acknowledges the support from Comunidad de Madrid through the Atracción de Talento grant 2017-T1/TIC-5213.

We analyzed Chandra X-ray observations of five galaxy clusters whose atmospheric cooling times, entropy parameters, and cooling time to free-fall time ratios within the central galaxies lie below 1 Gyr, below 30 keV cm^2, and between 20 < tcool/tff < 50, respectively. These thermodynamic properties are commonly associated with molecular clouds, bright H-alpha emission, and star formation in central galaxies. However, none of these clusters have detectable H-alpha indicated in the ACCEPT database, nor do they have significant star formation rates or detectable molecular gas. Among these, only RBS0533 has a detectable radio/X-ray bubble which are commonly observed in cooling atmospheres. Signatures of uplifted, high metallicity atmospheric gas are absent. Despite its prominent X-ray bubble, RBS0533 lacks significant levels of molecular gas. Cold gas is absent at appreciable levels in these systems perhaps because their radio sources have failed to lift low entropy atmospheric gas to an altitude where the ratio of the cooling time to the free-fall time falls below unity., Accepted for publication in ApJ. The draft includes significant changes to the Introduction intended to better motivate this work by placing the problem in the context of current theoretical understanding of thermally unstable cooling. The abstract and conclusion were edited for clarity

We present results of temporal and spectral analyses on four XMM-Newton EPIC pn observations of IRAS 17020+4544, a narrow-line Seyfert 1 galaxy with evidence of a radio jet. Analysis of the light curves reveals that this radio-loud source does not behave like the bulk population of its radio-quiet counterparts. A trend of spectral hardening with increased flux is found. Variability is found to increase with energy, though it decreases as the spectrum hardens. The first 40 ks of the most recent observation behave uniquely among the epochs, exhibiting a softer spectral state than at any other time. Possible non-stationarity at low energies is found, with no such effect present at higher energies, suggesting at least two distinct spectral components. A reverberation signature is confirmed, with the lag-frequency, lag-energy, and covariance spectra changing significantly during the soft-state epoch. The temporal analysis suggests a variable power-law in the presence of a reflection component, thus motivating such a fit for the 0.3-10 keV EPIC pn spectra from all epochs. We find an acceptable spectral fit using the timing-motivated parameters and report the detection of a broad Fe K emission line, requiring an additional model component beyond the reflection spectrum. We discuss links between this source and other narrow-line Seyfert 1 sources that show evidence of jet activity, finding similarities among this currently very limited sample of interesting objects., Comment: Accepted for publication in MNRAS. 19 pages, 14 figures

Axions/axion-like particles (ALPs) are a well motivated extension of the Standard Model and are generic within String Theory. The X-ray transparency of the intracluster medium (ICM) in galaxy clusters is a powerful probe of light ALPs (with mass $, Comment: Published in The Astrophysical Journal, 2020, vol 890, pp59. This posting corrects typos in equation (2) and the likelihood function just prior to equation (3). Typos were introduced into the draft at late stage and the analysis is correct