Source author record

Oscar Agertz

Oscar Agertz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

29works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

29 published item(s)

preprint2025arXiv

A new classification of ex-situ and in-situ Galactic globular clusters based on a method trained on Milky Way analogues in the TNG50 cosmological simulations

We present a novel method combining existing cosmological simulations and orbital integration to study the hierarchical assembly of globular cluster (GC) populations in the Milky Way (MW). Our method models the growth and evolution of GC populations across various galactic environments as well as the dynamical friction and mass-loss experienced by these objects. This allows us to follow the trajectory of $\sim$18,000 GCs over cosmic time in 198 MW-like galaxies from TNG50. This cosmological-scale tracking of the dynamics of in-situ and ex-situ GC populations with such a large statistical sample allows us to confirm the presence of an overlap between the two populations in MW-like galaxies, occurring below an energy threshold of $E < -0.7 |E_{\rm circ}(r_{\rm hm}^{*})|$ where $E_{\rm circ}(r_{\rm hm}^{*})$ is the energy of a circular orbit at the galaxy's stellar half-mass radius $r_{\rm hm}^{*}$. Our results challenge the validity of current classification schemes commonly adopted in the literature, which ultimately fail to provide a clear separation between the two populations. Instead, they tend to isolate only a subset of the ex-situ GCs. More precisely, we argue that it is highly unlikely to find in-situ clusters at $E > -0.7 |E_{\rm circ}(r_{\rm hm}^{*})|$, and that the real challenge lies in distinguishing the two populations below this energy threshold. In this context, we provide new predictions regarding the origins of the MW GCs observed with Gaia, as well as a comparison with existing literature. Additionally, we highlight that even if ex-situ clusters share a common origin, they inevitably lose their dynamical coherence in the $E$-$L_{z}$ space within MW-like galaxies. We observe a dispersion of GC groups as a function of $E$ and $L_{z}$, primarily driven by the evolution of the galactic potential over time and by dynamical friction, respectively.

preprint2025arXiv

The radial acceleration relation at the EDGE of galaxy formation: testing its universality in low-mass dwarf galaxies

A tight correlation between the baryonic and observed acceleration of galaxies has been reported over a wide range of mass ($10^8 < M_{\rm bar}/{\rm M}_\odot < 10^{11}$) - the Radial Acceleration Relation (RAR). This has been interpreted as evidence that dark matter is actually a manifestation of some modified weak-field gravity theory. In this paper, we study the radially resolved RAR of 12 nearby dwarf galaxies, with baryonic masses in the range $10^4 < M_{\rm bar}/{\rm M}_\odot < 10^{7.5}$, using a combination of literature data and data from the MUSE-Faint survey. We use stellar line-of-sight velocities and the Jeans modelling code GravSphere to infer the mass distributions of these galaxies, allowing us to compute the RAR. We compare the results with the EDGE simulations of isolated dwarf galaxies with similar stellar masses in a $Λ$CDM cosmology. We find that most of the observed dwarf galaxies lie systematically above the low-mass extrapolation of the RAR. Each galaxy traces a locus in the RAR space that can have a multi-valued observed acceleration for a given baryonic acceleration, while there is significant scatter from galaxy to galaxy. Our results indicate that the RAR does not apply to low-mass dwarf galaxies and that the inferred baryonic acceleration of these dwarfs does not contain enough information, on its own, to derive the observed acceleration. The simulated EDGE dwarfs behave similarly to the real data, lying systematically above the extrapolated RAR. We show that, in the context of modified weak-field gravity theories, these results cannot be explained by differential tidal forces from the Milky Way, nor by the galaxies being far from dynamical equilibrium, since none of the galaxies in our sample seems to experience strong tides. As such, our results provide further evidence for the need for invisible dark matter in the smallest dwarf galaxies.

preprint2023arXiv

Pre-supernova feedback sets the star cluster mass function to a power law and reduces the cluster formation efficiency

The star cluster initial mass function is observed to have an inverse power law exponent around 2, yet there is no consensus on what determines this distribution, and why some variation is observed in different galaxies. Furthermore, the cluster formation efficiency covers a range of values, particularly when considering different environments. These clusters are often used to empirically constrain star formation and as fundamental units for stellar feedback models. Detailed galaxy models must therefore accurately capture the basic properties of observed clusters to be considered predictive. We use hydrodynamical simulations of a dwarf galaxy as a laboratory to study star cluster formation. We test different combinations of stellar feedback mechanisms, including stellar winds, ionizing radiation, and supernovae. Each feedback mechanism affects the cluster formation efficiency and cluster mass function. Increasing the feedback budget by combining the different types of feedback decreases the cluster formation efficiency by reducing the number of massive clusters. Ionizing radiation is found to be especially influential. This effect depends on the timing of feedback initiation, as shown by comparing early and late feedback. Early feedback occurs from ionizing radiation and stellar winds with onset immediately after a massive star is formed. Late feedback occurs when energy injection only starts after the main-sequence lifetime of the most massive SN progenitor, a timing that is further influenced by the choice of the most massive SN progenitor. Late feedback alone results in a broad, flat mass function, approaching a log-normal shape in the complete absence of feedback. Early feedback, on the other hand, produces a power-law cluster mass function with lower formation efficiency, albeit with a steeper slope than that usually observed.

