-
Identification of Galaxy Protoclusters Based on the Spherical Top-hat Collapse Theory
Authors:
Jaehyun Lee,
Changbom Park,
Juhan Kim,
Christophe Pichon,
Brad K. Gibson,
Jihye Shin,
Yonghwi Kim,
Owain N. Snaith,
Yohan Dubois,
C. Gareth Few
Abstract:
We propose a new method for finding galaxy protoclusters that is motivated by structure formation theory and also directly applicable to observations. We adopt the conventional definition that a protocluster is a galaxy group whose virial mass $M_{\rm vir} < M_{\rm cl}$ at its epoch, where $M_{\rm cl}=10^{14}\,M_{\odot}$, but would exceed that limit when it evolves to $z=0$. We use the critical ov…
▽ More
We propose a new method for finding galaxy protoclusters that is motivated by structure formation theory and also directly applicable to observations. We adopt the conventional definition that a protocluster is a galaxy group whose virial mass $M_{\rm vir} < M_{\rm cl}$ at its epoch, where $M_{\rm cl}=10^{14}\,M_{\odot}$, but would exceed that limit when it evolves to $z=0$. We use the critical overdensity for complete collapse at $z = 0$ predicted by the spherical top-hat collapse model to find the radius and total mass of the regions that would collapse at $z=0$. If the mass of a region centered at a massive galaxy exceeds $M_{\rm cl}$, the galaxy is at the center of a protocluster. We define the outer boundary of protocluster as the zero-velocity surface at the turnaround radius so that the member galaxies are those sharing the same protocluster environment and showing some conformity in physical properties. We use the cosmological hydrodynamical simulation Horizon Run 5 (HR5) to calibrate this prescription and demonstrate its performance. We find that the protocluster identification method suggested in this study is quite successful. Its application to the high-redshift HR5 galaxies shows a tight correlation between the mass within the protocluster regions identified according to the spherical collapse model and the final mass to be found within the clusters at $z=0$, meaning that the regions can be regarded as the bona fide protoclusters with high reliability. We also confirm that the redshift-space distortion does not significantly affect the performance of the protocluster identification scheme.
△ Less
Submitted 29 November, 2023; v1 submitted 1 August, 2023;
originally announced August 2023.
-
Low-Surface-Brightness Galaxies are missing in the observed Stellar Mass Function
Authors:
Juhan Kim,
Jaehyun Lee,
Clotilde Laigle,
Yohan Dubois,
Yonghwi Kim,
Changbom Park,
Christophe Pichon,
Brad Gibson,
C. Gareth Few,
Jihye Shin,
Owain Snaith
Abstract:
We investigate the impact of the surface brightness (SB) limit on the galaxy stellar mass functions (GSMFs) using mock surveys generated from the Horizon Run 5 (HR5) simulation. We compare the stellar-to-halo-mass relation, GSMF, and size-stellar mass relation of the HR5 galaxies with empirical data and other cosmological simulations. The mean SB of simulated galaxies are computed using their effe…
▽ More
We investigate the impact of the surface brightness (SB) limit on the galaxy stellar mass functions (GSMFs) using mock surveys generated from the Horizon Run 5 (HR5) simulation. We compare the stellar-to-halo-mass relation, GSMF, and size-stellar mass relation of the HR5 galaxies with empirical data and other cosmological simulations. The mean SB of simulated galaxies are computed using their effective radii, luminosities, and colors. To examine the cosmic SB dimming effect, we compute $k$-corrections from the spectral energy distributions of individual simulated galaxy at each redshift, apply the $k$-corrections to the galaxies, and conduct mock surveys based on the various SB limits. We find that the GSMFs are significantly affected by the SB limits at a low-mass end. This approach can ease the discrepancy between the GSMFs obtained from simulations and observations at $0.625\le z\le 2$. We also find that a redshift survey with a SB selection limit of $\left<μ_r\right>^e =$ 28 mag arcsec${}^{-2}$ will miss 20% of galaxies with $M_\star^g=10^{9}~{\rm M_\odot}$ at $z=0.625$. The missing fraction of low-surface-brightness galaxies increases to 50%, 70%, and 98% at $z=0.9$, 1.1, and 1.9, respectively, at the SB limit.
△ Less
Submitted 2 May, 2023; v1 submitted 29 December, 2022;
originally announced December 2022.
-
Formation and Morphology of the First Galaxies in the Cosmic Morning
Authors:
Changbom Park,
Jaehyun Lee,
Juhan Kim,
Donghui Jeong,
Christophe Pichon,
Brad K. Gibson,
Owain N. Snaith,
Jihye Shin,
Yonghwi Kim,
Yohan Dubois,
C. Gareth Few
Abstract:
We investigate the formation and morphological evolution of the first galaxies in the cosmic morning ($10 \gtrsim z \gtrsim 4$) using the Horizon Run 5 (HR5) simulation. For galaxies above the stellar mass $M_{\star, {\rm min}} = 2\times 10^9\,M_{\odot}$, we classify them into disk, spheroid, and irregular types according to their asymmetry and stellar mass morphology. We find that about 2/3 of th…
▽ More
We investigate the formation and morphological evolution of the first galaxies in the cosmic morning ($10 \gtrsim z \gtrsim 4$) using the Horizon Run 5 (HR5) simulation. For galaxies above the stellar mass $M_{\star, {\rm min}} = 2\times 10^9\,M_{\odot}$, we classify them into disk, spheroid, and irregular types according to their asymmetry and stellar mass morphology. We find that about 2/3 of the galaxies have a Sérsic index $< 1.5$, reflecting the dominance of disk-type morphology in the cosmic morning. The rest are evenly distributed as incidental and transient irregulars or spheroids. These fractions are roughly independent of redshift and stellar mass up to $\sim10^{10}\,M_{\odot}$. Almost all the first galaxies with $M_{\star}> M_{\star, {\rm min}}$ at $z>4$ form at initial peaks of the matter density field. Large-scale structures in the universe emerge and grow like cosmic rhizomes as the underlying matter density fluctuations grow and form associations of galaxies in rare overdense regions and the realm of the galactic world is stretched into relatively lower-density regions along evolving filaments. The cosmic web of galaxies forms at lower redshifts when most rhizomes globally percolate. The primordial angular momentum produced by the induced tidal torques on protogalactic regions is correlated with the internal kinematics of galaxies and tightly aligned with the angular momentum of the total galaxy mass. The large-scale tidal field imprinted in the initial conditions seems responsible for the dominance of disk morphology, and for the tendency of galaxies to re-acquire a disk post-distortion.
△ Less
Submitted 29 July, 2022; v1 submitted 24 February, 2022;
originally announced February 2022.
-
The Horizon Run 5 Cosmological Hydrodynamic Simulation: Probing Galaxy Formation from Kilo- to Giga-parsec Scales
Authors:
Jaehyun Lee,
Jihye Shin,
Owain N. Snaith,
Yonghwi Kim,
C. Gareth Few,
Julien Devriendt,
Yohan Dubois,
Leah M. Cox,
Sungwook E. Hong,
Oh-Kyoung Kwon,
Chan Park,
Christophe Pichon,
Juhan Kim,
Brad K. Gibson,
Changbom Park
Abstract:
Horizon Run 5 (HR5) is a cosmological hydrodynamical simulation which captures the properties of the Universe on a Gpc scale while achieving a resolution of 1kpc. Inside the simulation box we zoom-in on a high-resolution cuboid region with a volume of $1049\times119\times127\,{\rm cMpc}^3$.The sub-grid physics chosen to model galaxy formation includes radiative heating/cooling, UV background, star…
▽ More
Horizon Run 5 (HR5) is a cosmological hydrodynamical simulation which captures the properties of the Universe on a Gpc scale while achieving a resolution of 1kpc. Inside the simulation box we zoom-in on a high-resolution cuboid region with a volume of $1049\times119\times127\,{\rm cMpc}^3$.The sub-grid physics chosen to model galaxy formation includes radiative heating/cooling, UV background, star formation, supernova feedback, chemical evolution tracking the enrichment of oxygen and iron, the growth of supermassive black holes and feedback from active galactic nuclei (AGN) in the form of a dual jet-heating mode. For this simulation we implemented a hybrid MPI-OMP version of RAMSES, specifically targeted for modern many-core many thread parallel architectures. In addition to the traditional simulation snapshots, light-cone data was generated on the fly. For the post-processing, we extended the Friends-of-Friend (FoF) algorithm and developed a new galaxy finder PGalF to analyse the outputs of HR5. The simulation successfully reproduces observations, such as the cosmic star formation history and connectivity of galaxy distribution, We identify cosmological structures at a wide range of scales, from filaments with a length of several cMpc, to voids with a radius of ~100 cMpc. The simulation also indicates that hydrodynamical effects on small scales impact galaxy clustering up to very large scales near and beyond the baryonic acoustic oscillation (BAO) scale. Hence, caution should be taken when using that scale as a cosmic standard ruler: one needs to carefully understand the corresponding biases. The simulation is expected to be an invaluable asset for the interpretation of upcoming deep surveys of the Universe.
△ Less
Submitted 5 October, 2023; v1 submitted 1 June, 2020;
originally announced June 2020.
-
High-resolution three-dimensional simulations of gas removal from ultrafaint dwarf galaxies. I. Stellar feedback
Authors:
Donatella Romano,
Francesco Calura,
Annibale D'Ercole,
C. Gareth Few
Abstract:
The faintest Local Group galaxies found lurking in and around the Milky Way halo provide a unique test bed for theories of structure formation and evolution on small scales. Deep Subaru and Hubble Space Telescope photometry demonstrates that their stellar populations are old, and that the star formation activity did not last longer than 2 Gyr in these systems. A few mechanisms that may lead to suc…
▽ More
The faintest Local Group galaxies found lurking in and around the Milky Way halo provide a unique test bed for theories of structure formation and evolution on small scales. Deep Subaru and Hubble Space Telescope photometry demonstrates that their stellar populations are old, and that the star formation activity did not last longer than 2 Gyr in these systems. A few mechanisms that may lead to such a rapid quenching have been investigated by means of hydrodynamic simulations, without providing any final assessment so far. This is the first in a series of papers aimed at analysing the roles of stellar feedback, ram pressure stripping, host-satellite tidal interactions and reionization in cleaning the lowest-mass Milky Way companions of their cold gas, by using high-resolution, three-dimensional hydrodynamic simulations. We simulate an isolated ultrafaint dwarf galaxy loosely modeled after Bootes I, and examine whether or not stellar feedback alone could drive a substantial fraction of the ambient gas out from the shallow potential well. In contrast to simple analytical estimates, but in agreement with previous hydrodynamical studies, we find that most of the cold gas reservoir is retained. Conversely, a significant fraction of the metal-enriched stellar ejecta crosses the boundaries of the computational box with velocities exceeding the local escape velocity and is, thus, likely lost from the system. Although the total energy output from multiple supernova explosions exceeds the binding energy of the gas, no galactic-scale outflow develops in our simulations and as such, most of the ambient medium remains trapped within the weak potential well of the model galaxy. It seems thus unavoidable that, in order to explain the dearth of gas in ultrafaint dwarf galaxies, we will have to resort to environmental effects. This will be the subject of a forthcoming paper.
△ Less
Submitted 5 September, 2019;
originally announced September 2019.
-
The time evolution of the Milky Way's oxygen abundance gradient
Authors:
M. Mollá,
Á. I. Díaz,
O. Cavichia,
B. K. Gibson,
W. J. Maciel,
R. D. D. Costa,
Y. Ascasibar,
C. G. Few
Abstract:
We study the evolution of oxygen abundance radial gradients as a function of time for the Milky Way Galaxy obtained with our {\sc Mulchem} chemical evolution model. We review the recent data of abundances for different objects observed in our Galactic disc. We analyse with our models the role of the growth of the stellar disc, as well as the effect of infall rate and star formation prescriptions,…
▽ More
We study the evolution of oxygen abundance radial gradients as a function of time for the Milky Way Galaxy obtained with our {\sc Mulchem} chemical evolution model. We review the recent data of abundances for different objects observed in our Galactic disc. We analyse with our models the role of the growth of the stellar disc, as well as the effect of infall rate and star formation prescriptions, or the pre-enrichment of the infall gas, on the time evolution of the oxygen abundance radial distribution. We compute the radial gradient of abundances within the {\sl disk}, and its corresponding evolution, taking into account the disk growth along time. We compare our predictions with the data compilation, showing a good agreement. Our models predict a very smooth evolution when the radial gradient is measured within the optical disc with a slight flattening of the gradient from $\sim -0.057$\,dex\,kpc$^{-1}$ at $z=4$ until values around $\sim -0.015$\,dex\,kpc$^{-1}$ at $z=1$ and basically the same gradient until the present, with small differences between models. Moreover, some models show a steepening at the last times, from $z=1$ until $z=0$ in agreement with data which give a variation of the gradient in a range from $-0.02$ to $-0.04$\,de\,kpc$^{-1}$ from $t=10$\,Gyr until now. The gradient measured as a function of the normalized radius $R/R_{\rm eff}$ is in good agreement with findings by CALIFA and MUSE, and its evolution with redshift falls within the error bars of cosmological simulations.
△ Less
Submitted 22 October, 2018;
originally announced October 2018.
-
Evolution of Starburst Galaxies in the Illustris Simulation
Authors:
C. L. Wilkinson,
K. A. Pimbblet,
J. P. Stott,
C. G. Few,
B. K. Gibson
Abstract:
There is a consensus in the literature that starburst galaxies are triggered by inter- action events. However, it remains an open question as to what extent both merging and non-merging interactions have in triggering starbursts? In this study, we make use of the Illustris simulation to test how different triggering mechanisms can effect starburst events. We examine star formation rate, colour and…
▽ More
There is a consensus in the literature that starburst galaxies are triggered by inter- action events. However, it remains an open question as to what extent both merging and non-merging interactions have in triggering starbursts? In this study, we make use of the Illustris simulation to test how different triggering mechanisms can effect starburst events. We examine star formation rate, colour and environment of starburst galaxies to determine if this could be why we witness a bimodality in post-starburst populations within observational studies. Further, we briefly test the extent of quenching due to AGN feedback. From Illustris, we select 196 starburst galaxies at z = 0.15 and split them into post-merger and pre-merger/harassment driven starburst samples. We find that 55% of this sample not undergone a merger in the past 2 Gyr. Both of our samples are located in low-density environments within the filament regions of the cosmic web, however we find that pre-merger/harassment driven starburst are in higher density environments than post-merger driven starbursts. We also find that pre-merger/harassment starbursts are redder than post-merger starbursts, this could be driven by environmental effects. Both however, produce nuclear starbursts of comparable strengths.
△ Less
Submitted 5 June, 2018;
originally announced June 2018.
-
The Role of Stellar Radial Motions in Shaping Galaxy Surface Brightness Profiles
Authors:
T. Ruiz-Lara,
C. G. Few,
E. Florido,
B. K. Gibson,
I. Pérez,
P. Sánchez-Blázquez
Abstract:
Aims: The physics driving features such as breaks observed in galaxy surface brightness (SB) profiles remains contentious. Here, we assess the importance of stellar radial motions in shaping their characteristics. Methods: We use the simulated Milky Way-mass, cosmological discs, from the Ramses Disc Environment Study (RaDES) to characterise the radial redistribution of stars in galaxies displaying…
▽ More
Aims: The physics driving features such as breaks observed in galaxy surface brightness (SB) profiles remains contentious. Here, we assess the importance of stellar radial motions in shaping their characteristics. Methods: We use the simulated Milky Way-mass, cosmological discs, from the Ramses Disc Environment Study (RaDES) to characterise the radial redistribution of stars in galaxies displaying type I (pure exponentials), II (downbending), and III (upbending) SB profiles. We compare radial profiles of the mass fractions and the velocity dispersions of different sub-populations of stars according to their birth and current locations. Results: Radial redistribution of stars is important in all galaxies regardless of their light profiles. Type II breaks seem to be a consequence of the combined effects of outward-moving and accreted stars. The former produces shallower inner profiles (lack of stars in the inner disc) and accumulate material around the break radius and beyond, strengthening the break; the latter can weaken or even convert the break into a pure exponential. Further accretion from satellites can concentrate material in the outermost parts, leading to type III breaks that can coexist with type II breaks, but situated further out. Type III galaxies would be the result of an important radial redistribution of material throughout the entire disc, as well as a concentration of accreted material in the outskirts. In addition, type III galaxies display the most efficient radial redistribution and the largest number of accreted stars, followed by type I and II systems, suggesting that type I galaxies may be an intermediate case between types II and III. In general, the velocity dispersion profiles of all galaxies tend to flatten or even ncrease around the locations where the breaks are found. The age and metallicity profiles are also affected, exhibiting...[abridged]
△ Less
Submitted 20 September, 2017;
originally announced September 2017.
-
The Gaia-ESO Survey: Matching Chemo-Dynamical Simulations to Observations of the Milky Way
Authors:
B. B. Thompson,
C. G. Few,
M. Bergemann,
B. K. Gibson,
B. A. MacFarlane,
A. Serenelli,
G. Gilmore,
S. Randich,
A. Vallenari,
E. J. Alfaro,
T. Bensby,
P. Francois,
A. J. Korn,
A. Bayo,
G. Carraro,
A. R. Casey,
M. T. Costado,
P Donati,
E. Franciosini,
A. Frasca,
A. Hourihane,
P. Jofre,
V. Hill,
U. Heiter,
S. E. Koposov
, et al. (17 additional authors not shown)
Abstract:
The typical methodology for comparing simulated galaxies with observational surveys is usually to apply a spatial selection to the simulation to mimic the region of interest covered by a comparable observational survey sample. In this work we compare this approach with a more sophisticated post-processing in which the observational uncertainties and selection effects (photometric, surface gravity…
▽ More
The typical methodology for comparing simulated galaxies with observational surveys is usually to apply a spatial selection to the simulation to mimic the region of interest covered by a comparable observational survey sample. In this work we compare this approach with a more sophisticated post-processing in which the observational uncertainties and selection effects (photometric, surface gravity and effective temperature) are taken into account. We compare a `solar neighbourhood analogue' region in a model Milky Way-like galaxy simulated with RAMSES-CH with fourth release Gaia-ESO survey data. We find that a simple spatial cut alone is insufficient and that observational uncertainties must be accounted for in the comparison. This is particularly true when the scale of uncertainty is large compared to the dynamic range of the data, e.g. in our comparison, the [Mg/Fe] distribution is affected much more than the more accurately determined [Fe/H] distribution. Despite clear differences in the underlying distributions of elemental abundances between simulation and observation, incorporating scatter to our simulation results to mimic observational uncertainty produces reasonable agreement. The quite complete nature of the Gaia-ESO survey means that the selection function has minimal impact on the distribution of observed age and metal abundances but this would become increasingly more important for surveys with narrower selection functions.
△ Less
Submitted 5 September, 2017;
originally announced September 2017.
-
The properties, origin and evolution of stellar clusters in galaxy simulations and observations
Authors:
C. L. Dobbs,
A. Adamo,
C. G. Few,
D. Calzetti,
D. A. Dale,
B. G. Elmegreen,
A. S. Evans,
D. A. Gouliermis,
K. Grasha,
E. K. Grebel,
K. E. Johnson,
H. Kim,
J. C. Lee,
M. Messa,
J. E. Ryon,
L. J. Smith,
D. A. Thilker,
L. Ubeda,
B. Whitmore
Abstract:
We investigate the properties and evolution of star particles in two simulations of isolated spiral galaxies, and two galaxies from cosmological simulations. Unlike previous numerical work, where typically each star particle represents one `cluster', for the isolated galaxies we are able to model features we term `clusters' with groups of particles. We compute the spatial distribution of stars wit…
▽ More
We investigate the properties and evolution of star particles in two simulations of isolated spiral galaxies, and two galaxies from cosmological simulations. Unlike previous numerical work, where typically each star particle represents one `cluster', for the isolated galaxies we are able to model features we term `clusters' with groups of particles. We compute the spatial distribution of stars with different ages, and cluster mass distributions, comparing our findings with observations including the recent LEGUS survey. We find that spiral structure tends to be present in older (100s Myrs) stars and clusters in the simulations compared to the observations. This likely reflects differences in the numbers of stars or clusters, the strength of spiral arms, and whether the clusters are allowed to evolve. Where we model clusters with multiple particles, we are able to study their evolution. The evolution of simulated clusters tends to follow that of their natal gas clouds. Massive, dense, long-lived clouds host massive clusters, whilst short-lived clouds host smaller clusters which readily disperse. Most clusters appear to disperse fairly quickly, in basic agreement with observational findings. We note that embedded clusters may be less inclined to disperse in simulations in a galactic environment with continuous accretion of gas onto the clouds than isolated clouds and correspondingly, massive young clusters which are no longer associated with gas tend not to occur in the simulations. Caveats of our models include that the cluster densities are lower than realistic clusters, and the simplistic implementation of stellar feedback.
△ Less
Submitted 1 September, 2016;
originally announced September 2016.
-
Testing Hydrodynamics Schemes in Galaxy Disc Simulations
Authors:
C. G. Few,
C. Dobbs,
A. Pettitt,
L. Konstandin
Abstract:
We examine how three fundamentally different numerical hydrodynamics codes follow the evolution of an isothermal galactic disc with an external spiral potential. We compare an adaptive mesh refinement code (RAMSES), a smoothed particle hydrodynamics code (sphNG), and a volume-discretised meshless code (GIZMO). Using standard refinement criteria, we find that RAMSES produces a disc that is less ver…
▽ More
We examine how three fundamentally different numerical hydrodynamics codes follow the evolution of an isothermal galactic disc with an external spiral potential. We compare an adaptive mesh refinement code (RAMSES), a smoothed particle hydrodynamics code (sphNG), and a volume-discretised meshless code (GIZMO). Using standard refinement criteria, we find that RAMSES produces a disc that is less vertically concentrated and does not reach such high densities as the sphNG or GIZMO runs. The gas surface density in the spiral arms increases at a lower rate for the RAMSES simulations compared to the other codes. There is also a greater degree of substructure in the sphNG and GIZMO runs and secondary spiral arms are more pronounced. By resolving the Jeans' length with a greater number of grid cells we achieve more similar results to the Lagrangian codes used in this study. Other alterations to the refinement scheme (adding extra levels of refinement and refining based on local density gradients) are less successful in reducing the disparity between RAMSES and sphNG/GIZMO. Although more similar, sphNG displays different density distributions and vertical mass profiles to all modes of GIZMO (including the smoothed particle hydrodynamics version). This suggests differences also arise which are not intrinsic to the particular method but rather due to its implementation. The discrepancies between codes (in particular, the densities reached in the spiral arms) could potentially result in differences in the locations and timescales for gravitational collapse, and therefore impact star formation activity in more complex galaxy disc simulations.
△ Less
Submitted 4 August, 2016; v1 submitted 31 May, 2016;
originally announced May 2016.
-
Origin of the Metallicity Distribution in the Thick Disc
Authors:
Maider S. Miranda,
Kate Pilkington,
Brad K. Gibson,
Christopher B. Brook,
Patricia Sánchez-Blázquez,
Ivan Minchev,
Christopher Gareth Few,
Rory Smith,
Rosa Domínguez-Tenreiro,
Aura Obreja,
Jeremy Bailin,
Greg S. Stinson
Abstract:
Aims. Using a suite of cosmological chemodynamical disc galaxy simulations, we assess how (a) radial metallicity gradients evolve with scaleheight; (b) the vertical metallicity gradients change through the thick disc; and (c) the vertical gradient of the stellar rotation velocity varies through the disc. We compare with the Milky Way to search for analogous trends. Methods. We analyse five simulat…
▽ More
Aims. Using a suite of cosmological chemodynamical disc galaxy simulations, we assess how (a) radial metallicity gradients evolve with scaleheight; (b) the vertical metallicity gradients change through the thick disc; and (c) the vertical gradient of the stellar rotation velocity varies through the disc. We compare with the Milky Way to search for analogous trends. Methods. We analyse five simulated spiral galaxies with masses comparable to the Milky Way. The simulations span a range of star formation and energy feedback strengths and prescriptions, particle- and grid-based hydrodynamical implementations, as well as initial conditions/assembly history. Results. Consistently, we find that the steeper, negative, radial metallicity gradients seen in the mid-plane flatten with increasing height away from the plane. In simulations with stronger (and/or more spatially-extended) feedback, the negative radial gradients invert, becoming positive for heights in excess of 1 kpc. Such behaviour is consistent with that inferred from recent observations. Our measurements of the vertical metallicity gradients show no clear correlation with galactocentric radius, and are in good agreement with those observed in the Milky Way's thick disc (locally). Conclusions. Simulations employing stronger/more extended feedback prescriptions possess radial and vertical metallicity and kinematic gradients more in line with recent observations. The inverted, positive, radial metallicity gradients seen in the simulated thick stellar discs originate from a population of younger, more metal-rich, stars formed in-situ, superimposed upon a background population of older migrators from the inner disc; the contrast provided by the former increases radially, due to the inside-out growth of the disc. A similar behaviour may be responsible for the same flattening seen in the radial gradients with scaleheight in the Milky Way.
△ Less
Submitted 14 December, 2015;
originally announced December 2015.
-
The imprint of satellite accretion on the chemical and dynamical properties of disc galaxies
Authors:
T. Ruiz-Lara,
C. G. Few,
B. K. Gibson,
I. Pérez,
E. Florido,
I. Minchev,
P. Sánchez-Blázquez
Abstract:
Aims: We study the effects of the cosmological assembly history on the chemical and dynamical properties of the discs of spiral galaxies as a function of radius. Methods: We make use of the simulated Milky-Way mass, fully-cosmological discs, from {\tt RaDES} (Ramses Disc Environment Study). We analyse their assembly history by examining the proximity of satellites to the galactic disc, instead of…
▽ More
Aims: We study the effects of the cosmological assembly history on the chemical and dynamical properties of the discs of spiral galaxies as a function of radius. Methods: We make use of the simulated Milky-Way mass, fully-cosmological discs, from {\tt RaDES} (Ramses Disc Environment Study). We analyse their assembly history by examining the proximity of satellites to the galactic disc, instead of their merger trees, to better gauge which satellites impact the disc. We present stellar age and metallicity profiles, Age-Metallicity Relation (AMR), Age-Velocity dispersion Relation (AVR), and Stellar Age Distribution (SAD) in several radial bins for the simulated galaxies. Results: Assembly histories can be divided into three different stages: i) a merger dominated phase, when a large number of mergers with mass ratios of $\sim$1:1 take place (lasting $\sim$3.2$\pm$0.4 Gyr on average); ii) a quieter phase, when $\sim$1:10 mergers take place (lasting $\sim$4.4$\pm$2.0 Gyr) - these two phases are able to kinematically heat the disc and produce a disc that is chemically mixed over its entire radial extension; iii) a "secular" phase where the few mergers that take place have mass ratios below 1:100, and which do not affect the disc properties (lasting $\sim$5.5$\pm$2.0 Gyr). Phase ii ends with a final merger event (at time $t_\mathrm{jump}$) marking the onset of important radial differences in the AMR, AVR, and SAD. Conclusions: Inverted AMR trends in the outer parts of discs, for stars younger than $t_\mathrm{jump}$, are found as the combined effect of radial motions and star formation in satellites temporarily located in these outer parts. "U-shaped" stellar age profiles change to an old plateau ($\sim$10 Gyr) in the outer discs for the entire {\tt RaDES} sample. This shape is a consequence of inside-out growth of the disc, radial motions of disc stars ... [abridged]
△ Less
Submitted 2 December, 2015;
originally announced December 2015.
-
Feedback from massive stars and gas expulsion from proto-Globular Clusters
Authors:
F. Calura,
C. G. Few,
D. Romano,
A. D'Ercole
Abstract:
Globular clusters are considerably more complex structures than previously thought, harbouring at least two stellar generations which present clearly distinct chemical abundances. Scenarios explaining the abundance patterns in globular clusters mostly assume that originally the clusters had to be much more massive than today, and that the second generation of stars originates from the gas shed by…
▽ More
Globular clusters are considerably more complex structures than previously thought, harbouring at least two stellar generations which present clearly distinct chemical abundances. Scenarios explaining the abundance patterns in globular clusters mostly assume that originally the clusters had to be much more massive than today, and that the second generation of stars originates from the gas shed by stars of the first generation (FG). The lack of metallicity spread in most globular clusters further requires that the supernova-enriched gas ejected by the FG is completely lost within ~30 Myr, a hypothesis never tested by means of three-dimensional hydrodynamic simulations. In this paper, we use 3D hydrodynamic simulations including stellar feedback from winds and supernovae, radiative cooling and self-gravity to study whether a realistic distribution of OB associations in a massive proto-GC of initial mass M_tot ~ 10^7 M_sun is sufficient to expel its entire gas content. Our numerical experiment shows that the coherence of different associations plays a fundamental role: as the bubbles interact, distort and merge, they carve narrow tunnels which reach deeper and deeper towards the innermost cluster regions, and through which the gas is able to escape. Our results indicate that after 3 Myr, the feedback from stellar winds is responsible for the removal of ~40% of the pristine gas, and that after 14 Myr, ~ 99% of the initial gas mass has been removed.
△ Less
Submitted 10 November, 2015;
originally announced November 2015.
-
Chemodynamics of a Simulated Disc Galaxy: Initial Mass Functions and SNIa Progenitors
Authors:
C. G. Few,
S. Courty,
B. K. Gibson,
L. Michel-Dansac,
F. Calura
Abstract:
We trace the formation and advection of several elements within a cosmological adaptive mesh refinement simulation of an L* galaxy. We use nine realisations of the same initial conditions with different stellar Initial Mass Functions (IMFs), mass limits for type-II and type-Ia supernovae (SNII, SNIa) and stellar lifetimes to constrain these sub-grid phenomena. Our code includes self-gravity, hydro…
▽ More
We trace the formation and advection of several elements within a cosmological adaptive mesh refinement simulation of an L* galaxy. We use nine realisations of the same initial conditions with different stellar Initial Mass Functions (IMFs), mass limits for type-II and type-Ia supernovae (SNII, SNIa) and stellar lifetimes to constrain these sub-grid phenomena. Our code includes self-gravity, hydrodynamics, star formation, radiative cooling and feedback from multiple sources within a cosmological framework. Under our assumptions of nucleosynthesis we find that SNII with progenitor masses of up to 100 Msun are required to match low metallicity gas oxygen abundances. Tardy SNIa are necessary to reproduce the classical chemical evolution knee in [O/Fe]-[Fe/H]: more prompt SNIa delayed time distributions do not reproduce this feature. Within our framework of hydrodynamical mixing of metals and galaxy mergers we find that chemical evolution is sensitive to the shape of the IMF and that there exists a degeneracy with the mass range of SNII. We look at the abundance plane and present the properties of different regions of the plot, noting the distinct chemical properties of satellites and a series of nested discs that have greater velocity dispersions, are more alpha-rich and metal poor with age.
△ Less
Submitted 21 August, 2014;
originally announced August 2014.
-
Properties of simulated Milky Way-mass galaxies in loose group and field environments
Authors:
C. G. Few,
B. K. Gibson,
S. Courty,
L. Michel-Dansac,
C. B. Brook,
G. S. Stinson
Abstract:
We test the validity of comparing simulated field disk galaxies with the empirical properties of systems situated within environments more comparable to loose groups, including the Milky Way's Local Group. Cosmological simulations of Milky Way-mass galaxies have been realised in two different environment samples: in the field and in environments with similar properties to the Local Group. Apart fr…
▽ More
We test the validity of comparing simulated field disk galaxies with the empirical properties of systems situated within environments more comparable to loose groups, including the Milky Way's Local Group. Cosmological simulations of Milky Way-mass galaxies have been realised in two different environment samples: in the field and in environments with similar properties to the Local Group. Apart from the environments of the galaxies, the samples are kept as homogeneous as possible with equivalent ranges in last major merger time, halo mass and halo spin. Comparison of these two samples allow for systematic differences in the simulations to be identified. Metallicity gradients, disk scale lengths, colours, magnitudes and age-velocity dispersion relations are studied for each galaxy in the suite and the strength of the link between these and environment of the galaxies is studied. The bulge-to-disk ratio of the galaxies show that these galaxies are less spheroid dominated than many other simulated galaxies in literature with the majority of both samples being disk dominated. We find that secular evolution and mergers dominate the spread of morphologies and metallicity gradients with no visible differences between the two environment samples. In contrast with this consistency in the two samples there is tentative evidence for a systematic difference in the velocity dispersion-age relations of galaxies in the different environments. Loose group galaxies appear to have more discrete steps in their velocity dispersion-age relations. We conclude that at the current resolution of cosmological galaxy simulations field environment galaxies are sufficiently similar to those in loose groups to be acceptable proxies for comparison with the Milky Way provided that a similar assembly history is considered.
△ Less
Submitted 3 October, 2012;
originally announced October 2012.
-
The stellar metallicity distribution of disc galaxies and bulges in cosmological simulations
Authors:
F. Calura,
B. K. Gibson,
L. Michel-Dansac,
G. S. Stinson,
M. Cignoni,
A. Dotter,
K. Pilkington,
E. L. House,
C. B. Brook,
C. G. Few,
J. Bailin,
H. M. P. Couchman,
J. Wadsley,
.
Abstract:
By means of high-resolution cosmological hydrodynamical simulations of Milky Way-like disc galaxies, we conduct an analysis of the associated stellar metallicity distribution functions (MDFs). After undertaking a kinematic decomposition of each simulation into spheroid and disc sub-components, we compare the predicted MDFs to those observed in the solar neighbourhood and the Galactic bulge. The ef…
▽ More
By means of high-resolution cosmological hydrodynamical simulations of Milky Way-like disc galaxies, we conduct an analysis of the associated stellar metallicity distribution functions (MDFs). After undertaking a kinematic decomposition of each simulation into spheroid and disc sub-components, we compare the predicted MDFs to those observed in the solar neighbourhood and the Galactic bulge. The effects of the star formation density threshold are visible in the star formation histories, which show a modulation in their behaviour driven by the threshold. The derived MDFs show median metallicities lower by 0.2-0.3 dex than the MDF observed locally in the disc and in the Galactic bulge. Possible reasons for this apparent discrepancy include the use of low stellar yields and/or centrally-concentrated star formation. The dispersions are larger than the one of the observed MDF; this could be due to simulated discs being kinematically hotter relative to the Milky Way. The fraction of low metallicity stars is largely overestimated, visible from the more negatively skewed MDF with respect to the observational sample. For our fiducial Milky Way analog, we study the metallicity distribution of the stars born "in situ" relative to those formed via accretion (from disrupted satellites), and demonstrate that this low-metallicity tail to the MDF is populated primarily by accreted stars. Enhanced supernova and stellar radiation energy feedback to the surrounding interstellar media of these pre-disrupted satellites is suggested as an important regulator of the MDF skewness.
△ Less
Submitted 6 September, 2012; v1 submitted 4 April, 2012;
originally announced April 2012.
-
RAMSES-CH: A New Chemodynamical Code for Cosmological Simulations
Authors:
C. Gareth Few,
Stephanie Courty,
Brad K. Gibson,
Daisuke Kawata,
Francesco Calura,
Romain Teyssier
Abstract:
We present a new chemodynamical code - Ramses-CH - for use in simulating the self-consistent evolution of chemical and hydrodynamical properties of galaxies within a fully cosmological framework. We build upon the adaptive mesh refinement code Ramses, which includes a treatment of self-gravity, hydrodynamics, star formation, radiative cooling, and supernovae feedback, to trace the dominant isotope…
▽ More
We present a new chemodynamical code - Ramses-CH - for use in simulating the self-consistent evolution of chemical and hydrodynamical properties of galaxies within a fully cosmological framework. We build upon the adaptive mesh refinement code Ramses, which includes a treatment of self-gravity, hydrodynamics, star formation, radiative cooling, and supernovae feedback, to trace the dominant isotopes of C, N, O, Ne, Mg, Si, and Fe. We include the contribution of Type Ia and II supernovae, in addition to low- and intermediate-mass asymptotic giant branch stars, relaxing the instantaneous recycling approximation. The new chemical evolution modules are highly flexible and portable, lending themselves to ready exploration of variations in the underpining stellar and nuclear physics. We apply Ramses-CH to the cosmological simulation of a typical L* galaxy, demonstrating the successful recovery of the basic empirical constraints regarding, [α/Fe]-[Fe/H] and Type Ia/II supernovae rates.
△ Less
Submitted 25 April, 2012; v1 submitted 28 February, 2012;
originally announced February 2012.
-
RAMSES-CH: A New Chemodynamics Code
Authors:
C. Gareth Few,
Stephanie Courty,
Brad K. Gibson
Abstract:
We present a new chemodynamical code based on the adaptive mesh refinement code RAMSES. The new code uses Eulerian hydrodynamics and N-body dynamics in a cosmological framework to trace the production and advection of several chemical species. It is the first such code to follow the self-consistent evolution of chemical elements in cosmological volumes while maintaining sub-kiloparsec resolution.…
▽ More
We present a new chemodynamical code based on the adaptive mesh refinement code RAMSES. The new code uses Eulerian hydrodynamics and N-body dynamics in a cosmological framework to trace the production and advection of several chemical species. It is the first such code to follow the self-consistent evolution of chemical elements in cosmological volumes while maintaining sub-kiloparsec resolution. The code will be used to simulate disk galaxies and explore the influence of chemical evolution models and star formation on galactic abundance ratios.
△ Less
Submitted 9 February, 2012;
originally announced February 2012.
-
Metallicity Gradients in Disks: Do Galaxies Form Inside-Out?
Authors:
K. Pilkington,
C. G. Few,
B. K. Gibson,
F. Calura,
L. Michel-Dansac,
R. J. Thacker,
M. Molla,
F. Matteucci,
A. Rahimi,
D. Kawata,
C. Kobayashi,
C. B. Brook,
G. S. Stinson,
H. M. P. Couchman,
J. Bailin,
J. Wadsley
Abstract:
We examine radial and vertical metallicity gradients using a suite of disk galaxy simulations, supplemented with two classic chemical evolution approaches. We determine the rate of change of gradient and reconcile differences between extant models and observations within the `inside-out' disk growth paradigm. A sample of 25 disks is used, consisting of 19 from our RaDES (Ramses Disk Environment St…
▽ More
We examine radial and vertical metallicity gradients using a suite of disk galaxy simulations, supplemented with two classic chemical evolution approaches. We determine the rate of change of gradient and reconcile differences between extant models and observations within the `inside-out' disk growth paradigm. A sample of 25 disks is used, consisting of 19 from our RaDES (Ramses Disk Environment Study) sample, realised with the adaptive mesh refinement code RAMSES. Four disks are selected from the MUGS (McMaster Unbiased Galaxy Simulations) sample, generated with the smoothed particle hydrodynamics (SPH) code GASOLINE, alongside disks from Rahimi et al. (GCD+) and Kobayashi & Nakasato (GRAPE-SPH). Two chemical evolution models of inside-out disk growth were employed to contrast the temporal evolution of their radial gradients with those of the simulations. We find that systematic differences exist between the predicted evolution of radial abundance gradients in the RaDES and chemical evolution models, compared with the MUGS sample; specifically, the MUGS simulations are systematically steeper at high-redshift, and present much more rapid evolution in their gradients. We find that the majority of the models predict radial gradients today which are consistent with those observed in late-type disks, but they evolve to this self-similarity in different fashions, despite each adhering to classical `inside-out' growth. We find that radial dependence of the efficiency with which stars form as a function of time drives the differences seen in the gradients; systematic differences in the sub-grid physics between the various codes are responsible for setting these gradients. Recent, albeit limited, data at redshift z=1.5 are consistent with the steeper gradients seen in our SPH sample, suggesting a modest revision of the classical chemical evolution models may be required.
△ Less
Submitted 13 February, 2012; v1 submitted 30 January, 2012;
originally announced January 2012.
-
The Cold Gas Content of Bulgeless Dwarf Galaxies
Authors:
K. Pilkington,
B. K. Gibson,
F. Calura,
A. M. Brooks,
L. Mayer,
C. B. Brook,
G. S. Stinson,
R. J. Thacker,
C. G. Few,
D. Cunnama,
J. Wadsley
Abstract:
We present an analysis of the neutral hydrogen (HI) properties of a fully cosmological hydrodynamical dwarf galaxy, run with varying simulation parameters. As reported by Governato et al. (2010), the high resolution, high star formation density threshold version of this galaxy is the first simulation to result in the successful reproduction of a (dwarf) spiral galaxy without any associated stellar…
▽ More
We present an analysis of the neutral hydrogen (HI) properties of a fully cosmological hydrodynamical dwarf galaxy, run with varying simulation parameters. As reported by Governato et al. (2010), the high resolution, high star formation density threshold version of this galaxy is the first simulation to result in the successful reproduction of a (dwarf) spiral galaxy without any associated stellar bulge. We have set out to compare in detail the HI distribution and kinematics of this simulated bulgeless disk with what is observed in a sample of nearby dwarfs. To do so, we extracted the radial gas density profiles, velocity dispersion (e.g., velocity ellipsoid, turbulence), and the power spectrum of structure within the cold interstellar medium from the simulations. The highest resolution dwarf, when using a high density star formation threshold comparable to densities of giant molecular clouds, possesses bulk characteristics consistent with those observed in nature, though the cold gas is not as radially extended as that observed in nearby dwarfs, resulting in somewhat excessive surface densities. The lines-of-sight velocity dispersion radial profiles have values that are in good agreement with observed dwarf galaxies, but due to the fact that only the streaming velocities of particles are tracked, a correction to include the thermal velocities can lead to profiles that are quite flat. The ISM power spectra of the simulations appear to possess more power on smaller spatial scales than that of the SMC. We conclude that unavoidable limitations remain due to the unresolved physics of star formation and feedback within pc-scale molecular clouds.
△ Less
Submitted 14 July, 2011;
originally announced July 2011.
-
Disc Heating: Comparing the Milky Way with Cosmological Simulations
Authors:
Elisa House,
Chris B. Brook,
Brad K. Gibson,
Patricia Sanchez-Blazquez,
Stephanie Courty,
Christopher G. Few,
Fabio Governato,
Daisuke Kawata,
Rok Roskar,
Matthias Steinmetz,
Gregory S. Stinson,
Romain Teyssier
Abstract:
We present the analysis of a suite of simulations run with different particle-and grid-based cosmological hydrodynamical codes and compare them with observational data of the Milky Way. This is the first study to make comparisons of properties of galaxies simulated with particle and grid-based codes. Our analysis indicates that there is broad agreement between these different modelling techniques.…
▽ More
We present the analysis of a suite of simulations run with different particle-and grid-based cosmological hydrodynamical codes and compare them with observational data of the Milky Way. This is the first study to make comparisons of properties of galaxies simulated with particle and grid-based codes. Our analysis indicates that there is broad agreement between these different modelling techniques. We study the velocity dispersion - age relation for disc stars at z=0 and find that four of the simulations are more consistent with observations by Holmberg et al. (2008) in which the stellar disc appears to undergo continual/secular heating. Two other simulations are in better agreement with the Quillen & Garnett (2001) observations that suggest a "saturation" in the heating profile for young stars in the disc. None of the simulations have thin discs as old as that of the Milky Way. We also analyse the kinematics of disc stars at the time of their birth for different epochs in the galaxies' evolution and find that in some simulations old stars are born cold within the disc and are subsequently heated, while other simulations possess old stellar populations which are born relatively hot. The models which are in better agreement with observations of the Milky Way's stellar disc undergo significantly lower minor-merger/assembly activity after the last major merger - i.e. once the disc has formed. All of the simulations are significantly "hotter" than the Milky Way disc; on top of the effects of mergers, we find a "floor" in the dispersion that is related to the underlying treatment of the heating and cooling of the interstellar medium, and the low density threshold which such codes use for star formation. This finding has important implications for all studies of disc heating that use hydrodynamical codes.
△ Less
Submitted 11 April, 2011;
originally announced April 2011.