Search

The Voids-within-Voids-within-Voids (VVV) project used dark-matter-only simulations to study the abundance and structure of dark matter haloes over the full mass range populated in the standard $\Lambda\mathrm{CDM}$ cosmology. Here we explore how baryonic effects modify these results for $z=0$ halo masses in the range $10^4$ to $10^7~\mathrm{M_\odot}$, below the threshold for galaxy formation. Our main study focuses on three simulations from identical initial conditions at $z=127$, one following dark matter only, one including non-radiative gas, and one additionally including the baryonic physics relevant in this halo mass range (cooling and photoheating). In the non-radiative simulation, above $10^{5.5}~\mathrm{M_\odot}$, halo abundance and internal structure are very similar to the dark-matter-only simulation, and the baryon to dark matter ratio is everywhere close to the cosmic value. At lower mass, this ratio drops and haloes are less concentrated and less massive in the non-radiative case. Test simulations at higher resolution show this to be mainly a resolution effect; the expected drop in baryon content due to residual pressure effects only becomes substantial for $z=0$ haloes below $\sim 10^{2.7}~\mathrm{M_\odot}$. However, gas is heated by reionization at $z=6$ in our ``full physics'' run, and this results in almost complete expulsion of gas from all haloes in our simulated mass range. This suppresses the halo mass function by $\sim 30 \%$, lowers halo concentration, and consequently weakens the dark matter annihilation signal by $\sim 40-60 \%$., Comment: 12+2 pages, 9 figures

The goal of this study is to conduct a timely analysis of the high-redshift star-forming galaxy populations, which will be informative in designing next-generation experiments and their extragalactic targets. We use the hydrodynamical simulation MillenniumTNG (MTNG) to model Lyman-alpha Emitting (LAE) galaxies to extract key properties such as their clustering and occupation statistics. We define LAEs through an empirical relation between star formation rate (SFR) and Lyman-alpha flux. We also explore two other definitions, finding that imposing an additional cut on the maximum stellar mass of the galaxy sample, which approximates the effect of a low escape fraction at high halo mass, leads to a 5-10\% decrease of the linear bias of the population. As expected, we find that the HOD mass parameters rapidly decrease with increasing number density. Additionally, the HOD parameter $\sigma$ also decreases with number density, implying that the SFR-halo mass relationship becomes tighter for low-luminosity objects. Surprisingly, the non-linear clustering, estimated by the parameter $r_0$, is fixed at fixed number density, whereas the linear bias parameter varies with redshift as $b(z) \propto (1 + z)$, suggesting that our LAE samples are relatively stable and long-lived. Finally, we study the amount of galaxy assembly bias present at $z = 2, \ 3$ and find that while at $z = 2$ it is roughly $\lesssim$10\%, at $z = 3$ it decreases significantly to $\lesssim$5\%. This suggests that assembly bias effects become less important at high $z$ likely due to the lower number of cumulative two-halo interactions (mergers, splashback, stripping, etc.). While our study is based on a single full-physics simulation, we expect our results to reflect the properties of LAEs in the Universe. We demonstrate that our findings are in good agreement with previous results using both observations and simulations., Comment: 19 pages, 10 figures, 1 table

Dwarf satellites in galaxy groups are distributed in an anisotropic and asymmetric manner, which is called the ``lopsided satellite distribution''. This lopsided signal has been observed not only in galaxy pairs but also in isolated systems. However, the physical origin of the lopsided signal in isolated systems is still unknown. In this work, we investigate this in the state-of-the-art hydrodynamical simulation of the MillenniumTNG Project by tracing each system back to high redshift. We find that the lopsided signal is dominated by satellites located in the outer regions of the halo and is also dominated by recently accreted satellites. The lopsided signal originates from the anisotropic accretion of galaxies from the surrounding large-scale structure and that, after accretion, the nonlinear evolution of satellites inside the dark-matter halo weakens the lopsidedness. The signal decreases as cosmic time passes because of a competition between anisotropic accretion and internal evolution within dark matter halos. Our findings provide a useful perspective for the study of galaxy evolution, especially for the origin of the spatial satellite galaxy distributions., Comment: 10 pages, 8 figures, 1 table. Submitted to MNRAS

For decades, the boundary of cosmic filaments have been a subject of debate. In this work, we determine the physically-motivated radii of filaments by constructing stacked galaxy number density profiles around the filament spines. We find that the slope of the profile changes with distance to the filament spine, reaching its minimum at approximately 1 Mpc at z = 0 in both state-of-the-art hydrodynamical simulations and observational data. This can be taken as the average value of the filament radius. Furthermore, we note that the average filament radius rapidly decreases from z = 4 to z = 1, and then slightly increases. Moreover, we find that the filament radius depends on the filament length, the distance from connected clusters, and the masses of the clusters. These results suggest a two-phase formation scenario of cosmic filaments. The filaments experience rapid contraction before z = 1, but their density distribution has remained roughly stable since then. The subsequent mass transport along the filaments to the connected clusters is likely to have contributed to the formation of the clusters themselves., Comment: 11 pages, 6 figures, submitted

The distribution of offsets between the brightest cluster galaxies of galaxy clusters and the centroid of their dark matter distributions is a promising probe of the underlying dark matter physics. In particular, since this distribution is sensitive to the shape of the potential in galaxy cluster cores, it constitutes a test of dark matter self-interaction on the largest mass scales in the universe. We examine these offsets in three suites of modern cosmological simulations; IllustrisTNG, MillenniumTNG and BAHAMAS. For clusters above $10^{14}\rm{M_\odot}$, we examine the dependence of the offset distribution on gravitational softening length, the method used to identify centroids, redshift, mass, baryonic physics, and establish the stability of our results with respect to various nuisance parameter choices. We find that offsets are overwhelmingly measured to be smaller than the minimum converged length scale in each simulation, with a median offset of $\sim 1\rm{kpc}$ in the highest resolution simulation considered, TNG300-1, which uses a gravitational softening length of $1.48\rm{kpc}$. We also find that centroids identified via source extraction on smoothed dark matter and stellar particle data are consistent with the potential minimum, but that observationally relevant methods sensitive to cluster strong gravitational lensing scales, or those using the the "light traces mass" approach, in this context meaning gas is used as a tracer for the potential, can overestimate offsets by factors of $\sim10$ and $\sim30$, respectively. This has the potential to reduce tensions with existing offset measurements which have served as evidence for a nonzero dark matter self-interaction cross section., Comment: 15 pages, 6 figures, comments welcome

