Milky Way-sized Galaxies Research Articles

Apostle–auriga: effects of stellar feedback subgrid models on the evolution of angular momentum in disc galaxies

ABSTRACT Utilizing the apostle–auriga simulations, which start from the same zoom-in initial conditions of Local Group-like systems, but run with different galaxy formation subgrid models and hydrodynamic solvers, we study the impact of stellar feedback models on the evolution of angular momentum in disc galaxies. At $z = 0$, auriga disc galaxies tend to exhibit higher specific angular momenta compared to their cross-matched apostle counterparts. By tracing the evolution history of the Lagrangian mass tracers of the in-situ star particles in the $z = 0$ galaxies, we find that the specific angular momentum distributions of the gas tracers from the two simulations at the halo accretion time are relatively similar. The present-day angular momentum difference is mainly driven by the physical processes occurring inside dark matter haloes, especially galactic fountains. Due to the different subgrid implementations of stellar feedback processes, auriga galaxies contain a high fraction of gas that has gone through recycled fountain (${\sim } 65$ per cent) which could acquire angular momentum through mixing with the high angular momentum circumgalactic medium (CGM). In apostle, however, the fraction of gas that has undergone the recycled fountain process is significantly lower (down to ${\sim } 20$ per cent for Milky Way-sized galaxies) and the angular momentum acquisition from the CGM is marginal. As a result, the present-day auriga galaxies overall have higher specific angular momenta.

The effect of local photoionization on the galaxy properties and the circumgalactic medium in simulations of Milky Way-sized galaxies

The effect of local photoionization on the galaxy properties and the circumgalactic medium in simulations of Milky Way-sized galaxies

Stellar cluster formation in a Milky Way-sized galaxy at z > 4 – I. The proto-globular cluster population and the imposter amongst us

ABSTRACT The formation history of globular clusters (GCs) at redshift z > 4 remains an unsolved problem. In this work, we use the cosmological, N-body hydrodynamical ‘zoom-in’ simulation GigaEris to study the properties and formation environment of proto-GC candidates in the region surrounding the progenitor of a Milky Way-sized galaxy. The simulation employs a modern implementation of smoothed-particle hydrodynamics, including metal-line cooling and metal and thermal diffusion. We define proto-GC candidate systems as gravitationally bound stellar systems with baryonic mass fraction Fb ≥ 0.75 and stellar velocity dispersion σ⋆ < 20 km s−1. At z = 4.4, we identify nine systems that satisfy our criteria, all of which form between 10 and 30 kpc from the centre of the main host. Their baryonic masses are in the range 105–107 M⊙. By the end of the simulation, they still have a relatively low stellar mass (M⋆ ∼ 104–105 M⊙) and a metallicity (−1.8 ≲ [Fe/H] ≲ −0.8) similar to the blue Galactic GCs. All of the identified systems except one appear to be associated with gas filaments accreting onto the main galaxy in the circum-galactic region and formed at z = 5–4. The exception is the oldest object, which appears to be a stripped compact dwarf galaxy that has interacted with the main halo between z = 5.8 and z = 5.2 and has lost its entire dark matter content due to tidal mass loss.

Spatial field reconstruction with INLA

Aims. Monte Carlo radiative transfer (MCRT) simulations are a powerful tool for understanding the role of dust in astrophysical systems and its influence on observations. However, due to the strong coupling of the radiation field and medium across the whole computational domain, the problem is non-local and non-linear, and such simulations are computationally expensive in the case of realistic 3D inhomogeneous dust distributions. We explore a novel technique for post-processing MCRT output to reduce the total computational run time by enhancing the output of computationally less expensive simulations of lower-quality. Methods. We combined principal component analysis (PCA) and non-negative matrix factorisation (NMF) as dimensionality reduction techniques together with Gaussian Markov random fields and the integrated nested Laplace approximation (INLA), an approximate method for Bayesian inference, to detect and reconstruct the non-random spatial structure in the images of lower signal-to-noise ratios or with missing data. Results. We tested our methodology using synthetic observations of a galaxy from the SKIRT Auriga project - a suite of high-resolution magnetohydrodynamic Milky Way-sized galaxies simulated in cosmological environment using a ‘zoom-in' technique. With this approach, we are able to reproduce high-photon-number reference images ~5 times faster with median residuals below ~20%.