preprint2022arXiv

EDGE: the puzzling ellipticity of Eridanus II's star cluster and its implications for dark matter at the heart of an ultra-faint dwarf

The Eridanus II (EriII) 'ultra-faint' dwarf has a large ($15\,\text{pc}$) and low mass ($4.3\times10^3\,\text{M}_\odot$) star cluster (SC) offset from its centre by $23\pm3\,\text{pc}$ in projection. Its size and offset are naturally explained if EriII has a central dark matter core, but such a core may be challenging to explain in a $Λ$CDM cosmology. In this paper, we revisit the survival and evolution of EriII's SC, focussing for the first time on its puzzlingly large ellipticity ($0.31^{+0.05}_{-0.06}$). We perform a suite of 960 direct $N$-body simulations of SCs, orbiting within a range of spherical background potentials fit to ultra-faint dwarf (UFD) galaxy simulations. We find only two scenarios that come close to explaining EriII's SC. In the first, EriII has a low density dark matter core (of size $\sim70\,\text{pc}$ and density $\lesssim2\times10^8\,\text{M}_{\odot}\,\text{kpc}^{-3}$). In this model, the high ellipticity of EriII's SC is set at birth, with the lack of tidal forces in the core allowing its ellipticity to remain frozen in for long times. In the second, EriII's SC orbits in a partial core, with its high ellipticity owing to its imminent tidal destruction. However, this latter model struggles to reproduce the large size of EriII's SC, and it predicts substantial tidal tails around EriII's SC that should have already been seen in the data. This leads us to favour the cored model. We discuss potential caveats to these findings, and the implications of the cored model for galaxy formation and the nature of dark matter.

preprint2022arXiv

EDGE: The sensitivity of ultra-faint dwarfs to a metallicity-dependent initial mass function

Motivated by the observed bottom-light initial mass function (IMF) in faint dwarfs, we study how a metallicity-dependent IMF affects the feedback budget and observables of an ultra-faint dwarf galaxy. We model the evolution of a low-mass ($\approx 8 \, \times \, 10^{8} \, \rm M_{\odot}$) dark matter halo with cosmological, zoomed hydrodynamical simulations capable of resolving individual supernovae explosions, which we complement with an empirically motivated subgrid prescription for systematic IMF variations. In this framework, at the low gas metallicities typical of faint dwarfs, the IMF of newborn stellar populations becomes top-heavy, increasing the efficiency of supernova and photoionization feedback in regulating star formation. This results in a 100-fold reduction of the final stellar mass of the dwarf compared to a canonical IMF, at fixed dynamical mass. The increase in the feedback budget is nonetheless met by increased metal production from more numerous massive stars, leading to nearly constant iron content at $z=0$. A metallicity-dependent IMF therefore provides a mechanism to produce low-mass ($\rm M_{\star}\sim 10^3 \rm M_{\odot}$), yet enriched ($\rm [Fe/H]\approx -2$) field dwarf galaxies, thus opening a self-consistent avenue to populate the plateau in $\rm [Fe/H]$ at the faintest end of the mass--metallicity relation.

preprint2022arXiv

EDGE: What shapes the relationship between HI and stellar observables in faint dwarf galaxies?

We show how the interplay between feedback and mass-growth histories introduces scatter in the relationship between stellar and neutral gas properties of field faint dwarf galaxies ($M_{\star} \lessapprox 10^{6} M_{\odot}$). Across a suite of cosmological, high-resolution zoomed simulations, we find that dwarf galaxies of stellar masses $10^5 \leq M_{\star} \leq 10^{6} M_{\odot}$ are bimodal in their cold gas content, being either HI-rich or HI-deficient. This bimodality is generated through the coupling between (i) the modulation of HI contents by the background of ultraviolet radiation (UVB) at late times and (ii) the significant scatter in the stellar-mass-halo-mass relationship induced by reionization. Furthermore, our HI-rich dwarfs exhibit disturbed and time-variable neutral gas distributions primarily due to stellar feedback. Over the last four billion years, we observe order-of-magnitude changes around the median $M_{HI}$, factor-of-a-few variations in HI spatial extents, and spatial offsets between HI and stellar components regularly exceeding the galaxies' optical sizes. Time variability introduces further scatter in the $M_{\star}-M_{HI}$ relation and affects a galaxy's detectability in HI at any given time. These effects will need to be accounted for when interpreting observations of the population of faint, HI-bearing dwarfs by the combination of optical and radio wide, deep surveys.

preprint2022arXiv

From giant clumps to clouds -- III. The connection between star formation and turbulence in the ISM

