Simulations
Cosmological hydrodynamical simulations, from large volumes to hundreds of massive clusters, including the new iCluster suite I ran with ~20 million CPU hours at LRZ.
- IllustrisTNG
- TNG-Cluster
- EAGLE
- SIMBA
- iCluster
Galaxies hold only part of their cosmic share of ordinary matter. I follow the rest: how black holes and stars push gas out, how environment strips it away, and how far it travels, from the hearts of galaxies to the cosmic web, and from the first two billion years to today.
At a glance
My research follows the gas from galaxy centers to the cosmic web, and from the first two billion years to today. The colors match the four chapters below.
Galaxies
star formation and quenching
Halo atmospheres
hot and cool gas in motion
Virial radius
satellites fall in and are stripped
Closure radius
where the baryon budget closes
Cosmic web
the IGM, and conformity over megaparsecs
The first two billion years: JWST’s early massive and quenched galaxies · Chapter 4
From z ~ 10 to today: galaxy scaling relations across cosmic time · Chapter 4
13 billion years of the baryon budget, in halos and the IGM · Chapter 1
Today: the closure radius, gas motions in clusters, galactic conformity · Chapters 1 to 3
Feedback and the open baryon cycle
Galaxies and their halos hold fewer baryons than their cosmic share. Energy from supermassive black holes and supernovae heats the gas and pushes it out, well beyond the traditional halo boundary. I measure how far it goes, in simulations and in X-ray observations, because the gas distribution today records the cumulative history of feedback.
I defined the closure radius: the distance from a halo’s center within which the enclosed baryon fraction reaches the cosmic value. Feedback pushes baryons well beyond the virial radius, out to this scale, from dwarfs to massive clusters (Ayromlou et al. 2023). In the TNG, EAGLE and SIMBA simulations, less than 30% of cosmic baryons sit inside galaxy-hosting halos today, against about 45% of the dark matter: over 70% of all baryons are in the intergalactic medium. The reach of this redistribution matches the scale on which the clustering of matter returns to its gravity-only form (Ayromlou et al. 2026).
In cluster cores, gas inflows and outflows are largely balanced; in the outskirts, the velocity distribution turns asymmetric, with a second peak from cosmic accretion (Ayromlou et al. 2024).
Every massive cluster today holds some cool gas, but most of it sits in the outskirts, carried in by infalling satellites (Staffehl et al. 2025).
Related work with collaborators: Introducing TNG-Cluster (Nelson et al. 2024); Hot atmospheres of cluster satellites (Rohr et al. 2024); The cooler past of the intracluster medium (Rohr et al. 2025); Core kinematics of Perseus-like clusters (Truong et al. 2024); Extreme AGN feedback in a fossil group (Eckert et al. 2025); X-ray clusters and the galaxy population (Comparat et al. 2025); ERGO-ML, X-ray clusters with contrastive learning (Chadayammuri et al. 2026).
Environment, stripping and conformity
Galaxies falling into groups and clusters lose their gas and stop forming stars. I develop ways to measure the local environment of every galaxy and to model ram-pressure stripping in semi-analytical models, to find out when, where and how fast environment quenches galaxies, and how far its influence reaches.
The gas around satellites shows no abrupt change at the virial radius, so stripping should begin before galaxies cross the halo boundary (Ayromlou et al. 2019).
Modeling ram-pressure stripping for all galaxies improves how well L-Galaxies reproduces observed quenched fractions and star formation rates, across environment, stellar mass and redshift (Ayromlou et al. 2021).
Galactic conformity (neighbors sharing the star-forming state of their central galaxy) reaches at least 5 Mpc in models and in SDSS and DESI data. Much, but not all, of it comes from galaxies near massive systems (Ayromlou et al. 2023).
Related work with collaborators: When and where ram pressure strips jellyfish galaxies (Rohr et al. 2023); Black-hole activity in stripped satellites (Kurinchi-Vendhan et al. 2025); The Hα luminosity function of the Virgo cluster (Boselli et al. 2023).
Models across techniques
Hydrodynamical simulations and semi-analytical models are the two workhorses of galaxy-formation theory, but they are rarely compared galaxy by galaxy. I ran L-Galaxies on the IllustrisTNG dark-matter merger trees to match the same galaxies in both, and to see where they agree and where they diverge. With the DECODE team, I also build semi-empirical models that reconstruct galaxy histories directly from data, and I run new cluster simulations (iCluster).
For the same galaxies, stellar masses agree to better than about 0.2 dex, but the transition from star-forming to quenched galaxies happens at a stellar mass about 0.5 dex lower in IllustrisTNG than in L-Galaxies (Ayromlou et al. 2021).
Semi-empirical models expose (in)consistencies among the stellar mass function, star-formation histories, satellite abundances and intracluster light (Fu et al. 2024).
Related work with collaborators: DECODE, a discrete semi-empirical model (Fu et al. 2022); Star formation and halo accretion in DECODE (Fu et al. 2025); Excursion set approach (Nikakhtar et al. 2018).
The early Universe
JWST has revealed galaxies that are massive, UV-bright, compact and, in some cases, already quenched within the first two billion years, and most galaxy-formation models struggle to make them. With my students and collaborators, I test which physics can: star formation and feedback in dense gas, mergers, and environmental quenching at the highest redshifts.
Spoiler: not yet, but star formation and feedback regulated by the local gas density close much of the gap.
Three versions of L-Galaxies broadly match observed scaling relations up to z ≈ 10, but they underpredict massive quenched galaxies, by up to a factor of 60 at z ≈ 3, and make them several times too large (Vani et al. 2025).
Making star formation efficient and stellar feedback ineffective in dense gas, with dissipative gas-rich mergers, makes massive UV-bright galaxies at z ≳ 9 and massive quenched galaxies at z ≈ 8 to 3 up to a hundred times more abundant than in the previous model, without changing the AGN physics or the initial mass function. Its compact galaxies follow the UV size-luminosity relation JWST observes at z ≈ 13 to 10 (Vani et al. 2026).
Among more than half a million simulated galaxies at z = 5, quenched galaxies are overwhelmingly satellites, although satellites are fewer than 10% of all galaxies. The low-mass quenched galaxies JWST finds may be quenched by their environment near massive halos, and often only temporarily (Döven et al. 2026).
A semi-empirical model built directly on observed UV luminosity functions predicts star formation efficiencies of 0.8 to 0.9 at z > 9, falling to 0.2 to 0.3 at later times; dust raises them further, while top-heavy initial mass functions lower them (Fu et al. 2026).
Looking back in time
z ~ 13 to 10
The first few hundred million years: JWST's bright, compact galaxies.
z ~ 8 to 3
Massive galaxies that have already stopped forming stars, and satellites quenched by their environment.
z = 0
Today: where we test the same models, from scaling relations to the gas around galaxies.
How I work
No single technique settles these questions, so I combine them: simulations to establish cause and effect, models to reach whole galaxy populations, and observations to test both.
Cosmological hydrodynamical simulations, from large volumes to hundreds of massive clusters, including the new iCluster suite I ran with ~20 million CPU hours at LRZ.
Galaxy-formation models, and measures and codes I developed to make theory directly comparable with data.
Surveys and X-ray data to test the models, and mock observations to design new tests.
Research with students. Students I have (co-)supervised led four of the papers on this page.
For studentsEvery paper, data release and public catalog is on the Publications page. If you would like to use the code, simulations or data, get in touch.