Tidal disruption of star clusters in galaxy formation simulations

ABSTRACT We investigate the evolution of the tidal field experienced by massive star clusters using cosmological simulations of Milky Way-sized galaxies. Clusters in our simulations experience the strongest tidal force in the first few hundred Myr after formation, when the maximum eigenvalue of the tidal tensor reaches several times 104 Gyr−2. After about 1 Gyr the tidal field plateaus at a lower value, with the median λm ∼ 3 × 103 Gyr−2. The fraction of time clusters spend in high tidal strength (λm > 3 × 104 Gyr−2) regions also decreases with their age from ∼20 per cent immediately after formation to less than 1 per cent after 1 Gyr. At early ages both the in situ and ex situ clusters experience similar tidal fields, while at older ages the in situ clusters in general experience stronger tidal field due to their lower orbits in host galaxy. This difference is reflected in the survival of clusters: we looked into cluster disruption calculated in simulation runtime and found that ex situ star clusters of the same initial mass typically end up with higher bound fraction at the last available simulation snapshot than the in situ ones.

The Dawn of Disk Formation in a Milky Way-sized Galaxy Halo: Thin Stellar Disks at z > 4

Abstract We present results from GigaEris, a cosmological, N-body hydrodynamical “zoom-in” simulation of the formation of a Milky Way-sized galaxy halo with unprecedented resolution, encompassing of order a billion particles within the refined region. The simulation employs a modern implementation of smoothed-particle hydrodynamics, including metal-line cooling and metal and thermal diffusion. We focus on the early assembly of the galaxy, down to redshift z = 4.4. The simulated galaxy has properties consistent with extrapolations of the main sequence of star-forming galaxies to higher redshifts and levels off to a star formation rate of ∼60 M ⊙ yr−1 at z = 4.4. A compact, thin rotating stellar disk with properties analogous to those of low-redshift systems arises already at z ∼ 8. The galaxy rapidly develops a multi-component structure, and the disk, at least at these early stages, does not grow “upside-down” as often reported in the literature. Rather, at any given time, newly born stars contribute to sustain a thin disk. The kinematics reflect the early, ubiquitous presence of a thin disk, as a stellar disk component with v ϕ /σ R larger than unity is already present at z ∼ 9–10. Our results suggest that high-resolution spectro-photometric observations of very high-redshift galaxies should find thin rotating disks, consistent with the recent discovery of cold rotating gas disks by ALMA. Finally, we present synthetic images for the James Webb Space Telescope NIRCam camera, showing how the early disk would be easily detectable already at those early times.

Searching signals of dark matter from unidentified Fermi-LAT objects with H.E.S.S

Milky Way-sized galaxies harbor a population of unmerged dark matter subhalos, as shown from cosmological N-body simulations. These subhalos could shine in gamma-rays and be eventually detected as unidentified sources in gamma-ray surveys. From a thorough selection of unidentified Fermi-LAT Objects (UFOs), we observe four UFOs with H.E.S.S. and we search for very high-energy (VHE, E ≥ 100 GeV) gamma-ray emission. Considering dark matter masses above a few hundred GeV, the observed UFOs could be identified as dark matter subhalos, given their hard gamma-ray spectra in the few-ten-to-hundred GeV energy range. Since no significant very-high-energy gamma-ray emission is detected in any of the four UFOs dataset nor in the combined one, we derive constraints on the product of the velocity-weighted annihilation cross-section 〈συ〉 by the J-factor for the dark matter models. We derive 95% CL upper limits on 〈συ〉 J in W + W − and τ+τ− annihilation channels for the TeV dark matter particles. Considering thermal WIMPs, we derive constraints on the J-factors from the H.E.S.S. observations. Assuming model-dependent predictions from cosmological N-body simulations on the J-factor distribution for Milky Way-sized galaxies, the dark matter models with masses greater than 0.3 TeV for the UFO emissions can be ruled out at high confidence level.

Phase spirals in cosmological simulations of Milky Way-sized galaxies