Supersonic gas turbulence is a ubiquitous property of the interstellar medium. The level of turbulence, quantified by the gas velocity dispersion ($σ_{\rm g}$), is observed to increase with the star formation rate (SFR) of a galaxy, but it is yet not established whether this trend is driven by stellar feedback or gravitational instabilities. In this work we carry out hydrodynamical simulations of entire disc galaxies, with different gas fractions, to understand the origins of the SFR-$σ_{\rm g}$ relation. We show that disc galaxies reach the same levels of turbulence regardless of the presence of stellar feedback processes, and argue that this is an outcome of the way disc galaxies regulate their gravitational stability. The simulations match the SFR-$σ_{\rm g}$ relation up to SFRs of the order of tens of M$_\odot$ yr$^{-1}$ and $σ_{\rm g}\sim 50$ km s$^{-1}$ in neutral hydrogen and molecular gas, but fail to reach the very large values ($>100$ km s$^{-1}$) reported in the literature for rapidly star-forming galaxies. We demonstrate that such high values of $σ_{\rm g}$ can be explained by 1) insufficient beam smearing corrections in observations, and 2) stellar feedback being coupled to the ionised gas phase traced by recombination lines. Given that the observed SFR-$σ_{\rm g}$ relation is composed of highly heterogeneous data, with $σ_{\rm g}$ at high SFRs almost exclusively being derived from H$α$ observations of high redshift galaxies with complex morphologies, we caution against analytical models that attempt to explain the SFR-$σ_{\rm g}$ relation without accounting for these effects.

preprint2021arXiv

From giant clumps to clouds -- II. The emergence of thick disc kinematics from the conditions of star formation in high redshift gas rich galaxies

