Researcher profile

Thorsten Naab

Thorsten Naab contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2026arXiv

FROST-CLUSTERS -- III. Metallicity-dependent intermediate mass black hole formation by runaway collisions in dense star clusters

We explore the formation of intermediate mass black holes (IMBHs), potential seeds for supermassive black holes (SMBHs), via runaway stellar collisions for a wide range of star cluster (surface) densities ($4\times10^3 M_\odot$ pc$^{-2} \lesssim Σ_\mathrm{h}$ $\lesssim 4\times10^6 M_\odot$ pc$^{-2}$) and metallicities ($0.01 Z_\odot \lesssim Z \lesssim 1.0 Z_\odot)$. Our sample of isolated $(>1400)$ and hierarchical ($30$) simulations of young, massive star clusters with up to $N=1.8\times10^6$ stars includes collisional stellar dynamics, stellar evolution, and post-Newtonian equations of motion for black holes using the BIFROST code. High stellar wind rates suppress IMBH formation at high metallicities $(Z \gtrsim 0.2 Z_\odot)$ and low collision rates prevent their formation at low densities ($Σ_\mathrm{h} \lesssim 3\times10^4 M_\odot$ pc$^{-2}$). The assumptions about stellar wind loss rates strongly affect the maximum final IMBH masses $(M_\bullet \sim 6000 M_\odot$ vs. $25000 M_\odot$). The total stellar mass loss from collisions and collisionally boosted winds before $t=3$ Myr can together reach up to 5-10% of the final cluster mass. We present fitting formulae for IMBH masses as a function of host star cluster $Σ_\mathrm{h}$ and Z, and formulate a model for the cosmic IMBH formation rate density. Depending on the cluster birth densities, the IMBH formation rates peak at $z\sim2$-$4$ at up to $\sim10^{-7}$ yr$^{-1}$cMpc$^{-3}$. As more than 50% form below $z\lesssim1.5$-$3$, the model challenges a view in which all local IMBHs are failed early Universe SMBH seeds.

preprint2025arXiv

NOEMA$^{\rm 3D}$: A first kpc resolution study of a $z\sim1.5$ main sequence barred galaxy channeling gas into a growing bulge

We present a very deep CO(3-2) observation of a massive, gas-rich, main sequence, barred spiral galaxy at $z\approx1.52$. Our data were taken with the IRAM-NOEMA interferometer for a 12-antenna equivalent on-source integration time of $\sim$ 50 hours. We fit the major axis kinematics using forward modelling of a rotating disk, and then subtract the two-dimensional beam convolved best-fit model revealing signatures of planar non-circular motions in the residuals. The inferred in-plane radial velocities are remarkably large, of the order of $\approx60$ km/s. Direct comparisons with a high-resolution, simulated, gas-rich, barred galaxy, obtained with the moving mesh code AREPO and the TNG sub-grid model, show that the observed non-circular gas flows can be explained as radial flows driven by the central bar, with an inferred net inflow rate of the order of the SFR. Given the recent evidence for a higher-than-expected fraction of barred disk galaxies at cosmic noon, our results suggest that rapid gas inflows due to bars could be important evolutionary drivers for the dominant population of star-forming galaxies at the peak epoch of star and galaxy formation.

preprint2025arXiv

Rapid formation of a very massive star >50000 $M_\odot$ and subsequently an IMBH from runaway collisions. Direct N-body and Monte Carlo simulations of dense star clusters