ABSTRACT We study the vertical perturbations in the galactic disc of the Milky Way-size high-resolution hydrodynamical cosmological simulation named GARROTXA. We detect phase spirals in the vertical projection Z − VZ of disc’s stellar particles for the first time in this type of simulations. Qualitatively similar structures were detected in the recent Gaia data, and their origin is still under study. In our model, the spiral-like structures in the phase space are present in a wide range of times and locations across the disc. By accounting for an evolving mix of stellar populations, we observe that, as seen in the data, the phase spirals are better observed in the range of younger-intermediate star particles. We measure the intensity of the spiral with a Fourier decomposition and find that these structures appear stronger near satellite pericentres. Current dynamical models of the phase spiral considering a single perturber required a mass at least of the order of 1010 M⊙, but all three of our satellites have masses of the order of ∼108 M⊙. We suggest that there are other mechanisms at play which appear naturally in our model such as the physics of gas, collective effect of multiple perturbers, and a dynamically cold population that is continuously renovated by the star formation. Complementing collisionless isolated N-body models with the use of fully cosmological simulations with enough resolution can provide new insights into the nature/origin of the phase spiral.

Evolving beyond z=0: insights about the future of stars and the intergalactic medium

ABSTRACT We present results from seven cosmological simulations that have been extended beyond the present era as far as redshift z ≈ −0.99 or $t\approx 96.0\, {\rm Gyr}$, using the Enzo simulation code. We adopt the calibrated star formation and feedback prescriptions from our previous work on reproducing the Milky Way with Enzo, with modifications to the simulation code, chemistry, and cooling library. We then consider the future behaviour of the halo mass function (HMF), the equation of state (EOS) of the IGM, and the cosmic star formation history (SFH). Consistent with previous work, we find a freeze out in the HMF at z ≈ −0.6 or $t\approx 28.1\, {\rm Gyr}$. The evolution of the EOS of the IGM presents an interesting case study of the cosmological coincidence problem, where there is a sharp decline in the IGM temperature immediately after z = 0. For the SFH, the simulations produce a peak and a subsequent decline into the future. However, we do find a turnaround in the SFH after z ≈ −0.98 or $t\approx 82.4\, {\rm Gyr}$ in some simulations, which we attribute to limitations of the criteria used for star formation. By integrating the SFH in time up to z ≈ −0.92 or $t\approx 55.1\, {\rm Gyr}$, the simulation with the best spatial resolution predicts an asymptotic total stellar mass that is very close to that obtained from extrapolating the fit of the observed SFR. Lastly, we investigate the future evolution of the partition of baryons within a Milky Way-sized galaxy, using both a zoom and a box simulation. Despite vastly different resolutions, these simulations predict individual haloes containing an equal fraction of baryons in stars and gas at the time of freeze out.

Evolution of disc thickness in simulated high-redshift galaxies

ABSTRACT We study the growth of stellar discs of Milky Way-sized galaxies using a suite of cosmological simulations. We calculate the half-mass axis lengths and axis ratios of stellar populations split by age in galaxies with stellar mass $M_{*}=10^7\!-\!10^{10}\, \mathrm{M}_{\odot }$ at redshifts z > 1.5. We find that in our simulations stars always form in relatively thin discs, and at ages below 100 Myr are contained within half-mass height z1/2 ∼ 0.1 kpc and short-to-long axial ratio z1/2/x1/2 ∼ 0.15. Disc thickness increases with the age of stellar population, reaching median z1/2 ∼ 0.8 kpc and z1/2/x1/2 ∼ 0.6 for stars older than 500 Myr. We trace the same group of stars over the simulation snapshots and show explicitly that their intrinsic shape grows more spheroidal over time. We identify a new mechanism that contributes to the observed disc thickness: rapid changes in the orientation of the galactic plane mix the configuration of young stars. The frequently mentioned ‘upside-down’ formation scenario of galactic discs, which posits that young stars form in already thick discs at high redshift, may be missing this additional mechanism of quick disc inflation. The actual formation of stars within a fairly thin plane is consistent with the correspondingly flat configuration of dense molecular gas that fuels star formation.

Self-gravitating strange dark matter halos around galaxies