The statistics of thermal gas pressure are a new and promising probe of cosmology and astrophysics. The large-scale cross-correlation between galaxies and the thermal Sunyaev-Zeldovich effect gives the bias-weighted mean electron pressure, $\langle b_\mathrm{h}P_e\rangle$. In this paper, we show that $\langle b_\mathrm{h}P_e\rangle$ is sensitive to the amplitude of fluctuations in matter density, for example $\langle b_\mathrm{h}P_e\rangle\propto \left(\sigma_8\Omega_\mathrm{m}^{0.81}h^{0.67}\right)^{3.14}$ at redshift $z=0$. We find that at $z<0.5$ the observed $\langle b_\mathrm{h}P_e\rangle$ is smaller than that predicted by the state-of-the-art hydrodynamical simulations of galaxy formation, MillenniumTNG, by a factor of $0.93$. This can be explained by a lower value of $\sigma_8$ and $\Omega_\mathrm{m}$, similar to the so-called "$S_8$ tension'' seen in the gravitational lensing effect, although the influence of astrophysics cannot be completely excluded. The difference between Magneticum and MillenniumTNG at $z<2$ is small, indicating that the difference in the galaxy formation models used by these simulations has little impact on $\langle b_\mathrm{h}P_e\rangle$ at this redshift range. At higher $z$, we find that both simulations are in a modest tension with the existing upper bounds on $\langle b_\mathrm{h}P_e\rangle$. We also find a significant difference between these simulations there, which we attribute to a larger sensitivity to the galaxy formation models in the high redshift regime. Therefore, more precise measurements of $\langle b_\mathrm{h}P_e\rangle$ at all redshifts will provide a new test of our understanding of cosmology and galaxy formation., Comment: 16 pages, 10 figures

We present a study of the evolution of cosmic filaments across redshift with an emphasis on some important properties: filament lengths, growth rates, and radial profiles of galaxy densities. Following an observation-driven approach, we build cosmic filament catalogues at z=0,1,2,3, and 4 from the galaxy distributions of the large hydro-dynamical run of the MilleniumTNG project. We employ the extensively used DisPerSE cosmic web finder code, for which we provide a user-friendly guide, including the details of a physics-driven calibration procedure, with the hope of helping future users. We perform the first statistical measurements of the evolution of connectivity in a large-scale simulation, finding that the connectivity of cosmic nodes (defined as the number of filaments attached) globally decreases from early to late times. The study of cosmic filaments in proper coordinates reveals that filaments grow in length and radial extent, as expected from large-scale structures in an expanding Universe. But the most interesting results arise once the Hubble flow is factored out. We find remarkably stable comoving filament length functions and over-density profiles, showing only little evolution of the total population of filaments in the past ~12.25 Gyrs. However, by tracking the spatial evolution of individual structures, we demonstrate that filaments of different lengths actually follow different evolutionary paths. While short filaments preferentially contract, long filaments expand along their longitudinal direction with growth rates that are the highest in the early, matter-dominated Universe. Filament diversity at fixed redshift is also shown by the different (~$5 \sigma$) density values between the shortest and longest filaments. Our results hint that cosmic filaments can be used as additional probes for dark energy, but further theoretical work is still needed., Comment: 17 pages, accepted in Astronomy & Astrophysics

Wide-area deep imaging surveys have discovered large numbers of extremely low surface brightness dwarf galaxies, which challenge galaxy formation theory and, potentially, offer new constraints on the nature of dark matter. Here we discuss one as-yet unexplored formation mechanism that may account for a fraction of low surface brightness dwarfs. We call this the `ghost galaxy' scenario. In this scenario, inefficient radiative cooling prevents star formation in the `main branch' of the merger tree of a low mass dark matter halo, such that almost all its stellar mass is acquired through mergers with less massive (but nevertheless star-forming) progenitors. Present-day systems formed in this way would be `ghostly' isolated stellar halos with no central galaxy. We use merger trees based on the Extended Press-Schechter formalism and the COCO cosmological N-body simulation to demonstrate that mass assembly histories of this kind can occur for low-mass halos in Lambda-CDM, but they are rare. They are most probable in isolated halos of present-day mass ~4x10^9 M_sun, occurring for ~5 per cent of all halos of that mass under standard assumptions about the timing and effect of cosmic reionization. The stellar masses of star-forming progenitors in these systems are highly uncertain; abundance-matching arguments imply a bimodal present-day mass function having a brighter population (median M_star ~3x10^6 M_sun) consistent with the tail of the observed luminosity function of ultra-diffuse galaxies. This suggests observable analogues of these systems may await discovery. We find that a stronger ionizing background (globally or locally) produces brighter and more extended ghost galaxies., Comment: 19 pages, 13 figures, ApJ in press

We study weak gravitational lensing convergence maps produced from the MillenniumTNG (MTNG) simulations by direct projection of the mass distribution on the past backwards lightcone of a fiducial observer. We explore the lensing maps over a large dynamic range in simulation mass and angular resolution, allowing us to establish a clear assessment of numerical convergence. By comparing full physics hydrodynamical simulations with corresponding dark-matter-only runs we quantify the impact of baryonic physics on the most important weak lensing statistics. Likewise, we predict the impact of massive neutrinos reliably far into the non-linear regime. We also demonstrate that the "fixed & paired" variance suppression technique increases the statistical robustness of the simulation predictions on large scales not only for time slices but also for continuously output lightcone data. We find that both baryonic and neutrino effects substantially impact weak lensing shear measurements, with the latter dominating over the former on large angular scales. Thus, both effects must explicitly be included to obtain sufficiently accurate predictions for stage IV lensing surveys. Reassuringly, our results agree accurately with other simulation results where available, supporting the promise of simulation modelling for precision cosmology far into the non-linear regime., Comment: 14 pages, 10 figures, comments welcome

In the coming years, Sunyaev-Zel'dovich (SZ) measurements can dramatically improve our understanding of the Intergalactic Medium (IGM) and the role of feedback processes on galaxy formation, allowing us to calibrate important astrophysical systematics in cosmological constraints from weak lensing galaxy clustering surveys. However, the signal is only measured in a two-dimensional projection, and its correct interpretation relies on understanding the connection between observable quantities and the underlying intrinsic properties of the gas, in addition to the relation between the gas and the underlying matter distribution. One way to address these challenges is through the use of hydrodynamical simulations such as the high-resolution, large-volume MillenniumTNG suite. We find that measurements of the optical depth, $\tau$, and the Compton-y parameter, $Y$, receive large line-of-sight contributions which can be removed effectively by applying a Compensated Aperture Photometry (CAP) filter. In contrast with other $\tau$ probes (e.g., X-rays and Fast Radio Bursts), the kSZ-inferred $\tau$ receives most of its signal from a confined cylindrical region around the halo due to the velocity decorrelation along the line-of-sight. Additionally, we perform fits to the $Y-M$ and $\tau-M$ scaling relations and report best-fit parameters adopting the smoothly broken power law (SBPL) formalism. We note that subgrid physics modeling can broaden the error bar on these by 30\% for intermediate-mass halos ($\sim$$10^{13} \, {\rm M}_{\odot}$). The scatter of the scaling relations can be captured by an intrinsic dependence on concentration, and an extrinsic dependence on tidal shear. Finally, we comment on the effect of using galaxies rather than halos in real observations, which can bias the inferred SZ profiles by $\sim$20\% for $L_\ast$-galaxies., Comment: 14 pages, 6 figures

