-
The role of prompt cusps in driving the core collapse of SIDM halos
Authors:
Vinh Tran,
Daniel Gilman,
M. Sten Delos,
Xuejian Shen,
Oliver Zier,
Mark Vogelsberger,
David Xu
Abstract:
Prompt cusps (PCs) form from the direct collapse of overdensities in the early Universe, reside at the center of every dark matter halo, and have density profiles steeper than $r^{-1}$ NFW cusps. Using a suite of high-resolution N-body simulations, we study the evolution of isolated halos in self-interacting dark matter (SIDM) with massive PCs embedded at their centers, a scenario that could be re…
▽ More
Prompt cusps (PCs) form from the direct collapse of overdensities in the early Universe, reside at the center of every dark matter halo, and have density profiles steeper than $r^{-1}$ NFW cusps. Using a suite of high-resolution N-body simulations, we study the evolution of isolated halos in self-interacting dark matter (SIDM) with massive PCs embedded at their centers, a scenario that could be realized in certain SIDM models with light mediators that predict a small-scale suppression of the linear matter power spectrum. We track the evolution of three equally concentrated $10^7\,{\rm{M}}_\odot$ halos, hosting PCs of various total masses, and quantify how the presence of a PC affects the processes of core formation and collapse. Early in the core-formation phase, halos with more prominent PCs exhibit a delayed evolution by a factor of $\sim 2$ due to smaller velocity dispersion gradients in the inner region. During most of the core-collapse phase, the halo evolution becomes closely aligned in physical time, with appropriate rescaling of densities, radii, and velocity dispersions. The scale densities and radii preserve the virial mass of the original halos, but with increased concentration. Deviations occur at the late phase of core-collapse at the level of $\sim 5\%$ relative to the reference collapse track of an NFW halo. These deviations depend non-trivially on both the increased concentration incurred by the PCs, as well as the velocity dispersion (temperature) of the outer halo regions, which can inhibit or enhance the heat transfer process. Our simulations illustrate the complex interplay between the inner and outer halo profiles in determining the onset of core collapse and motivate future studies in the full cosmological context.
△ Less
Submitted 28 November, 2025;
originally announced December 2025.
-
Metallicity Gradients in Modern Cosmological Simulations II: The Role of Bursty Versus Smooth Feedback at High-Redshift
Authors:
Alex M. Garcia,
Paul Torrey,
Aniket Bhagwat,
Xuejian Shen,
Mark Vogelsberger,
William McClymont,
Jaya Nagarajan-Swenson,
Sophia G. Ridolfo,
Peixin Zhu,
Dhruv T. Zimmerman,
Oliver Zier,
Sarah Biddle,
Arnab Sarkar,
Priyanka Chakraborty,
Ruby J. Wright,
Kathryn Grasha,
Tiago Costa,
Laura Keating,
Rahul Kannan,
Aaron Smith,
Enrico Garaldi,
Ewald Puchwein,
Benedetta Ciardi,
Lars Hernquist,
Lisa J. Kewley
Abstract:
The distribution of gas-phase metals within galaxies encodes the impact of stellar feedback on galactic evolution. At high-redshift, when galaxies are rapidly assembling, feedback-driven outflows and turbulence can strongly reshape radial metallicity gradients. In this work, we use the FIRE-2, SPICE, Thesan and Thesan Zoom cosmological simulations -- spanning a range of stellar feedback from burst…
▽ More
The distribution of gas-phase metals within galaxies encodes the impact of stellar feedback on galactic evolution. At high-redshift, when galaxies are rapidly assembling, feedback-driven outflows and turbulence can strongly reshape radial metallicity gradients. In this work, we use the FIRE-2, SPICE, Thesan and Thesan Zoom cosmological simulations -- spanning a range of stellar feedback from bursty (time-variable) to smooth (steady) -- to investigate how these feedback modes shape gas-phase metallicity gradients at $3<z\lesssim11$. Across all models, we find that galaxies with bursty feedback (FIRE-2, SPICE Bursty, and Thesan Zoom) develop systematically flatter (factors of $\sim2-10$) metallicity gradients than those with smooth feedback (SPICE Smooth and Thesan Box), particularly at stellar masses $M_\star > 10^{9}~{\rm M_\odot}$. These results demonstrate that bursty stellar feedback provides sufficient turbulence to prevent strong negative gradients from forming, while smooth stellar feedback does not generically allow for efficient radial redistribution of metals thereby keeping gradients steep. Finally, we compare with recent observations, finding that the majority -- but, notably, not all -- of the observed gradients may favor a bursty stellar feedback scenario. In all, these results highlight the utility of high-resolution observations of gas-phase metallicity at high-redshift as a key discriminator of these qualitatively different feedback types.
△ Less
Submitted 30 October, 2025;
originally announced October 2025.
-
The THESAN-ZOOM project: The Hidden Neighbours of OI Absorbers during Reionization
Authors:
Giulia Pruto,
Laura Keating,
Rahul Kannan,
Ewald Puchwein,
Aaron Smith,
Josh Borrow,
Enrico Garaldi,
Mark Vogelsberger,
Oliver Zier,
William McClymont,
Xuejian Shen,
Sandro Tacchella
Abstract:
Metal absorbers represent a powerful probe of galaxy feedback and reionization, as highlighted by both observational and theoretical results showing an increased abundance of low-ionised metal species at higher redshifts. The origin of such absorbers is currently largely unknown because of the low number of galaxy counterparts detected, suggesting that they might be surrounded by low-mass faint so…
▽ More
Metal absorbers represent a powerful probe of galaxy feedback and reionization, as highlighted by both observational and theoretical results showing an increased abundance of low-ionised metal species at higher redshifts. The origin of such absorbers is currently largely unknown because of the low number of galaxy counterparts detected, suggesting that they might be surrounded by low-mass faint sources that fall below the detection threshold of current instruments. We use the THESAN-ZOOM radiation hydrodynamic simulations to investigate the connection between properties of neutral oxygen (OI) absorbers and galaxies within the redshift range $z = 5 - 8$. We find that the circumgalactic medium of galaxies becomes progressively ionised with cosmic time, leading to a decrease of $\approx 0.2$ in the covering fraction of neutral oxygen, while the total oxygen covering fraction remains constant. The observable absorbers ($N_{\rm OI} \gtrsim 10^{13}\,\text{cm}^{-2}$) are not confined to haloes: at $z \geq 5$ the majority ($\gtrsim 60\%$) arise beyond $R_{\rm{vir}}$, and including these systems is essential to reproduce the observed increase in absorber incidence with redshift. The simulated absorbers preferentially reside in overdensities rich in low-mass galaxies ($M_\star \leq 10^8\,\rm{M}_\odot$), explaining the scarcity of detected counterparts, while not excluding the possibility of nearby star-forming sources ($\geq 5\,\text{M}_\odot\,\text{yr}^{-1}$) similar to those suggested by the latest ALMA observations and, at larger distances, by the JWST. These results establish OI absorbers as sensitive tracers of the evolving ionisation structure around faint galaxies to be probed by forthcoming deep spectroscopic surveys.
△ Less
Submitted 15 October, 2025;
originally announced October 2025.
-
Formation of protostars and the launching of stellar core outflows with moving-mesh radiation non-ideal magnetohydrodynamics
Authors:
Alexander C. Mayer,
Rüdiger Pakmor,
Thorsten Naab,
Oliver Zier,
Alexei V. Ivlev,
Tommaso Grassi,
Paola Caselli,
Volker Springel
Abstract:
We present an implementation of radiative transfer with flux-limited diffusion (FLD) for the moving-mesh code {\small AREPO} and use the method in a physical model for the formation of protostars with non-ideal radiation-magnetohydrodynamics (RMHD). We follow previous work in splitting the additional terms to the hydrodynamical equations arising from the inclusion of radiation into terms to be int…
▽ More
We present an implementation of radiative transfer with flux-limited diffusion (FLD) for the moving-mesh code {\small AREPO} and use the method in a physical model for the formation of protostars with non-ideal radiation-magnetohydrodynamics (RMHD). We follow previous work in splitting the additional terms to the hydrodynamical equations arising from the inclusion of radiation into terms to be integrated explicitly and implicitly, as the diffusion and coupling terms would impose very restrictive timestep criteria. We validate the scheme with standard test problems for radiation diffusion, matter-gas coupling, and radiative shocks from the literature. Our implementation is compatible with local timestepping, which often presents problems for implicit schemes, and we found very good agreement with results obtained with global timesteps. We present an example application of the new implementation to the collapse of a $1\,{\rm M}_\odot$ molecular cloud core to a second Larson core modelled with radiation non-ideal magnetohydrodynamics. A high-velocity jet with v$_{\rm rad}> 10\, {\rm km\,s^{-1}}$ is self-consistently launched from the second core, nested within the first core, which produces a lower-velocity magnetorotational outflow. We observe magnetic field amplification up to more than $\vert \mathbf{B}\vert_{\rm max}>10^5$ G in the second core, which is surrounded by a small (<0.5 au) disk. This application demonstrates the robustness of our scheme in multi-scale and high-resolution simulations on arbitrary meshes and, as such, the model can be readily used for further simulations of protostar formation at high resolution.
△ Less
Submitted 20 December, 2025; v1 submitted 14 October, 2025;
originally announced October 2025.
-
The Cosmic Rush Hour: Rapid Formation of Bright, Massive, Disky, Star-Forming Galaxies as Signatures of Early-Universe Physics
Authors:
Xuejian Shen,
Oliver Zier,
Mark Vogelsberger,
Michael Boylan-Kolchin,
Lars Hernquist,
Sandro Tacchella,
Rohan P. Naidu
Abstract:
Early JWST observations have revealed a high-redshift universe more vibrant than predicted by canonical galaxy-formation models within $Λ$CDM, showing an excess of ultraviolet(UV)-bright, massive, and morphologically mature galaxies. Departures from $Λ$CDM prior to recombination can imprint signatures on non-linear structure formation at high redshift. In this paper, we investigate one such scenar…
▽ More
Early JWST observations have revealed a high-redshift universe more vibrant than predicted by canonical galaxy-formation models within $Λ$CDM, showing an excess of ultraviolet(UV)-bright, massive, and morphologically mature galaxies. Departures from $Λ$CDM prior to recombination can imprint signatures on non-linear structure formation at high redshift. In this paper, we investigate one such scenario - Early Dark Energy, originally proposed to resolve the Hubble tension - and its implications for these high-redshift challenges. We present the first large-scale cosmological hydrodynamic simulations of these models. Modifications to the pre-recombination expansion history accelerate early structure formation and produce UV luminosity and stellar mass functions in excellent agreement with JWST measurements, requiring essentially no additional calibrations. Predictions converge to $Λ$CDM at lower redshifts ($z \lesssim 3$), thereby preserving all successes of $Λ$CDM. This model also accelerates the emergence of stellar and gaseous disks, increasing their number densities by $\sim 0.5$ dex at $z\simeq 6$-7, primarily due to the higher abundance of massive galaxies. Taken together, these results demonstrate how early-universe physics can simultaneously reconcile multiple high-redshift challenges and the Hubble tension while retaining the core achievements of $Λ$CDM. This opens a pathway to constraining a broad class of beyond-$Λ$CDM models with forthcoming observations.
△ Less
Submitted 23 September, 2025;
originally announced September 2025.
-
High-redshift AGN population in radiation-hydrodynamics simulations
Authors:
Teodora-Elena Bulichi,
Oliver Zier,
Aaron Smith,
Mark Vogelsberger,
Anna-Christina Eilers,
Rahul Kannan,
Xuejian Shen,
Ewald Puchwein,
Enrico Garaldi,
Josh Borrow
Abstract:
High-redshift active galactic nuclei (AGN) have long been recognized as key probes of early black hole growth and galaxy evolution. However, modeling this population remains difficult due to the wide range of luminosities and black hole masses involved, and the high computational costs of capturing the hydrodynamic response of gas and evolving radiation fields on-the-fly. In this study, we present…
▽ More
High-redshift active galactic nuclei (AGN) have long been recognized as key probes of early black hole growth and galaxy evolution. However, modeling this population remains difficult due to the wide range of luminosities and black hole masses involved, and the high computational costs of capturing the hydrodynamic response of gas and evolving radiation fields on-the-fly. In this study, we present a new suite of simulations based on the IllustrisTNG galaxy formation framework, enhanced with on-the-fly radiative transfer, to examine AGN at high redshift (z > 5) in a protocluster environment extracted from the MillenniumTNG simulation. We focus on the co-evolution of black holes and their host galaxies, as well as the radiative impact on surrounding intergalactic gas. The model predicts that black holes form in overdense regions and lie below the local black hole-stellar mass relation, with stellar mass assembly preceding significant black hole accretion. Ionizing photons are primarily produced by stars, which shape the morphology of ionized regions and drive reionization. Given the restrictive black hole growth in the original IllustrisTNG model, we reduce the radiative efficiency from 0.2 to 0.1, resulting in higher accretion rates for massive black holes, more bursty growth, and earlier AGN-driven quenching. However, the resulting AGN remain significantly fainter than observed high-redshift quasars. As such, to incorporate this missing population, we introduce a quasar boosted model, in which we artificially boost the AGN luminosity. This results in strong effects on the surrounding gas, most notably a proximity effect, and large contributions to He ionization.
△ Less
Submitted 20 October, 2025; v1 submitted 15 July, 2025;
originally announced July 2025.
-
The THESAN-ZOOM project: Mystery N/O more -- uncovering the origin of peculiar chemical abundances and a not-so-fundamental metallicity relation at $3<z<12$
Authors:
William McClymont,
Sandro Tacchella,
Aaron Smith,
Rahul Kannan,
Enrico Garaldi,
Ewald Puchwein,
Yuki Isobe,
Xihan Ji,
Xuejian Shen,
Zihao Wang,
Vasily Belokurov,
Josh Borrow,
Francesco D'Eugenio,
Laura Keating,
Roberto Maiolino,
Stephanie Monty,
Mark Vogelsberger,
Oliver Zier
Abstract:
We present an analysis of metallicities and chemical abundances at $3<z<12$ in the THESAN-ZOOM simulations. We find that smoothly curved gas-phase and stellar mass-metallicity relations (MZR) are already in place at $z\approx12$ and evolve slowly ($\sim$0.2 dex increase for gas, $\sim$0.4 dex increase for stars at a fixed stellar mass) down to $z=3$, governed largely by the efficiency with which g…
▽ More
We present an analysis of metallicities and chemical abundances at $3<z<12$ in the THESAN-ZOOM simulations. We find that smoothly curved gas-phase and stellar mass-metallicity relations (MZR) are already in place at $z\approx12$ and evolve slowly ($\sim$0.2 dex increase for gas, $\sim$0.4 dex increase for stars at a fixed stellar mass) down to $z=3$, governed largely by the efficiency with which galaxies retain their metals, rather than gas fraction. The canonical fundamental metallicity relation (FMR) survives in stars but breaks down and inverts for gas in low-mass galaxies ($M_\ast\lesssim10^{9}\mathrm{M_\odot}$) due to regular dilution by low-metallicity gas inflow. We find broad agreement of gas-phase N/O, Fe/O, and C/O with high-redshift observations, including the presence of nitrogen-rich galaxies (NRGs; $\log(\mathrm{N/O})>-0.6$) without the need for exotic yields in our chemical network. Instead, bursty star formation naturally generates order-of-magnitude excursions in N/O on $\lesssim$100 Myr timescales due to temporally differential galactic winds; after a starburst, stellar feedback expels gas, leaving a large population of asymptotic-giant-branch stars to dominate the enrichment of the relatively low-mass interstellar medium. NRGs lie below the main sequence and typically exhibit $\mathrm{EW}[H$β$]\lesssim40$ Å, in apparent tension with observed high-EW NRGs. This tension is reconciled if observed NRGs are in the initial stages of a subsequent starburst, illuminating previously enriched gas, which is supported by the finding of high SFR surface density nitrogen-rich giant molecular clouds.
△ Less
Submitted 11 July, 2025;
originally announced July 2025.
-
Protostellar disks in their natural habitat -- the formation of protostars and their accretion disks in the turbulent and magnetized interstellar medium
Authors:
Alexander C. Mayer,
Thorsten Naab,
Paola Caselli,
Alexei V. Ivlev,
Tommaso Grassi,
Oliver Zier,
Rüdiger Pakmor,
Stefanie Walch,
Volker Springel
Abstract:
We present simulations of the supernova-driven turbulent interstellar medium (ISM) in a simulation domain of volume $(256\,{\rm pc})^3$ within which we resolve the formation of protostellar accretion disks and their stellar cores to spatial scales of $\sim 10^{-4}$ au, using the moving-mesh code {\small AREPO}. We perform simulations with no magnetic fields, ideal magnetohydrodynamics (MHD) and am…
▽ More
We present simulations of the supernova-driven turbulent interstellar medium (ISM) in a simulation domain of volume $(256\,{\rm pc})^3$ within which we resolve the formation of protostellar accretion disks and their stellar cores to spatial scales of $\sim 10^{-4}$ au, using the moving-mesh code {\small AREPO}. We perform simulations with no magnetic fields, ideal magnetohydrodynamics (MHD) and ambipolar diffusion, and compare the resulting first Larson cores and their associated structures, including the accretion disks, their location within the larger-scale structure and the streamers connecting these. We find that disks of sizes $10-100\,{\rm au}$ form early in the simulations without magnetic fields, while there are no disks larger than 10 au with ideal MHD. Ambipolar diffusion causes large disks to form in a subset of cases (two out of six cores), and generally reduces the strength of outflows, which are seen to play a central role. When they are able to carry away significant angular momentum, they prevent the formation of a rotationally supported disk. Magnetic fields strengths grow from $0.1 - 1$ mG in the protostellar core to more than 10 G in the first Larson core in all simulations with ideal MHD. The rotationally supported disks which form can have rotation speeds $> 1$ km s$^{-1}$ even out to further than 100 au from the centre, become gravitationally unstable and form complex spiral substructures with Toomre $Q < 1$. We conclude that the impact of magnetic fields and non-ideal MHD on the formation of protostellar disks is substantial in realistic formation scenarios from the turbulent ISM.
△ Less
Submitted 3 September, 2025; v1 submitted 17 June, 2025;
originally announced June 2025.
-
Overmassive black holes in the early Universe can be explained by gas-rich, dark matter-dominated galaxies
Authors:
William McClymont,
Sandro Tacchella,
Xihan Ji,
Rahul Kannan,
Roberto Maiolino,
Charlotte Simmonds,
Aaron Smith,
Ewald Puchwein,
Enrico Garaldi,
Mark Vogelsberger,
Francesco D'Eugenio,
Laura Keating,
Xuejian Shen,
Bartolomeo Trefoloni,
Oliver Zier
Abstract:
JWST has revealed the apparent evolution of the black hole (BH)-stellar mass ($M_\mathrm{BH}$-$M_\rm{\ast}$) relation in the early Universe, while remaining consistent the BH-dynamical mass ($M_\mathrm{BH}$-$M_\mathrm{dyn}$) relation. We predict BH masses for $z>3$ galaxies in the high-resolution THESAN-ZOOM simulations by assuming the $M_\mathrm{BH}$-$M_\mathrm{dyn}$ relation is fundamental. Even…
▽ More
JWST has revealed the apparent evolution of the black hole (BH)-stellar mass ($M_\mathrm{BH}$-$M_\rm{\ast}$) relation in the early Universe, while remaining consistent the BH-dynamical mass ($M_\mathrm{BH}$-$M_\mathrm{dyn}$) relation. We predict BH masses for $z>3$ galaxies in the high-resolution THESAN-ZOOM simulations by assuming the $M_\mathrm{BH}$-$M_\mathrm{dyn}$ relation is fundamental. Even without live BH modelling, our approach reproduces the JWST-observed $M_\mathrm{BH}$ distribution, including overmassive BHs relative to the local $M_\mathrm{BH}$-$M_\mathrm{\ast}$ relation. We find that $M_\mathrm{BH}/M_\mathrm{\ast}$ declines with $M_\mathrm{\ast}$, evolving from $\sim$0.1 at $M_\mathrm{\ast}=10^6\,\mathrm{M_\odot}$ to $\sim$0.01 at $M_\mathrm{\ast}=10^{10.5}\,\mathrm{M_\odot}$. This trend reflects the dark matter ($f_\mathrm{DM}$) and gas fractions ($f_\mathrm{gas}$), which decrease with $M_\mathrm{\ast}$ but show little redshift evolution down to $z=3$, resulting in small $M_\mathrm{\ast}/M_\mathrm{dyn}$ ratios and thus overmassive BHs in low-mass galaxies. We use $\texttt{Prospector}$-derived stellar masses and star-formation rates to infer $f_\mathrm{gas}$ across 48,022 galaxies in JADES at $3<z<9$, finding excellent agreement with our simulation. Our results demonstrate that overmassive BHs would naturally result from a fundamental $M_\mathrm{BH}$-$M_\mathrm{dyn}$ relation and be typical of the gas-rich, dark matter-dominated nature of low-mass, high-redshift galaxies. Such overmassive, rapidly growing BHs may strongly influence the earliest stages of galaxy formation.
△ Less
Submitted 16 June, 2025;
originally announced June 2025.
-
The THESAN-ZOOM project: Star formation efficiency from giant molecular clouds to galactic scale in high-redshift starbursts
Authors:
Zihao Wang,
Xuejian Shen,
Mark Vogelsberger,
Hui Li,
Rahul Kannan,
Ewald Puchwein,
Aaron Smith,
Josh Borrow,
Enrico Garaldi,
Laura Keating,
Oliver Zier,
William McClymont,
Sandro Tacchella,
Yang Ni,
Lars Hernquist
Abstract:
Star formation in galaxies is inherently complex, involving the interplay of physical processes over a hierarchy of spatial scales. In this work, we investigate the connection between global (galaxy-scale) and local (cloud-scale) star formation efficiencies (SFEs) at high redshifts ($z\gtrsim 3$), using the state-of-the-art cosmological zoom-in simulation suite THESAN-ZOOM. We find that the galaxy…
▽ More
Star formation in galaxies is inherently complex, involving the interplay of physical processes over a hierarchy of spatial scales. In this work, we investigate the connection between global (galaxy-scale) and local (cloud-scale) star formation efficiencies (SFEs) at high redshifts ($z\gtrsim 3$), using the state-of-the-art cosmological zoom-in simulation suite THESAN-ZOOM. We find that the galaxy-scale average SFE, $\langle ε^{\rm gal}_{\rm ff} \rangle$, scales with $M_{\rm halo}^{1/3}\,(1+z)^{1/2} \sim V_{\rm vir}$, consistent with expectations from feedback-regulated models. On cloud scales, we identify giant molecular clouds (GMCs) in a broad sample of high-redshift starbursts spanning a wide range of halo masses and redshifts. Star formation in these systems is predominantly hosted by filamentary GMCs embedded in a dense and highly turbulent interstellar medium (ISM). GMCs exhibit remarkably universal properties, including mass function, size, turbulence, and surface density, regardless of the environment in which they are identified. The global gas depletion time (and the Kennicutt-Schmidt relation) is determined by the GMC mass fraction in the ISM, while the cloud-scale SFE shows little variation. In particular, we find a nearly constant gas surface density of $Σ_{\rm GMC} \approx 70\,{\rm M}_{\odot}\,{\rm pc}^{-2}$ across different host galaxies. Nevertheless, we identify two regimes where phases with high SFE can arise. First, stars may form efficiently in the shock fronts generated by feedback from a preceding starburst. Second, the increasing background dark matter surface density with redshift may contribute to the gravitational potential of clouds at $z \gtrsim 8$ and confine them in high-SFE phases over extended periods.
△ Less
Submitted 8 May, 2025;
originally announced May 2025.
-
The THESAN-ZOOM project: central starbursts and inside-out quenching govern galaxy sizes in the early Universe
Authors:
William McClymont,
Sandro Tacchella,
Aaron Smith,
Rahul Kannan,
Ewald Puchwein,
Josh Borrow,
Enrico Garaldi,
Laura Keating,
Mark Vogelsberger,
Oliver Zier,
Xuejian Shen,
Filip Popovic
Abstract:
We explore the evolution of galaxy sizes at high redshift ($3 < z < 13$) using the high-resolution THESAN-ZOOM radiation-hydrodynamics simulations, focusing on the mass range of $10^6\,\mathrm{M}_{\odot} < \mathrm{M}_{\ast} < 10^{10}\,\mathrm{M}_{\odot}$. Our analysis reveals that galaxy size growth is tightly coupled to bursty star formation. Galaxies above the star-forming main sequence experien…
▽ More
We explore the evolution of galaxy sizes at high redshift ($3 < z < 13$) using the high-resolution THESAN-ZOOM radiation-hydrodynamics simulations, focusing on the mass range of $10^6\,\mathrm{M}_{\odot} < \mathrm{M}_{\ast} < 10^{10}\,\mathrm{M}_{\odot}$. Our analysis reveals that galaxy size growth is tightly coupled to bursty star formation. Galaxies above the star-forming main sequence experience rapid central compaction during starbursts, followed by inside-out quenching and spatially extended star formation that leads to expansion, causing oscillatory behavior around the size-mass relation. Notably, we find a positive intrinsic size-mass relation at high redshift, consistent with observations but in tension with large-volume simulations. We attribute this discrepancy to the bursty star formation captured by our multi-phase interstellar medium framework, but missing from simulations using the effective equation-of-state approach with hydrodynamically decoupled feedback. We also find that the normalization of the size-mass relation follows a double power law as a function of redshift, with a break at $z\approx6$, because the majority of galaxies at $z > 6$ show rising star-formation histories, and therefore are in a compaction phase. We demonstrate that H$α$ emission is systematically extended relative to the UV continuum by a median factor of 1.7, consistent with recent JWST studies. However, in contrast to previous interpretations that link extended H$α$ sizes to inside-out growth, we find that Lyman-continuum (LyC) emission is spatially disconnected from H$α$. Instead, a simple Strömgren sphere argument reproduces observed trends, suggesting that extreme LyC production during central starbursts is the primary driver of extended nebular emission.
△ Less
Submitted 11 March, 2025; v1 submitted 6 March, 2025;
originally announced March 2025.
-
The THESAN-ZOOM project: Population III star formation continues until the end of reionization
Authors:
Oliver Zier,
Rahul Kannan,
Aaron Smith,
Ewald Puchwein,
Mark Vogelsberger,
Josh Borrow,
Enrico Garaldi,
Laura Keating,
William McClymont,
Xuejian Shen,
Lars Hernquist
Abstract:
Population III (Pop III) stars are the first stars in the Universe, forming from pristine, metal-free gas and marking the end of the cosmic dark ages. Their formation rate is expected to sharply decline after redshift $z \approx 15$ due to metal enrichment from previous generations of stars. In this paper, we analyze 14 zoom-in simulations from the THESAN-ZOOM project, which evolves different halo…
▽ More
Population III (Pop III) stars are the first stars in the Universe, forming from pristine, metal-free gas and marking the end of the cosmic dark ages. Their formation rate is expected to sharply decline after redshift $z \approx 15$ due to metal enrichment from previous generations of stars. In this paper, we analyze 14 zoom-in simulations from the THESAN-ZOOM project, which evolves different haloes from the THESAN-1 cosmological box down to redshift $z=3$. The high mass resolution of up to $142 M_\odot$ per cell in the gas phase combined with a multiphase model of the interstellar medium (ISM), radiative transfer including Lyman-Werner radiation, dust physics, and a non-equilibrium chemistry network that tracks molecular hydrogen, allows for a realistic but still approximate description of Pop III star formation in pristine gas. Our results show that Pop III stars continue to form in low-mass haloes ranging from $10^6 M_\odot$ to $10^9 M_\odot$ until the end of reionization at around $z=5$. At this stage, photoevaporation suppresses further star formation in these minihaloes, which subsequently merge into larger central haloes. Hence, the remnants of Pop III stars primarily reside in the satellite galaxies of larger haloes at lower redshifts. While direct detection of Pop III stars remains elusive, these results hint that lingering primordial star formation could leave observable imprints or indirectly affect the properties of high-redshift galaxies. Explicit Pop III feedback and specialized initial mass function modelling within the THESAN-ZOOM framework would further help interpreting emerging constraints from the James Webb Space Telescope.
△ Less
Submitted 20 June, 2025; v1 submitted 5 March, 2025;
originally announced March 2025.
-
The THESAN-ZOOM project: Long-term imprints of external reionization on galaxy evolution
Authors:
Oliver Zier,
Rahul Kannan,
Aaron Smith,
Ewald Puchwein,
Mark Vogelsberger,
Josh Borrow,
Enrico Garaldi,
Laura Keating,
William McClymont,
Xuejian Shen,
Lars Hernquist
Abstract:
We investigate the impact of ionizing external ultraviolet (UV) radiation on low-mass haloes ($M_{h}<10^{10}M_\odot$) at high redshift using $1140M_\odot$ baryonic resolution zoom-in simulations of seven regions from the THESAN-ZOOM project. We compare three simulation sets that differ in the treatment of external UV radiation: one employing a uniform UV background initiated at z=10.6 in addition…
▽ More
We investigate the impact of ionizing external ultraviolet (UV) radiation on low-mass haloes ($M_{h}<10^{10}M_\odot$) at high redshift using $1140M_\odot$ baryonic resolution zoom-in simulations of seven regions from the THESAN-ZOOM project. We compare three simulation sets that differ in the treatment of external UV radiation: one employing a uniform UV background initiated at z=10.6 in addition to radiation transport for local sources, another with the same background starting at z=5.5, and the default configuration in which the large-scale radiation field from the parent THESAN-1 simulation box acts as a boundary condition. The multi-phase interstellar medium (ISM) model, combined with its high mass resolution, allows us to resolve all star-forming haloes and capture the back-reaction of ionizing radiation on galaxy properties during the epoch of reionization. When present, external UV radiation efficiently unbinds gas in haloes with masses below $10^9M_\odot$ and suppresses subsequent star formation. As a result, in simulations with early reionization, minihaloes fail to form stars from pristine gas, leading to reduced metal enrichment of gas later accreted by more massive haloes. Consequently, haloes with masses below $10^{10}M_\odot$ at all simulated epochs (z>3) exhibit lower metallicities and altered metallicity distributions. The more accurate and realistic shielding from external UV radiation, achieved through self-consistent radiative transfer, permits the existence of a cold but low-density gas phase down to z=3. These findings highlight the importance of capturing a patchy reionization history in high-resolution simulations targeting high-redshift galaxy formation. We conclude that at minimum, a semi-numerical model that incorporates spatially inhomogeneous reionization and a non-uniform metallicity floor is necessary to accurately emulate metal enrichment in minihaloes.
△ Less
Submitted 25 June, 2025; v1 submitted 4 March, 2025;
originally announced March 2025.
-
The THESAN-ZOOM project: Star-formation efficiencies in high-redshift galaxies
Authors:
Xuejian Shen,
Rahul Kannan,
Ewald Puchwein,
Aaron Smith,
Mark Vogelsberger,
Josh Borrow,
Enrico Garaldi,
Laura Keating,
Oliver Zier,
William McClymont,
Sandro Tacchella,
Zihao Wang,
Lars Hernquist
Abstract:
Recent JWST observations hint at unexpectedly intense cosmic star-formation in the early Universe, often attributed to enhanced star-formation efficiencies (SFEs). Here, we analyze the SFE in THESAN-ZOOM, a novel zoom-in radiation-hydrodynamic simulation campaign of high-redshift ($z \gtrsim 3$) galaxies employing a state-of-the-art galaxy formation model resolving the multiphase interstellar medi…
▽ More
Recent JWST observations hint at unexpectedly intense cosmic star-formation in the early Universe, often attributed to enhanced star-formation efficiencies (SFEs). Here, we analyze the SFE in THESAN-ZOOM, a novel zoom-in radiation-hydrodynamic simulation campaign of high-redshift ($z \gtrsim 3$) galaxies employing a state-of-the-art galaxy formation model resolving the multiphase interstellar medium (ISM). The halo-scale SFE ($ε^{\ast}_{\rm halo}$) - the fraction of baryons accreted by a halo that are converted to stars - follows a double power-law dependence on halo mass, with a mild redshift evolution above $M_{\rm halo} \gtrsim 10^{9.5}\,{\rm M}_{\odot}$. The power-law slope is roughly $1/3$ at large halo masses, consistent with expectations when gas outflows are momentum-driven. At lower masses, the slope is roughly $2/3$ and is more aligned with the energy-driven outflow scenario. $ε^{\ast}_{\rm halo}$ is a factor of $2-3$ larger than commonly assumed in empirical galaxy-formation models at $M_{\rm halo} \lesssim 10^{11}\,{\rm M}_{\odot}$. On galactic (kpc) scales, the Kennicutt-Schmidt (KS) relation of neutral gas is universal in THESAN-ZOOM, following $Σ_{\rm SFR} \propto Σ_{\rm gas}^2$, indicative of a turbulent energy balance in the ISM maintained by stellar feedback. The rise of $ε^{\ast}_{\rm halo}$ with halo mass can be traced primarily to increasing gas surface densities in massive galaxies, while the underlying KS relation and neutral, star-forming gas fraction remain unchanged. Although the increase in $ε^{\ast}_{\rm halo}$ with redshift is relatively modest, it is sufficient to explain the large observed number density of UV-bright galaxies at $z \gtrsim 12$. However, reproducing the brightest sources at $M_{\rm UV} \lesssim -21$ may require extrapolating the SFE beyond the halo mass range directly covered by THESAN-ZOOM.
△ Less
Submitted 3 March, 2025;
originally announced March 2025.
-
The THESAN-ZOOM project: Burst, quench, repeat -- unveiling the evolution of high-redshift galaxies along the star-forming main sequence
Authors:
William McClymont,
Sandro Tacchella,
Aaron Smith,
Rahul Kannan,
Ewald Puchwein,
Josh Borrow,
Enrico Garaldi,
Laura Keating,
Mark Vogelsberger,
Oliver Zier,
Xuejian Shen,
Filip Popovic,
Charlotte Simmonds
Abstract:
Characterizing the evolution of the star-forming main sequence (SFMS) at high redshift is crucial to contextualize the observed extreme properties of galaxies in the early Universe. We present an analysis of the SFMS and its scatter in the THESAN-ZOOM simulations, where we find a redshift evolution of the SFMS normalization scaling as $\propto (1+z)^{2.64\pm0.03}$, significantly stronger than is t…
▽ More
Characterizing the evolution of the star-forming main sequence (SFMS) at high redshift is crucial to contextualize the observed extreme properties of galaxies in the early Universe. We present an analysis of the SFMS and its scatter in the THESAN-ZOOM simulations, where we find a redshift evolution of the SFMS normalization scaling as $\propto (1+z)^{2.64\pm0.03}$, significantly stronger than is typically inferred from observations. We can reproduce the flatter observed evolution by filtering out weakly star-forming galaxies, implying that current observational fits are biased due to a missing population of lulling galaxies or overestimated star-formation rates. We also explore star-formation variability using the scatter of galaxies around the SFMS ($σ_{\mathrm{MS}}$). At the population level, the scatter around the SFMS increases with cosmic time, driven by the increased importance of long-term environmental effects in regulating star formation at later times. To study short-term star-formation variability, or ''burstiness'', we isolate the scatter on timescales shorter than 50 Myr. The short-term scatter is larger at higher redshift, indicating that star formation is indeed more bursty in the early Universe. We identify two starburst modes: (i) externally driven, where rapid large-scale inflows trigger and fuel prolonged, extreme star formation episodes, and (ii) internally driven, where cyclical ejection and re-accretion of the interstellar medium in low-mass galaxies drive bursts, even under relatively steady large-scale inflow. Both modes occur at all redshifts, but the increased burstiness of galaxies at higher redshift is due to the increasing prevalence of the more extreme external mode of star formation.
△ Less
Submitted 11 March, 2025; v1 submitted 28 February, 2025;
originally announced March 2025.
-
Introducing the THESAN-ZOOM project: radiation-hydrodynamic simulations of high-redshift galaxies with a multi-phase interstellar medium
Authors:
Rahul Kannan,
Ewald Puchwein,
Aaron Smith,
Josh Borrow,
Enrico Garaldi,
Laura Keating,
Mark Vogelsberger,
Oliver Zier,
William McClymont,
Xuejian Shen,
Filip Popovic,
Sandro Tacchella,
Lars Hernquist,
Volker Springel
Abstract:
We introduce the THESAN-ZOOM project, a comprehensive suite of high-resolution zoom-in simulations of $14$ high-redshift ($z>3$) galaxies selected from the THESAN simulation volume. This sample encompasses a diverse range of halo masses, with $M_\mathrm{halo} \approx 10^8 - 10^{13}~\mathrm{M}_\odot$ at $z=3$. At the highest-resolution, the simulations achieve a baryonic mass of…
▽ More
We introduce the THESAN-ZOOM project, a comprehensive suite of high-resolution zoom-in simulations of $14$ high-redshift ($z>3$) galaxies selected from the THESAN simulation volume. This sample encompasses a diverse range of halo masses, with $M_\mathrm{halo} \approx 10^8 - 10^{13}~\mathrm{M}_\odot$ at $z=3$. At the highest-resolution, the simulations achieve a baryonic mass of $142~\mathrm{M}_\odot$ and a gravitational softening length of $17~\mathrm{cpc}$. We employ a state-of-the-art multi-phase interstellar medium (ISM) model that self-consistently includes stellar feedback, radiation fields, dust physics, and low-temperature cooling through a non-equilibrium thermochemical network. Our unique framework incorporates the impact of patchy reionization by adopting the large-scale radiation field topology from the parent THESAN simulation box rather than assuming a spatially uniform UV background. In total, THESAN-ZOOM comprises $60$ simulations, including both fiducial runs and complementary variations designed to investigate the impact of numerical and physical parameters on galaxy properties. The fiducial simulation set reproduces a wealth of high-redshift observational data such as the stellar-to-halo-mass relation, the star-forming main sequence, the Kennicutt-Schmidt relation, and the mass-metallicity relation. While our simulations slightly overestimate the abundance of low-mass and low-luminosity galaxies they agree well with observed stellar and UV luminosity functions at the higher mass end. Moreover, the star-formation rate density closely matches the observational estimates from $z=3-14$. These results indicate that the simulations effectively reproduce many of the essential characteristics of high-redshift galaxies, providing a realistic framework to interpret the exciting new observations from JWST.
△ Less
Submitted 8 October, 2025; v1 submitted 27 February, 2025;
originally announced February 2025.
-
A Novel Density Profile for Isothermal Cores of Dark Matter Halos
Authors:
Vinh Tran,
Xuejian Shen,
Mark Vogelsberger,
Daniel Gilman,
Stephanie O'Neil,
Cian Roche,
Oliver Zier,
Jiarun Gao
Abstract:
We present a novel analytic density profile for halos in self-interacting dark matter (SIDM) models, which accurately captures the isothermal-core configuration, i.e. where both the density and velocity dispersion profiles exhibit central plateaus in the halo innermost region. Importantly, the profile retains a simple and tractable functional form. We demonstrate analytically how our density profi…
▽ More
We present a novel analytic density profile for halos in self-interacting dark matter (SIDM) models, which accurately captures the isothermal-core configuration, i.e. where both the density and velocity dispersion profiles exhibit central plateaus in the halo innermost region. Importantly, the profile retains a simple and tractable functional form. We demonstrate analytically how our density profile satisfies the aforementioned conditions, with comparisons to other contemporary functional choices. We further validate the profile using idealized N-body simulations, showing that it provides excellent representations of both the density and velocity dispersion profiles across a broad range of evolutionary stages, from the early thermalization phase to the late core-collapse regime. As a result of its accuracy and simplicity, the proposed profile offers a robust framework for analyzing halo evolution in a variety of SIDM scenarios. It also holds practical utility in reducing simulation needs and in generating initial conditions for simulations targeting the deep core-collapse regime.
△ Less
Submitted 18 November, 2025; v1 submitted 18 November, 2024;
originally announced November 2024.
-
The THESAN project: tracking the expansion and merger histories of ionized bubbles during the Epoch of Reionization
Authors:
Nathan Jamieson,
Aaron Smith,
Meredith Neyer,
Rahul Kannan,
Enrico Garaldi,
Mark Vogelsberger,
Lars Hernquist,
Oliver Zier,
Xuejian Shen,
Koki Kakiichi
Abstract:
The growth of ionized hydrogen bubbles in the intergalactic medium around early luminous objects is a fundamental process during the Epoch of Reionization (EoR). In this study, we analyze bubble sizes and their evolution using the state-of-the-art THESAN radiation-hydrodynamics simulation suite, which self-consistently models radiation transport and realistic galaxy formation throughout a large (9…
▽ More
The growth of ionized hydrogen bubbles in the intergalactic medium around early luminous objects is a fundamental process during the Epoch of Reionization (EoR). In this study, we analyze bubble sizes and their evolution using the state-of-the-art THESAN radiation-hydrodynamics simulation suite, which self-consistently models radiation transport and realistic galaxy formation throughout a large (95.5 cMpc)^3 volume of the Universe. Analogous to the accretion and merger tree histories employed in galaxy formation simulations, we characterize the growth and merger rates of ionized bubbles by focusing on the spatially-resolved redshift of reionization. By tracing the chronological expansion of bubbles, we partition the simulation volume and construct a natural ionization history. We identify three distinct stages of ionized bubble growth: (1) initial slow expansion around the earliest ionizing sources seeding formation sites, (2) accelerated growth through percolation as bubbles begin to merge, and (3) rapid expansion dominated by the largest bubble. Notably, we find that the largest bubble emerges by z=9-10, well before the midpoint of reionization. This bubble becomes dominant during the second growth stage, and defines the third stage by rapidly expanding to eventually encompass the remainder of the simulation volume and becoming one of the few bubbles actively growing. Additionally, we observe a sharp decline in the number of bubbles with radii around ~10 cMpc compared to smaller sizes, indicating a characteristic scale in the final segmented bubble size distribution. Overall, these chronologically sequenced spatial reconstructions offer crucial insights into the physical mechanisms driving ionized bubble growth during the EoR and provide a framework for interpreting the structure and evolution of reionization itself.
△ Less
Submitted 13 November, 2024;
originally announced November 2024.
-
Moving-mesh non-ideal magnetohydrodynamical simulations of the collapse of cloud cores to protostars
Authors:
Alexander C. Mayer,
Oliver Zier,
Thorsten Naab,
Rüdiger Pakmor,
Paola Caselli,
Alexei V. Ivlev,
Volker Springel,
Stefanie Walch
Abstract:
Magnetic fields have been shown both observationally and through theoretical work to be an important factor in the formation of protostars and their accretion disks. Accurate modelling of the evolution of the magnetic field in low-ionization molecular cloud cores requires the inclusion of non-ideal magnetohydrodynamics (MHD) processes, specifically Ohmic and ambipolar diffusion and the Hall effect…
▽ More
Magnetic fields have been shown both observationally and through theoretical work to be an important factor in the formation of protostars and their accretion disks. Accurate modelling of the evolution of the magnetic field in low-ionization molecular cloud cores requires the inclusion of non-ideal magnetohydrodynamics (MHD) processes, specifically Ohmic and ambipolar diffusion and the Hall effect. These have a profound influence on the efficiency of magnetic removal of angular momentum from protostellar disks and simulations that include them can avoid the `magnetic-braking catastrophe' in which disks are not able to form. However, the impact of the Hall effect, in particular, is complex and remains poorly studied. In this work, we perform a large suite of simulations of the collapse of cloud cores to protostars with several non-ideal MHD chemistry models and initial core geometries using the moving-mesh code {\small AREPO}. We find that the efficiency of angular momentum removal is significantly reduced with respect to ideal MHD, in line with previous results. The Hall effect has a varied influence on the evolution of the disk which depends on the initial orientation of the magnetic field. This extends to the outflows seen in a subset of the models, where this effect can act to enhance or suppress them and open up new outflow channels. We conclude, in agreement with a subset of the previous literature, that the Hall effect is the dominant non-ideal MHD process in some collapse scenarios and thus should be included in simulations of protostellar disk formation.
△ Less
Submitted 10 January, 2025; v1 submitted 9 October, 2024;
originally announced October 2024.
-
Adapting AREPO-RT for Exascale Computing: GPU Acceleration and Efficient Communication
Authors:
Oliver Zier,
Rahul Kannan,
Aaron Smith,
Mark Vogelsberger,
Erkin Verbeek
Abstract:
Radiative transfer (RT) is a crucial ingredient for self-consistent modelling of numerous astrophysical phenomena across cosmic history. However, on-the-fly integration into radiation-hydrodynamics (RHD) simulations is computationally demanding, particularly due to the stringent time-stepping conditions and increased dimensionality inherent in multi-frequency collisionless Boltzmann physics. The e…
▽ More
Radiative transfer (RT) is a crucial ingredient for self-consistent modelling of numerous astrophysical phenomena across cosmic history. However, on-the-fly integration into radiation-hydrodynamics (RHD) simulations is computationally demanding, particularly due to the stringent time-stepping conditions and increased dimensionality inherent in multi-frequency collisionless Boltzmann physics. The emergence of exascale supercomputers, equipped with extensive CPU cores and GPU accelerators, offers new opportunities for enhancing RHD simulations. We present a novel optimization of AREPO-RT explicitly tailored for such high-performance computing environments. We implement a novel node-to-node communication strategy that utilizes shared memory to substitute intra-node communication with direct memory access. Furthermore, combining multiple inter-node messages into a single message substantially enhances network bandwidth utilization and performance for large-scale simulations on modern supercomputers. The single-message node-to-node approach also improves performance on smaller-scale machines with less optimized networks. Furthermore, by transitioning all RT-related calculations to GPUs, we achieve a significant computational speedup of around 15 for standard benchmarks compared to the original CPU implementation. As a case study, we perform cosmological RHD simulations of the Epoch of Reionization, employing a similar setup as the THESAN project. In this context, RT becomes sub-dominant such that even without modifying the core AREPO codebase, there is an overall threefold improvement in efficiency. The advancements presented here have broad implications, potentially transforming the complexity and scalability of future simulations for a wide variety of astrophysical studies.
△ Less
Submitted 26 April, 2024;
originally announced April 2024.
-
The THESAN project: galaxy sizes during the epoch of reionization
Authors:
Xuejian Shen,
Mark Vogelsberger,
Josh Borrow,
Yongao Hu,
Evan Erickson,
Rahul Kannan,
Aaron Smith,
Enrico Garaldi,
Lars Hernquist,
Takahiro Morishita,
Sandro Tacchella,
Oliver Zier,
Guochao Sun,
Anna-Christina Eilers,
Hui Wang
Abstract:
We investigate galaxy sizes at redshift $z\gtrsim 6$ with the cosmological radiation-magneto-hydrodynamic simulation suite THESAN(-HR). These simulations simultaneously capture the reionization of the large-scale intergalactic medium and resolved galaxy properties. The intrinsic size ($r^{\ast}_{1/2}$) of simulated galaxies increases moderately with stellar mass at…
▽ More
We investigate galaxy sizes at redshift $z\gtrsim 6$ with the cosmological radiation-magneto-hydrodynamic simulation suite THESAN(-HR). These simulations simultaneously capture the reionization of the large-scale intergalactic medium and resolved galaxy properties. The intrinsic size ($r^{\ast}_{1/2}$) of simulated galaxies increases moderately with stellar mass at $M_{\ast} \lesssim 10^{8}\,{\rm M}_{\odot}$ and decreases fast at larger masses, resulting in a hump feature at $M_{\ast}\sim 10^{8}\,{\rm M}_{\odot}$ that is insensitive to redshift. Low-mass galaxies are in the initial phase of size growth and are better described by a spherical shell model with feedback-driven gas outflows competing with the cold inflows. In contrast, massive galaxies fit better with the disk formation model. They generally experience a phase of rapid compaction and gas depletion, likely driven by internal disk instability rather than external processes. We identify four compact quenched galaxies in the $(95.5\,{\rm cMpc})^{3}$ volume of THESAN-1 at $z\simeq 6$, and their quenching follows reaching a characteristic stellar surface density akin to the massive compact galaxies at cosmic noon. Compared to observations, we find that the median UV effective radius ($R^{\rm UV}_{\rm eff}$) of simulated galaxies is at least three times larger than the observed ones at $M_{\ast}\lesssim 10^{9}\,{\rm M}_{\odot}$ or $M_{\rm UV}\gtrsim -20$ at $6 \lesssim z \lesssim 10$. This inconsistency, related to the hump feature of the intrinsic size--mass relation, persists across many other cosmological simulations with different galaxy formation models and demonstrates the potential of using galaxy morphology to constrain the physics of galaxy formation at high redshifts.
△ Less
Submitted 16 September, 2024; v1 submitted 13 February, 2024;
originally announced February 2024.
-
Non-ideal magnetohydrodynamics on a moving mesh II: Hall effect
Authors:
Oliver Zier,
Alexander C. Mayer,
Volker Springel
Abstract:
In this work we extend the non-ideal magnetohydrodynamics (MHD) solver in the moving mesh code AREPO to include the Hall effect. The core of our algorithm is based on an estimation of the magnetic field gradients by a least-square reconstruction on the unstructured mesh, which we also used in the companion paper for Ohmic and ambipolar diffusion. In an extensive study of simulations of a magnetic…
▽ More
In this work we extend the non-ideal magnetohydrodynamics (MHD) solver in the moving mesh code AREPO to include the Hall effect. The core of our algorithm is based on an estimation of the magnetic field gradients by a least-square reconstruction on the unstructured mesh, which we also used in the companion paper for Ohmic and ambipolar diffusion. In an extensive study of simulations of a magnetic shock, we show that without additional magnetic diffusion our algorithm for the Hall effect becomes unstable at high resolution. We can however stabilise it by artificially increasing the Ohmic resistivity, $η_{\rm OR}$, so that it satisfies the condition $η_{\rm OR} \geq η_{\rm H} /5$, where $η_{\rm H}$ is the Hall diffusion coefficient. Adopting this solution we find second order convergence for the C-shock and are also able to accurately reproduce the dispersion relation of the whistler waves. As a first application of the new scheme, we simulate the collapse of a magnetised cloud with constant Hall parameter $η_{\rm H}$ and show that, depending on the sign of $η_{\rm H}$, the magnetic braking can either be weakened or strengthened by the Hall effect. The quasi-Lagrangian nature of the moving mesh method used here automatically increases the resolution in the forming core, making it well suited for more realistic studies with non-constant magnetic diffusivities in the future.
△ Less
Submitted 27 September, 2023;
originally announced September 2023.
-
The Accretion Mode in Sub-Eddington Supermassive Black Holes: Getting into the Central Parsecs of Andromeda
Authors:
C. Alig,
A. Prieto,
M. Blaña,
M. Frischman,
C. Metzl,
A. Burkert,
O. Zier,
A. Streblyanska
Abstract:
The inner kiloparsec regions surrounding sub-Eddington (luminosity less than 10$^{-3}$ in Eddington units, L$_{Edd}$) supermassive black holes (BHs) often show a characteristic network of dust filaments that terminate in a nuclear spiral in the central parsecs. Here we study the role and fate of these filaments in one of the least accreting BHs known, M31 (10$^{-7}$ L$_{Edd}$) using hydrodynamical…
▽ More
The inner kiloparsec regions surrounding sub-Eddington (luminosity less than 10$^{-3}$ in Eddington units, L$_{Edd}$) supermassive black holes (BHs) often show a characteristic network of dust filaments that terminate in a nuclear spiral in the central parsecs. Here we study the role and fate of these filaments in one of the least accreting BHs known, M31 (10$^{-7}$ L$_{Edd}$) using hydrodynamical simulations. The evolution of a streamer of gas particles moving under the barred potential of M31 is followed from kiloparsec distance to the central parsecs. After an exploratory study of initial conditions, a compelling fit to the observed dust/ionized gas morphologies and line-of-sight velocities in the inner hundreds of parsecs is produced. After several million years of streamer evolution, during which friction, thermal dissipation, and self-collisions have taken place, the gas settles into a disk tens of parsecs wide. This is fed by numerous filaments that arise from an outer circumnuclear ring and spiral toward the center. The final configuration is tightly constrained by a critical input mass in the streamer of several 10$^3$ M$_{\odot}$ (at an injection rate of 10$^{-4}$ M$_{\odot}$ yr$^{-1}$); values above or below this lead to filament fragmentation or dispersion respectively, which are not observed. The creation of a hot gas atmosphere in the region of $\sim$10$^6$ K is key to the development of a nuclear spiral during the simulation. The final inflow rate at 1pc from the center is $\sim$1.7 $\times$ 10$^{-7}$ M$_{\odot}$ yr$^{-1}$, consistent with the quiescent state of the M31 BH.
△ Less
Submitted 9 August, 2023;
originally announced August 2023.
-
Non-ideal magnetohydrodynamics on a moving mesh I: Ohmic and ambipolar diffusion
Authors:
Oliver Zier,
Volker Springel,
Alexander C. Mayer
Abstract:
Especially in cold and high-density regions, the assumptions of ideal magnetohydrodynamics (MHD) can break down, making first order non-ideal terms such as Ohmic and ambipolar diffusion as well as the Hall effect important. In this study we present a new numerical scheme for the first two resistive terms, which we implement in the moving-mesh code AREPO using the single-fluid approximation combine…
▽ More
Especially in cold and high-density regions, the assumptions of ideal magnetohydrodynamics (MHD) can break down, making first order non-ideal terms such as Ohmic and ambipolar diffusion as well as the Hall effect important. In this study we present a new numerical scheme for the first two resistive terms, which we implement in the moving-mesh code AREPO using the single-fluid approximation combined with a new gradient estimation technique based on a least-squares fit per interface. Through various test calculations including the diffusion of a magnetic peak, the structure of a magnetic C-shock, and the damping of an Alfvén wave, we show that we can achieve an accuracy comparable to the state-of-the-art code ATHENA++. We apply the scheme to the linear growth of the magnetorotational instability and find good agreement with the analytical growth rates. By simulating the collapse of a magnetised cloud with constant magnetic diffusion, we show that the new scheme is stable even for large density contrasts. Thanks to the Lagrangian nature of the moving mesh method the new scheme is thus well suited for intended future applications where a high resolution in the dense cores of collapsing protostellar clouds needs to be achieved. In a forthcoming work we will extend the scheme to the Hall effect.
△ Less
Submitted 21 July, 2023;
originally announced July 2023.
-
Gravito-turbulence in local disk simulations with an adaptive moving mesh
Authors:
Oliver Zier,
Volker Springel
Abstract:
Self-gravity plays an important role in the evolution of rotationally supported systems such as protoplanetary disks, accretion disks around black holes, or galactic disks, as it can both feed turbulence or lead to gravitational fragmentation. While such systems can be studied in the shearing box approximation with high local resolution, the large density contrasts that are possible in the case of…
▽ More
Self-gravity plays an important role in the evolution of rotationally supported systems such as protoplanetary disks, accretion disks around black holes, or galactic disks, as it can both feed turbulence or lead to gravitational fragmentation. While such systems can be studied in the shearing box approximation with high local resolution, the large density contrasts that are possible in the case of fragmentation still limit the utility of Eulerian codes with constant spatial resolution. In this paper, we present a novel self-gravity solver for the shearing box based on the TreePM method of the moving-mesh code AREPO. The spatial gravitational resolution is adaptive which is important to make full use of the quasi-Lagrangian hydrodynamical resolution of the code. We apply our new implementation to two- and three-dimensional, self-gravitating disks combined with a simple $β$-cooling prescription. For weak cooling we find a steady, gravito-turbulent state, while for strong cooling the formation of fragments is inevitable. To reach convergence for the critical cooling efficiency above which fragmentation occurs, we require a smoothing of the gravitational force in the two dimensional case that mimics the stratification of the three-dimensional simulations. The critical cooling efficiency we find, $β\approx 3$, as well as box-averaged quantities characterizing the gravito-turbulent state, agree well with various previous results in the literature. Interestingly, we observe stochastic fragmentation for $β> 3$, which slightly decreases the cooling efficiency required to observe fragmentation over the lifetime of a protoplanetary disk. The numerical method outlined here appears well suited to study the problem of galactic disks as well as magnetized, self-gravitating disks.
△ Less
Submitted 27 September, 2023; v1 submitted 5 December, 2022;
originally announced December 2022.
-
Simulating the magnetorotational instability on a moving mesh with the shearing box approximation
Authors:
Oliver Zier,
Volker Springel
Abstract:
The magnetorotational instability (MRI) is an important process in sufficiently ionized accretion disks, as it can create turbulence that acts as an effective viscosity, mediating angular momentum transport. Due to its local nature, it is often analyzed in the shearing box approximation with Eulerian methods, which otherwise would suffer from large advection errors in global disk simulations. In t…
▽ More
The magnetorotational instability (MRI) is an important process in sufficiently ionized accretion disks, as it can create turbulence that acts as an effective viscosity, mediating angular momentum transport. Due to its local nature, it is often analyzed in the shearing box approximation with Eulerian methods, which otherwise would suffer from large advection errors in global disk simulations. In this work, we report on an extensive study that applies the quasi-Lagrangian, moving-mesh code AREPO, combined with the Dedner cleaning scheme to control deviations from $\nabla\cdot B=0$, to the problem of magnetized flows in shearing boxes. We find that we can resolve the analytical linear growth rate of the MRI with mean background magnetic field well. In the zero net flux case, there is a threshold value for the strength of the divergence cleaning above which the turbulence eventually dies out, and in contrast to previous Eulerian simulations, the strength of the MRI does not decrease with increasing resolution. In boxes with larger vertical aspect ratio we find a mean-field dynamo, as well as an active shear current effect that can sustain MRI turbulence for at least 200 orbits. In stratified simulations, we obtain an active $αω$ dynamo and the characteristic butterfly diagram. Our results compare well with previous results obtained with static grid codes such as ATHENA. We thus conclude that AREPO represents an attractive approach for global disk simulations due to its quasi-Lagrangian nature, and for shearing box simulations with large density variations due to its continuously adaptive resolution.
△ Less
Submitted 27 September, 2023; v1 submitted 1 August, 2022;
originally announced August 2022.
-
Simulating cold shear flows on a moving mesh
Authors:
Oliver Zier,
Volker Springel
Abstract:
Rotationally supported, cold, gaseous disks are ubiquitous in astrophysics and appear in a diverse set of systems, such as protoplanetary disks, accretion disks around black holes, or large spiral galaxies. Capturing the gas dynamics accurately in these systems is challenging in numerical simulations due to the low sound speed compared to the bulk velocity of the gas, the resolution limitations of…
▽ More
Rotationally supported, cold, gaseous disks are ubiquitous in astrophysics and appear in a diverse set of systems, such as protoplanetary disks, accretion disks around black holes, or large spiral galaxies. Capturing the gas dynamics accurately in these systems is challenging in numerical simulations due to the low sound speed compared to the bulk velocity of the gas, the resolution limitations of full disk models, and the fact that numerical noise can easily source spurious growth of fluid instabilities if not suppressed sufficiently well, negatively interfering with real physical instabilities present in such disks (like the magneto-rotational instability). Here we implement the so-called shearing-box approximation in the moving-mesh code AREPO in order to facilitate achieving high resolution in local regions of differentially rotating disks and to address these problems. While our new approach offers manifest translational invariance across the shearing-box boundaries and offers continuous local adaptivity, we demonstrate that the unstructured mesh of AREPO introduces unwanted levels of "grid-noise" in the default version of the code. We show that this can be rectified by high-order integrations of the flux over mesh boundaries. With our new techniques we obtain highly accurate results for shearing-box calculations of the magneto-rotational instability that are superior to other Lagrangian techniques. These improvements are also of value for other applications of the code that feature strong shear flows.
△ Less
Submitted 27 September, 2023; v1 submitted 16 May, 2022;
originally announced May 2022.
-
On the interaction of a Bonnor-Ebert sphere with a stellar wind
Authors:
Oliver Zier,
Andreas Burkert,
Christian Alig
Abstract:
The structure of protostellar cores can often be approximated by isothermal Bonnor-Ebert spheres (BES) which are stabilized by an external pressure. For the typical pressure of $10^4k_B\,\mathrm{K\,cm^{-3}}$ to $10^5k_B\,\mathrm{K\,cm^{-3}}$ found in molecular clouds, cores with masses below $1.5\,{\rm M_\odot}$ are stable against gravitational collapse. In this paper, we analyze the efficiency of…
▽ More
The structure of protostellar cores can often be approximated by isothermal Bonnor-Ebert spheres (BES) which are stabilized by an external pressure. For the typical pressure of $10^4k_B\,\mathrm{K\,cm^{-3}}$ to $10^5k_B\,\mathrm{K\,cm^{-3}}$ found in molecular clouds, cores with masses below $1.5\,{\rm M_\odot}$ are stable against gravitational collapse. In this paper, we analyze the efficiency of triggering a gravitational collapse by a nearby stellar wind, which represents an interesting scenario for triggered low-mass star formation. We derive analytically a new stability criterion for a BES compressed by a stellar wind, which depends on its initial nondimensional radius $ξ_{max}$. If the stability limit is violated the wind triggers a core collapse. Otherwise, the core is destroyed by the wind. We estimate its validity range to $2.5<ξ_{max}<4.2$ and confirm this in simulations with the SPH Code GADGET-3. The efficiency to trigger a gravitational collapse strongly decreases for $ξ_{max}<2.5$ since in this case destruction and acceleration of the whole sphere begin to dominate. We were unable to trigger a collapse for $ξ_{max}<2$, which leads to the conclusion that a stellar wind can move the smallest unstable stellar mass to $0.5\,\mathrm{M_\odot}$ and destabilizing even smaller cores would require an external pressure larger than $10^5k_B\,\mathrm{K\,cm^{-3}}$. For $ξ_{max}>4.2$ the expected wind strength according to our criterion is small enough so that the compression is slower than the sound speed of the BES and sound waves can be triggered. In this case our criterion underestimates somewhat the onset of collapse and detailed numerical analyses are required.
△ Less
Submitted 10 May, 2021;
originally announced May 2021.
-
Simulating cosmic structure formation with the GADGET-4 code
Authors:
Volker Springel,
Rüdiger Pakmor,
Oliver Zier,
Martin Reinecke
Abstract:
Numerical methods have become a powerful tool for research in astrophysics, but their utility depends critically on the availability of suitable simulation codes. This calls for continuous efforts in code development, which is necessitated also by the rapidly evolving technology underlying today's computing hardware. Here we discuss recent methodological progress in the GADGET code, which has been…
▽ More
Numerical methods have become a powerful tool for research in astrophysics, but their utility depends critically on the availability of suitable simulation codes. This calls for continuous efforts in code development, which is necessitated also by the rapidly evolving technology underlying today's computing hardware. Here we discuss recent methodological progress in the GADGET code, which has been widely applied in cosmic structure formation over the past two decades. The new version offers improvements in force accuracy, in time-stepping, in adaptivity to a large dynamic range in timescales, in computational efficiency, and in parallel scalability through a special MPI/shared-memory parallelization and communication strategy, and a more-sophisticated domain decomposition algorithm. A manifestly momentum conserving fast multipole method (FMM) can be employed as an alternative to the one-sided TreePM gravity solver introduced in earlier versions. Two different flavours of smoothed particle hydrodynamics, a classic entropy-conserving formulation and a pressure-based approach, are supported for dealing with gaseous flows. The code is able to cope with very large problem sizes, thus allowing accurate predictions for cosmic structure formation in support of future precision tests of cosmology, and at the same time is well adapted to high dynamic range zoom-calculations with extreme variability of the particle number density in the simulated volume. The GADGET-4 code is publicly released to the community and contains infrastructure for on-the-fly group and substructure finding and tracking, as well as merger tree building, a simple model for radiative cooling and star formation, a high dynamic range power spectrum estimator, and an initial conditions generator based on second-order Lagrangian perturbation theory.
△ Less
Submitted 18 August, 2021; v1 submitted 7 October, 2020;
originally announced October 2020.