A new family of nonrelativistic, Newtonian, non-quantum equilibrium configurations describing galactic halos is introduced, by considering strange quark matter conglomerates with masses larger than about 8 GeV as new possible components of the dark matter. Originally introduced to explain the state of matter in neutron stars, such conglomerates may also form in the high-density and temperature conditions of the primordial Universe and then decouple from ordinary baryonic matter, providing the fundamental components of dark matter for the formation of pristine gravitational potential wells and the subsequent evolution of cosmic structures. The obtained results for halo mass and radius are consistent with the rotational velocity curve observed in the Galaxy. Additionally, the average density of such dark matter halos is similar to that derived for halos of dwarf spheroidal galaxies, which can therefore be interpreted as downscaled versions of larger dark matter distributions around Milky Way-sized galaxies and hint for a common origin of the two families of cosmic structures.

A dark matter profile to model diverse feedback-induced core sizes of ΛCDM haloes

ABSTRACT We analyse the cold dark matter density profiles of 54 galaxy haloes simulated with Feedback In Realistic Environments (FIRE)-2 galaxy formation physics, each resolved within $0.5{{\ \rm per\ cent}}$ of the halo virial radius. These haloes contain galaxies with masses that range from ultrafaint dwarfs ($M_\star \simeq 10^{4.5}\, \mathrm{M}_{\odot }$) to the largest spirals ($M_\star \simeq 10^{11}\, \mathrm{M}_{\odot }$) and have density profiles that are both cored and cuspy. We characterize our results using a new, analytic density profile that extends the standard two-parameter Einasto form to allow for a pronounced constant density core in the resolved innermost radius. With one additional core-radius parameter, rc, this three-parameter core-Einasto profile is able to characterize our feedback-impacted dark matter haloes more accurately than other three-parameter profiles proposed in the literature. To enable comparisons with observations, we provide fitting functions for rc and other profile parameters as a function of both M⋆ and M⋆/Mhalo. In agreement with past studies, we find that dark matter core formation is most efficient at the characteristic stellar-to-halo mass ratio M⋆/Mhalo ≃ 5 × 10−3, or $M_{\star } \sim 10^9 \, \mathrm{M}_{\odot }$, with cores that are roughly the size of the galaxy half-light radius, rc ≃ 1−5 kpc. Furthermore, we find no evidence for core formation at radii $\gtrsim 100\ \rm pc$ in galaxies with M⋆/Mhalo < 5 × 10−4 or $M_\star \lesssim 10^6 \, \mathrm{M}_{\odot }$. For Milky Way-size galaxies, baryonic contraction often makes haloes significantly more concentrated and dense at the stellar half-light radius than DMO runs. However, even at the Milky Way scale, FIRE-2 galaxy formation still produces small dark matter cores of ≃ 0.5−2 kpc in size. Recent evidence for a ∼2 kpc core in the Milky Way’s dark matter halo is consistent with this expectation.

Artificial versus biological intelligence in the Cosmos: clues from a stochastic analysis of the Drake equation

AbstractThe Drake equation has been used many times to estimate the number of observable civilizations in the galaxy. However, the uncertainty of the outcome is so great that any individual result is of limited use, as predictions can range from a handful of observable civilizations in the observable universe to tens of millions per Milky Way-sized galaxy. A statistical investigation shows that the Drake equation, despite its uncertainties, delivers robust predictions of the likelihood that the prevalent form of intelligence in the universe is artificial rather than biological. The likelihood of artificial intelligence far exceeds the likelihood of biological intelligence in all cases investigated. This conclusion is contingent upon a limited number of plausible assumptions. The significance of this outcome for the Fermi paradox is discussed.

Origin of giant stellar clumps in high-redshift galaxies

ABSTRACT We examine the nature of kpc-scale clumps seen in high-redshift galaxies using a suite of cosmological simulations of galaxy formation. We identify rest-frame UV clumps in mock HST images smoothed to 500 pc resolution, and compare them with the intrinsic 3D clumps of young stars identified in the simulations with 100 pc resolution. According to this comparison for the progenitors of Milky Way-sized galaxies probed by our simulations, we expect that the stellar masses of the observed clumps are overestimated by as much as an order of magnitude, and that the sizes of these clumps are also overestimated by factor of several, due to a combination of spatial resolution and projection. The masses of young stars contributing most of the UV emission can also be overestimated by factor of a few. We find that most clumps of young stars present in a simulation at one time dissolve on a timescale shorter than ∼150 Myr. Some clumps with dense cores can last longer but eventually disperse. Most of the clumps are not bound structures, with virial parameter αvir > 1. We find similar results for clumps identified in mock maps of H α emission measure. We examine the predictions for effective clump sizes from the linear theory of gravitational perturbations and conclude that they are inconsistent with being formed by global disc instabilities. Instead, the observed clumps represent random projections of multiple compact star-forming regions.