The intrinsic alignment (IA) of observed galaxy shapes with the underlying cosmic web is a source of contamination in weak lensing surveys. Sensitive methods to identify the IA signal will therefore need to be included in the upcoming weak lensing analysis pipelines. Hydrodynamical cosmological simulations allow us to directly measure the intrinsic ellipticities of galaxies and thus provide a powerful approach to predict and understand the IA signal. Here we employ the novel, large-volume hydrodynamical simulation MTNG740, a product of the MillenniumTNG (MTNG) project, to study the IA of galaxies. We measure the projected correlation functions between the intrinsic shape/shear of galaxies and various tracers of large-scale structure, $w_{+g},\ w_{+m},\ w_{++}$ over the radial range $r_{\rm p} \in [0.02 , 200]\,h^{-1}{\rm Mpc}$ and at redshifts $z=0.0$, $0.5$ and $1.0$. We detect significant signal-to-noise IA signals with the density field for both elliptical and spiral galaxies. We also find significant intrinsic shear-shear correlations for ellipticals. We further examine correlations of the intrinsic shape of galaxies with the local tidal field. Here we find a significant IA signal for elliptical galaxies assuming a linear model. We also detect a weak IA signal for spiral galaxies under a quadratic tidal torquing model. Lastly, we measure the alignment between central galaxies and their host dark-matter halos, finding small to moderate misalignments between their principal axes that decline with halo mass., Comment: 16 pages, 14 figures

We take advantage of the unprecedented dynamical range provided by the "Cosmic-Zoom" project to study the mass accretion history (MAH) of present-day dark matter haloes over the entire mass range present in the $\Lambda$CDM paradigm when the dark matter is made of weakly interacting massive particles of mass $100\ \mathrm{GeV}$. In particular, we complement previous studies by exploring the MAHs of haloes with mass from $10^8\ h^{-1}\mathrm{M_{\odot}}$ down to Earth mass, $10^{-6}\ h^{-1}\mathrm{M_{\odot}}$. The formation redshift of low-mass haloes anti-correlates weakly with mass, peaking at $z=3$ for haloes of mass $10^{-4}\ h^{-1}\mathrm{M_{\odot}}$. Even lower masses are affected by the free-streaming cutoff in the primordial spectrum of density fluctuations and form at lower redshift. We compare MAHs in our simulations with predictions from two analytical models based on the extended Press-Schechter theory (EPS), and three empirical models derived by fitting and extrapolating either results from cosmological $N$-body simulations or Monte Carlo realizations of halo growth. All models fit our simulations reasonably well over the mass range for which they were calibrated. While the empirical models match better for more massive haloes, $M>10^{10}\ h^{-1}\mathrm{M_{\odot}}$, the analytical models do better when extrapolated down to Earth mass. At the higher masses, we explore the correlation between local environment density and MAH, finding that biases are relatively weak, with typical MAHs for haloes in extremely low-density and in typical regions differing by less than $20$ percent at high redshift. We conclude that EPS theory predicts the hierarchical build-up of dark matter haloes quite well over the entire mass range.

We investigate cosmological structure formation in Fuzzy Dark Matter (FDM) with an attractive self-interaction (SI) with numerical simulations. Such a SI would arise if the FDM boson were an ultra-light axion, which has a strong CP symmetry-breaking scale (decay constant). Although weak, the attractive SI may be strong enough to counteract the quantum 'pressure' and alter structure formation. We find in our simulations that the SI can enhance small-scale structure formation, and soliton cores above a critical mass undergo a phase transition, transforming from dilute to dense solitons., Comment: 10 pages, 3 figures, submitted to mnras

In this series of papers we present an emulator-based halo model for the non-linear clustering of galaxies in modified gravity cosmologies. In the first paper, we present emulators for the following halo properties: the halo mass function, concentration-mass relation and halo-matter cross-correlation function. The emulators are trained on data extracted from the \textsc{FORGE} and \textsc{BRIDGE} suites of $N$-body simulations, respectively for two modified gravity (MG) theories: $f(R)$ gravity and the DGP model, varying three standard cosmological parameters $\Omega_{\mathrm{m0}}, H_0, \sigma_8$, and one MG parameter, either $\bar{f}_{R0}$ or $r_{\mathrm{c}}$. Our halo property emulators achieve an accuracy of $\lesssim 1\%$ on independent test data sets. We demonstrate that the emulators can be combined with a galaxy-halo connection prescription to accurately predict the galaxy-galaxy and galaxy-matter correlation functions using the halo model framework., Comment: 18 pages, 15 figures

We use the Voids-within-Voids-within-Voids (VVV) simulations, a suite of successive nested N-body simulations with extremely high resolution (denoted, from low to high resolution, by L0 to L7), to test the Press-Schechter (PS), Sheth-Tormen (ST), and extended Press-Schechter (EPS) formulae for the halo abundance over the entire mass range, from mini-haloes of $10^{-6}\ \mathrm{M_\odot}$, to cluster haloes of $10^{15}\ \mathrm{M_\odot}$, at different redshifts, from $z=30$ to the present. We find that at $z=0$ and $z=2$, ST best reproduces the results of L0, which has the mean cosmic density (overdensity $\delta=0$), at $10^{11-15} ~\mathrm{M_\odot}$. The higher resolution levels (L1-L7) are biased underdense regions ($\delta<-0.6$). The EPS formalism takes this into account since it gives the mass function of a region conditioned, in this case, on having a given underdensity. EPS provides good matches to these higher levels, with deviations $\lesssim 20\%$, at $10^{-6-12.5} ~\mathrm{M_\odot}$. At $z \sim 7-15$, the ST predictions for L0 and the EPS for L1-L7 show somewhat larger deviations from the simulation results. However, at even higher redshifts, $z \sim 30$, EPS fits the simulations well again. We confirm our results by picking more subvolumes from the L0 simulation, finding that our conclusions depend only weakly on the size and overdensity of the region. The good agreement of EPS with the higher-level simulations implies that PS (or ST) gives an accurate description of the total halo mass function in representative regions of the universe., Comment: 10 pages, 5 figures (additional 2 figures in the appendix)