Context. We present simulations of a massive young star cluster using \textsc{Nbody6++GPU} and \textsc{MOCCA}. The cluster is initially more compact than previously published models, with one million stars, a total mass of $5.86 \times 10^5~\mathrm{M}_{\odot}$, and a half-mass radius of $0.1~\mathrm{pc}$. Aims. We analyse the formation and growth of a very massive star (VMS) through successive stellar collisions and investigate the subsequent formation of an intermediate-mass black hole (IMBH) in the core of a dense star cluster. Methods. We use both direct \textit{N}-body and Monte Carlo simulations, incorporating updated stellar evolution prescriptions (SSE/BSE) tailored to massive stars and VMSs. These include revised treatments of stellar radii, rejuvenation, and mass loss during collisions. While the prescriptions represent reasonable extrapolations into the VMS regime, the internal structure and thermal state of VMSs formed through stellar collisions remain uncertain, and future work may require further refinement. Results. We find that runaway stellar collisions in the cluster core produce a VMS exceeding $5 \times 10^4~\mathrm{M}_{\odot}$ within 5 Myr, which subsequently collapses into an IMBH. Conclusions. Our model suggests that dense stellar environments may enable the formation of very massive stars and massive black hole seeds through runaway stellar collisions. These results provide a potential pathway for early black hole growth in star clusters and offer theoretical context for interpreting recent JWST observations of young, compact clusters at high redshift.

preprint2025arXiv

The Noctua Suite of Simulations -- The Difficulty of Growing Massive Black Holes in Low-Mass Dwarf Galaxies

Aims. We study the individual and cumulative impact of stellar feedback processes on massive black hole (MBH) growth in a simulated low-mass dwarf galaxy. Methods. A suite of high-resolution radiation-hydrodynamic simulations called Noctua is performed, using the ArepoNoctua numerical framework for BHs in galaxy simulations. The chemical evolution of the gas is explicitly modelled in a time-dependent non-equilibrium way. Two types of stellar feedback are considered: individually-traced type II supernova (SNII) explosions, and radiatively transferred (on-the-fly) ionising stellar radiation (ISR) from OB stars. As part of the numerical framework, we develop and apply a novel physically-motivated model for MBH gas accretion, taking into account the angular momentum of the gas in the radiatively efficient regime, to estimate the gas accretion rate from the sub-grid accretion disc. Results. Without any stellar feedback, an initial $10^4~\mathrm{M}_\odot$ MBH is able to steadily grow over time, roughly doubling its mass after 800 Myr. Surprisingly, the growth of the MBH is more than doubled when only ISR feedback is considered, compared to the no stellar feedback run. This is due to the star formation rate (SFR) being highly suppressed (to a similar level or slightly above that when SNII feedback is considered), enabling a higher cumulative net gas inflow onto the MBH from not only the cold neutral- and molecular medium phases, but also the unstable- and warm neutral medium phases. With SNII feedback included, the gas accretion onto the MBH is episodic over time, and is suppressed by more than an order of magnitude already during the first 150 Myr. When combining SNII with ISR feedback, the growth of the MBH remains suppressed due to SNII feedback, but to a lesser extent compared to the SNII-only feedback run, due to a slightly lower SFR, and hence a reduced number of SNII events.

preprint2024arXiv

Evolution of eccentric stellar disks around supermassive black holes: the complex disk disruption dynamics and the milliparsec stars

We study the 10 Myr evolution of parsec-scale stellar disks with initial masses of $M_{\mathrm{disk}} = 1.0$ - $7.5 \times 10^4 M_\odot$ and eccentricities $e_\mathrm{init}=0.1$-$0.9$ around supermassive black holes (SMBHs). Our disk models are embedded in a spherical background potential and have top-heavy single and binary star initial mass functions (IMF) with slopes of $0.25$-$1.7$. The systems are evolved with the N-body code $\texttt{BIFROST}$ including post-Newtonian (PN) equations of motion and simplified stellar evolution. All disks are unstable and evolve on Myr timescales towards similar eccentricity distributions peaking at $e_\star \sim 0.3$-$0.4$. Models with high $e_\mathrm{init}$ also develop a very eccentric $(e_\star\gtrsim0.9)$ stellar population. For higher disk masses $M_\mathrm{disk} \gtrsim3 \times10^4\;\mathrm{M_\odot}$, the disk disruption dynamics is more complex than the standard secular eccentric disk instability with opposite precession directions at different disk radii - a precession direction instability. We present an analytical model describing this behavior. A milliparsec population of $N\sim10$-$100$ stars forms around the SMBH in all models. For low $e_\mathrm{init}$ stars migrate inward while for $e_\mathrm{init}\gtrsim0.6$ stars are captured by the Hills mechanism. Without PN, after $6$ Myr the captured stars have a sub-thermal eccentricity distribution. We show that including PN effects prevents this thermalization by suppressing resonant relaxation effects and cannot be ignored. The number of tidally disrupted stars is similar or larger than the number of milliparsec stars. None of the simulated models can simultaneously reproduce the kinematic and stellar population properties of the Milky Way center clockwise disk and the S-cluster.

