Data and analysis scripts for: Solving the cooling flow problem with combined jet-wind AGN feedback
Data files
May 27, 2026 version files 447.89 MB
-
He_2025.zip
447.88 MB
-
README.md
5.48 KB
Abstract
This repository contains the primary simulation data and analysis scripts associated with "Solving the cooling flow problem with combined jet-wind AGN feedback". The original data were generated using the Athena++ hydrodynamics code. Active galactic nucleus (AGN) feedback is widely viewed as the most promising solution to the long-standing cooling flow problem in galaxy clusters, yet previous models prescribe jet properties inconsistent with accretion physics. We perform an idealized hydrodynamic simulation of a galaxy cluster with no merger history and a relaxed state, with its other properties similar to the Perseus cluster, using the MACER framework, incorporating both jets and winds whose properties are constrained by general relativistic magnetohydrodynamic simulations of black hole accretion and observations. The combined feedback reproduces key observables, including cold gas mass, star formation rate, thermodynamic radial profiles, and black hole growth, while jet-only or wind-only models fail. The success arises from turbulence driven by jet-wind shear that enhances kinetic-to-thermal energy conversion, boosting heating efficiency by factors of three and six relative to wind-only and jet-only cases, respectively.
Project Title
Simulation Data and Plotting Scripts for AGN Feedback in Galaxy Clusters (He et al. 2025)
Overview
This repository contains the primary simulation data and analysis scripts associated with "Solving the cooling flow problem with combined jet-wind AGN feedback" (He et al. 2025). The original data were generated using the Athena++ hydrodynamics code.
Repository Structure
The data are organized into three main directories, representing the different physics models explored in the paper:
JetWind/: JetWind Model.JetOnly/: JetOnly Model.WindOnly/: WindOnly Model.
The root directory contains the primary visualization script:
he+2025.ipynb: Jupyter Notebook (Python 3) for reading the data and generating all manuscript figures.
Data Description
Each model directory (JetWind/, WindOnly/, JetOnly/) contains the following files:
1. Temporal Evolution Data (ASCII Text Files)
| File Name | Description | Columns / Units |
|---|---|---|
He_2025.zip |
The zip file of all the Data and analysis scripts. | None |
cal_avg.txt |
Integrated AGN and SFR properties | 1: Time (Gyr), 2: $M_{\text{dot, Bondi}}$ at inner boundary ($M_\odot$/yr), 3: Star Formation Rate ($M_\odot$/yr), 4: Hot wind power (erg/s), 5: Cold wind power (erg/s), 6: Jet power (erg/s), 7: Total AGN power (erg/s), 8: Total AGN kinetic power (erg/s) |
E_tur.txt |
Turbulence dissipation rate | 1: Time (Gyr), 2: 100 Myr averaged turbulence dissipation rate (erg/s) |
m_cold.txt |
Cold gas mass evolution | 1: Time (Gyr), 2: Instantaneous cold gas mass ($M_\odot$), 3: 15 Myr averaged cold gas mass ($M_\odot$) |
R_bondi.txt |
Bondi radius evolution | 1: Time (Myr), 2: Mass-flux weighted Bondi radius (pc) |
shock.txt |
Shock heating evolution | 1: Time (Myr), 2: Shock heating rate (erg/s) |
yuan18.hst |
Athena++ history file | 1: Time (Gyr), 20: Black Hole mass ($M_\odot$), 25: BH Bolometric luminosity ($L_{\text{BH}}$ in erg/s) |
2. Radial Profiles & 3D Grid Data
| File Name | Format | Description |
|---|---|---|
profiles.h5 |
HDF5 | Radial profile data for gas Entropy, Temperature, and Density. |
t_cool_ff.pkl |
Pickle | Radial profiles of the ratio between cooling time and free-fall time ($t_{\text{cool}}/t_{\text{ff}}$). |
vr_mean.npy |
NumPy | Time-averaged radial velocity ($v_r$) data on the grid during the hot mode phase. |
r_all.npy |
NumPy | Radial coordinate values ($r$) for the 3D grid setup. |
Requirements
To execute the plotting script he+2025.ipynbThe following Python libraries are required:
numpymatplotlibh5pyscipypandaspickle(standard library)
Usage
- Keep the
JetWind,WindOnly, andJetOnlyfolders in the same directory as the.ipynbfile. - Open
he+2025.ipynbin a Jupyter environment. - Run the cells to reproduce the figures. The script is configured to automatically handle paths and unit conversions for different models.
Contact
For further information regarding the simulation setup or the MACER3D framework, please contact the corresponding author as listed in the publication.