Observations of the high redshift universe provide a promising avenue for constraining the nature of the dark matter (DM). This will be even more true following the now successful launch of the James Webb Space Telescope (JWST). We run cosmological simulations of galaxy formation as a part of the Effective Theory of Structure Formation (ETHOS) project to compare the properties of high redshift galaxies in Cold (CDM) and alternative DM models which have varying relativistic coupling and self-interaction strengths. Phenomenologically, the interacting DM scenarios result in a cutoff in the linear power spectrum on small-scales, followed by a series of "dark acoustic oscillations" (DAOs). We find that DM interactions suppress the abundance of galaxies below $M_\star \sim 10^8\, M_\odot$ for the models considered. The cutoff in the linear power spectrum generally causes a delay in structure formation relative to CDM. Objects in ETHOS that end up at the same final masses as their CDM counterparts are characterised by a more vigorous phase of early star formation. While galaxies with $M_\star \lesssim 10^6\,{M_\odot}$ make up more than 60 per cent of star formation in CDM at $z\approx 10$, they contribute only about half the star formation density in ETHOS. These differences in star formation diminish with decreasing redshift. We find that the effects of DM self-interactions are negligible compared to effects of relativistic coupling (i.e. the effective initial conditions for galaxy formation) in all properties of the galaxy population we examine. Finally, we show that the clustering strength of galaxies at high redshifts depends sensitively on DM physics, although these differences are manifest on scales that may be too small to be measurable by JWST., Comment: 15 pages, 11 figures, 3 Tables. Comments are welcome. Submitted to MNRAS; fixed a bug in figure 4 and type in the explanation, general results are unchanged

We employ the hydrodynamical simulation IllustrisTNG to inform the galaxy-halo connection of the Luminous Red Galaxy (LRG) and Emission Line Galaxy (ELG) samples of the Dark Energy Spectroscopic Instrument (DESI) survey at redshift z ~ 0.8. Specifically, we model the galaxy colors of IllustrisTNG and apply sliding DESI color-magnitude cuts, matching the DESI target densities. We study the halo occupation distribution model (HOD) of the selected samples by matching them to their corresponding dark matter halos in the IllustrisTNG dark matter run. We find the HOD of both the LRG and ELG samples to be consistent with their respective baseline models, but also we find important deviations from common assumptions about the satellite distribution, velocity bias, and galaxy secondary biases. We identify strong evidence for concentration-based and environment-based occupational variance in both samples, an effect known as "galaxy assembly bias". The central and satellite galaxies have distinct dependencies on secondary halo properties, showing that centrals and satellites have distinct evolutionary trajectories and should be modelled separately. These results serve to inform the necessary complexities in modeling galaxy-halo connection for DESI analyses and also prepare for building high-fidelity mock galaxies. Finally, we present a shuffling-based clustering analysis that reveals a 10-15% excess in the LRG clustering of modest statistical significance due to secondary galaxy biases. We also find a similar excess signature for the ELGs, but with much lower statistical significance. When a larger hydrodynamical simulation volume becomes available, we expect our analysis pipeline to pinpoint the exact sources of such excess clustering signatures., Comment: Accepted to MNRAS, comments welcome

The detection of Intermediate-Mass Black Holes (IMBHs) in dwarf galaxies is crucial to closing the gap in the wide mass distribution of black holes ($\sim 3 \, \rm M_{\odot}$ to $\sim 5 \times 10^{10} \, \rm M_{\odot}$). IMBHs originally located at the center of dwarfs that later collide with the Milky Way (MW) could be wandering, undetected, in our Galaxy. We used TNG50, the highest-resolution run of the IllustrisTNG project, to study the kinematics and dynamics of star clusters, in the appropriate mass range, acting as IMBH proxies in a MW analog galaxy. We showed that $\sim 87\%$ of our studied IMBHs drift inward. The radial velocity of these sinking IMBHs has a median magnitude of $\sim 0.44 \, \rm ckpc \, h^{-1} \, Gyr^{-1}$ and no dependence on the black hole mass. The central $1 \, \rm ckpc \, h^{-1}$ has the highest number density of IMBHs in the galaxy. A physical toy model with linear drag forces was developed to explain the orbital circularization with time. These findings constrain the spatial distribution of IMBHs, suggesting that future searches should focus on the central regions of the Galaxy. Additionally, we found that the 3D velocity distribution of IMBHs with respect to the galactic center has a mean of $\sim 180 \, \rm km \, s^{-1}$ and larger variance with decreasing radius. Remarkably, the velocity distribution relative to the local gas shows significantly lower values, with a mean of $\sim 88 \, \rm km \, s^{-1}$. These results are instrumental for predicting the accretion and radiation properties of IMBHs, facilitating their detection with future surveys., Comment: Accepted for publication in MNRAS. This is the final version of the manuscript. 10 pages, 12 figures

We describe best practices for early career researchers to develop an effective and fair approach to refereeing manuscripts.

We show that gravitational waves have the potential to unravel the microphysical properties of dark matter due to the dependence of the binary black hole merger rate on cosmic structure formation, which is itself highly dependent on the dark matter scenario. In particular, we demonstrate that suppression of small-scale structure -- such as that caused by interacting, warm, or fuzzy dark matter -- leads to a significant reduction in the rate of binary black hole mergers at redshifts $z\gtrsim5$. This shows that future gravitational-wave observations will provide a new probe of physics beyond the $\Lambda$CDM model., Comment: 15 pages, 12 figures, matches version published in PRD

We describe best practices for early career researchers to develop an effective and fair approach to refereeing manuscripts.