preprint2022arXiv

A panchromatic view of star cluster formation in a simulated dwarf galaxy starburst

We present a photometric analysis of star and star cluster (SC) formation in a high-resolution simulation of a dwarf galaxy starburst that allows the formation of individual stars to be followed. Previous work demonstrated that the properties of the SCs formed in the simulation are in good agreement with observations. In this paper, we create mock spectral energy distributions and broad-band photometric images using the radiative transfer code SKIRT 9. We test several observational star formation rate (SFR) tracers and find that $24$ $μ$m, total infrared and H$α$ trace the underlying SFR during the (post)starburst phase, while UV tracers yield a more accurate picture of star formation during quiescent phases prior to and after the merger. We then place the simulated galaxy at distances of $10$ and $50$ Mpc and use aperture photometry at Hubble Space Telescope resolution to analyse the simulated SC population. During the starburst phase, a hierarchically forming set of SCs leads inaccurate source separation because of crowding. This results in estimated SC mass function slopes that are up to $\sim0.3$ shallower than the true slope of $\sim-1.9$ to $-2$ found for the bound clusters identified from the particle data in the simulation. The masses of the largest clusters are overestimated by a factor of up to $2.9$ due to unresolved clusters within the apertures. The aperture-based analysis also produces a relation between cluster formation efficiency and SFR surface density that is slightly flatter than that recovered from bound clusters. The differences are strongest in quiescent SF environments.

preprint2022arXiv

LYRA II: Cosmological dwarf galaxy formation with inhomogeneous Population III enrichment

We present the simulation of a $2\times10^9 M_\odot$ halo mass cosmological dwarf galaxy run to $z=0$ at 4 solar mass gas resolution with resolved supernova feedback. We compare three simple subgrid implementations for the inhomogeneous chemical enrichment from Population III stars and compare them to constraints from Local Group dwarf galaxies. The employed model, LYRA, is a novel high resolution galaxy formation model built for the moving mesh code AREPO, which is marked by a resolved multi-phase interstellar medium, single stars and individual supernova events. The resulting reionization relic is characterized by a short ($<1.5$ Gyr) star formation history that is repeatedly brought to a standstill by violent bursts of feedback. Star formation is reignited for a short duration due to a merger at $z\approx4$ and then again at $z\approx0.2-0$ after sustained gas accretion. Our model $z=0$ galaxy matches the stellar mass, size, stellar kinematics and metallicity relations of Local Group dwarf galaxies well. The dark matter profile does not exhibit a core in any version of the model. We show that the host halo masses of Population III stars affect the assembly history of dwarf galaxies. This manifests itself through the initial gaseous collapse in the progenitor halos, affecting the central density of the stellar component and through the accretion of luminous substructure.

preprint2022arXiv

Signatures of the Many Supermassive Black Hole Mergers in a Cosmologically Forming Massive Early-Type Galaxy

We model here the merger histories of the supermassive black hole (SMBH) population in the late stages of a cosmological simulation of a $\sim 2 \times 10^{13} M_\odot$ galaxy group. The gravitational dynamics around the several tens of SMBHs ($M_{\bullet} \gtrsim 7.5\times 10^7 M_\odot$) hosted by the galaxies in the group is computed at high accuracy using regularized integration with the KETJU code. The 11 SMBHs which form binaries and hierarchical triplets eventually merge after hardening through dynamical friction, stellar scattering, and gravitational wave (GW) emission. The binaries form at eccentricities of $e \sim 0.3$-$0.9$, with one system evolving to a very high eccentricity of $e = 0.998$, and merge on timescales of a few tens to several hundred megayears. During the simulation the merger-induced GW recoil kicks eject one SMBH remnant from the central host galaxy. This temporarily drives the galaxy off the $M_{\bullet}$-$σ_{\star}$ relation, however the galaxy returns to the relation due to subsequent galaxy mergers, which bring in new SMBHs. This showcases a possible mechanism contributing to the observed scatter of the $M_{\bullet}$-$σ_{\star}$ relation. Finally, we show that Pulsar Timing Arrays and LISA would be able to detect parts of the GW signals from the SMBH mergers that occur during the $\sim 4\,\mathrm{Gyr}$ time span simulated with KETJU.

