Search

The AGORA Cosmorun (arXiv:2106.09738) is a set of hydrodynamical cosmological zoom-in simulations carried out within the AGORA High-resolution Galaxy Simulations Comparison Project (arXiv:1308.2669,arXiv:1610.03066). These simulations show the formation and evolution of a Milky Way-sized galaxy using eight of the most widely used numerical codes in the community (Art-I, Enzo, Ramses, Changa, Gadget-3, Gear, Gizmo, and Arepo). In this short report, we describe the public release of the raw output data from all of these simulations at z = 8, 7, 6, 5, 4, 3, 2 (plus at z=1, 0 when available), and several metadata files containing the halo centers, virial quantities, and merger trees. The data from even thinner timesteps will be released as soon as the upcoming collaboration papers (VII-IX) are submitted and accepted., Comment: ArXiv Data Release from the AGORA Collaboration. arXiv admin note: text overlap with arXiv:2001.04354

Black holes are believed to be crucial in regulating star formation in massive galaxies, which makes it essential to faithfully represent the physics of these objects in cosmological hydrodynamics simulations. Limited spatial and mass resolution and the associated discreteness noise make following the dynamics of black holes especially challenging. In particular, dynamical friction, which is responsible for driving massive black holes towards the centres of galaxies, cannot be accurately modelled with softened $N$-body interactions. A number of subgrid models have been proposed to mimic dynamical friction or directly include its full effects in simulations. Each of these methods has its individual benefits and shortcomings, while all suffer from a common issue of being unable to represent black holes with masses below a few times the simulated dark matter particle mass. In this paper, we propose a correction for unresolved dynamical friction, which has been calibrated on simulations run with the code KETJU, in which gravitational interactions of black holes are not softened. We demonstrate that our correction is able to sink black holes with masses greater than the dark matter particle mass at the correct rate. We show that the impact of stochasticity is significant for low-mass black holes ($M_{\rm BH} \leq 5 M_{\rm DM}$) and propose a correction for stochastic heating. Combined, this approach is applicable to next generation cosmological hydrodynamics simulations that jointly track galaxy and black hole growth with realistic black hole orbits., Comment: 20 pages, 14 figures, accepted in MNRAS

In this fourth paper from the AGORA Collaboration, we study the evolution down to redshift $z=2$ and below of a set of cosmological zoom-in simulations of a Milky Way mass galaxy by eight of the leading hydrodynamic simulation codes. We also compare this CosmoRun suite of simulations with dark matter-only simulations by the same eight codes. We analyze general properties of the halo and galaxy at $z=4$ and 3, and before the last major merger, focusing on the formation of well-defined rotationally-supported disks, the mass-metallicity relation, the specific star formation rate, the gas metallicity gradients, and the non-axisymmetric structures in the stellar disks. Codes generally converge well to the stellar-to-halo mass ratios predicted by semi-analytic models at $z\sim$2. We see that almost all the hydro codes develop rotationally-supported structures at low redshifts. Most agree within 0.5 dex with the observed MZR at high and intermediate redshifts, and reproduce the gas metallicity gradients obtained from analytical models and low-redshift observations. We confirm that the inter-code differences in the halo assembly history reported in the first paper of the collaboration also exist in CosmoRun, making the code-to-code comparison more difficult. We show that such differences are mainly due to variations in code-dependent parameters that control the time-stepping strategy of the gravity solver. We find that variations in the early stellar feedback can also result in differences in the timing of the low-redshift mergers. All the simulation data down to $z=2$ and the auxiliary data will be made publicly available., Comment: Submitted to ApJ. Visit the AGORA Collaboration website (www.agorasimulations.org ) for more information