Cosmological simulations are an important theoretical pillar for understanding nonlinear structure formation in our Universe and for relating it to observations on large scales. In several papers, we introduce our MillenniumTNG (MTNG) project that provides a comprehensive set of high-resolution, large volume simulations of cosmic structure formation aiming to better understand physical processes on large scales and to help interpreting upcoming large-scale galaxy surveys. We here focus on the full physics box MTNG740 that computes a volume of $(740\,\mathrm{Mpc})^3$ with a baryonic mass resolution of $3.1\times~10^7\,\mathrm{M_\odot}$ using \textsc{arepo} with $80.6$~billion cells and the IllustrisTNG galaxy formation model. We verify that the galaxy properties produced by MTNG740 are consistent with the TNG simulations, including more recent observations. We focus on galaxy clusters and analyse cluster scaling relations and radial profiles. We show that both are broadly consistent with various observational constraints. We demonstrate that the SZ-signal on a deep lightcone is consistent with Planck limits. Finally, we compare MTNG740 clusters with galaxy clusters found in Planck and the SDSS-8 RedMaPPer richness catalogue in observational space, finding very good agreement as well. However, {\it simultaneously} matching cluster masses, richness, and Compton-$y$ requires us to assume that the SZ mass estimates for Planck clusters are underestimated by $0.2$~dex on average. Thanks to its unprecedented volume for a high-resolution hydrodynamical calculation, the MTNG740 simulation offers rich possibilities to study baryons in galaxies, galaxy clusters, and in large scale structure, and in particular their impact on upcoming large cosmological surveys., Comment: 18 pages, 14 figures, accepted for publication by MNRAS, comments welcome

Luminous red galaxies (LRGs) and blue star-forming emission-line galaxies (ELGs) are key tracers of large-scale structure used by cosmological surveys. Theoretical predictions for such data are often done via simplistic models for the galaxy-halo connection. In this work, we use the large, high-fidelity hydrodynamical simulation of the MillenniumTNG project (MTNG) to inform a new phenomenological approach for obtaining an accurate and flexible galaxy-halo model on small scales. Our aim is to study LRGs and ELGs at two distinct epochs, $z = 1$ and $z = 0$, and recover their clustering down to very small scales, $r \sim 0.1 \ {\rm Mpc}/h$, i.e. the one-halo regime, while a companion paper extends this to a two-halo model for larger distances. The occupation statistics of ELGs in MTNG inform us that: (1) the satellite occupations exhibit a slightly super-Poisson distribution, contrary to commonly made assumptions, and (2) that haloes containing at least one ELG satellite are twice as likely to host a central ELG. We propose simple recipes for modeling these effects, each of which calls for the addition of a single free parameter to simpler halo occupation models. To construct a reliable satellite population model, we explore the LRG and ELG satellite radial and velocity distributions and compare them with those of subhalos and particles in the simulation. We find that ELGs are anisotropically distributed within halos, which together with our occupation results provides strong evidence for cooperative galaxy formation (manifesting itself as one-halo galaxy conformity); i.e.~galaxies with similar properties form in close proximity to each other. Our refined galaxy-halo model represents a useful improvement of commonly used analysis tools and thus can be of help to increase the constraining power of large-scale structure surveys., Comment: 14 pages, 10 figures; submitted to MNRAS

Modern redshift surveys are tasked with mapping out the galaxy distribution over enormous distance scales. Existing hydrodynamical simulations, however, do not reach the volumes needed to match upcoming surveys. We present results for the clustering of galaxies using a new, large volume hydrodynamical simulation as part of the MillenniumTNG (MTNG) project. With a computational volume that is $\approx15$ times larger than the next largest such simulation currently available, we show that MTNG is able to accurately reproduce the observed clustering of galaxies as a function of stellar mass. When separated by colour, there are some discrepancies with respect to the observed population, which can be attributed to the quenching of satellite galaxies in our model. We combine MTNG galaxies with those generated using a semi-analytic model to emulate the sample selection of luminous red galaxies (LRGs) and emission line galaxies (ELGs), and show that although the bias of these populations is approximately (but not exactly) constant on scales larger than $\approx10$ Mpc, there is significant scale-dependent bias on smaller scales. The amplitude of this effect varies between the two galaxy types, and also between the semi-analytic model and MTNG. We show that this is related to the distribution of haloes hosting LRGs and ELGs. Using mock SDSS-like catalogues generated on MTNG lightcones, we demonstrate the existence of prominent baryonic acoustic features in the large-scale galaxy clustering. We also demonstrate the presence of realistic redshift space distortions in our mocks, finding excellent agreement with the multipoles of the redshift-space clustering measured in SDSS data., Comment: 15 pages, 11 figures. Accepted for publication in MNRAS

We introduce a novel technique for constraining cosmological parameters and galaxy assembly bias using non-linear redshift-space clustering of galaxies. We scale cosmological N-body simulations and insert galaxies with the SubHalo Abundance Matching extended (SHAMe) empirical model to generate over 175,000 clustering measurements spanning all relevant cosmological and SHAMe parameter values. We then build an emulator capable of reproducing the projected galaxy correlation function at the monopole, quadrupole and hexadecapole level for separations between $0.1\,h^{-1}{\rm Mpc}$ and $25\,h^{-1}{\rm Mpc}$. We test this approach by using the emulator and Monte Carlo Markov Chain (MCMC) inference to jointly estimate cosmology and assembly bias parameters both for the MTNG740 hydrodynamic simulation and for a semi-analytical galaxy formation model (SAM) built on the MTNG740-DM dark matter-only simulation, obtaining unbiased results for all cosmological parameters. For instance, for MTNG740 and a galaxy number density of $n\sim 0.01 h^{3}{\rm Mpc}^{-3}$, we obtain $\sigma_{8}=0.799^{+0.039}_{-0.044}$ ($\sigma_{8,{\rm MTNG}} =$ 0.8159), and $\Omega_\mathrm{M}h^2= 0.138^{+ 0.025}_{- 0.018}$ ($\Omega_{\mathrm{M}} h^2_{\rm MTNG} =$ 0.142). For fixed Hubble parameter ($h$), the constraint becomes $\Omega_\mathrm{M}h^2= 0.137^{+ 0.011}_{- 0.012}$. Our method performs similarly well for the SAM and for other tested sample densities. We almost always recover the true amount of galaxy assembly bias within one sigma. The best constraints are obtained when scales smaller than $2\,h^{-1}{\rm Mpc}$ are included, as well as when at least the projected correlation function and the monopole are incorporated. These methods offer a powerful way to constrain cosmological parameters using galaxy surveys., Comment: 18 pages, 12 figures, accepted by MNRAS