preprint2022arXiv

The Metallicity Distribution Function in Outer Halo Fields of Simulated Elliptical Galaxies Compared to Observations of NGC 5128

Stellar metallicity distribution functions (MDF) have been measured for resolved stellar populations in the outer halos of many galaxies in nearby groups. Among them, the MDF of NGC 5128, the central giant elliptical in the Centaurus group, provides essential constraints for theories of massive galaxy formation and hierarchical assembly. To investigate the formation and chemical evolution history of the outer halo of giant elliptical galaxies, we examine the chemical properties of three zoom-in high resolution cosmological hydrodynamical simulations of an NGC 5128-like giant elliptical galaxy and compare their outer halo MDFs to the observed one of NGC 5128. Even though the simulated galaxies have different merging histories and age distributions, all predicted MDFs are in good qualitative agreement with the observed one. The median metallicity of the simulated galaxies is on average $\rm [M/H]=-0.41 \pm 0.06$ compared to the observed value of $\rm [M/H]=-0.38 \pm 0.02$ for NGC 5128, and the dispersion in metallicity is $\sim 0.77$ dex for both observed and simulated galaxies. We investigate the origin of the stars ending up in the outer halo field of simulated galaxies and show that most have an `accreted&#39; origin, formed in other small galaxies and later accreted in mergers. Only $\sim 15$ percent of the stars are formed `in situ&#39; within the main progenitor of galaxy and radially migrate outwards. We show that the contribution of metal-rich in situ stars is sub-dominant in the outer halos of our simulated galaxies, but can be prominent in the inner regions.

preprint2022arXiv

The origin of the [CII]-deficit in a simulated dwarf galaxies starburst

We present [CII] synthetic observations of smoothed particle hydrodynamics (SPH) simulations of a dwarf galaxy merger. The merging process varies the star-formation rate by more than three orders of magnitude. Several star clusters are formed, the feedback of which disperses and unbinds the dense gas through expanding HII regions and supernova (SN) explosions. For galaxies with properties similar to the modelled ones, we find that the [CII] emission remains optically thin throughout the merging process. We identify the Warm Neutral Medium ($3<\log T_{\rm gas}<4$ with $χ_{\rm HI}>2χ_{\rm H2}$) to be the primary source of [CII] emission ($\sim58\%$ contribution), although at stages when the HII regions are young and dense (during star cluster formation or SNe in the form of ionized bubbles) they can contribute $\gtrsim50\%$ to the total [CII] emission. We find that the [CII]/FIR ratio decreases due to thermal saturation of the [CII] emission caused by strong FUV radiation fields emitted by the massive star clusters, leading to a [CII]-deficit medium. We investigate the [CII]-SFR relation and find an approximately linear correlation which agrees well with observations, particularly those from the Dwarf Galaxy Survey. Our simulation reproduces the observed trends of [CII]/FIR versus $Σ_{\rm SFR}$ and $Σ_{\rm FIR}$, and it agrees well with the Kennicutt relation of SFR-FIR luminosity. We propose that local peaks of [CII] in resolved observations may provide evidence for ongoing massive cluster formation.

preprint2021arXiv

Black hole mergers in compact star clusters and massive black hole formation beyond the mass-gap