High redshift disc galaxies are more gas rich, clumpier, and more turbulent than local Universe galaxies. This early era of galaxy formation imprints the distribution and kinematics of the stars that we observe today, but it is not yet well established how. In this work, we use simulations of isolated Milky Way-mass disc galaxies to study how kinematic properties of stars change when varying the gas fraction. This allows us to quantify the roles played by internal processes, e.g. gas turbulence and gravitational scattering off massive gas clumps, in establishing the observed stellar velocity dispersions and orbital eccentricities. We find that models with gas fractions $>20$ per cent feature a turbulent and clumpy interstellar medium (ISM), leading to zero-age stellar velocity dispersions $\sim 20-30~{\rm km\, s}^{-1}$ and high mean orbital eccentricities. Low eccentricities cannot arise from these physical conditions. For gas fractions below $20$ per cent, the ISM becomes less turbulent, with stellar velocity dispersions $<10~{\rm km\, s}^{-1}$, and nearly circular orbits for young stars. The turbulence present in gas-rich high redshift galaxies hence acts as a `barrier' against the formation of thin discs. We compare our findings to the Milky Way's age-velocity dispersion relation and argue that velocity dispersions imprinted already at star formation by the ISM contribute significantly at all times. Finally, we show that observed orbital eccentricities in the Milky Way's thick and thin discs can be explained entirely as imprints by the star-forming ISM, rather than by mergers or secular processes.

preprint2020arXiv

EDGE: From quiescent to gas-rich to star-forming low-mass dwarf galaxies

We study how star formation is regulated in low-mass field dwarf galaxies ($10^5 \leq M_{\star} \leq 10^6 \, \text{M}_{\odot}$), using cosmological high-resolution ($3 \, \text{pc}$) hydrodynamical simulations. Cosmic reionization quenches star formation in all our simulated dwarfs, but three galaxies with final dynamical masses of $3 \times 10^{9} \,\text{M}_{\odot}$ are subsequently able to replenish their interstellar medium by slowly accreting gas. Two of these galaxies re-ignite and sustain star formation until the present day at an average rate of $10^{-5} \, \text{M}_{\odot} \, \text{yr}^{-1}$, highly reminiscent of observed low-mass star-forming dwarf irregulars such as Leo T. The resumption of star formation is delayed by several billion years due to residual feedback from stellar winds and Type Ia supernovae; even at $z=0$, the third galaxy remains in a temporary equilibrium with a large gas content but without any ongoing star formation. Using the "genetic modification'' approach, we create an alternative mass growth history for this gas-rich quiescent dwarf and show how a small $(0.2\,\mathrm{dex})$ increase in dynamical mass can overcome residual stellar feedback, re-igniting star formation. The interaction between feedback and mass build-up produces a diversity in the stellar ages and gas content of low-mass dwarfs, which will be probed by combining next-generation HI and imaging surveys.

preprint2020arXiv

How runaway stars boost galactic outflows

Roughly ten per cent of OB stars are kicked out of their natal clusters before ending their life as supernovae. These so-called runaway stars can travel hundreds of parsecs into the low-density interstellar medium, where momentum and energy from stellar feedback is efficiently deposited. In this work we explore how this mechanism affects large scale properties of the galaxy, such as outflows. To do so we use a new model which treats OB stars and their associated feedback processes on a star-by-star basis. With this model we compare two hydrodynamical simulations of Milky Way-like galaxies, one where we include runaways, and one where we ignore them. Including runaway stars leads to twice as many supernovae explosions in regions with gas densities ranging from 1e-5 cm^-3 to 1e-3 cm^-3. This results in more efficient heating of the inter-arm regions, and drives strong galactic winds with mass loading factors boosted by up to one order of magnitude. These outflows produce a more massive and extended multi-phase circumgalactic medium, as well as a population of dense clouds in the halo. Conversely, since less energy and momentum is released in the dense star forming regions, the cold phase of the interstellar medium is less disturbed by feedback effects.

preprint2020arXiv

Rapid Filamentary Accretion as the Origin of Extended Thin Discs

Galactic outflows driven by stellar feedback are crucial for explaining the inefficiency of star formation in galaxies. Although strong feedback can promote the formation of galactic discs by limiting star formation at early times and removing low angular momentum gas, it is not understood how the same feedback can result in diverse objects such as elliptical galaxies or razor thin spiral galaxies. We investigate this problem using cosmological zoom-in simulations of two galaxies forming within $10^{12}~\mathrm{M_\odot}$ halos with almost identical mass accretion histories and halo spin parameters. However, the two resulting galaxies end up with very different bulge-to-disc ratios at $z = 0$. At $z>1.5$, the two galaxies feature a surface density of star formation $Σ_{\rm SFR}\simeq 10~\mathrm{M_\odot}~{\rm yr}^{-1}~{\rm kpc}^{-2}$, leading to strong outflows. After the last starburst episode, both galaxies feature a dramatic gaseous disc growth from 1~kpc to 5~kpc during 1~Gyr, a decisive event we dub "the Grand Twirl". After this event, the evolutionary tracks diverge strongly, with one galaxy ending up as a bulge-dominated galaxy, whereas the other ends up as a disc-dominated galaxy. The origins of this dichotomy are the angular momentum of the accreted gas, and whether it adds constructively to the initial disc angular momentum. The build-up of this extended disc leads to a rapid lowering of $Σ_{\rm SFR}$ by over two orders of magnitude with $Σ_{\rm SFR} \lesssim 0.1~\mathrm{M_\odot}~{\rm yr}^{-1}~{\rm kpc}^{-2}$, in remarkable agreement with what is derived from Milky Way stellar populations. As a consequence, supernovae explosions are spread out and cannot launch galactic outflows anymore, allowing for the persistence of a thin, gently star forming, extended disc.

preprint2019arXiv

Constraining churning and blurring in the Milky Way using large spectroscopic surveys -- an exploratory study

We have investigated the possibilities to quantify how much stars move in the Milky Way stellar disk due to diffuse processes (i.e. so called blurring) and due to influences from spiral arms and the bar (i.e. so called churning). To this end we assume that it is possible to infer the formation radius of a star if we know their elemental abundances and age as well as the metallicity profile of the interstellar medium at the time of the formation of the star. Using this information, coupled with orbital information derived from Gaia DR2 data and radial velocities from large spectroscopic surveys, we show that it is possible to isolate stellar samples such that we can start to quantify how important the role of churning is. We use data from APOGEE DR14, parallaxes from Gaia and stellar ages based on C and N elemental abundances in the stars. In our sample, we find that about half of the stars have experienced some sort of radial migration (based solely on their orbital properties), 10 % have likely have suffered only from churning, whilst a modest 5-7 % of stars have never experienced either churning or blurring making them ideal tracers of the original properties of the cool stellar disk. Our investigation shows that it is possible to put up a framework where we can begin to quantify churning and blurring an important. Important aspects for future work would include to investigate how selection effects should be accounted for.

preprint2016arXiv

Column Density Profiles of Multi-Phase Gaseous Halos

We analyze circumgalactic medium (CGM) in a suite of high-resolution cosmological re-simulations of a Milky-Way size galaxy and show that CGM properties are quite sensitive to details of star formation--feedback loop modelling. The simulation that produces a realistic late-type galaxy, fails to reproduce existing observations of the CGM. In contrast, simulation that does not produce a realistic galaxy has the predicted CGM in better agreement with observations. This illustrates that properties of galaxies and properties of their CGM provide strong ${\it complementary}$ constraints on the processes governing galaxy formation. Our simulations predict that column density profiles of ions are well described by an exponential function of projected distance $d$: $N \propto e^{-d/h_s}$. Simulations thus indicate that the sharp drop in absorber detections at larger distances in observations does not correspond to a "boundary" of an ion, but reflects the underlying steep exponential column density profile. Furthermore, we find that ionization energy of ions is tightly correlated with the scale height $h_s$: $h_s \propto E_{\rm ion}^{0.74}$. At $z \approx 0$, warm gas traced by low-ionization species (e.g., Mg II and C IV) has $ h_s \approx 0.03-0.07 R_{\rm vir}$, while higher ionization species (O VI and Ne VIII) have $h_s \approx 0.32-0.45R_{\rm vir}$. Finally, the scale heights of ions in our simulations evolve slower than the virial radius for $z\leq 2$, but similarly to the halo scale radius, $r_s$. Thus, we suggest that the column density profiles of galaxies at different redshifts should be scaled by $r_s$ rather than the halo virial radius.

preprint2016arXiv

The origin of the Milky Way globular clusters

We present a cosmological zoom-in simulation of a Milky Way-like galaxy used to explore the formation and evolution of star clusters. We investigate in particular the origin of the bimodality observed in the colour and metallicity of globular clusters, and the environmental evolution through cosmic times in the form of tidal tensors. Our results self-consistently confirm previous findings that the blue, metal-poor clusters form in satellite galaxies which are accreted onto the Milky Way, while the red, metal-rich clusters form mostly in situ or, to a lower extent in massive, self-enriched galaxies merging with the Milky Way. By monitoring the tidal fields these populations experience, we find that clusters formed in situ (generally centrally concentrated) feel significantly stronger tides than the accreted ones, both in the present-day, and when averaged over their entire life. Furthermore, we note that the tidal field experienced by Milky Way clusters is significantly weaker in the past than at present-day, confirming that it is unlikely that a power-law cluster initial mass function like that of young massive clusters, is transformed into the observed peaked distribution in the Milky Way with relaxation-driven evaporation in a tidal field.

preprint2016arXiv

The roles of stellar feedback and galactic environment in star forming molecular clouds

Feedback from massive stars is thought to play an important role in the evolution of molecular clouds. In this work we analyse the effects of stellar winds and supernovae (SNe) in the evolution of two massive ($\sim 10^6\,M_\odot$) giant molecular clouds (GMCs): one gravitationally bound collapsing cloud and one unbound cloud undergoing disruption by galactic shear. These two clouds have been extracted from a large scale galaxy model and are re-simulated at a spatial resolution of $\sim 0.01$ pc, including feedback from winds, SNe, and the combined effect of both. We find that stellar winds stop accretion of gas onto sink particles, and can also trigger star formation in the shells formed by the winds, although the overall effect is to reduce the global star formation rate of both clouds. Furthermore, we observe that winds tend to escape through the corridors of diffuse gas. The effect of SNe is not so prominent and the star formation rate is similar to models neglecting stellar feedback. We find that most of the energy injected by the SNe is radiated away, but overdense areas are created by multiple and concurrent SN events especially in the most virialised cloud. Our results suggest that the impact of stellar feedback is sensitive to the morphology of star forming clouds, which is set by large scale galactic flows, being of greater importance in clouds undergoing gravitational collapse.

preprint2015arXiv

Characterizing gravitational instability in turbulent multi-component galactic discs

Gravitational instabilities play an important role in galaxy evolution and in shaping the interstellar medium (ISM). The ISM is observed to be highly turbulent, meaning that observables like the gas surface density and velocity dispersion depend on the size of the region over which they are measured. In this work we investigate, using simulations of Milky Way-like disc galaxies with a resolution of $\sim 9$ pc, the nature of turbulence in the ISM and how this affects the gravitational stability of galaxies. By accounting for the measured average turbulent scalings of the density and velocity fields in the stability analysis, we can more robustly characterize the average level of stability of the galaxies as a function of scale, and in a straightforward manner identify scales prone to fragmentation. Furthermore, we find that the stability of a disc with feedback-driven turbulence can be well described by a "Toomre-like" $Q$ stability criterion on all scales, whereas the classical $Q$ can formally lose its meaning on small scales if violent disc instabilities occur in models lacking pressure support from stellar feedback.

preprint2015arXiv

Galaxies that Shine: radiation-hydrodynamical simulations of disk galaxies

Radiation feedback is typically implemented using subgrid recipes in hydrodynamical simulations of galaxies. Very little work has so far been performed using radiation-hydrodynamics (RHD), and there is no consensus on the importance of radiation feedback in galaxy evolution. We present RHD simulations of isolated galaxy disks of different masses with a resolution of 18 pc. Besides accounting for supernova feedback, our simulations are the first galaxy-scale simulations to include RHD treatments of photo-ionisation heating and radiation pressure, from both direct optical/UV radiation and multi-scattered, re-processed infrared (IR) radiation. Photo-heating smooths and thickens the disks and suppresses star formation about as much as the inclusion of ("thermal dump") supernova feedback does. These effects decrease with galaxy mass and are mainly due to the prevention of the formation of dense clouds, as opposed to their destruction. Radiation pressure, whether from direct or IR radiation, has little effect, but for the IR radiation we show that its impact is limited by our inability to resolve the high optical depths for which multi-scattering becomes important. While artificially boosting the IR optical depths does reduce the star formation, it does so by smoothing the gas rather than by generating stronger outflows. We conclude that although higher-resolution simulations, and potentially also different supernova implementations, are needed for confirmation, our findings suggest that radiation feedback is more gentle and less effective than is often assumed in subgrid prescriptions.

preprint2015arXiv

Multiple populations in globular clusters: the distinct kinematic imprints of different formation scenarios

Several scenarios have been proposed to explain the presence of multiple stellar populations in globular clusters. Many of them invoke multiple generations of stars to explain the observed chemical abundance anomalies, but it has also been suggested that self-enrichment could occur via accretion of ejecta from massive stars onto the circumstellar disc of low-mass pre-main sequence stars. These scenarios imply different initial conditions for the kinematics of the various stellar populations. Given some net angular momentum initially, models for which a second generation forms from gas that collects in a cooling flow into the core of the cluster predict an initially larger rotational amplitude for the polluted stars compared to the pristine stars. This is opposite to what is expected from the accretion model, where the polluted stars are the ones crossing the core and are on preferentially radial (low-angular momentum) orbits, such that their rotational amplitude is lower. Here we present the results of a suite of $N$-body simulations with initial conditions chosen to capture the distinct kinematic properties of these pollution scenarios. We show that initial differences in the kinematics of polluted and pristine stars can survive to the present epoch in the outer parts of a large fraction of Galactic globular clusters. The differential rotation of pristine and polluted stars is identified as a unique kinematic signature that could allow us to distinguish between various scenarios, while other kinematic imprints are generally very similar from one scenario to the other.

preprint2015arXiv

Novel Adaptive softening for collisionless N-body simulations: Eliminating spurious halos

We describe a NOVel form of Adaptive softening (NovA) for collisionless $N$-body simulations, implemented in the Ramses adaptive mesh refinement code. We introduce a refinement criterion that the particle distribution within each cell be sufficiently isotropic, as measured by its moment of inertia tensor. In this way, collapse is only refined if it occurs along all three axes, ensuring that the softening $ε$ is always of order twice the largest inter-particle spacing in a cell. This more conservative force softening criterion is designed to minimise spurious two-body effects, while maintaining high force resolution in collapsed regions of the flow. We test NovA using an antisymmetric perturbed plane wave collapse (`Valinia' test) before applying it to warm dark matter (WDM) simulations. For the Valinia test, we show that -- unlike the standard $N$-body method -- NovA produces no numerical fragmentation while still being able to correctly capture fine caustics and shells around the collapsing regions. For the WDM simulations, we find that NovA converges significantly more rapidly than standard $N$-body, producing little or no spurious halos on small scales. We show, however, that determining whether or not halos exist below the free streaming mass $M_{\rm fs}$ is complicated by the fact that our halo finder (AHF) likely incorrectly labels some caustics and criss-crossing filaments as halos, while one or two particularly massive filaments appear to fragment in any version of NovA where refinement is allowed. Such massive filaments may be physically unstable to collapse, as is the case for infinite, static, self-gravitating cylinders. We will use NovA in forthcoming papers to study the issue of halo formation below $M_{\rm fs}$; filament stability; and to obtain new constraints on the temperature of dark matter.