Accelerated orbital decay of supermassive black hole binaries in merging nuclear star clusters

ABSTRACT The coalescence of supermassive black holes (SMBHs) should generate the strongest sources of gravitational waves (GWs) in the Universe. However, the dynamics of their coalescence is the subject of much debate. In this study, we use a suite of N-body simulations to follow the merger of two nuclear star clusters (NSCs), each hosting an SMBH in their centre. We find that the presence of distinct star clusters around each SMBH has important consequences for the dynamical evolution of the SMBH binary: (i) The separation between the SMBHs decreases by a few orders of magnitude in the first few Myrs by the combined effects of dynamical friction and a drag force caused by tidally stripped stars. In fact, this is a significant speedup for equal mass ratio binaries, and becomes extreme for unequal mass ratios, e.g. 1:10 or 1:100, which traditional dynamical friction alone would not permit to bind. (ii) The subsequent binary hardening is driven by the gravitational slingshots between the SMBH binary and stars, and also depends on the mass ratio between the SMBHs. Thus, with this additional drag force, we find that all SMBHs in our suite coalesce within a Hubble time. Given that about 50 per cent of Milky Way-sized galaxies host NSCs, our results are encouraging for upcoming GW observations with the Laser Interferometer Space Antenna – LISA – which will detect SMBH coalescence in the 104–107 M⊙ mass range.

Formation and evolution of globular clusters in cosmological simulations

AbstractIn a series of three papers, we introduced a novel cluster formation model that describes the formation, growth, and disruption of star clusters in high-resolution cosmological simulations. We tested this model on a Milky Way-sized galaxy and found that various properties of young massive clusters, such as the mass function and formation efficiency, are consistent with observations in the local universe. Interestingly, most massive clusters – globular cluster candidates – are preferentially formed during major merger events. We follow the dynamical evolution of clusters in the galactic tidal field. Due to tidal disruption, the cluster mass function evolves from initial power law to a peaked shape. The surviving clusters at z = 0 show a broad range of metallicity [Fe/H] from -3 to -0.5. A robust prediction of the model is the age–metallicity relation, in which metal-rich clusters are systematically younger than metal-poor clusters by up to 3 Gyr.

Star cluster formation in cosmological simulations – III. Dynamical and chemical evolution

In previous papers of this series, we developed a new algorithm for modelling the formation of star clusters in galaxy formation simulations. Here we investigate how dissolution of bound star clusters affects the shape of the cluster mass function and the metallicity distribution of surviving clusters. Cluster evolution includes the loss of stars that become unbound due to tidal disruption as well as mass-loss due to stellar evolution. We calculate the tidal tensor along cluster trajectories and use it to estimate the instantaneous mass-loss rate. The typical tidal tensor exhibits large variations on a time-scale of ∼100 Myr, with maximum eigenvalue of 107 Gyr−2, and median value of 104 Gyr−2 for the first Gyr after cluster formation. As a result of dynamical disruption, at the final available output of our simulations at redshift z ≈ 1.5, the cluster mass function has an approximately lognormal shape peaked at |${\sim }10^{4.3}\, \mathrm{ M}_{\odot }$|⁠. Extrapolation of the disruption to z = 0 results in too many low-mass clusters compared to the observed Galactic globular clusters (GCs). Over 70 per cent of GC candidates are completely disrupted before the present; only 10 per cent of the total GC candidate mass remains in surviving clusters. The total mass of surviving clusters at z = 0 varies from run to run in the range |$(2\!-\!6)\times 10^7\, \mathrm{ M}_{\odot }$|⁠, consistent with the observed mass of GC systems in Milky Way-sized galaxies. The metallicity distributions of all massive star clusters and of the surviving GCs have similar shapes but different normalization because of cluster disruption. The model produces a larger fraction of very metal-poor clusters than observed. A robust prediction of the model is the age–metallicity relation, in which metal-rich clusters are systematically younger than metal-poor clusters by up to 3 Gyr.

Optical integral field spectroscopy observations applied to simulated galaxies: testing the fossil record method