We present direct N-body simulations, carried out with Nbody6++GPU, of young and compact low metallicity star clusters with $1.1\times 10^5$ stars, a velocity dispersion of $\sim$ 10 $\mathrm{km\,s^{-1}}$, a half mass radius $R_h=0.6$ pc, and a binary fraction of $10\%$ including updated evolution models for stellar winds and pair-instability supernovae (PISNe). Within the first tens of megayears of evolution, each cluster hosts several black hole (BH) merger events which nearly cover the complete mass range of primary and secondary BH masses for current LIGO/Virgo/Kagra gravitational wave detections. The importance of gravitational recoil is estimated statistically. We present several possible formation paths of massive BHs above the assumed lower PISNe mass-gap limit ($45 M_\odot$) into the intermediate-mass BH (IMBH) regime ($> 100 M_\odot$) which include collisions of stars and BHs as well as the direct collapse of stellar merger remnants with low mass cores. The stellar evolution updates result in the early formation of higher mass stellar BHs than for the previous model. The resulting higher collision rates with massive stars support the rapid formation of massive BHs. For models assuming a high accretion efficiency for star-BH mergers, we present a first-generation formation scenario for GW190521-like events, a merger of two BHs in the PISN mass-gap, which is dominated by star-BH mergers. This IMBH formation path is independent of gravitational recoil and therefore conceivable in dense stellar systems with low escape velocities. One simulated cluster even forms an IMBH binary (153$M_\odot$,173$M_\odot$) which is expected to merge within a Hubble time.

preprint2021arXiv

EMERGE: Empirical predictions of galaxy merger rates since $z\sim6$

We explore the galaxy-galaxy merger rate with the empirical model for galaxy formation, Emerge. On average, we find that between $2$ per cent and $20$ per cent of massive galaxies ($\log_{10}(m_{*}/M_{\odot}) \geq 10.3$) will experience a major merger per Gyr. Our model predicts galaxy merger rates that do not scale as a power-law with redshift when selected by descendant stellar mass, and exhibit a clear stellar mass and mass-ratio dependence. Specifically, major mergers are more frequent at high masses and at low redshift. We show mergers are significant for the stellar mass growth of galaxies $\log_{10}(m_{*}/M_{\odot}) \gtrsim 11.0$. For the most massive galaxies major mergers dominate the accreted mass fraction, contributing as much as $90$ per cent of the total accreted stellar mass. We reinforce that these phenomena are a direct result of the stellar-to-halo mass relation, which results in massive galaxies having a higher likelihood of experiencing major mergers than low mass galaxies. Our model produces a galaxy pair fraction consistent with recent observations, exhibiting a form best described by a power-law exponential function. Translating these pair fractions into merger rates results in an inaccurate prediction compared to the model intrinsic values when using published observation timescales. We find the pair fraction can be well mapped to the intrinsic merger rate by adopting an observation timescale that decreases linearly with redshift as $T_{\mathrm{obs}} = -0.36(1+z)+2.39$ [Gyr], assuming all observed pairs merge by $z=0$.

preprint2021arXiv

LYRA I: Simulating the multi-phase ISM of a dwarf galaxy with variable energy supernovae from individual stars

We introduce the LYRA project, a new high resolution galaxy formation model built within the framework of the cosmological hydro-dynamical moving mesh code AREPO. The model resolves the multi-phase interstellar medium down to 10 K. It forms individual stars sampled from the initial mass function (IMF), and tracks their lifetimes and death pathways individually. Single supernova (SN) blast waves with variable energy are followed within the hydrodynamic calculation to interact with the surrounding interstellar medium (ISM). In this paper, we present the methods and apply the model to a $10^{10} M_{\odot}$ isolated halo. We demonstrate that the majority of supernovae are Sedov-resolved at our fiducial gas mass resolution of $4 M_{\odot}$. We show that our SN feedback prescription self-consistently produces a hot phase within the ISM that drives significant outflows, reduces the gas density and suppresses star formation. Clustered SN play a major role in enhancing the effectiveness of feedback, because the majority of explosions occur in low density material. Accounting for variable SN energy allows the feedback to respond directly to stellar evolution. We show that the ISM is sensitive to the spatially distributed energy deposition. It strongly affects the outflow behaviour, reducing the mass loading by a factor of 2-3, thus allowing the galaxy to retain a higher fraction of mass and metals. LYRA makes it possible to use a comprehensive multi-physics ISM model directly in cosmological (zoom) simulations of dwarf and higher mass galaxies.