Cosmological inference with large galaxy surveys requires theoretical models that combine precise predictions for large-scale structure with robust and flexible galaxy formation modelling throughout a sufficiently large cosmic volume. Here, we introduce the MillenniumTNG (MTNG) project which combines the hydrodynamical galaxy formation model of IllustrisTNG with the large volume of the Millennium simulation. Our largest hydrodynamic simulation, covering (500 Mpc/h)^3 = (740 Mpc)^3, is complemented by a suite of dark-matter-only simulations with up to 4320^3 dark matter particles (a mass resolution of 1.32 x 10^8 Msun/h) using the fixed-and-paired technique to reduce large-scale cosmic variance. The hydro simulation adds 4320^3 gas cells, achieving a baryonic mass resolution of 2 x 10^7 Msun/h. High time-resolution merger trees and direct lightcone outputs facilitate the construction of a new generation of semi-analytic galaxy formation models that can be calibrated against both the hydro simulation and observation, and then applied to even larger volumes - MTNG includes a flagship simulation with 1.1 trillion dark matter particles and massive neutrinos in a volume of (3000 Mpc)^3. In this introductory analysis we carry out convergence tests on basic measures of non-linear clustering such as the matter power spectrum, the halo mass function and halo clustering, and we compare simulation predictions to those from current cosmological emulators. We also use our simulations to study matter and halo statistics, such as halo bias and clustering at the baryonic acoustic oscillation scale. Finally we measure the impact of baryonic physics on the matter and halo distributions., Comment: submitted to MNRAS, 23 pages, 19 figures, for future public data release, see https://www.mtng-project.org

The early release science results from $\textit{JWST}$ have yielded an unexpected abundance of high-redshift luminous galaxies that seems to be in tension with current theories of galaxy formation. However, it is currently difficult to draw definitive conclusions form these results as the sources have not yet been spectroscopically confirmed. It is in any case important to establish baseline predictions from current state-of-the-art galaxy formation models that can be compared and contrasted with these new measurements. In this work, we use the new large-volume ($L_\mathrm{box}\sim 740 \, \mathrm{cMpc}$) hydrodynamic simulation of the MillenniumTNG project, suitably scaled to match results from higher resolution - smaller volume simulations, to make predictions for the high-redshift ($z\gtrsim8$) galaxy population and compare them to recent $\textit{JWST}$ observations. We show that the simulated galaxy population is broadly consistent with observations until $z\sim10$. From $z\approx10-12$, the observations indicate a preference for a galaxy population that is largely dust-free, but is still consistent with the simulations. Beyond $z\gtrsim12$, however, our simulation results underpredict the abundance of luminous galaxies and their star-formation rates by almost an order of magnitude. This indicates either an incomplete understanding of the new $\textit{JWST}$ data or a need for more sophisticated galaxy formation models that account for additional physical processes such as Population~III stars, variable stellar initial mass functions, or even deviations from the standard $\Lambda$CDM model. We emphasise that any new process invoked to explain this tension should only significantly influence the galaxy population beyond $z\gtrsim10$, while leaving the successful galaxy formation predictions of the fiducial model intact below this redshift., Comment: Accepted for publication in MNRAS -- Part of the initial set of papers introducing the MillenniumTNG project. Visit www.mtng-project.org for more details

Approximate methods to populate dark matter halos with galaxies are of great utility to large galaxy surveys. However, the limitations of simple halo occupation models (HODs) preclude a full use of small-scale galaxy clustering data and call for more sophisticated models. We study two galaxy populations, luminous red galaxies (LRGs) and star-forming emission-line galaxies (ELGs), at two epochs, $z=1$ and $z=0$, in the large volume, high-resolution hydrodynamical simulation of the MillenniumTNG project. In a partner study we concentrated on the small-scale, one-halo regime down to $r\sim 0.1{\rm Mpc}/h$, while here we focus on modeling galaxy assembly bias in the two-halo regime, $r\gtrsim 1{\rm Mpc}/h$. Interestingly, the ELG signal exhibits scale dependence out to relatively large scales ($r\sim 20{\rm Mpc}/h$), implying that the linear bias approximation for this tracer is invalid on these scales, contrary to common assumptions. The 10-15\% discrepancy present in the standard halo model prescription is only reconciled when we augment our halo occupation model with a dependence on extrinsic halo properties ("shear" being the best-performing one) rather than intrinsic ones (e.g., concentration, peak mass). We argue that this fact constitutes evidence for two-halo galaxy conformity. Including tertiary assembly bias (i.e. a property beyond mass and "shear") is not an essential requirement for reconciling the galaxy assembly bias signal of LRGs, but the combination of external and internal properties is beneficial for recovering ELG the clustering. We find that centrals in low-mass haloes dominate the assembly bias signal of both populations. Finally, we explore the predictions of our model for higher-order statistics such as nearest-neighbor counts. The latter supplies additional information about galaxy assembly bias and can be used to break degeneracies between halo model parameters., Comment: 16 pages, 8 figures; submitted to MNRAS

Upcoming large galaxy surveys will subject the standard cosmological model, $\Lambda$CDM, to new precision tests. These can be tightened considerably if theoretical models of galaxy formation are available that can predict galaxy clustering and galaxy-galaxy lensing on the full range of measurable scales throughout volumes as large as those of the surveys and with sufficient flexibility that uncertain aspects of the underlying astrophysics can be marginalised over. This, in particular, requires mock galaxy catalogues in large cosmological volumes that can be directly compared to observation, and can be optimised empirically by Monte Carlo Markov Chains or other similar schemes to eliminate or estimate astrophysical parameters related to galaxy formation when constraining cosmology. Semi-analytic galaxy formation methods implemented on top of cosmological dark matter simulations offer a computationally efficient approach to construct physically based and flexibly parametrised galaxy formation models, and as such they are more potent than still faster, but purely empirical models. Here we introduce an updated methodology for the semi-analytic L-GALAXIES code, allowing it to be applied to simulations of the new MillenniumTNG project, producing galaxies directly on fully continuous past lightcones, potentially over the full sky, out to high redshift, and for all galaxies more massive than $\sim 10^8\,{\rm M}_\odot$. We investigate the numerical convergence of the resulting predictions, and study the projected galaxy clustering signals of different samples. The new methodology can be viewed as an important step towards more faithful forward-modelling of observational data, helping to reduce systematic distortions in the comparison of theory to observations., Comment: 24 pages, 18 figures, published in MNRAS

