-
Prospects for Detecting Gaps in Globular Cluster Stellar Streams in External Galaxies with the Nancy Grace Roman Space Telescope
Authors:
Christian Aganze,
Sarah Pearson,
Tjitske Starkenburg,
Gabriella Contardo,
Kathryn V. Johnston,
Kiyan Tavangar,
Adrian M. Price-Whelan,
Adam J. Burgasser
Abstract:
Stellar streams form through the tidal disruption of satellite galaxies or globular clusters orbiting a host galaxy. Globular cluster streams are exciting since they are thin (dynamically cold) and, therefore sensitive to perturbations from low-mass subhalos. Since the subhalo mass function differs depending on the dark matter composition, these gaps can provide unique constraints on dark matter m…
▽ More
Stellar streams form through the tidal disruption of satellite galaxies or globular clusters orbiting a host galaxy. Globular cluster streams are exciting since they are thin (dynamically cold) and, therefore sensitive to perturbations from low-mass subhalos. Since the subhalo mass function differs depending on the dark matter composition, these gaps can provide unique constraints on dark matter models. However, current samples are limited to the Milky Way. With its large field of view, deep imaging sensitivity, and high angular resolution, the upcoming Nancy Grace Roman Space Telescope (Roman) presents a unique opportunity to increase the number of observed streams and gaps significantly. This paper presents a first exploration of the prospects for detecting gaps in streams in M31 and other nearby galaxies with resolved stars. We simulate the formation of gaps in a Palomar-5-like stream and generate mock observations of these gaps with background stars in M31 and the foreground Milky Way stellar fields. We assess Roman's ability to detect gaps out to 10 Mpc through visual inspection and with the gap-finding tool ${\texttt{FindTheGap}}$. We conclude that gaps of $\approx 1.5$ kpc in streams that are created from subhalos of masses $\geq5 \times 10^6$M$_{\odot}$ are detectable within a 2-3 Mpc volume in exposures of 1000s to 1 hour. This volume contains $\approx 150$ galaxies, including $\approx 8$ galaxies with luminosities $>10^{9}~$L$_{\odot}$. Large samples of stream gaps in external galaxies will open up a new era of statistical analyses of gap characteristics in stellar streams and help constrain dark matter models.
△ Less
Submitted 15 February, 2024; v1 submitted 19 May, 2023;
originally announced May 2023.
-
JWST constraints on the UV luminosity density at cosmic dawn: implications for 21-cm cosmology
Authors:
Sultan Hassan,
Christopher C. Lovell,
Piero Madau,
Marc Huertas-Company,
Rachel S. Somerville,
Blakesley Burkhart,
Keri L. Dixon,
Robert Feldmann,
Tjitske K. Starkenburg,
John F. Wu,
Christian Kragh Jespersen,
Joseph D. Gelfand,
Ankita Bera
Abstract:
An unprecedented array of new observational capabilities are starting to yield key constraints on models of the epoch of first light in the Universe. In this Letter we discuss the implications of the UV radiation background at cosmic dawn inferred by recent JWST observations for radio experiments aimed at detecting the redshifted 21-cm hyperfine transition of diffuse neutral hydrogen. Under the ba…
▽ More
An unprecedented array of new observational capabilities are starting to yield key constraints on models of the epoch of first light in the Universe. In this Letter we discuss the implications of the UV radiation background at cosmic dawn inferred by recent JWST observations for radio experiments aimed at detecting the redshifted 21-cm hyperfine transition of diffuse neutral hydrogen. Under the basic assumption that the 21-cm signal is activated by the Ly$α$ photon field produced by metal-poor stellar systems, we show that a detection at the low frequencies of the EDGES and SARAS3 experiments may be expected from a simple extrapolation of the declining UV luminosity density inferred at $z\lesssim 14$ from JWST early galaxy data. Accounting for an early radiation excess above the CMB suggests a shallower or flat evolution to simultaneously reproduce low and high-$z$ current UV luminosity density constraints, which cannot be entirely ruled out, given the large uncertainties from cosmic variance and the faint-end slope of the galaxy luminosity function at cosmic dawn. Our findings raise the intriguing possibility that a high star formation efficiency at early times may trigger the onset of intense Ly$α$ emission at redshift $z\lesssim 20$ and produce a cosmic 21-cm absorption signal 200 Myr after the Big Bang.
△ Less
Submitted 11 October, 2023; v1 submitted 4 May, 2023;
originally announced May 2023.
-
VINTERGATAN-GM: The cosmological imprints of early mergers on Milky-Way-mass galaxies
Authors:
Martin P. Rey,
Oscar Agertz,
Tjitske K. Starkenburg,
Florent Renaud,
Gandhali D. Joshi,
Andrew Pontzen,
Nicolas F. Martin,
Diane K. Feuillet,
Justin I. Read
Abstract:
We present a new suite of cosmological zoom-in hydrodynamical ($\approx 20\, \mathrm{pc}$ spatial resolution) simulations of Milky-Way mass galaxies to study how a varying mass ratio for a Gaia-Sausage-Enceladus (GSE) progenitor impacts the $z=0$ chemodynamics of halo stars. Using the genetic modification approach, we create five cosmological histories for a Milky-Way-mass dark matter halo (…
▽ More
We present a new suite of cosmological zoom-in hydrodynamical ($\approx 20\, \mathrm{pc}$ spatial resolution) simulations of Milky-Way mass galaxies to study how a varying mass ratio for a Gaia-Sausage-Enceladus (GSE) progenitor impacts the $z=0$ chemodynamics of halo stars. Using the genetic modification approach, we create five cosmological histories for a Milky-Way-mass dark matter halo ($M_{200} \approx 10^{12} \, M_\mathrm{\odot}$), incrementally increasing the stellar mass ratio of a $z\approx2$ merger from 1:25 to 1:2, while fixing the galaxy's final dynamical, stellar mass and large-scale environment. We find markedly different morphologies at $z=0$ following this change in early history, with a growing merger resulting in increasingly compact and bulge-dominated galaxies. Despite this structural diversity, all galaxies show a radially-biased population of inner halo stars like the Milky-Way's GSE which, surprisingly, has a similar magnitude, age, $\rm [Fe/H]$ and $\rm [α/Fe]$ distribution whether the $z\approx2$ merger is more minor or major. This arises because a smaller ex-situ population at $z\approx2$ is compensated by a larger population formed in an earlier merger-driven starburst whose contribution to the GES can grow dynamically over time, with both populations strongly overlapping in the $\rm [Fe/H]-\rm [α/Fe]$ plane. Our study demonstrates that multiple high-redshift histories can lead to similar $z=0$ chemodynamical features in the halo, highlighting the need for additional constraints to distinguish them, and the importance of considering the full spectrum of progenitors when interpreting $z=0$ data to reconstruct our Galaxy's past.
△ Less
Submitted 15 March, 2023; v1 submitted 28 November, 2022;
originally announced November 2022.
-
Code Comparison in Galaxy Scale Simulations with Resolved Supernova Feedback: Lagrangian vs. Eulerian Methods
Authors:
Chia-Yu Hu,
Matthew C. Smith,
Romain Teyssier,
Greg L. Bryan,
Robbert Verbeke,
Andrew Emerick,
Rachel S. Somerville,
Blakesley Burkhart,
Yuan Li,
John C. Forbes,
Tjitske Starkenburg
Abstract:
We present a suite of high-resolution simulations of an isolated dwarf galaxy using four different hydrodynamical codes: {\sc Gizmo}, {\sc Arepo}, {\sc Gadget}, and {\sc Ramses}. All codes adopt the same physical model which includes radiative cooling, photoelectric heating, star formation, and supernova (SN) feedback. Individual SN explosions are directly resolved without resorting to sub-grid mo…
▽ More
We present a suite of high-resolution simulations of an isolated dwarf galaxy using four different hydrodynamical codes: {\sc Gizmo}, {\sc Arepo}, {\sc Gadget}, and {\sc Ramses}. All codes adopt the same physical model which includes radiative cooling, photoelectric heating, star formation, and supernova (SN) feedback. Individual SN explosions are directly resolved without resorting to sub-grid models, eliminating one of the major uncertainties in cosmological simulations. We find reasonable agreement on the time-averaged star formation rates as well as the joint density-temperature distributions between all codes. However, the Lagrangian codes show significantly burstier star formation, larger supernova-driven bubbles, and stronger galactic outflows compared to the Eulerian code. This is caused by the behavior in the dense, collapsing gas clouds when the Jeans length becomes unresolved: gas in Lagrangian codes collapses to much higher densities than in Eulerian codes, as the latter is stabilized by the minimal cell size. Therefore, more of the gas cloud is converted to stars and SNe are much more clustered in the Lagrangian models, amplifying their dynamical impact. The differences between Lagrangian and Eulerian codes can be reduced by adopting a higher star formation efficiency in Eulerian codes, which significantly enhances SN clustering in the latter. Adopting a zero SN delay time reduces burstiness in all codes, resulting in vanishing outflows as SN clustering is suppressed.
△ Less
Submitted 3 May, 2023; v1 submitted 22 August, 2022;
originally announced August 2022.
-
walter: A Tool for Predicting Resolved Stellar Population Observations with Applications to the Roman Space Telescope
Authors:
Lachlan Lancaster,
Sarah Pearson,
Benjamin F. Williams,
Kathryn V. Johnston,
Tjitske K. Starkenburg,
Erin Kado-Fong,
Anil C. Seth,
Eric F. Bell
Abstract:
Studies of resolved stellar populations in the Milky Way and nearby galaxies reveal an amazingly detailed and clear picture of galaxy evolution. Within the Local Group, the ability to probe the stellar populations of small and large galaxies opens up the possibility of exploring key questions such as the nature of dark matter, the detailed formation history of different galaxy components, and the…
▽ More
Studies of resolved stellar populations in the Milky Way and nearby galaxies reveal an amazingly detailed and clear picture of galaxy evolution. Within the Local Group, the ability to probe the stellar populations of small and large galaxies opens up the possibility of exploring key questions such as the nature of dark matter, the detailed formation history of different galaxy components, and the role of accretion in galactic formation. Upcoming wide-field surveys promise to extend this ability to all galaxies within 10~Mpc, drastically increasing our capability to decipher galaxy evolution and enabling statistical studies of galaxies' stellar populations. To facilitate the optimum use of these upcoming capabilities we develop a simple formalism to predict the density of resolved stars for an observation of a stellar population at fixed surface brightness and population parameters. We provide an interface to calculate all quantities of interest to this formalism via a public release of the code: \texttt{walter}. This code enables calculation of (i) the expected number density of detected stars, (ii) the exposure time needed to reach certain population features, such as the horizontal branch, and (iii) an estimate of the crowding limit, among other features. These calculations will be very useful for planning surveys with NASA's upcoming Nancy Grace Roman Space Telescope (Roman, formerly WFIRST), which we use for example calculations throughout this work.
△ Less
Submitted 5 July, 2022;
originally announced July 2022.
-
Preparing for low surface brightness science with the Vera C. Rubin Observatory: characterisation of tidal features from mock images
Authors:
G. Martin,
A. E. Bazkiaei,
M. Spavone,
E. Iodice,
J. C. Mihos,
M. Montes,
J. A. Benavides,
S. Brough,
J. L. Carlin,
C. A. Collins,
P. A. Duc,
F. A. Gómez,
G. Galaz,
H. M. Hernández-Toledo,
R. A. Jackson,
S. Kaviraj,
J. H. Knapen,
C. Martínez-Lombilla,
S. McGee,
D. O'Ryan,
D. J. Prole,
R. M. Rich,
J. Román,
E. A. Shah,
T. K. Starkenburg
, et al. (28 additional authors not shown)
Abstract:
Tidal features in the outskirts of galaxies yield unique information about their past interactions and are a key prediction of the hierarchical structure formation paradigm. The Vera C. Rubin Observatory is poised to deliver deep observations for potentially of millions of objects with visible tidal features, but the inference of galaxy interaction histories from such features is not straightforwa…
▽ More
Tidal features in the outskirts of galaxies yield unique information about their past interactions and are a key prediction of the hierarchical structure formation paradigm. The Vera C. Rubin Observatory is poised to deliver deep observations for potentially of millions of objects with visible tidal features, but the inference of galaxy interaction histories from such features is not straightforward. Utilising automated techniques and human visual classification in conjunction with realistic mock images produced using the NEWHORIZON cosmological simulation, we investigate the nature, frequency and visibility of tidal features and debris across a range of environments and stellar masses. In our simulated sample, around 80 per cent of the flux in the tidal features around Milky Way or greater mass galaxies is detected at the 10-year depth of the Legacy Survey of Space and Time (30-31 mag / sq. arcsec), falling to 60 per cent assuming a shallower final depth of 29.5 mag / sq. arcsec. The fraction of total flux found in tidal features increases towards higher masses, rising to 10 per cent for the most massive objects in our sample (M*~10^{11.5} Msun). When observed at sufficient depth, such objects frequently exhibit many distinct tidal features with complex shapes. The interpretation and characterisation of such features varies significantly with image depth and object orientation, introducing significant biases in their classification. Assuming the data reduction pipeline is properly optimised, we expect the Rubin Observatory to be capable of recovering much of the flux found in the outskirts of Milky Way mass galaxies, even at intermediate redshifts (z<0.2).
△ Less
Submitted 7 May, 2022; v1 submitted 15 March, 2022;
originally announced March 2022.
-
Hot-mode accretion and the physics of thin-disk galaxy formation
Authors:
Zachary Hafen,
Jonathan Stern,
James Bullock,
Alex B. Gurvich,
Sijie Yu,
Claude-Andre Faucher-Giguere,
Drummond B. Fielding,
Daniel Angles-Alcazar,
Eliot Quataert,
Andrew Wetzel,
Tjitske Starkenburg,
Michael Boylan-Kolchin,
Jorge Moreno,
Robert Feldmann,
Kareem El-Badry,
T. K. Chan,
Cameron Trapp,
Dusan Keres,
Philip F. Hopkins
Abstract:
We use FIRE simulations to study disk formation in z~0, Milky Way-mass galaxies, and conclude that a key ingredient for the formation of thin stellar disks is the ability for accreting gas to develop an aligned angular momentum distribution via internal cancellation *prior* to joining the galaxy. Among galaxies with a high fraction (>70%) of their young stars in a thin disk (h/R~0.1) we find that:…
▽ More
We use FIRE simulations to study disk formation in z~0, Milky Way-mass galaxies, and conclude that a key ingredient for the formation of thin stellar disks is the ability for accreting gas to develop an aligned angular momentum distribution via internal cancellation *prior* to joining the galaxy. Among galaxies with a high fraction (>70%) of their young stars in a thin disk (h/R~0.1) we find that: (i) hot, virial-temperature gas dominates the inflowing gas mass on halo scales (>~20 kpc), with radiative losses offset by compression heating; (ii) this hot accretion proceeds until angular momentum support slows inward motion, at which point the gas cools to T~10^4 K or less; (iii) prior to cooling, the accreting gas develops an angular momentum distribution that is aligned with the galaxy disk, and while cooling transitions from a quasi-spherical spatial configuration to a more flattened, disk-like configuration. We show that the existence of this "rotating cooling flow" accretion mode is strongly correlated with the fraction of stars forming in a thin disk among a sample of 17 z~0 galaxies spanning a halo mass range of 10^10.5 solar masses to 10^12 solar masses, or a stellar mass range 10^8 solar masses to 10^11 solar masses. Notably, galaxies with a thick disk or irregular morphology do not undergo significant angular momentum alignment of gas prior to accretion and show no correspondence between halo gas cooling and flattening. Our results suggest that rotating cooling flows (or, more generally, rotating subsonic flows) that become coherent and angular momentum-supported prior to accretion onto the galaxy are likely a necessary condition for the formation of thin, star-forming disk galaxies in a LambdaCDM universe.
△ Less
Submitted 6 June, 2022; v1 submitted 18 January, 2022;
originally announced January 2022.
-
Deep Hubble Space Telescope Observations of GW170817: Complete Light Curves and the Properties of the Galaxy Merger of NGC 4993
Authors:
Charles D. Kilpatrick,
Wen-fai Fong,
Peter K. Blanchard,
Joel Leja,
Anya E. Nugent,
Antonella Palmese,
Kerry Paterson,
Tjitske Starkenburg,
Kate D. Alexander,
Edo Berger,
Ryan Chornock,
Aprajita Hajela,
Raffaella Margutti
Abstract:
We present the complete set of {\it Hubble Space Telescope} imaging of the binary neutron star merger GW170817 and its optical counterpart AT 2017gfo. Including deep template imaging in F814W, F110W, F140W, and F160W at 3.4 years post-merger, we re-analyze the full light curve of AT 2017gfo across 12 bands from 5--1273 rest-frame days after merger. We obtain four new detections of the short $γ$-ra…
▽ More
We present the complete set of {\it Hubble Space Telescope} imaging of the binary neutron star merger GW170817 and its optical counterpart AT 2017gfo. Including deep template imaging in F814W, F110W, F140W, and F160W at 3.4 years post-merger, we re-analyze the full light curve of AT 2017gfo across 12 bands from 5--1273 rest-frame days after merger. We obtain four new detections of the short $γ$-ray burst (GRB) 170817A afterglow from 109--170 rest-frame days post-merger. These detections are consistent with the previously observed $β=-0.6$ spectral index in the afterglow light curve with no evidence for spectral evolution. We also analyze our limits in the case of novel late-time optical and IR emission signatures, such as a kilonova afterglow or infrared dust echo, but find our limits are not constraining in these contexts. We use the new data to construct deep optical and infrared stacks, reaching limits of $M=-6.3$ to $-4.6$ mag, to analyze the local environment around AT 2017gfo and low surface brightness features in its host galaxy NGC 4993. We rule out the presence of any globular cluster at the position of AT 2017gfo to $2.3 \times 10^{4} L_{\odot}$, including those with the reddest $V-H$ colors. Finally, we analyze the substructure of NGC 4993 in deep residual imaging, and find shell features which extend up to 71.8\arcsec\ (14.2 kpc) from the center of the galaxy. We find that the shells have a cumulative stellar mass of $6.3\times10^{8} M_{\odot}$, roughly 2% the total stellar mass of NGC 4993, and mass-weighted ages of $>$3 Gyr. We conclude that it was unlikely the GW170817 progenitor system formed in the galaxy merger, which based on dynamical signatures and the stellar population in the shells mostly likely occurred 220--685 Myr ago.
△ Less
Submitted 13 September, 2021;
originally announced September 2021.
-
The In-situ Origins of Dwarf Stellar Outskirts in FIRE-2
Authors:
Erin Kado-Fong,
Robyn E. Sanderson,
Jenny E. Greene,
Emily C. Cunningham,
Coral Wheeler,
T. K. Chan,
Kareem El-Badry,
Philip F. Hopkins,
Andrew Wetzel,
Michael Boylan-Kolchin,
Claude-André Faucher-Giguère,
Song Huang,
Eliot Quataert,
Tjitske Starkenburg
Abstract:
Extended, old, and round stellar halos appear to be ubiquitous around high-mass dwarf galaxies ($10^{8.5}<M_\star/M_\odot<10^{9.6}$) in the observed universe. However, it is unlikely that these dwarfs have undergone a sufficient number of minor mergers to form stellar halos that are composed of predominantly accreted stars. Here, we demonstrate that FIRE-2 (Feedback in Realistic Environments) cosm…
▽ More
Extended, old, and round stellar halos appear to be ubiquitous around high-mass dwarf galaxies ($10^{8.5}<M_\star/M_\odot<10^{9.6}$) in the observed universe. However, it is unlikely that these dwarfs have undergone a sufficient number of minor mergers to form stellar halos that are composed of predominantly accreted stars. Here, we demonstrate that FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulations are capable of producing dwarf galaxies with realistic structure, including both a thick disk and round stellar halo. Crucially, these stellar halos are formed in-situ, largely via the outward migration of disk stars. However, there also exists a large population of "non-disky" dwarfs in FIRE-2 that lack a well-defined disk/halo and do not resemble the observed dwarf population. These non-disky dwarfs tend to be either more gas poor or to have burstier recent star formation histories than the disky dwarfs, suggesting that star formation feedback may be preventing disk formation. Both classes of dwarfs underscore the power of a galaxy's intrinsic shape -- which is a direct quantification of the distribution of the galaxy's stellar content -- to interrogate the feedback implementation in simulated galaxies.
△ Less
Submitted 2 May, 2022; v1 submitted 10 September, 2021;
originally announced September 2021.
-
The Hough Stream Spotter: A New Method for Detecting Linear Structure in Resolved Stars and Application to the Stellar Halo of M31
Authors:
Sarah Pearson,
Susan E. Clark,
Alexis J. Demirjian,
Kathryn V. Johnston,
Melissa K. Ness,
Tjitske K. Starkenburg,
Benjamin F. Williams,
Rodrigo A. Ibata
Abstract:
Stellar streams from globular clusters (GCs) offer constraints on the nature of dark matter and have been used to explore the dark matter halo structure and substructure of our Galaxy. Detection of GC streams in other galaxies would broaden this endeavor to a cosmological context, yet no such streams have been detected to date. To enable such exploration, we develop the Hough Stream Spotter (HSS),…
▽ More
Stellar streams from globular clusters (GCs) offer constraints on the nature of dark matter and have been used to explore the dark matter halo structure and substructure of our Galaxy. Detection of GC streams in other galaxies would broaden this endeavor to a cosmological context, yet no such streams have been detected to date. To enable such exploration, we develop the Hough Stream Spotter (HSS), and apply it to the Pan-Andromeda Archaeological Survey (PAndAS) photometric data of resolved stars in M31's stellar halo. We first demonstrate that our code can re-discover known dwarf streams in M31. We then use the HSS to blindly identify 27 linear GC stream-like structures in the PAndAS data. For each HSS GC stream candidate, we investigate the morphologies of the streams and the colors and magnitudes of all stars in the candidate streams. We find that the five most significant detections show a stronger signal along the red giant branch in color-magnitude diagrams (CMDs) than spurious non-stream detections. Lastly, we demonstrate that the HSS will easily detect globular cluster streams in future Nancy Grace Roman Space Telescope data of nearby galaxies. This has the potential to open up a new discovery space for GC stream studies, GC stream gap searches, and for GC stream-based constraints on the nature of dark matter.
△ Less
Submitted 3 January, 2022; v1 submitted 30 June, 2021;
originally announced July 2021.
-
IQ Collaboratory III: The Empirical Dust Attenuation Framework -- Taking Hydrodynamical Simulations with a Grain of Dust
Authors:
ChangHoon Hahn,
Tjitske K. Starkenburg,
Daniel Anglés-Alcázar,
Ena Choi,
Romeel Davé,
Claire Dickey,
Kartheik G. Iyer,
Ariyeh H. Maller,
Rachel S. Somerville,
Jeremy L. Tinker,
L. Y. Aaron Yung
Abstract:
We present the Empirical Dust Attenuation (EDA) framework -- a flexible prescription for assigning realistic dust attenuation to simulated galaxies based on their physical properties. We use the EDA to forward model synthetic observations for three state-of-the-art large-scale cosmological hydrodynamical simulations: SIMBA, IllustrisTNG, and EAGLE. We then compare the optical and UV color-magnitud…
▽ More
We present the Empirical Dust Attenuation (EDA) framework -- a flexible prescription for assigning realistic dust attenuation to simulated galaxies based on their physical properties. We use the EDA to forward model synthetic observations for three state-of-the-art large-scale cosmological hydrodynamical simulations: SIMBA, IllustrisTNG, and EAGLE. We then compare the optical and UV color-magnitude relations, $(g-r) - M_r$ and $(FUV-NUV)-M_r$, of the simulations to a $M_r < -20$ and UV complete SDSS galaxy sample using likelihood-free inference. Without dust, none of the simulations match observations, as expected. With the EDA, however, we can reproduce the observed color-magnitude with all three simulations. Furthermore, the attenuation curves predicted by our dust prescription are in good agreement with the observed attenuation-slope relations and attenuation curves of star-forming galaxies. However, the EDA does not predict star-forming galaxies with low $A_V$ since simulated star-forming galaxies are intrinsically much brighter than observations. Additionally, the EDA provides, for the first time, predictions on the attenuation curves of quiescent galaxies, which are challenging to measure observationally. Simulated quiescent galaxies require shallower attenuation curves with lower amplitude than star-forming galaxies. The EDA, combined with forward modeling, provides an effective approach for shedding light on dust in galaxies and probing hydrodynamical simulations. This work also illustrates a major limitation in comparing galaxy formation models: by adjusting dust attenuation, simulations that predict significantly different galaxy populations can reproduce the same UV and optical observations.
△ Less
Submitted 17 June, 2021;
originally announced June 2021.
-
How cosmological merger histories shape the diversity of stellar haloes
Authors:
Martin P. Rey,
Tjitske K. Starkenburg
Abstract:
We introduce and apply a new approach to probe the response of galactic stellar haloes to the interplay between cosmological merger histories and galaxy formation physics. We perform dark matter-only, zoomed simulations of two Milky Way-mass hosts and make targeted, controlled changes to their cosmological histories using the genetic modification technique. Populating each history's stellar halo w…
▽ More
We introduce and apply a new approach to probe the response of galactic stellar haloes to the interplay between cosmological merger histories and galaxy formation physics. We perform dark matter-only, zoomed simulations of two Milky Way-mass hosts and make targeted, controlled changes to their cosmological histories using the genetic modification technique. Populating each history's stellar halo with a semi-empirical, particle tagging approach then enables a controlled study, with all instances converging to the same large-scale structure, dynamical and stellar mass at $z=0$ as their reference. These related merger scenarios alone generate an extended spread in stellar halo mass fractions (1.5 dex) comparable to the observed population, with the largest scatter achieved by growing late ($z\leq1$) major mergers that spread out existing stars to create massive, in-situ dominated stellar haloes. Increasing a last major merger at $z\sim2$ brings more accreted stars into the inner regions, resulting in smaller scatter in the outskirts which are predominantly built by subsequent minor events. Exploiting the flexibility of our semi-empirical approach, we show that the diversity of stellar halo masses across scenarios is reduced by allowing shallower slopes in the stellar mass--halo mass relation for dwarf galaxies, while it remains conserved when central stars are born with hotter kinematics across cosmic time. The merger-dependent diversity of stellar haloes thus responds distinctly to assumptions in modelling the central and dwarf galaxies respectively, opening exciting prospects to constrain star formation and feedback at different galactic mass-scales with the coming generation of deep, photometric observatories.
△ Less
Submitted 20 January, 2022; v1 submitted 17 June, 2021;
originally announced June 2021.
-
IQ Collaboratory II: The Quiescent Fraction of Isolated, Low Mass Galaxies Across Simulations and Observations
Authors:
Claire M Dickey,
Tjitske K Starkenburg,
Marla Geha,
ChangHoon Hahn,
Daniel Anglés-Alcázar,
Ena Choi,
Romeel Davé,
Shy Genel,
Kartheik G Iyer,
Ariyeh H Maller,
Nir Mandelker,
Rachel S Somerville,
L Y Aaron Yung
Abstract:
We compare three major large-scale hydrodynamical galaxy simulations (EAGLE, Illustris-TNG, and SIMBA) by forward modeling simulated galaxies into observational space and computing the fraction of isolated and quiescent low mass galaxies as a function of stellar mass. Using SDSS as our observational template, we create mock surveys and synthetic spectroscopic and photometric observations of each s…
▽ More
We compare three major large-scale hydrodynamical galaxy simulations (EAGLE, Illustris-TNG, and SIMBA) by forward modeling simulated galaxies into observational space and computing the fraction of isolated and quiescent low mass galaxies as a function of stellar mass. Using SDSS as our observational template, we create mock surveys and synthetic spectroscopic and photometric observations of each simulation, adding realistic noise and observational limits. All three simulations show a decrease in the number of quiescent, isolated galaxies in the mass range $\mathrm{M}_* = 10^{9-10} \ \mathrm{M}_\odot$, in broad agreement with observations. However, even after accounting for observational and selection biases, none of the simulations reproduce the observed absence of quiescent field galaxies below $\mathrm{M}_*=10^{9} \ \mathrm{M}_\odot$. We find that the low mass quiescent populations selected via synthetic observations have consistent quenching timescales, despite apparent variation in the late time star formation histories. The effect of increased numerical resolution is not uniform across simulations and cannot fully mitigate the differences between the simulations and the observations. The framework presented here demonstrates a path towards more robust and accurate comparisons between theoretical simulations and galaxy survey observations, while the quenching threshold serves as a sensitive probe of feedback implementations.
△ Less
Submitted 2 October, 2020;
originally announced October 2020.
-
The breakBRD Breakdown: Using IllustrisTNG to Track the Quenching of an Observationally-Motivated Sample of Centrally Star-Forming Galaxies
Authors:
Claire Kopenhafer,
Tjitske K. Starkenburg,
Stephanie Tonnesen,
Sarah Tuttle
Abstract:
The observed breakBRD ("break bulges in red disks") galaxies are a nearby sample of face-on disk galaxies with particularly centrally concentrated star formation: they have red disks but recent star formation in their centers as measured by the D$_n$4000 spectral index (Tuttle & Tonnesen 2020). In this paper, we search for breakBRD analogues in the IllustrisTNG simulation and describe their histor…
▽ More
The observed breakBRD ("break bulges in red disks") galaxies are a nearby sample of face-on disk galaxies with particularly centrally concentrated star formation: they have red disks but recent star formation in their centers as measured by the D$_n$4000 spectral index (Tuttle & Tonnesen 2020). In this paper, we search for breakBRD analogues in the IllustrisTNG simulation and describe their history and future. We find that a small fraction (${\sim}4\%$ at $z=0$; ${\sim}1\%$ at $z=0.5$) of galaxies fulfill the breakBRD criteria, in agreement with observations. In comparison with the mass-weighted parent IllustrisTNG sample, these galaxies tend to consist of a higher fraction of satellite and splashback galaxies. However, the central, non-splashback breakBRD galaxies show similar environments, black hole masses, and merger rates, indicating that there is not a single formation trigger for inner star formation and outer quenching. We determine that breakBRD analogue galaxies as a whole are in the process of quenching. The breakBRD state - with its highly centrally concentrated star formation - is uncommon in the history of either currently quiescent or star-forming galaxies; however, approximately 10% of $10^{10} < M_\ast/M_{\odot} < 10^{11}$ quiescent galaxies at $z=0$ have experienced SFR concentrations comparable to those of the breakBRDs in their past. Additionally, the breakBRD state is short-lived, lasting a few hundred Myr up to ${\sim}2$ Gyr. The observed breakBRD galaxies may therefore be a unique sample of outside-in quenching galaxies.
△ Less
Submitted 29 September, 2020;
originally announced September 2020.
-
The time-scales probed by star formation rate indicators for realistic, bursty star formation histories from the FIRE simulations
Authors:
José A. Flores Velázquez,
Alexander B. Gurvich,
Claude-André Faucher-Giguère,
James S. Bullock,
Tjitske K. Starkenburg,
Jorge Moreno,
Alexandres Lazar,
Francisco J. Mercado,
Jonathan Stern,
Martin Sparre,
Christopher C. Hayward,
Andrew Wetzel,
Kareem El-Badry
Abstract:
Understanding the rate at which stars form is central to studies of galaxy formation. Observationally, the star formation rates (SFRs) of galaxies are measured using the luminosity in different frequency bands, often under the assumption of a time-steady SFR in the recent past. We use star formation histories (SFHs) extracted from cosmological simulations of star-forming galaxies from the FIRE pro…
▽ More
Understanding the rate at which stars form is central to studies of galaxy formation. Observationally, the star formation rates (SFRs) of galaxies are measured using the luminosity in different frequency bands, often under the assumption of a time-steady SFR in the recent past. We use star formation histories (SFHs) extracted from cosmological simulations of star-forming galaxies from the FIRE project to analyze the time-scales to which the H$α$ and far-ultraviolet (FUV) continuum SFR indicators are sensitive. In these simulations, the SFRs are highly time variable for all galaxies at high redshift, and continue to be bursty to z=0 in dwarf galaxies. When FIRE SFHs are partitioned into their bursty and time-steady phases, the best-fitting FUV time-scale fluctuates from its ~10 Myr value when the SFR is time-steady to >~100 Myr immediately following particularly extreme bursts of star formation during the bursty phase. On the other hand, the best-fitting averaging time-scale for H$α$ is generally insensitive to the SFR variability in the FIRE simulations and remains ~5 Myr at all times. These time-scales are shorter than the 100 Myr and 10 Myr time-scales sometimes assumed in the literature for FUV and H$α$, respectively, because while the FUV emission persists for stellar populations older than 100 Myr, the time-dependent luminosities are strongly dominated by younger stars. Our results confirm that the ratio of SFRs inferred using H$α$ vs. FUV can be used to probe the burstiness of star formation in galaxies.
△ Less
Submitted 15 December, 2020; v1 submitted 19 August, 2020;
originally announced August 2020.
-
The Diversity and Variability of Star Formation Histories in Models of Galaxy Evolution
Authors:
Kartheik G. Iyer,
Sandro Tacchella,
Shy Genel,
Christopher C. Hayward,
Lars Hernquist,
Alyson M. Brooks,
Neven Caplar,
Romeel Davé,
Benedikt Diemer,
John C. Forbes,
Eric Gawiser,
Rachel S. Somerville,
Tjitske K. Starkenburg
Abstract:
Understanding the variability of galaxy star formation histories (SFHs) across a range of timescales provides insight into the underlying physical processes that regulate star formation within galaxies. We compile the SFHs of galaxies at $z=0$ from an extensive set of models, ranging from cosmological hydrodynamical simulations (Illustris, IllustrisTNG, Mufasa, Simba, EAGLE), zoom simulations (FIR…
▽ More
Understanding the variability of galaxy star formation histories (SFHs) across a range of timescales provides insight into the underlying physical processes that regulate star formation within galaxies. We compile the SFHs of galaxies at $z=0$ from an extensive set of models, ranging from cosmological hydrodynamical simulations (Illustris, IllustrisTNG, Mufasa, Simba, EAGLE), zoom simulations (FIRE-2, g14, and Marvel/Justice League), semi-analytic models (Santa Cruz SAM) and empirical models (UniverseMachine), and quantify the variability of these SFHs on different timescales using the power spectral density (PSD) formalism. We find that the PSDs are well described by broken power-laws, and variability on long timescales ($\gtrsim1$ Gyr) accounts for most of the power in galaxy SFHs. Most hydrodynamical models show increased variability on shorter timescales ($\lesssim300$ Myr) with decreasing stellar mass. Quenching can induce $\sim0.4-1$ dex of additional power on timescales $>1$ Gyr. The dark matter accretion histories of galaxies have remarkably self-similar PSDs and are coherent with the in-situ star formation on timescales $>3$ Gyr. There is considerable diversity among the different models in their (i) power due to SFR variability at a given timescale, (ii) amount of correlation with adjacent timescales (PSD slope), (iii) evolution of median PSDs with stellar mass, and (iv) presence and locations of breaks in the PSDs. The PSD framework is a useful space to study the SFHs of galaxies since model predictions vary widely. Observational constraints in this space will help constrain the relative strengths of the physical processes responsible for this variability.
△ Less
Submitted 15 July, 2020;
originally announced July 2020.
-
First results from SMAUG: The need for preventative stellar feedback and improved baryon cycling in semi-analytic models of galaxy formation
Authors:
Viraj Pandya,
Rachel S. Somerville,
Daniel Anglés-Alcázar,
Christopher C. Hayward,
Greg L. Bryan,
Drummond B. Fielding,
John C. Forbes,
Blakesley Burkhart,
Shy Genel,
Lars Hernquist,
Chang-Goo Kim,
Stephanie Tonnesen,
Tjitske Starkenburg
Abstract:
Semi-analytic models (SAMs) are a promising means of tracking the physical processes associated with galaxy formation, but many of their approximations have not been rigorously tested. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, we compare predictions from the FIRE-2 hydrodynamical "zoom-in" simulations to those from the Santa Cruz SAM run on the same…
▽ More
Semi-analytic models (SAMs) are a promising means of tracking the physical processes associated with galaxy formation, but many of their approximations have not been rigorously tested. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, we compare predictions from the FIRE-2 hydrodynamical "zoom-in" simulations to those from the Santa Cruz SAM run on the same halo merger trees, with an emphasis on the global mass flow cycle. Our study includes 13 halos spanning low-mass dwarfs (M_vir~10^10 M_sun at z=0), intermediate-mass dwarfs (M_vir~10^11 M_sun) and Milky Way-mass galaxies (M_vir~10^12 M_sun). The SAM and FIRE-2 predictions agree relatively well with each other in terms of stellar and interstellar mass, but differ dramatically on circumgalactic mass (the SAM is lower than FIRE-2 by ~3 orders of magnitude for dwarfs). Strikingly, the SAM predicts higher gas accretion rates for dwarfs compared to FIRE-2 by factors of ~10-100, and this is compensated for with higher mass outflow rates in the SAM. We argue that the most severe model discrepancies are caused by the lack of preventative stellar feedback and the assumptions for halo gas cooling and recycling in the SAM. As a first step towards resolving these model tensions, we present a simple yet promising new preventative stellar feedback model in which the energy carried by supernova-driven winds is allowed to heat some fraction of gas outside of halos to at least the virial temperature such that accretion is suppressed.
△ Less
Submitted 30 October, 2020; v1 submitted 29 June, 2020;
originally announced June 2020.
-
Decoupling the rotation of stars and gas -- II: the link between black hole activity and MaNGA kinematics in TNG
Authors:
Christopher Duckworth,
Tjitske K. Starkenburg,
Shy Genel,
Timothy Davis,
Melanie Habouzit,
Katarina Kraljic,
Rita Tojeiro
Abstract:
We study the relationship between supermassive black hole (BH) feedback, BH luminosity and the kinematics of stars and gas for galaxies in IllustrisTNG. We use a sample of galaxies with mock MaNGA observations to identify kinematic misalignment at $z=0$ (difference in rotation of stars and gas), for which we follow the evolutionary history of BH activity and gas properties over the last 8 Gyrs. Mi…
▽ More
We study the relationship between supermassive black hole (BH) feedback, BH luminosity and the kinematics of stars and gas for galaxies in IllustrisTNG. We use a sample of galaxies with mock MaNGA observations to identify kinematic misalignment at $z=0$ (difference in rotation of stars and gas), for which we follow the evolutionary history of BH activity and gas properties over the last 8 Gyrs. Misaligned low mass galaxies typically have boosted BH luminosity, BH growth and have had more energy injected into the gas over the last 8 Gyr in comparison to aligned galaxies. These properties likely lead to outflows and gas loss, in agreement with active low mass galaxies in observations. We show that splitting on BH luminosity at $z=0$ produces statistically consistent distributions of kinematic misalignment at $z=0$, however, splitting on the maximum BH luminosity over the last 8 Gyrs does not. While instantaneous correlation at $z=0$ is difficult due to misalignment persisting on longer timescales, the relationship between BH activity and misalignment is clear. High mass quenched galaxies with misalignment typically have similar BH luminosities, show no overall gas loss, and have typically lower gas phase metallicity over the last 8 Gyrs in comparison to those aligned; suggesting external origin.
△ Less
Submitted 12 November, 2019;
originally announced November 2019.
-
Multiple retrograde substructures in the Galactic halo: A shattered view of Galactic history
Authors:
Helmer H. Koppelman,
Amina Helmi,
Davide Massari,
Adrian M. Price-Whelan,
Tjitske K. Starkenburg
Abstract:
Aims. Several kinematic and chemical substructures have been recently found amongst Milky Way halo stars with retrograde motions. It is currently unclear how these various structures are related to each other. This Letter aims to shed light on this issue. Methods. We explore the retrograde halo with an augmented version of the Gaia DR2 RVS sample, extended with data from three large spectroscopic…
▽ More
Aims. Several kinematic and chemical substructures have been recently found amongst Milky Way halo stars with retrograde motions. It is currently unclear how these various structures are related to each other. This Letter aims to shed light on this issue. Methods. We explore the retrograde halo with an augmented version of the Gaia DR2 RVS sample, extended with data from three large spectroscopic surveys, namely RAVE, APOGEE and LAMOST. In this dataset, we identify several structures using the HDBSCAN clustering algorithm. We discuss their properties and possible links using all the available chemical and dynamical information. Results. In concordance with previous work, we find that stars with [Fe/H] $<-1$ have more retrograde motions than those with [Fe/H] $>-1$. The retrograde halo contains a mixture of debris from objects like Gaia-Enceladus, Sequoia, and even the chemically defined thick-disc. We find that the Sequoia has a smaller range in orbital energies than previously suggested and is confined to high-energy. Sequoia could be a small galaxy in itself, but since it overlaps both in integrals-of-motion space and chemical abundance space with the less bound debris of Gaia-Enceladus, its nature cannot be fully settled yet. In the low-energy part of the halo we find evidence for at least one more distinct structure: Thamnos. Stars in Thamnos are on low inclination, mildly eccentric retrograde orbits, moving at $v_φ\approx-150$ km/s, and are chemically distinct from the other structures. Conclusions. Even with the excellent Gaia DR2 data it remains challenging to piece together all the fragments found in the retrograde halo. At this point, we are very much in need of large datasets with high-quality high-resolution spectra and tailored high-resolution hydrodynamical simulations of galaxy mergers.
△ Less
Submitted 19 September, 2019;
originally announced September 2019.
-
Detecting Thin Stellar Streams in External Galaxies: Resolved Stars & Integrated Light
Authors:
Sarah Pearson,
Tjitske K. Starkenburg,
Kathryn V. Johnston,
Benjamin F. Williams,
Rodrigo A. Ibata
Abstract:
The morphology of thin stellar streams can be used to test the nature of dark matter. It is therefore crucial to extend searches for globular cluster streams to other galaxies than the Milky Way. In this paper, we investigate the current and future prospects of detecting globular cluster streams in external galaxies in resolved stars (e.g. with WFIRST) and using integrated light (e.g. with HSC, LS…
▽ More
The morphology of thin stellar streams can be used to test the nature of dark matter. It is therefore crucial to extend searches for globular cluster streams to other galaxies than the Milky Way. In this paper, we investigate the current and future prospects of detecting globular cluster streams in external galaxies in resolved stars (e.g. with WFIRST) and using integrated light (e.g. with HSC, LSST and Euclid). In particular, we inject mock-streams to data from the PAndAS M31 survey, and produce simulated M31 backgrounds mimicking what WFIRST will observe in M31. Additionally, we estimate the distance limit to which globular cluster streams will be observable. Our results demonstrate that for a 1 hour (1000 sec.) exposure, using conservative estimates, WFIRST should detect globular cluster streams in resolved stars in galaxies out to distances of ~3.5 Mpc (~2 Mpc). This volume contains 199 (122) galaxies of which >90% are dwarfs. With integrated light, thin streams can be resolved out to ~100 Mpc with HSC and LSST and to ~600 Mpc with WFIRST and Euclid. The low surface brightness of the streams (typically >30 mag/arcsec$^2$), however, will make them difficult to detect, unless the streams originate from very young clusters. We emphasize that if the external galaxies do not host spiral arms or galactic bars, gaps in their stellar streams provide an ideal test case for evidence of interactions with dark matter subhalos. Furthermore, obtaining a large samples of thin stellar streams can help constrain the orbital structure and hence the potentials of external halos.
△ Less
Submitted 7 June, 2019;
originally announced June 2019.
-
On the origin of star-gas counterrotation in low-mass galaxies
Authors:
Tjitske K. Starkenburg,
Laura V. Sales,
Shy Genel,
Christina Manzano-King,
Gabriela Canalizo,
Lars Hernquist
Abstract:
Stars in galaxies form from the cold rotationally supported gaseous disks that settle at the center of dark matter halos. In the simplest models, such angular momentum is acquired early on at the time of collapse of the halo and preserved thereafter, implying a well-aligned spin for the stellar and gaseous component. Observations however have shown the presence of gaseous disks in counterrotation…
▽ More
Stars in galaxies form from the cold rotationally supported gaseous disks that settle at the center of dark matter halos. In the simplest models, such angular momentum is acquired early on at the time of collapse of the halo and preserved thereafter, implying a well-aligned spin for the stellar and gaseous component. Observations however have shown the presence of gaseous disks in counterrotation with the stars. We use the Illustris numerical simulations to study the origin of such counterrotation in low mass galaxies ($M_\star = 2 \times 10^9$ - $5 \times 10^{10}\; \rm M_\odot$), a sample where mergers have not played a significant role. Only ${\sim}1\%$ of our sample shows a counterrotating gaseous disk at $z=0$. These counterrotating disks arise in galaxies that have had a significant episode of gas removal followed by the acquisition of new gas with misaligned angular momentum. In our simulations, we identify two main channels responsible for the gas loss: a strong feedback burst and gas stripping during a fly-by passage through a more massive group environment. Once settled, counterrotation can be long-lived with several galaxies in our sample displaying misaligned components consistently for more than $2$ Gyr. As a result, no major correlation with the present day environment or structural properties might remain, except for a slight preference for early type morphologies and a lower than average gas content at a given stellar mass.
△ Less
Submitted 8 March, 2019;
originally announced March 2019.
-
The host galaxy of GRB 980425 / SN1998bw: a collisional ring galaxy
Authors:
M. Arabsalmani,
S. Roychowdhury,
T. K. Starkenburg,
L. Christensen,
E. Le Floc'h,
N. Kanekar,
F. Bournaud,
M. A. Zwaan,
J. P. U. Fynbo,
P. Møller,
E. Pian
Abstract:
We report Giant Metrewave Radio Telescope (GMRT) , Very Large Telescope (VLT) and Spitzer Space Telescope observations of ESO 184$-$G82, the host galaxy of GRB 980425/SN 1998bw, that yield evidence of a companion dwarf galaxy at a projected distance of 13 kpc. The companion, hereafter GALJ193510-524947, is a gas-rich, star-forming galaxy with a star formation rate of…
▽ More
We report Giant Metrewave Radio Telescope (GMRT) , Very Large Telescope (VLT) and Spitzer Space Telescope observations of ESO 184$-$G82, the host galaxy of GRB 980425/SN 1998bw, that yield evidence of a companion dwarf galaxy at a projected distance of 13 kpc. The companion, hereafter GALJ193510-524947, is a gas-rich, star-forming galaxy with a star formation rate of $\rm0.004\,M_{\odot}\, yr^{-1}$, a gas mass of $10^{7.1\pm0.1} M_{\odot}$, and a stellar mass of $10^{7.0\pm0.3} M_{\odot}$. The interaction between ESO 184$-$G82 and GALJ193510-524947 is evident from the extended gaseous structure between the two galaxies in the GMRT HI 21 cm map. We find a ring of high column density HI gas, passing through the actively star forming regions of ESO 184$-$G82 and the GRB location. This ring lends support to the picture in which ESO 184$-$G82 is interacting with GALJ193510-524947. The massive stars in GALJ193510-524947 have similar ages to those in star-forming regions in ESO 184$-$G82, also suggesting that the interaction may have triggered star formation in both galaxies. The gas and star formation properties of ESO 184$-$G82 favour a head-on collision with GALJ193510-524947 rather than a classical tidal interaction. We perform state-of-the art simulations of dwarf--dwarf mergers and confirm that the observed properties of ESO 184$-$G82 can be reproduced by collision with a small companion galaxy. This is a very clear case of interaction in a gamma ray burst host galaxy, and of interaction-driven star formation giving rise to a gamma ray burst in a dense environment.
△ Less
Submitted 1 March, 2019;
originally announced March 2019.
-
What Is Inside Matters: Simulated Green Valley Galaxies Have Centrally Concentrated Star Formation
Authors:
Tjitske K. Starkenburg,
Stephanie Tonnesen,
Claire Kopenhafer
Abstract:
In spatially resolved galaxy observations, star formation rate radial profiles are found to correlate with total specific star formation rates. A central depletion in star formation is thought to correlate with the globally depressed star formation rates of, for example, galaxies within the Green Valley. We present, for the first time, radial specific star formation rate profiles for a statistical…
▽ More
In spatially resolved galaxy observations, star formation rate radial profiles are found to correlate with total specific star formation rates. A central depletion in star formation is thought to correlate with the globally depressed star formation rates of, for example, galaxies within the Green Valley. We present, for the first time, radial specific star formation rate profiles for a statistical sample of simulated galaxies from the Illustris and EAGLE large cosmological simulations. For galaxies on the star-forming sequence, simulated specific star formation rate profiles are in loose agreement with observations, although galaxies from the EAGLE simulation are too centrally peaked, and galaxies from Illustris have a steeper decline at large radii. However, both galaxy samples show centrally concentrated star formation for galaxies in the Green Valley at all galaxy stellar masses, indicating that quenching occurs from the outside-in, in strong conflict with observations of inside-out quenching. These results appear in spite of the different feedback models in these two simulations. We conclude that the distribution of star formation within galaxies is a strong additional constraint for star formation and feedback models in simulations, in particular related to the quenching of star formation.
△ Less
Submitted 18 March, 2019; v1 submitted 3 December, 2018;
originally announced December 2018.
-
IQ-Collaboratory 1.1: the Star-Forming Sequence of Simulated Central Galaxies
Authors:
ChangHoon Hahn,
Tjitske K. Starkenburg,
Ena Choi,
Romeel Davé,
Claire M. Dickey,
Marla C. Geha,
Shy Genel,
Christopher C. Hayward,
Ariyeh H. Maller,
Nityasri Mandyam,
Viraj Pandya,
Gergö Popping,
Mika Rafieferantsoa,
Rachel S. Somerville,
Jeremy L. Tinker
Abstract:
A tightly correlated star formation rate-stellar mass relation of star forming galaxies, or star-forming sequence (SFS), is a key feature in galaxy property-space that is predicted by modern galaxy formation models. We present a flexible data-driven approach for identifying this SFS over a wide range of star formation rates and stellar masses using Gaussian mixture modeling (GMM). Using this metho…
▽ More
A tightly correlated star formation rate-stellar mass relation of star forming galaxies, or star-forming sequence (SFS), is a key feature in galaxy property-space that is predicted by modern galaxy formation models. We present a flexible data-driven approach for identifying this SFS over a wide range of star formation rates and stellar masses using Gaussian mixture modeling (GMM). Using this method, we present a consistent comparison of the $z=0$ SFSs of central galaxies in the Illustris, EAGLE, and Mufasa hydrodynamic simulations and the Santa Cruz semi-analytic model (SC-SAM), alongside data from the Sloan Digital Sky Survey. We find, surprisingly, that the amplitude of the SFS varies by up to ${\sim} 0.7\,\mathrm{dex}$ (factor of ${\sim} 5$) among the simulations with power-law slopes range from $0.7$ to $1.2$. In addition to the SFS, our GMM method also identifies sub-components in the star formation rate-stellar mass relation corresponding to star-burst, transitioning, and quiescent sub-populations. The hydrodynamic simulations are similarly dominated by SFS and quiescent sub-populations unlike the SC-SAM, which predicts substantial fractions of transitioning and star-burst galaxies at stellar masses above and below $10^{10} M_\odot$, respectively. All of the simulations also produce an abundance of low-mass quiescent central galaxies in apparent tension with observations. These results illustrate that, even among models that well reproduce many observables of the galaxy population, the $z=0$ SFS and other sub-populations still show marked differences that can provide strong constraints on galaxy formation models.
△ Less
Submitted 5 September, 2018;
originally announced September 2018.
-
Dark influences III. Structural characterization of minor mergers of dwarf galaxies with dark satellites
Authors:
Tjitske K. Starkenburg,
Amina Helmi,
Laura V. Sales
Abstract:
In the current concordance cosmology small halos are expected to be completely dark and can significantly perturb low-mass galaxies during minor merger interactions. These interactions may well contribute to the diversity of the dwarf galaxy population. Dwarf galaxies in the field are often observed to have peculiarities in their structure, morphology, and kinematics, as well as strong bursts of s…
▽ More
In the current concordance cosmology small halos are expected to be completely dark and can significantly perturb low-mass galaxies during minor merger interactions. These interactions may well contribute to the diversity of the dwarf galaxy population. Dwarf galaxies in the field are often observed to have peculiarities in their structure, morphology, and kinematics, as well as strong bursts of star formation without apparent cause. We aim to characterize the signatures of minor mergers of dwarf galaxies with dark satellites to aid their observational identification. We explore and quantify a variety of structural, morphological, and kinematic indicators of merging dwarf galaxies and their remnants using a suite of hydrodynamical simulations. The most sensitive indicators of mergers with dark satellites are large asymmetries in the gaseous and stellar distributions, enhanced central surface brightness and starbursts, and velocity offsets and misalignments between the cold gas and stellar components. In general, merging systems span a wide range of values of the most commonly used indicators, while isolated objects tend to have more confined values. Interestingly, we find in our simulations that a significantly off-centered burst of star formation can pinpoint the location of the dark satellite. Observational systems with such characteristics are perhaps the most promising for unveiling the presence of the hitherto, missing satellites.
△ Less
Submitted 12 September, 2016; v1 submitted 31 December, 2015;
originally announced January 2016.
-
Dark influences II: gas and star formation in minor mergers of dwarf galaxies with dark satellites
Authors:
Tjitske K. Starkenburg,
Amina Helmi,
Laura V. Sales
Abstract:
It has been proposed that mergers induce starbursts and lead to important morphological changes in galaxies. Most studies so far have focused on large galaxies, but dwarfs might also experience such events, since the halo mass function is scale-free in the concordance cosmological model. Notably, because of their low mass, most of their interactions will be with dark satellites. In this paper we f…
▽ More
It has been proposed that mergers induce starbursts and lead to important morphological changes in galaxies. Most studies so far have focused on large galaxies, but dwarfs might also experience such events, since the halo mass function is scale-free in the concordance cosmological model. Notably, because of their low mass, most of their interactions will be with dark satellites. In this paper we follow the evolution of gas-rich disky dwarf galaxies as they experience a minor merger with a dark satellite. We aim to characterize the effects of such an interaction on the dwarf's star formation, morphology, and kinematical properties. We performed a suite of carefully set-up hydrodynamical simulations of dwarf galaxies that include dark matter, gas, and stars merging with a satellite consisting solely of dark matter. For the host system we vary the gas fraction, disk size and thickness, halo mass, and concentration, while we explore different masses, concentrations, and orbits for the satellite. We find that the interactions cause strong starbursts of both short and long duration in the dwarfs. Their star formation rates increase by factors of a few to 10 or more. They are strongest for systems with extended gas disks and high gas fractions merging with a high-concentration satellite on a planar, radial orbit. In contrast to analogous simulations of Milky Way-mass galaxies, many of the systems experience strong morphological changes and become spheroidal even in the presence of significant amounts of gas. The simulated systems compare remarkably well with the observational properties of a large selection of irregular dwarf galaxies and blue compact dwarfs. This implies that mergers with dark satellites might well be happening but not be fully evident, and may thus play a role in the diversity of the dwarf galaxy population.
△ Less
Submitted 2 November, 2015; v1 submitted 25 August, 2015;
originally announced August 2015.
-
Dark influences: imprints of dark satellites on dwarf galaxies
Authors:
Tjitske K. Starkenburg,
Amina Helmi
Abstract:
In the context of the current $Λ$CDM cosmological model small dark matter haloes are abundant and satellites of dwarf galaxies are expected to be predominantly dark. Since low mass galaxies have smaller baryon fractions interactions with these satellites may leave particularly dramatic imprints. We uncover the influence of the most massive of these dark satellites on disky dwarf galaxies and the p…
▽ More
In the context of the current $Λ$CDM cosmological model small dark matter haloes are abundant and satellites of dwarf galaxies are expected to be predominantly dark. Since low mass galaxies have smaller baryon fractions interactions with these satellites may leave particularly dramatic imprints. We uncover the influence of the most massive of these dark satellites on disky dwarf galaxies and the possible dynamical and morphological transformations that result from these interactions. We use a suite of carefully set-up, controlled simulations of isolated dwarf galaxies. The primary dwarf galaxies have solely a stellar disk in the dark matter halo and the secundaries are completely devoid of baryons. We vary the disk mass, halo concentration, initial disk thickness and inclination of the satellite orbit. The disky dwarf galaxies are heated and disrupted due to the minor merger event, more extremely for higher satellite over disk mass ratios, and the morphology and kinematics are significantly altered. Moreover, for less concentrated haloes the minor merger can completely destroy the disk leaving a low-luminosity spheroidal-like galaxy instead. We conclude that dwarf galaxies are very much susceptible to being disturbed by dark galaxies and that even a minor merger event can significantly disrupt and alter the structure and kinematics of a dwarf galaxy. This process may be seen as a new channel for the formation of dwarf spheroidal galaxies.
△ Less
Submitted 27 November, 2014; v1 submitted 5 November, 2014;
originally announced November 2014.
-
Dark satellites and the morphology of dwarf galaxies
Authors:
Amina Helmi,
L. V. Sales,
E. Starkenburg,
T. K. Starkenburg,
C. A. Vera-Ciro,
G. De Lucia,
Y. -S. Li
Abstract:
One of the strongest predictions of the LambdaCDM cosmological model is the presence of dark satellites orbiting all types of galaxies. We focus here on the dynamical effects of such satellites on disky dwarf galaxies, and demonstrate that these encounters can be dramatic. Although mergers with M_sat > M_d are not very common, because of the lower baryonic content they occur much more frequently o…
▽ More
One of the strongest predictions of the LambdaCDM cosmological model is the presence of dark satellites orbiting all types of galaxies. We focus here on the dynamical effects of such satellites on disky dwarf galaxies, and demonstrate that these encounters can be dramatic. Although mergers with M_sat > M_d are not very common, because of the lower baryonic content they occur much more frequently on the dwarf scale than for L_*-galaxies. As an example, we present a numerical simulation of a 20% (virial) mass ratio merger between a dark satellite and a disky dwarf (akin to the Fornax dwarf galaxy in luminosity) that shows that the merger remnant has a spheroidal morphology. We conclude that perturbations by dark satellites provide a plausible path for the formation of dSph systems and also could trigger starbursts in gas rich dwarf galaxies. Therefore the transition from disky to the often amorphous, irregular, or spheroidal morphologies of dwarfs could be a natural consequence of the dynamical heating of hitherto unobservable dark satellites.
△ Less
Submitted 11 June, 2012;
originally announced June 2012.