Constraints on turbulent pressure in the X-ray halos of giant elliptical galaxies from resonant scattering

The dense cores of X-ray emitting gaseous halos of large elliptical galaxies with temperatures below about 0.8 keV show two prominent Fe XVII emission features, which provide a sensitive diagnostic tool to measure the effects of resonant scattering. We present here high-resolution spectra of five bright nearby elliptical galaxies, obtained with the Reflection Grating Spectrometers (RGS) on the XMM-Newton satellite. The spectra for the cores of four of the galaxies show the Fe XVII line at 15.01 Angstrom being suppressed by resonant scattering. The data for NGC 4636 in particular allow the effects of resonant scattering to be studied in detail and to prove that the 15.01 Angstrom line is suppressed only in the dense core and not in the surrounding regions. Using deprojected density and temperature profiles for this galaxy obtained with the Chandra satellite, we model the radial intensity profiles of the strongest resonance lines, accounting for the effects of resonant scattering, for different values of the characteristic turbulent velocity. Comparing the model to the data, we find that the isotropic turbulent velocities on spatial scales smaller than about 1 kpc are less than 100 km/s and the turbulent pressure support in the galaxy core is smaller than 5% of the thermal pressure at the 90% confidence level, and less than 20% at 95% confidence. Neglecting the effects of resonant scattering in spectral fitting of the inner 2 kpc core of NGC 4636 will lead to underestimates of the chemical abundances of Fe and O by ~10-20%.

💡 Research Summary

The paper investigates turbulent motions in the hot X‑ray emitting halos of massive elliptical galaxies by exploiting resonant scattering of Fe XVII lines. The authors selected five nearby, X‑ray bright ellipticals (NGC 4636, NGC 5813, NGC 1404, NGC 4649, NGC 4472) whose core temperatures lie between 0.5 and 0.8 keV, a regime where a substantial fraction of iron is in the Fe XVII ionization state. High‑resolution spectra were obtained with the Reflection Grating Spectrometers (RGS) on XMM‑Newton.

Resonant scattering becomes important when the optical depth τ of a strong resonance transition exceeds unity. The Fe XVII 15.01 Å line (2p–3d) has a large oscillator strength (f ≈ 2.73) and therefore a high τ, while the blended 17.05/17.10 Å lines (2p–3s) have f ≈ 0.12 and are essentially optically thin. Consequently, the ratio (I₁₇₀₅ + I₁₇₁₀)/I₁₅₀₁ is expected to be nearly constant for an optically thin plasma, with only a weak dependence on temperature (≈0.6–0.8 keV). Any spatial variation of this ratio therefore directly signals resonant scattering.

The RGS data were processed with SAS v8.0.0, soft‑proton flares were filtered, and background was modeled using the standard extended‑source background routine because the galaxies fill the entire field of view. Since the RGS is slitless, source extent broadens each line by Δλ = 0.138 Δθ Å; the authors accounted for this by convolving the instrument response with the surface‑brightness profile derived from EPIC/MOS images along the dispersion direction. A scaling factor s was introduced for each line to allow for differences between the line emissivity distribution and the broadband X‑ray surface brightness. Spectral fitting was performed in the 10–28 Å band with the SPEX package, assuming a single‑temperature collisional ionisation equilibrium plasma absorbed by the Galactic column (fixed to LAB values). Free parameters included temperature, normalisation, and the abundances of N, O, Ne, and Fe; C‑statistics were used.

The 13.8–15.5 Å region, which contains the strongest Fe XVII and Fe XVIII resonance lines, was initially excluded from the fit because these lines may be suppressed. After fitting the continuum and non‑resonant lines, the authors measured the Fe XVII line fluxes directly. Four of the five galaxies (NGC 4636, NGC 5813, NGC 1404, NGC 4472) show a clear suppression of the 15.01 Å line in their central 0.5′ extraction (corresponding to 1.3–2.2 kpc). NGC 4649 does not show a statistically significant suppression, likely due to its slightly higher core temperature and lower Fe XVII fraction.

NGC 4636, the brightest and best‑studied object, was examined in greater detail. Spectra were extracted not only from the central 0.5′ region but also from two adjacent 2.25′ strips, allowing a radial profile of the line ratio to be constructed. The ratio is low (~0.6) in the core and rises to the optically‑thin value (~1.0) at ≳2 kpc, confirming that resonant scattering is confined to the dense central region.

To interpret these measurements, the authors used deep Chandra ACIS observations of NGC 4636 to derive deprojected electron density and temperature profiles. These profiles were fed into a Monte‑Carlo radiative transfer model that computes the emergent intensity of the 15.01 Å line for different assumed turbulent velocities σₜ (0–200 km s⁻¹). Turbulent motions broaden the line, reducing the effective optical depth and thus weakening the scattering effect. By comparing the modeled radial intensity profiles with the observed RGS profiles, they constrained σₜ < 100 km s⁻¹ (on spatial scales ≲1 kpc) at the 90 % confidence level. This translates into a turbulent pressure fraction Pₜ/Pₜₕ < 5 % in the core; at the 95 % confidence level the limit relaxes to Pₜ/Pₜₕ < 20 %.

The authors also quantified the bias introduced by neglecting resonant scattering in standard spectral analyses. Fitting the core spectrum with an optically‑thin model underestimates the Fe and O abundances by roughly 10–20 %, because the suppressed 15.01 Å line is interpreted as a lower elemental abundance rather than scattering.

Overall, the study demonstrates that the hot gas in the cores of massive ellipticals is largely quiescent, with turbulent motions contributing only a small fraction of the total pressure support. The method of using Fe XVII resonant scattering provides a powerful, independent probe of sub‑kiloparsec turbulence, complementing other techniques such as line broadening measurements with future calorimeter missions (e.g., XRISM, Athena). The results have implications for our understanding of heat transport, mixing, and the stability of cooling flows in galaxy‑scale hot atmospheres.