preprint2015arXiv

On the interplay between star formation and feedback in galaxy formation simulations

We investigate the star formation-feedback cycle in cosmological galaxy formation simulations, focusing on progenitors of Milky Way (MW)-sized galaxies. We find that in order to reproduce key properties of the MW progenitors, such as semi-empirically derived star formation histories and the shape of rotation curves, our implementation of star formation and stellar feedback requires 1) a combination of local early momentum feedback via radiation pressure and stellar winds and subsequent efficient supernovae feedback, and 2) efficacy of feedback that results in self-regulation of the global star formation rate on kiloparsec scales. We show that such feedback-driven self-regulation is achieved globally for a local star formation efficiency per free fall time of $ε_{\rm ff}\approx 10\%$. Although this value is larger that the $ε_{\rm ff}\sim 1\%$ value usually inferred from the Kennicutt-Schmidt (KS) relation, we show that it is consistent with direct observational estimates of $ε_{\rm ff}$ in molecular clouds. Moreover, we show that simulations with local efficiency of $ε_{\rm ff}\approx 10\%$ reproduce the global observed KS relation. Such simulations also reproduce the cosmic star formation history of the Milky Way sized galaxies and satisfy a number of other observational constraints. Conversely, we find that simulations that a priori assume an inefficient mode of star formation, instead of achieving it via stellar feedback regulation, fail to produce sufficiently vigorous outflows and do not reproduce observations. This illustrates the importance of understanding the complex interplay between star formation and feedback and the detailed processes that contribute to the feedback-regulated formation of galaxies.