preprint2020arXiv

Formation channels of slowly rotating early-type galaxies

We study the evidence for a diversity of formation processes in early-type galaxies by presenting the first complete volume-limited sample of slow rotators with both integral-field kinematics from the ATLAS3D Project and high spatial resolution photometry from the Hubble Space Telescope. Analysing the nuclear surface brightness profiles of 12 newly imaged slow rotators, we classify their light profiles as core-less, and place an upper limit to the core size of ~10 pc. Considering the full magnitude and volume-limited ATLAS3D sample, we correlate the presence or lack of cores with stellar kinematics, including the proxy for the stellar angular momentum and the velocity dispersion within one half-light radius, stellar mass, stellar age, $α$-element abundance, and age and metallicity gradients. More than half of the slow rotators have core-less light profiles, and they are all less massive than $10^{11}$ Msun. Core-less slow rotators show evidence for counter-rotating flattened structures, have steeper metallicity gradients, and a larger dispersion of gradient values than core slow rotators. Our results suggest that core and core-less slow rotators have different assembly processes, where the former are the relics of massive dissipation-less merging in the presence of central supermassive black holes. Formation processes of core-less slow rotators are consistent with accretion of counter-rotating gas or gas-rich mergers of special orbital configurations, which lower the final net angular momentum of stars, but support star formation. We also highlight core fast rotators as galaxies that share properties of core slow rotators and core-less slow rotators. Formation processes similar to those for core-less slow rotators can be invoked to explain the assembly of core fast rotators, with the distinction that these processes form or preserve cores.[Abridged]

preprint2020arXiv

Hot phase generation by supernovae: resolution, chemistry and thermal conduction

Supernovae (SN) generate hot gas in the interstellar medium (ISM), help setting the ISM structure and support the driving of outflows. It is important to resolve the hot gas generation for galaxy formation simulations at solar mass and sub-parsec resolution which realise individual supernova (SN) explosions with ambient densities varying by several orders of magnitude in a realistic multi-phase ISM. We test resolution requirements by simulating SN blast waves at three metallicities ($Z = 0.01, 0.1$ and $1 Z_{\odot}$), six densities and their respective equilibrium chemical compositions ($n=0.001$ cm$^{-3}$ - $100$ cm$^{-3}$), and four mass resolutions ($0.1$ - $100$ M$_{\odot}$), in three dimensions. We include non-equilibrium cooling and chemistry, a homogenous interstellar radiation field, and shielding with a modern pressure-energy smoothed particle hydrodynamics (SPH) method including isotropic thermal conduction and a meshless-finite-mass (MFM) solver. We find stronger resolution requirements for chemistry and hot phase generation than for momentum generation. While at $10$ M$_{\odot}$ the radial momenta at the end of the Sedov phase start converging, the hot phase generation and chemistry require higher resolutions to represent the neutral to ionised hydrogen fraction at the end of the Sedov phase correctly. Thermal conduction typically reduces the hot phase by $0.2$ dex and has little impact on the chemical composition. In general, our $1$, and $0.1$ M$_{\odot}$ results agree well with previous numerical and analytic estimates. We conclude that for the thermal energy injection SN model presented here resolutions higher than $10$ M$_{\odot}$ are required to model the chemistry, momentum and hot phase generation in a multi-phase ISM.

preprint2020arXiv

Lyman-αabsorption beyond the disk of simulated spiral galaxies