We analyze and compare the satellite halo populations at $z\sim2$ in the high-resolution cosmological zoom-in simulations of a $10^{12}\,{\rm M}_{\odot}$ target halo ($z=0$ mass) carried out on eight widely-used astrophysical simulation codes ({\sc Art-I}, {\sc Enzo}, {\sc Ramses}, {\sc Changa}, {\sc Gadget-3}, {\sc Gear}, {\sc Arepo-t}, and {\sc Gizmo}) for the {\it AGORA} High-resolution Galaxy Simulations Comparison Project. We use slightly different redshift epochs near $z=2$ for each code (hereafter ``$z\sim2$') at which the eight simulations are in the same stage in the target halo's merger history. After identifying the matched pairs of halos between the {\it CosmoRun} simulations and the DMO simulations, we discover that each {\it CosmoRun} halo tends to be less massive than its DMO counterpart. When we consider only the halos containing stellar particles at $z\sim2$, the number of satellite {\it galaxies} is significantly fewer than that of dark matter halos in all participating {\it AGORA} simulations, and is comparable to the number of present-day satellites near the Milky Way or M31. The so-called ``missing satellite problem' is fully resolved across all participating codes simply by implementing the common baryonic physics adopted in {\it AGORA} and the stellar feedback prescription commonly used in each code, with sufficient numerical resolution ($\lesssim100$ proper pc at $z=2$). We also compare other properties such as the stellar mass$-$halo mass relation and the mass$-$metallicity relation. Our work highlights the value of comparison studies such as {\it AGORA}, where outstanding problems in galaxy formation theory are studied simultaneously on multiple numerical platforms., Comment: Accepted for publication to the ApJ, 19 pages, 10 figures. Visit the AGORA Collaboration website (www.agorasimulations.org) for more information. For summary video, please see https://www.youtube.com/watch?v=HShRJYvPe1k

We analyze the circumgalactic medium (CGM) for eight commonly-used cosmological codes in the AGORA collaboration. The codes are calibrated to use identical initial conditions, cosmology, heating and cooling, and star formation thresholds, but each evolves with its own unique code architecture and stellar feedback implementation. Here, we analyze the results of these simulations in terms of the structure, composition, and phase dynamics of the CGM. We show properties such as metal distribution, ionization levels, and kinematics are effective tracers of the effects of the different code feedback and implementation methods, and as such they can be highly divergent between simulations. This is merely a fiducial set of models, against which we will in the future compare multiple feedback recipes for each code. Nevertheless, we find that the large parameter space these simulations establish can help disentangle the different variables that affect observable quantities in the CGM, e.g., showing that abundances for ions with higher ionization energy are more strongly determined by the simulation's metallicity, while abundances for ions with lower ionization energy are more strongly determined by the gas density and temperature., Comment: 32 pages, 17 figures. Published in ApJ

We present a high resolution numerical study of the sinking and merging of massive black holes (MBHs) with masses in the range of $10^3 - 10^7 \, \mathrm{M}_\odot$ in multiple minor mergers of low mass dark matter halos without and with galaxies ($4\times 10^8 \, \mathrm{M}_\odot \lesssim \mathrm{M}_{\mathrm{halo}} \lesssim 2\times 10^{10} \, \mathrm{M}_\odot)$. The Ketju simulation code, a combination of the Gadget tree solver with accurate regularised integration, uses unsoftened forces between the star/dark matter components and the MBHs for an accurate treatment of dynamical friction and scattering of dark matter/stars by MBH binaries or multiples. Post-Newtonian corrections up to order 3.5 for MBH interactions allow for coalescence by gravitational wave emission and gravitational recoil kicks. Low mass MBHs ($\lesssim 10^5 \, \mathrm{M}_\odot$) hardly sink to the centre or merge. Sinking MBHs have various complex evolution paths - binaries, triplets, free-floating MBHs, and dynamically or recoil ejected MBHs. Collisional interactions with dark matter alone can drive MBHs to coalescence. The highest mass MBHs of $\gtrsim 10^6 M_\odot$ mostly sink to the centre and trigger the scouring of dark matter and stellar cores. The scouring can transform a centrally baryon dominated system to a dark matter dominated system. Our idealized high-resolution study highlights the difficulty to bring in and keep low mass MBHs in the centres of low mass halos/galaxies - a remaining challenge for merger assisted MBH seed growth mechanisms., Comment: 22 pages, 15 figures

We investigate the build-up of the accreted stellar and dark matter haloes of Milky Way-like galaxies in APOSTLE suite of cosmological hydrodynamics simulations. We show that the stellar halo is made up primarily of stars stripped from a small number of massive dwarfs, most of which are disrupted by the present day. The dark matter halo, on the other hand, is made up primarily of small unresolved subhaloes ($\lesssim 10^6$ M$_{\odot}$) and a ``smooth'' component consisting of particles which were never bound to a subhalo. Despite these differences, the massive dwarfs that make up the majority of the stellar halo also contribute a significant fraction of the dark matter. The stars and dark matter stripped from these dwarfs are related through their kinematics and this leaves imprints in the phase-space structure of the haloes. We examine the relation between the location of features, such as caustics, in the phase space of the stars and dark halo properties. We show that the ``edge'' of the stellar halo is a probe of dark matter halo mass and assembly history. The edges of Milky Way-mass galaxies should be visible at a surface brightness of 31-36 mag arcsec$^{-2}$., Comment: 20 pages, 14 figures, Accepted for publication in MNRAS

The low dark matter density in the Fornax dwarf galaxy is often interpreted as being due to the presence of a constant density `core', but it could also be explained by the effects of Galactic tides. The latter interpretation has been disfavoured because it is apparently inconsistent with the orbital parameters and star formation history of Fornax. We revisit these arguments with the help of the APOSTLE cosmological hydrodynamics simulations. We show that simulated dwarfs with similar properties to Fornax are able to form stars after infall, so that star formation is not necessarily a good tracer of infall time. We also examine the constraints on the pericentre of Fornax and point out that small pericentres (<50 kpc) are not currently ruled out by the data, allowing for Fornax to be tidally influenced on its current orbit. Furthermore, we find that some dwarfs with large orbital pericentres can be stripped prior to infall due to interactions with more massive galaxies. Tidal effects lead to a reduction in the dark matter density, while the profile remains cuspy. Navarro-Frenk-White profiles are consistent with the kinematic data within 3$\sigma$ in the innermost regions, while profiles with shallow cusps or cores provide a better fit. We predict that if the reduction of the dark matter density in Fornax occurs, at least in part, because of the action of Galactic tides, then tidal tails should be visible with a surface brightness limit of $\sim35-36$ mag arcsec$^{-2}$ over a survey area of $\gtrsim$100 deg$^2$., Comment: 23 pages, 12 figures, accepted in MNRAS

We present an updated analysis of the gamma-ray flux from the directions of classical dwarf spheroidal galaxies, deriving new constraints on WIMP dark matter (DM) annihilation using a decade of Fermi-LAT data. Among the major novelties, we infer the dwarfs' J-factors by including new observations without imposing any a priori parametric profile for the DM distribution. While statistically compatible with results obtained from more conventional parameterisations, this procedure reduces the theoretical bias imposed on the data. Furthermore, we retain the full data-driven shape of the J-factors' empirical probability distributions when setting limits on DM, without imposing log-normality as is typically done. In conjunction with the data-driven J-factors, we use a new method for estimating the probability distribution function of the astrophysical background at the dwarf position, fully profiling over background uncertainties. We show that, for most classical dwarfs, the background systematic uncertainty dominates over the uncertainty on their J-factors. Raw distributions of J- and D-factors (the latter being the analogous of J-factors for decaying DM) are available upon request., Comment: 14 pages, 7 figures, 3 tables. Matches published version in JCAP. Raw distributions of J- and D-factors are available upon request

We test a non-parametric higher-order Jeans analysis method, GravSphere, on 32 simulated dwarf galaxies comparable to classical Local Group dwarfs like Fornax. The galaxies are selected from the APOSTLE suite of cosmological hydrodynamics simulations with Cold Dark Matter (CDM) and Self-Interacting Dark Matter (SIDM) models, allowing us to investigate cusps and cores in density distributions. We find that, for CDM dwarfs, the recovered enclosed mass profiles have a bias of no more than 10 per cent, with a 50 per cent scatter in the inner regions and a 20 per cent scatter near the half-light radius, consistent with standard mass estimators. The density profiles are also recovered with a bias of no more than 10 per cent and a scatter of 30 per cent in the inner regions. For SIDM dwarfs, the mass and density profiles are recovered within our 95 per cent confidence intervals, but are biased towards cuspy dark matter distributions. This is mainly due to a lack of sufficient constraints from the data. We explore the sources of scatter in the accuracy of the recovered profiles and suggest a $\chi^2$ statistic to separate successful models from biased ones. Finally, we show that the uncertainties on the mass profiles obtained with GravSphere are smaller than those for comparable Jeans methods, and that they can be further improved if stronger priors, motivated by cosmological simulations, are placed on the velocity anisotropy. We conclude that GravSphere is a promising Jeans-based approach for modelling dark matter distributions in dwarf galaxies., Comment: 22 pages, 13 figures. Accepted in MNRAS

A number of Local Group dwarf galaxies are known to have two spatially segregated stellar metallicity populations, a centrally concentrated metal-rich population and a more extended metal-poor population. In this work we discuss mechanisms that lead to the formation of two spatially segregated metallicity populations. Using a set of high-resolution hydrodynamical simulations of Local Group-like environments, we select a sample of satellite and field galaxies, spanning the stellar mass range $10^6-10^9$M$_{\odot}$, that exhibit bimodality in their metallicity distributions. Among those, we identify a subsample with a strong spatial segregation in the two populations. We find three distinct mechanisms for their formation. In field dwarfs and in a small fraction of satellites, a merger causes the metal-poor stars to migrate to larger radii and encourages the available gas to sink to the centre of the dwarf. Most of the gas is subsequently blown out of the halo through star formation feedback, but the remaining gas is consumed in the formation of a metal-rich population. In the exclusive case of satellites that have retained some of their gas at infall, it is the compression of this gas by ram pressure near pericentre that triggers the formation of metal-rich stars, whilst simultaneously preventing star formation at larger radii through stripping. Additionally, in a small number of field and satellite dwarfs, interactions with gaseous filaments and other galaxies can result in the formation of a metal-rich population. Regardless of the formation mechanism, a history of mergers typically enhances the spatial segregation., Comment: 20 pages, 19 figures, accepted in MNRAS

The existence of two kinematically and chemically distinct stellar subpopulations in the Sculptor and Fornax dwarf galaxies offers the opportunity to constrain the density profile of their matter haloes by measuring the mass contained within the well-separated half-light radii of the two metallicity subpopulations. Walker and Penarrubia have used this approach to argue that data for these galaxies are consistent with constant-density `cores' in their inner regions and rule out `cuspy' Navarro-Frenk-White (NFW) profiles with high statistical significance, particularly in the case of Sculptor. We test the validity of these claims using dwarf galaxies in the APOSTLE (A Project Of Simulating The Local Environment) Lambda cold dark matter cosmological hydrodynamics simulations of analogues of the Local Group. These galaxies all have NFW dark matter density profiles and a subset of them develop two distinct metallicity subpopulations reminiscent of Sculptor and Fornax. We apply a method analogous to that of Walker and Penarrubia to a sample of 53 simulated dwarfs and find that this procedure often leads to a statistically significant detection of a core in the profile when in reality there is a cusp. Although multiple factors contribute to these failures, the main cause is a violation of the assumption of spherical symmetry upon which the mass estimators are based. The stellar populations of the simulated dwarfs tend to be significantly elongated and, in several cases, the two metallicity populations have different asphericity and are misaligned. As a result, a wide range of slopes of the density profile are inferred depending on the angle from which the galaxy is viewed., Comment: 15 pages, 13 figures, accepted for publication in MNRAS

The recently discovered object Triangulum II appears to be an ultra faint dwarf spheroidal galaxy which may be one of the most dark matter dominated objects yet known. In this work we try to estimate the potential of this object for studies of the indirect detection of self-annihilating dark matter by obtaining its astrophysical J-factor. We perform a basic estimate of the velocity gradient to look for signs of the halo being tidally disrupted but show that the observed value is statistically compatible with zero velocity gradient. We solve the spherical Jeans equation using Markov Chain Monte Carlo (MCMC) engine GreAT and the Jeans analysis part of the CLUMPY package. We find the results point towards a very large J-factor, appearing to make Triangulum II one of the best targets in the search for dark matter. However we stress that the very small number of line of sight velocities currently available for this object make follow up studies essential., Comment: More stars used and text improved, results essentially unchanged. Matches version accepted in MNRAS, 7 figures, 7 pages

We extend the Generalised Langevin Equation (GLE) method [Phys. Rev. B 89, 134303 (2014)] to model a central classical region connected to two realistic thermal baths at two different temperatures. In such nonequilibrium conditions a heat flow is established, via the central system, in between the two baths. The GLE-2B (GLE two baths) scheme permits us to have a realistic description of both the dissipative central system and its surrounding baths. Following the original GLE approach, the extended Langevin dynamics scheme is modified to take into account two sets of auxiliary degrees of freedom corresponding to the mapping of the vibrational properties of each bath. These auxiliary variables are then used to solve the non-Markovian dissipative dynamics of the central region. The resulting algorithm is used to study a model of a short Al nanowire connected to two baths. The results of the simulations using the GLE-2B approach are compared to the results of other simulations that were carried out using standard thermostatting approaches (based on Markovian Langevin and Nose-Hoover thermostats). We concentrate on the steady state regime and study the establishment of a local temperature profile within the system. The conditions for obtaining a flat profile or a temperature gradient are examined in detail, in agreement with earlier studies. The results show that the GLE-2B approach is able to treat, within a single scheme, two widely different thermal transport regimes, i.e. ballistic systems, with no temperature gradient, and diffusive systems with a temperature gradient., Comment: present version accepted for publication in Phys. Rev. B (Apr 2016)