preprint2014arXiv

A Systematic Look at the Effects of Radiative Feedback on Disc Galaxy Formation

Galaxy formation models and simulations rely on various feedback mechanisms to reproduce the observed baryonic scaling relations and galaxy morphologies. Although dwarf galaxy and giant elliptical properties can be explained using feedback from supernova and active galactic nuclei, Milky Way-sized galaxies still represent a challenge to current theories of galaxy formation. In this paper, we explore the possible role of feedback from stellar radiation in regulating the main properties of disk galaxies such as our own Milky Way. We have performed a suite of cosmological simulations of the same $\sim10^{12} {\rm M}_{\odot}$ halo selected based on its rather typical mass accretion history. We have implemented radiative feedback from young stars using a crude model of radiative transfer for ultraviolet (UV) and infrared (IR) radiation. However, the model is realistic enough such that the dust opacity plays a direct role in regulating the efficiency of our feedback mechanism. We have explored various models for the dust opacity, assuming different constant dust temperatures, as well as a varying dust temperature model. We find that while strong radiative feedback appears as a viable mechanism to regulate the stellar mass fraction in massive galaxies, it also prevents the formation of discs with reasonable morphologies. In models with strong stellar radiation feedback, stellar discs are systematically too thick while the gas disc morphology is completely destroyed due to the efficient mixing between the feedback-affected gas and its surroundings. At the resolution of our simulation suite, we find it impossible to preserve spiral disc morphology while at the same time expelling enough baryons to satisfy the abundance matching constraints.