The stellar halo of the Milky Way records the history of its interactions with dwarf galaxies, whose subsequent destruction results in the formation of an extended stellar component. Recent works have suggested that galaxies with masses comparable to the Large Magellanic Cloud (LMC, $M_\star \sim 10^9\,{\rm M}_\odot$) may be the primary building blocks of the stellar halo of our Galaxy. We use cosmological simulations of the $\Lambda$ Cold Dark Matter model to investigate LMC-mass galaxies at $z=1-2$ using a semi-analytic model of galaxy formation. We find that LMC analogues at $z=2$ evolve until the present day along three distinct pathways: (1) those that are destroyed in Milky Way-mass hosts; (2) those that are themselves the main progenitors of Milky Way-mass galaxies; and (3) those that survive until $z=0$, with stellar mass $\sim$1.0 dex lower than typical Milky Ways. We predict that the properties of these galaxies at $z=2$ (stellar metallicities, sizes, gas content etc.) are largely indistinguishable, irrespective of which of these pathways is eventually taken; a survey targeting such galaxies in this redshift range would struggle to tell apart a 'destroyed' stellar halo progenitor from a 'surviving' LMC analogue. The only factor that determines the eventual fate of these galaxies is their proximity to a neighbouring Milky Way main progenitor at $z=2$: while the mean separation to a 'surviving' galaxy is around 7 Mpc, it is only 670 kpc to a 'destroyed' galaxy. This suggests that old stellar populations in the Milky Way may share intrinsic (i.e. non-dynamical) properties that are essentially indistinguishable from progenitors of its stellar halo., Comment: 9 pages, 5 figures, 1 table. Accepted for publication in MNRAS

In this series of papers, we present a simulation-based model for the non-linear clustering of galaxies based on separate modelling of clustering in real space and velocity statistics. In the first paper, we present an emulator for the real-space correlation function of galaxies, whereas the emulator of the real-to-redshift space mapping based on velocity statistics is presented in the second paper. Here, we show that a neural network emulator for real-space galaxy clustering trained on data extracted from the Dark Quest suite of N-body simulations achieves sub-per cent accuracies on scales $1 < r < 30 $ $h^{-1} \,\mathrm{Mpc}$, and better than $3\%$ on scales $r < 1$ $h^{-1}\mathrm{Mpc}$ in predicting the clustering of dark-matter haloes with number density $10^{-3.5}$ $(h^{-1}\mathrm{Mpc})^{-3}$, close to that of SDSS LOWZ-like galaxies. The halo emulator can be combined with a galaxy-halo connection model to predict the galaxy correlation function through the halo model. We demonstrate that we accurately recover the cosmological and galaxy-halo connection parameters when galaxy clustering depends only on the mass of the galaxies' host halos. Furthermore, the constraining power in $\sigma_8$ increases by about a factor of $2$ when including scales smaller than $5$ $h^{-1} \,\mathrm{Mpc}$. However, when mass is not the only property responsible for galaxy clustering, as observed in hydrodynamical or semi-analytic models of galaxy formation, our emulator gives biased constraints on $\sigma_8$. This bias disappears when small scales ($r < 10$ $h^{-1}\mathrm{Mpc}$) are excluded from the analysis. This shows that a vanilla halo model could introduce biases into the analysis of future datasets.

$JWST$'s first glimpse of the $z>10$ Universe has yielded a surprising abundance of luminous galaxy candidates. Here we present the most extreme of these systems: CEERS-1749. Based on $0.6-5\mu$m photometry, this strikingly luminous ($\approx$26 mag) galaxy appears to lie at $z\approx17$. This would make it an $M_{\rm{UV}}\approx-22$, $M_{\rm{\star}}\approx5\times10^{9}M_{\rm{\odot}}$ system that formed a mere $\sim220$ Myrs after the Big Bang. The implied number density of this galaxy and its analogues challenges virtually every early galaxy evolution model that assumes $\Lambda$CDM cosmology. However, there is strong environmental evidence supporting a secondary redshift solution of $z\approx5$: all three of the galaxy's nearest neighbors at $<2.5$" have photometric redshifts of $z\approx5$. Further, we show that CEERS-1749 may lie in a $z\approx5$ protocluster that is $\gtrsim5\times$ overdense compared to the field. Intense line emission at $z\approx5$ from a quiescent galaxy harboring ionized gas, or from a dusty starburst, may provide satisfactory explanations for CEERS-1749's photometry. The emission lines at $z\approx5$ conspire to boost the $>2\mu$m photometry, producing an apparent blue slope as well as a strong break in the SED. Such a perfectly disguised contaminant is possible only in a narrow redshift window ($\Delta z\lesssim0.1$), implying that the permitted volume for such interlopers may not be a major concern for $z>10$ searches, particularly when medium-bands are deployed. If CEERS-1749 is confirmed to lie at $z\approx5$, it will be the highest-redshift quiescent galaxy, or one of the lowest mass dusty galaxies of the early Universe detected to-date. Both redshift solutions of this intriguing galaxy hold the potential to challenge existing models of early galaxy evolution, making spectroscopic follow-up of this source critical., Comment: Submitted to ApJL. Fig. 1 shows the "superposed" state, Fig. 2 presents the potential z~5 protocluster, Fig. 3 shows some atypical z~5 SEDs, Fig. 4 studies the morphology, Fig. 5 summarizes the extraordinary situation at z~17, and Fig. 6 argues for why such interlopers may not be a major issue for z>10 searches. Comments warmly welcomed

© 2019 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society We use the IllustrisTNG (TNG) simulations to explore the galaxy-halo connection as inferred from state-of-the-art cosmological, magnetohydrodynamical simulations. With the high-mass resolution and large volume achieved by combining the 100 Mpc (TNG100) and 300 Mpc (TNG300) volumes, we establish the mean occupancy of central and satellite galaxies and their dependence on the properties of the dark matter haloes hosting them. We derive best-fitting HOD parameters from TNG100 and TNG300 for target galaxy number densities of n¯ g = 0.032 and n¯ g = 0.016 h3 Mpc−3, respectively, corresponding to a minimum galaxy stellar mass of M* ∼ 1.9 × 109 and M* ∼ 3.5 × 109 M☉, respectively, in hosts more massive than 1011 M☉. Consistent with previous work, we find that haloes located in dense environments, with low concentrations, later formation times, and high angular momenta are richest in their satellite population. At low mass, highly concentrated haloes and those located in overdense regions are more likely to contain a central galaxy. The degree of environmental dependence is sensitive to the definition adopted for the physical boundary of the host halo. We examine the extent to which correlations between galaxy occupancy and halo properties are independent and demonstrate that HODs predicted by halo mass and present-day concentration capture the qualitative dependence on the remaining halo properties. At fixed halo mass, concentration is a strong predictor of the stellar mass of the central galaxy, which may play a defining role in the fate of the satellite population. The radial distribution of satellite galaxies, which exhibits a universal form across a wide range of host halo mass, is described accurately by the best-fitting NFW density profile of their host haloes.

We show that gravitational waves have the potential to unravel the microphysical properties of dark matter due to the dependence of the binary black hole merger rate on cosmic structure formation, which is itself highly dependent on the dark matter scenario. In particular, we demonstrate that suppression of small-scale structure -- such as that caused by interacting, warm, or fuzzy dark matter -- leads to a significant reduction in the rate of binary black hole mergers at redshifts $z\gtrsim5$. This shows that future gravitational-wave observations will provide a new probe of physics beyond the $\Lambda$CDM model., Comment: 15 pages, 12 figures, matches version published in PRD