We present an analysis of the origin and properties of the circum-galactic medium (CGM) in a suite of 11 cosmological zoom simulations resembling present day spiral galaxies. On average the galaxies retain about 50\% of the cosmic fraction in baryons, almost equally divided into disc (interstellar medium) gas, cool CGM gas and warm-hot CGM gas. At radii smaller than 50 kpc the CGM is dominated by recycled warm-hot gas injected from the central galaxy, while at larger radii it is dominated by cool gas accreted onto the halo. The recycled gas typically accounts for one-third of the CGM mass. We introduce the novel publicly available analysis tool \textsc{pygad} to compute ion abundances and mock absorption spectra. For Lyman-$α$ absorption we find good agreement of the simulated equivalent width (EW) distribution and observations out to large radii. Disc galaxies with quiescent assembly histories show significantly more absorption along the disc major axis. By comparing the EW and HI column densities we find that CGM Lyman-$α$ absorbers are best represented by an effective line-width $b\approx 50 - 70$ km s$^{-1}$ that increases mildly with halo mass, larger than typically assumed.

preprint2020arXiv

MSTAR -- a fast parallelised algorithmically regularised integrator with minimum spanning tree coordinates

We present the novel algorithmically regularised integration method MSTAR for high accuracy ($|ΔE/E| \gtrsim 10^{-14}$) integrations of N-body systems using minimum spanning tree coordinates. The two-fold parallelisation of the $\mathcal{O}(N_\mathrm{part}^2)$ force loops and the substep divisions of the extrapolation method allows for a parallel scaling up to $N_\mathrm{CPU} = 0.2 \times N_\mathrm{part}$. The efficient parallel scaling of MSTAR makes the accurate integration of much larger particle numbers possible compared to the traditional algorithmic regularisation chain (AR-CHAIN) methods, e.g. $N_\mathrm{part} = 5000$ particles on $400$ CPUs for $1$ Gyr in a few weeks of wall-clock time. We present applications of MSTAR on few particle systems, studying the Kozai mechanism and N-body systems like star clusters with up to $N_\mathrm{part} =10^4$ particles. Combined with a tree or a fast multipole based integrator the high performance of MSTAR removes a major computational bottleneck in simulations with regularised subsystems. It will enable the next generation galactic-scale simulations with up to $10^9$ stellar particles (e.g. $m_\star = 100 M_\odot$ for a $M_\star = 10^{11} M_\odot$ galaxy) including accurate collisional dynamics in the vicinity of nuclear supermassive black holes.

preprint2020arXiv

The GRIFFIN project -- Formation of star clusters with individual massive stars in a simulated dwarf galaxy starburst

We describe a population of young star clusters (SCs) formed in a hydrodynamical simulation of a gas-rich dwarf galaxy merger resolved with individual massive stars at sub-parsec spatial resolution. The simulation is part of the GRIFFIN (Galaxy Realizations Including Feedback From INdividual massive stars) project. The star formation environment during the simulation spans seven orders of magnitude in gas surface density and thermal pressure, and the global star formation rate surface density ($Σ_\mathrm{SFR}$) varies by more than three orders of magnitude during the simulation. Young SCs more massive than $M_{\mathrm{*,cl}}\sim 10^{2.5}\,M_{\odot}$ form along a mass function with a power-law index $α\sim-1.7$ ($α\sim-2$ for $M_{\mathrm{*,cl}}\gtrsim10^{3}\,M_{\odot}$) at all merger phases, while the normalization and the highest SC masses (up to $\sim 10^6 M_{\odot}$) correlate with $Σ_\mathrm{SFR}$. The cluster formation efficiency varies from $Γ\sim20\%$ in early merger phases to $Γ\sim80\%$ at the peak of the starburst and is compared to observations and model predictions. The massive SCs ($\gtrsim10^4\,M_{\odot}$) have sizes and mean surface densities similar to observed young massive SCs. Simulated lower mass clusters appear slightly more concentrated than observed. All SCs form on timescales of a few Myr and lose their gas rapidly resulting in typical stellar age spreads between $σ\sim0.1-2$ Myr ($1σ$), consistent with observations. The age spreads increase with cluster mass, with the most massive cluster ($\sim10^6\, M_{\odot}$) reaching a spread of $5\, \mathrm{Myr}$ once its hierarchical formation finishes. Our study shows that it is now feasible to investigate the SC population of entire galaxies with novel high-resolution numerical simulations.