preprint2014arXiv

Larson's scaling laws, and the gravitational instability of clumpy discs at high redshift

Gravitational instabilities play a primary role in shaping the clumpy structure and powering the star formation activity of gas-rich high-redshift galaxies. Here we analyse the stability of such systems, focusing on the size and mass ranges of unstable regions in the disc. Our analysis takes into account the mass-size and linewidth-size scaling relations observed in molecular gas, originally discovered by Larson. We show that such relations can have a strong impact on the size and mass of star-forming clumps, as well as on the stability properties of the disc at all observable scales, making the classical Toomre parameter a highly unreliable indicator of gravitational instability. For instance, a disc with Q=1 can be far from marginal instability, while a disc with Q<<1 can be marginally unstable. Our work raises an important caveat: if clumpy discs at high redshift have scale-dependent surface densities and velocity dispersions, as implied by the observed clump scaling relations, then we cannot thoroughly understand their stability and star formation properties unless we perform multi-scale observations. This will soon be possible thanks to dedicated ALMA surveys, which will explore the physical properties of super-giant molecular clouds at the peak of cosmic star formation and beyond.

preprint2014arXiv

The Smith Cloud and its dark matter halo: Survival of a Galactic disc passage

The current velocity of the Smith Cloud indicates that it has undergone at least one passage of the Galactic disc. Using hydrodynamic simulations we examine the present day structure of the Smith Cloud. We find that a dark matter supported cloud is able to reproduce the observed present day neutral hydrogen mass, column density distribution and morphology. In this case the dark matter halo becomes elongated, owing to the tidal interaction with the Galactic disc. Clouds in models neglecting dark matter confinement are destroyed upon disc passage, unless the initial cloud mass is well in excess of what is observed today. We then determine integrated flux upper limits to the gamma-ray emission around such a hypothesised dark matter core in the Smith Cloud. No statistically significant core or extended gamma-ray emission are detected down to a 95% confidence level upper limit of $1.4\times10^{-10}$ ph cm$^{-2}$ s$^{-1}$ in the 1-300 GeV energy range. For the derived distance of 12.4 kpc, the Fermi upper limits set the first tentative constraints on the dark matter cross sections annihilating into $τ^{+}τ^{-}$ and $b\bar{b}$ for a high-velocity cloud.

preprint2013arXiv

Simulations of disk galaxies with cosmic ray driven galactic winds

We present results from high-resolution hydrodynamic simulations of isolated SMC- and Milky Way-sized galaxies that include a model for feedback from galactic cosmic rays (CRs). We find that CRs are naturally able to drive winds with mass loading factors of up to ~10 in dwarf systems. The scaling of the mass loading factor with circular velocity between the two simulated systems is consistent with \propto v_c^{1-2} required to reproduce the faint end of the galaxy luminosity function. In addition, simulations with CR feedback reproduce both the normalization and the slope of the observed trend of wind velocity with galaxy circular velocity. We find that winds in simulations with CR feedback exhibit qualitatively different properties compared to SN driven winds, where most of the acceleration happens violently in situ near star forming sites. In contrast, the CR-driven winds are accelerated gently by the large-scale pressure gradient established by CRs diffusing from the star-forming galaxy disk out into the halo. The CR-driven winds also exhibit much cooler temperatures and, in the SMC-sized system, warm (T~10^4 K) gas dominates the outflow. The prevalence of warm gas in such outflows may provide a clue as to the origin of ubiquitous warm gas in the gaseous halos of galaxies detected via absorption lines in quasar spectra.

preprint2012arXiv

Towards a complete accounting of energy and momentum from stellar feedback in galaxy formation simulations

Stellar feedback plays a key role in galaxy formation by regulating star formation, driving interstellar turbulence and generating galactic scale outflows. Although modern simulations of galaxy formation can resolve scales of 10-100 pc, star formation and feedback operate on smaller, "subgrid" scales. Great care should therefore be taken in order to properly account for the effect of feedback on global galaxy evolution. We investigate the momentum and energy budget of feedback during different stages of stellar evolution, and study its impact on the interstellar medium using simulations of local star forming regions and galactic disks at the resolution affordable in modern cosmological zoom-in simulations. In particular, we present a novel subgrid model for the momentum injection due to radiation pressure and stellar winds from massive stars during early, pre-supernova evolutionary stages of young star clusters. Early injection of momentum acts to clear out dense gas in star forming regions, hence limiting star formation. The reduced gas density mitigates radiative losses of thermal feedback energy from subsequent supernova explosions, leading to an increased overall efficiency of stellar feedback. The detailed impact of stellar feedback depends sensitively on the implementation and choice of parameters. Somewhat encouragingly, we find that implementations in which feedback is efficient lead to approximate self-regulation of global star formation efficiency. We compare simulation results using our feedback implementation to other phenomenological feedback methods, where thermal feedback energy is allowed to dissipate over time scales longer than the formal gas cooling time. We find that simulations with maximal momentum injection suppress star formation to a similar degree as is found in simulations adopting adiabatic thermal feedback.