Abstract We present a comprehensive analysis of the shape of dark matter (DM) halos in a sample of 25 Milky Way-like galaxies in TNG50 simulation. Using an enclosed volume iterative method, we infer an oblate-to-triaxial shape for the DM halo with median T ≃ 0.24. We group DM halos into three different categories. Simple halos (32% of the population) establish principal axes whose ordering in magnitude does not change with radius and whose orientations are almost fixed throughout the halo. Twisted halos (32%) experience levels of gradual rotations throughout their radial profiles. Finally, stretched halos (36%) demonstrate a stretching in the lengths of their principal axes where the ordering of different eigenvalues changes with radius. Subsequently, the halo experiences a “rotation” of ∼90° where the stretching occurs. Visualizing the 3D ellipsoid of each halo, for the first time, we report signs of a reorienting ellipsoid in twisted and stretched halos. We examine the impact of baryonic physics on DM halo shape through a comparison to dark matter only (DMO) simulations. This suggests a triaxial (prolate) halo. We analyze the impacts of substructure on DM halo shape in both hydrodynamical and DMO simulations and confirm that they are subdominant. We study the distribution of satellites in our sample. In simple and twisted halos, the angle between satellites’ angular momentum and the galaxy’s angular momentum grows with radius. However, stretched halos show a flat distribution of angles. Overlaying our theoretical outcome on the observational results presented in the literature establishes a fair agreement.

Refereeing is a crucial component of publishing astronomical research, but few professional astronomers receive formal training on how to effectively referee a manuscript. In this article, we lay out considerations and best practices for referees. This document is intended as a tool for early career researchers to develop a fair, effective, and efficient approach to refereeing., Comment: Submitted to the Bulletin of the AAS

Observations of the high redshift universe provide a promising avenue for constraining the nature of the dark matter (DM). This will be even more true with the advent of the James Webb Space Telescope (JWST). We run cosmological simulations of galaxy formation as part of the Effective Theory of Structure Formation (ETHOS) project to compare high redshift galaxies in Cold (CDM) and alternative DM models which have varying relativistic coupling and self-interaction strengths. The interacting DM scenarios produce a cutoff in the linear power spectrum on small-scales, followed by a series of "dark acoustic oscillations". We find that DM interactions suppress the abundance of galaxies below $M_\star \sim 10^8\,M_\odot$ for the models considered. The cutoff in the power spectrum delays structure formation relative to CDM. Objects in ETHOS that end up at the same final masses as their CDM counterparts are characterised by a more vigorous phase of early star formation. While galaxies with $M_\star \lesssim 10^6\,M_\odot$ make up more than 60 per cent of star formation in CDM at $z\approx 10$, they contribute only about half the star formation density in ETHOS. These differences diminish with decreasing redshift. We find that the effects of DM self-interactions are negligible compared to effects of relativistic coupling (i.e. the effective initial conditions for galaxy formation) in all properties of the galaxy population we examine. Finally, we show that the clustering strength of galaxies at high redshifts depends sensitively on DM physics, although these differences are manifest on scales that may be too small to be measurable by JWST., Comment: 15 pages, 11 figures, 3 Tables. Comments are welcome. Updated based on referee's comments. Accepted to MNRAS on July 29, 2022

We employ the hydrodynamical simulation IllustrisTNG to inform the galaxy-halo connection of the Luminous Red Galaxy (LRG) and Emission Line Galaxy (ELG) samples of the Dark Energy Spectroscopic Instrument (DESI) survey at redshift z ~ 0.8. Specifically, we model the galaxy colors of IllustrisTNG and apply sliding DESI color-magnitude cuts, matching the DESI target densities. We study the halo occupation distribution model (HOD) of the selected samples by matching them to their corresponding dark matter halos in the IllustrisTNG dark matter run. We find the HOD of both the LRG and ELG samples to be consistent with their respective baseline models, but also we find important deviations from common assumptions about the satellite distribution, velocity bias, and galaxy secondary biases. We identify strong evidence for concentration-based and environment-based occupational variance in both samples, an effect known as "galaxy assembly bias". The central and satellite galaxies have distinct dependencies on secondary halo properties, showing that centrals and satellites have distinct evolutionary trajectories and should be modelled separately. These results serve to inform the necessary complexities in modeling galaxy-halo connection for DESI analyses and also prepare for building high-fidelity mock galaxies. Finally, we present a shuffling-based clustering analysis that reveals a 10-15% excess in the LRG clustering of modest statistical significance due to secondary galaxy biases. We also find a similar excess signature for the ELGs, but with much lower statistical significance. When a larger hydrodynamical simulation volume becomes available, we expect our analysis pipeline to pinpoint the exact sources of such excess clustering signatures., Comment: Accepted to MNRAS, comments welcome

The detection of Intermediate-Mass Black Holes (IMBHs) in dwarf galaxies is crucial to closing the gap in the wide mass distribution of black holes ($\sim 3 \, \rm M_{\odot}$ to $\sim 5 \times 10^{10} \, \rm M_{\odot}$). IMBHs originally located at the center of dwarfs that later collide with the Milky Way (MW) could be wandering, undetected, in our Galaxy. We used TNG50, the highest-resolution run of the IllustrisTNG project, to study the kinematics and dynamics of star clusters, in the appropriate mass range, acting as IMBH proxies in a MW analog galaxy. We showed that $\sim 87\%$ of our studied IMBHs drift inward. The radial velocity of these sinking IMBHs has a median magnitude of $\sim 0.44 \, \rm ckpc \, h^{-1} \, Gyr^{-1}$ and no dependence on the black hole mass. The central $1 \, \rm ckpc \, h^{-1}$ has the highest number density of IMBHs in the galaxy. A physical toy model with linear drag forces was developed to explain the orbital circularization with time. These findings constrain the spatial distribution of IMBHs, suggesting that future searches should focus on the central regions of the Galaxy. Additionally, we found that the 3D velocity distribution of IMBHs with respect to the galactic center has a mean of $\sim 180 \, \rm km \, s^{-1}$ and larger variance with decreasing radius. Remarkably, the velocity distribution relative to the local gas shows significantly lower values, with a mean of $\sim 88 \, \rm km \, s^{-1}$. These results are instrumental for predicting the accretion and radiation properties of IMBHs, facilitating their detection with future surveys., Comment: Accepted for publication in MNRAS. This is the final version of the manuscript. 10 pages, 12 figures