By means of post-processed cosmological Hydrodynamics simulations we explore the ability of the fossil record method to recover the stellar mass, age gradients, and global/radial star formation and mass growth histories of galaxies observed with an optical Integral Field Unity. We use two simulations of representative Milky Way-sized galaxies, one disk dominated and another bulge-dominated. We generate data cubes emulating MaNGA/SDSS-IV instrumental and observational conditions, and analyze them using the Pipe3D pipeline. For optimal MaNGA instrumental/observational setups, the global masses, ages, and dust extinctions of both galaxies are well recovered within the uncertainties, while for the radial distributions, the trend is to get slightly flatter age gradients and weaker inside-out growth modes than the true ones. The main bias is in recovering younger stellar populations in the inner, older regions, and slightly older stellar populations in the outer, younger regions, as well as lower global stellar masses. These trends are enhanced as lower is the spatial sampling/resolution (less number of fibers), more inclined and dust-attenuated is the observed galaxy, lower is the signal-to-noise ratio, and worse is the seeing. The recovered properties and histories are more affected for the disk-dominated galaxy with a prominent inside-out growth mode than for the bulge-dominated galaxy assembled earlier and more coherently. We make publicly available the code and the mock IFU datacubes used in this paper.

Imprints of Chameleon f(R) Gravity on Galaxy Rotation Curves

Current constraints on gravity are relatively weak on galactic and intergalactic scales. Screened modified gravity models can exhibit complex behaviour there without violating stringent tests of gravity within our Solar System. They might hence provide viable extensions of the theory of gravity. Here, we use galaxy kinematics to constrain screened modified gravity models. We focus on chameleon $f(R)$ gravity and predict its impact on galaxy rotation curves and radial acceleration relations. This is achieved by post-processing state-of-the-art galaxy formation simulations from the \textsc{auriga project}, using the \textsc{mg-gadget} code. For a given galaxy, the surface dividing screened and un-screened regions adopts an oblate shape, reflecting the disc morphology of the galaxy's mass distribution. At the `screening radius'---the radius at which screening is triggered in the disc plane---characteristic `upturns' are present in both rotation curves and radial acceleration relations. The locations of these features depend on various factors, such as the galaxy mass, the concentration of the density profile and the value of the background field amplitude $f_{R0}$. Self-screening of stars and environmental screening also play a role. For Milky Way-size galaxies, we find that a model with $|f_{R0}|=10^{-7}$ results in rotation curves that are indistinguishable from $\Lambda$CDM, while for $|f_{R0}| \geq 2 \times 10^{-6}$ the simulated galaxies are entirely unscreened, violating Solar System constraints. For intermediate values, distinct upturns are present. With a careful statistical analysis of existing samples of observed rotation curves, including lower mass objects, constraints on $f(R)$ gravity with a sensitivity down to $|f_{R0}|\sim10^{-7}$ should be possible.

Preferential Accretion in the Supermassive Black Holes of Milky Way-size Galaxies Due to Direct Feeding by Satellites

Abstract Using a new, high-resolution cosmological hydrodynamic simulation of a Milky Way-type (MW-type) galaxy, we explore how a merger-rich assembly history affects the mass budget of the central supermassive black hole (SMBH). We examine a MW-mass halo at the present epoch whose evolution is characterized by several major mergers to isolate the importance of merger history on black hole (BH) accretion. This study is an extension of Bellovary et al. (2013), which analyzed the accretion of high mass, high-redshift galaxies and their central BHs, and found that the gas content of the central BH reflects what is accreted by the host galaxy halo. In this study, we find that a merger-rich galaxy will have a central SMBH preferentially fed by gas accreted through mergers. Moreover, we find that the gas composition of the inner ∼10 kpc of the galaxy can account for the increase of merger-accreted gas fueling the SMBH. Through an investigation of the angular momentum of the gas entering the host and its SMBH, we determine that gas accreted through mergers enters the galaxy halo with lower angular momentum compared to smooth accretion, partially accounting for the preferential fueling witnessed in the SMBH. In addition, the presence of mergers, particularly major mergers, also helps funnel low angular momentum gas more readily to the center of the galaxy. Our results imply that galaxy mergers play an important role in feeding the SMBH in MW-type galaxies with merger-rich histories.