preprint2010arXiv

A Toomre-like stability criterion for the clumpy and turbulent interstellar medium

We explore the gravitational instability of clumpy and turbulent gas discs, taking into account the Larson-type scaling laws observed in giant molecular clouds (GMCs) and HI, as well as more general scaling relations. This degree of freedom is of special interest in view of the coming high-z ISM surveys, and is thus potentially important for understanding the dynamical effects of turbulence at all epochs of galaxy evolution. Our analysis shows that turbulence has a deep impact on the gravitational instability of the disc. It excites a rich variety of stability regimes, several of which have no classical counterpart. Among other diagnostics, we provide two useful tools for observers and simulators: (1) the stability map of turbulence, which illustrates our stability scenario and relates it to the phenomenology of interstellar turbulence: GMC/HI observations, simulations and models; (2) a Toomre-like stability criterion, $Q\geq\bar{Q}$, which applies to a large class of clumpy/turbulent discs. We make specific predictions about GMC and cold-HI turbulence, and point out the implications of our analysis for high-z galaxy surveys.

preprint2010arXiv

Systematic uncertainties in the determination of the local dark matter density

A precise determination of the local dark matter density and an accurate control over the corresponding uncertainties are of paramount importance for Dark Matter (DM) searches. Using very recent high-resolution numerical simulations of a Milky Way like object, we study the systematic uncertainties that affect the determination of the local dark matter density based on dynamical measurements in the Galaxy. In particular, extracting from the simulation with baryons the orientation of the Galactic stellar disk with respect to the DM distribution, we study the DM density for an observer located at $\sim$8 kpc from the Galactic center {\it on the stellar disk}, $ρ_0$. This quantity is found to be always larger than the average density in a spherical shell of same radius $\barρ_0$, which is the quantity inferred from dynamical measurements in the Galaxy, and to vary in the range $ρ_0/\barρ_0=1.01-1.41$. This suggests that the actual dark matter density in the solar neighbourhood is on average 21\% larger than the value inferred from most dynamical measurements, and that the associated systematic errors are larger than the statistical errors recently discussed in the literature.

preprint2010arXiv

The formation of disc galaxies in a LCDM universe

We study the formation of disc galaxies in a fully cosmological framework using adaptive mesh refinement simulations. We perform an extensive parameter study of the main subgrid processes that control how gas is converted into stars and the coupled effect of supernovae feedback. We argue that previous attempts to form disc galaxies have been unsuccessful because of the universal adoption of strong feedback combined with high star formation efficiencies. Unless extreme amounts of energy are injected into the interstellar medium during supernovae events, these star formation parameters result in bulge dominated S0/Sa galaxies as star formation is too efficient at z~3. We show that a low efficiency of star-formation more closely models the subparsec physical processes, especially at high redshift. We highlight the successful formation of extended disc galaxies with scale lengths r_d=4-5 kpc, flat rotation curves and bulge to disc ratios of B/D~1/4. Not only do we resolve the formation of a Milky Way-like spiral galaxy, we also observe the secular evolution of the disc as it forms a pseudo-bulge. The disc properties agree well with observations and are compatible with the photometric and baryonic Tully-Fisher relations, the Kennicutt-Schmidt relation and the observed angular momentum content of spiral galaxies. We conclude that underlying small-scale star formation physics plays a larger role than previously considered in simulations of galaxy formation.

preprint2009arXiv

Disk formation and the origin of clumpy galaxies at high redshift

Observations of high redshift galaxies have revealed a multitude of large clumpy rapidly star-forming galaxies. Their formation scenario and their link to present day spirals is still unknown. In this Letter we perform adaptive mesh refinement simulations of disk formation in a cosmological context that are unrivalled in terms of mass and spatial resolution. We find that the so called "chain-galaxies" and "clump-clusters" are a natural outcome of early epochs of enhanced gas accretion from cold dense streams as well as tidally and ram-pressured stripped material from minor mergers and satellites. Through interaction with the hot halo gas, this freshly accreted cold gas settles into a large disk-like system, not necessarily aligned to an older stellar component, that undergoes fragmentation and subsequent star formation, forming large clumps in the mass range 10^7-10^9 M_sun. Galaxy formation is a complex process at this important epoch when most of the central baryons are being acquired through a range of different mechanisms - we highlight that a rapid mass loading epoch is required to fuel the fragmentation taking place in the massive arms in the outskirts of extended disks, an accretion mode that occurs naturally in the hierarchical assembly process at early epochs.