Showing 1–50 of 129 results for author: Gnedin, O

GD-1 Stellar Stream and Cocoon in the DESI Early Data Release

Authors: Monica Valluri, Parker Fagrelius, Sergey. E. Koposov, Ting S. Li, Oleg Y. Gnedin, Eric F. Bell, Raymond G. Carlberg, Andrew P. Cooper, Jessia N. Aguilar, Carlos Allende Prieto, Vasily Belokurov, Leandro Beraldo e Silva, David Brooks, Amanda Byström, Todd Claybaugh, Kyle Dawson, Arjun Dey, Peter Doel, Jaime E. Forero-Romero, Enrique Gaztañaga, Satya Gontcho A Gontcho, Klaus Honscheid, T . Kisner, Anthony Kremin, A. Lambert , et al. (27 additional authors not shown)

Abstract: We present ~ 126 new spectroscopically identified members of the GD-1 tidal stream obtained with the 5000-fiber Dark Energy Spectroscopic Instrument (DESI). We confirm the existence of a ``cocoon'' which is broad (FWHM~2.932deg~460pc) and kinematically hot (velocity dispersion, sigma~5-8km/s) component that surrounds a narrower (FWHM~0.353deg~55pc) and colder (sigma~ 2.2-2.6km/s) thin stream compo… ▽ More We present ~ 126 new spectroscopically identified members of the GD-1 tidal stream obtained with the 5000-fiber Dark Energy Spectroscopic Instrument (DESI). We confirm the existence of a ``cocoon'' which is broad (FWHM~2.932deg~460pc) and kinematically hot (velocity dispersion, sigma~5-8km/s) component that surrounds a narrower (FWHM~0.353deg~55pc) and colder (sigma~ 2.2-2.6km/s) thin stream component (based on a median per star velocity precision of 2.7km/s). The cocoon extends over at least a ~ 20deg segment of the stream observed by DESI. The thin and cocoon components have similar mean values of [Fe/H]: -2.54+/- 0.04dex and -2.45+/-0.06dex suggestive of a common origin. The data are consistent with the following scenarios for the origin of the cocoon. The progenitor of the GD-1 stream was an accreted globular cluster (GC) and: (a) the cocoon was produced by pre-accretion tidal stripping of the GC while it was still inside its parent dwarf galaxy; (b) the cocoon is debris from the parent dwarf galaxy; (c) an initially thin GC tidal stream was heated by impacts from dark subhalos in the Milky Way; (d) an initially thin GC stream was heated by a massive Sagittarius dwarf galaxy; or a combination of some these. In the first two cases the velocity dispersion and mean metallicity are consistent with the parent dwarf galaxy having a halo mass of ~0^9\msun. Future DESI spectroscopy and detailed modeling may enable us to distinguish between these possible origins. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: Submitted to ApJ, 23 pages, 13 figures 4 tables

AuriDESI: Mock Catalogues for the DESI Milky Way Survey

Authors: Namitha Kizhuprakkat, Andrew P. Cooper, Alexander H. Riley, Sergey E. Koposov, Jessica Nicole Aguilar, Steven Ahlen, Carlos Allende Prieto, David Brooks, Todd Claybaugh, Kyle Dawson, Axel de la Macorra, Peter Doel, Jaime E. Forero-Romero, Carlos Frenk, Enrique Gaztañaga, Oleg Y. Gnedin, Robert J. J. Grand, Satya Gontcho A Gontcho, Klaus Honscheid, Robert Kehoe, Martin Landriau, Marc Manera, Aaron Meisner, Ramon Miquel, Jundan Nie , et al. (9 additional authors not shown)

Abstract: The Dark Energy Spectroscopic Instrument Milky Way Survey (DESI MWS) will explore the assembly history of the Milky Way by characterising remnants of ancient dwarf galaxy accretion events and improving constraints on the distribution of dark matter in the outer halo. We present mock catalogues that reproduce the selection criteria of MWS and the format of the final MWS data set. These catalogues c… ▽ More The Dark Energy Spectroscopic Instrument Milky Way Survey (DESI MWS) will explore the assembly history of the Milky Way by characterising remnants of ancient dwarf galaxy accretion events and improving constraints on the distribution of dark matter in the outer halo. We present mock catalogues that reproduce the selection criteria of MWS and the format of the final MWS data set. These catalogues can be used to test methods for quantifying the properties of stellar halo substructure and reconstructing the Milky Way's accretion history with the MWS data, including the effects of halo-to-halo variance. The mock catalogues are based on a phase-space kernel expansion technique applied to star particles in the Auriga suite of six high-resolution $Λ$CDM magneto-hydrodynamic zoom-in simulations. They include photometric properties (and associated errors) used in DESI target selection and the outputs of the MWS spectral analysis pipeline (radial velocity, metallicity, surface gravity, and temperature). They also include information from the underlying simulation, such as the total gravitational potential and information on the progenitors of accreted halo stars. We discuss how the subset of halo stars observable by MWS in these simulations corresponds to their true content and properties. These mock Milky Ways have rich accretion histories, resulting in a large number of substructures that span the whole stellar halo out to large distances and have substantial overlap in the space of orbital energy and angular momentum. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: Accepted for publication in MNRAS, 31 pages, 27 figues, 7 tables. The mock catalogues are available at https://data.desi.lbl.gov/public/papers/mws/auridesi/v1

Forecasting the Population of Globular Cluster Streams in Milky Way-type Galaxies

Authors: Sarah Pearson, Ana Bonaca, Yingtian Chen, Oleg Y. Gnedin

Abstract: Thin stellar streams originating from globular clusters are among the most sensitive tracers of low-mass dark-matter subhalos. Joint analysis of the entire population of stellar streams will place the most robust constraints on the dark-matter subhalo mass function, and therefore the nature of dark matter. Here we use a hierarchical model of globular cluster formation to forecast the total number,… ▽ More Thin stellar streams originating from globular clusters are among the most sensitive tracers of low-mass dark-matter subhalos. Joint analysis of the entire population of stellar streams will place the most robust constraints on the dark-matter subhalo mass function, and therefore the nature of dark matter. Here we use a hierarchical model of globular cluster formation to forecast the total number, masses and radial distribution of dissolved globular cluster in Milky Way-like galaxies. Furthermore, we generate mock stellar streams from these progenitors' orbital histories taking into account the clusters' formation and accretion time, mass, and metallicity. Out of $\sim$10,000 clusters more massive than $10^4$ M$_{\odot}$, $\sim$9000 dissolved in the central bulge and are fully phase-mixed at the present, while the remaining $\sim$1000 survive as coherent stellar streams. This suggests that the current census of $\sim$80 globular cluster streams in the Milky Way is severely incomplete. Beyond 15 kpc from the Galactic center we are missing hundreds of streams, of which the vast majority are from accreted GCs. Deep Rubin photometry $(g\lesssim27.5)$ would be able to detect these streams, even the most distant ones beyond $> 75$ kpc. We also find that M31 will have an abundance of streams at galactocentric radii of 30-100 kpc. We conclude that future surveys will find a multitude of stellar streams from globular clusters which can be used for dark matter subhalo searches. △ Less

Submitted 24 May, 2024; originally announced May 2024.

Comments: 21 pages, 13 figures, submitted to AAS journals, comments welcome

A SN Ia Near a Globular Cluster in the Early-Type Galaxy NGC 5353

Authors: Joel N. Bregman, Oleg Y. Gnedin, Patrick O. Seitzer, Zhijie Qu

Abstract: No progenitor of a Type Ia supernova is known, but in old population early-type galaxies, one may find SN Ia associated with globular clusters, yielding a population age and metallicity. It also provides insight into the formation path and the SN enhancement rate in globular clusters. We sought to find such associations and identified SN 2019ein to be within the ground-based optical positional unc… ▽ More No progenitor of a Type Ia supernova is known, but in old population early-type galaxies, one may find SN Ia associated with globular clusters, yielding a population age and metallicity. It also provides insight into the formation path and the SN enhancement rate in globular clusters. We sought to find such associations and identified SN 2019ein to be within the ground-based optical positional uncertainty of a globular cluster candidate within the early-type galaxy NGC 5353 at about 30 Mpc distance. We reduced the positional uncertainties by obtaining Hubble Space Telescope images with the Advanced Camera for Surveys, using filters F475W and F814W and obtained in June 2020. We find that the globular cluster candidate has a magnitude, color, and angular extent that are consistent with it being a typical globular cluster. The separation between the globular cluster and SN 2019ein is 0.43'', or 59 pc in projection. The chance occurrence with a random globular cluster is about 3%, favoring but not proving an association. If the SN progenitor originated in the globular cluster, one scenario is that SN 2019ein was previously a double degenerate white dwarf binary that was dynamically ejected from the globular cluster and exploded within 10 Myr; models do not predict this to be common. Another, but less likely scenario is where the progenitor remained bound to the globular cluster, allowing the double degenerate binary to inspiral on a much longer timescale before producing a SN. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: 7 pages with 1 figure

Galaxy assembly revealed by globular clusters

Authors: Yingtian Chen, Oleg Y. Gnedin

Abstract: Many observable properties of globular clusters (GCs) provide valuable insights for unveiling the hierarchical assembly of their host galaxy. For the Milky Way (MW) in particular, GCs from different accreted satellite galaxies show distinct chemical, spatial, kinematic, and age distributions. Here we examine such clustering features for model GC populations in simulated galaxies, which are careful… ▽ More Many observable properties of globular clusters (GCs) provide valuable insights for unveiling the hierarchical assembly of their host galaxy. For the Milky Way (MW) in particular, GCs from different accreted satellite galaxies show distinct chemical, spatial, kinematic, and age distributions. Here we examine such clustering features for model GC populations in simulated galaxies, which are carefully selected to match various observational constraints of the MW assembly. We evaluate several widely used clustering, dimensionality reduction, and supervised classification methods on these model GCs, using 10 properties that are observable in the MW. We can categorize in-situ and ex-situ formed GCs with about 90% accuracy, based solely on their clustering features in these 10 variables. The methods are also effective in distinguishing the last major merger in MW-analogs with similar accuracy. Although challenging, we still find it possible to identify one, and only one, additional smaller satellite. We develop a new technique to classify the progenitors of MW GCs by combining several methods and weighting them by the validated accuracy. According to this technique, about 60% of GCs belong to the in-situ group, 20% are associated with the Gaia-Sausage/Enceladus event, and 10% are associated with the Sagittarius dwarf galaxy. The remaining 10% of GCs cannot be reliably associated with any single accretion event. △ Less

Submitted 19 March, 2024; v1 submitted 30 January, 2024; originally announced January 2024.

Comments: 24 pages, 13 figures, accepted for publication in The Open Journal of Astrophysics. Classification of the Galactic globular clusters available at https://websites.umich.edu/~ognedin/mw_gc_classification.txt

Catalogue of model star clusters in the Milky Way and M31 galaxies

Authors: Yingtian Chen, Oleg Y. Gnedin

Abstract: Detailed understanding of the formation and evolution of globular clusters (GCs) has been recently advanced through a combination of numerical simulations and analytical models. We employ a state-of-the-art model to create a comprehensive catalogue of simulated clusters in three Milky Way (MW) and three Andromeda (M31) analogue galaxies. Our catalogue aims to connect the chemical and kinematic pro… ▽ More Detailed understanding of the formation and evolution of globular clusters (GCs) has been recently advanced through a combination of numerical simulations and analytical models. We employ a state-of-the-art model to create a comprehensive catalogue of simulated clusters in three Milky Way (MW) and three Andromeda (M31) analogue galaxies. Our catalogue aims to connect the chemical and kinematic properties of GCs to the assembly histories of their host galaxies. We apply the model to a selected sample of simulated galaxies that closely match the virial mass, circular velocity profile, and defining assembly events of the MW and M31. The resulting catalogue has been calibrated to successfully reproduce key characteristics of the observed GC systems, including total cluster mass, mass function, metallicity distribution, radial profile, and velocity dispersion. We find that clusters in M31 span a wider range of age and metallicity, relative to the MW, possibly due to M31's recent major merger. Such a merger also heated up the in-situ GC population to higher orbital energy and introduced a large number of ex-situ clusters at large radii. Understanding the impacts of galaxy mergers and accretion on the GC populations is crucial for uncovering the galaxy assembly histories. △ Less

Submitted 30 October, 2023; v1 submitted 23 September, 2023; originally announced September 2023.

Comments: 16 pages, 8 figures, accepted for publication in MNRAS. The catalogue of model globular clusters available at https://github.com/ognedin/gc_model_mw

RomAndromeda: The Roman Survey of the Andromeda Halo

Authors: Arjun Dey, Joan Najita, Carrie Filion, Jiwon Jesse Han, Sarah Pearson, Rosemary Wyse, Adrien C. R. Thob, Borja Anguiano, Miranda Apfel, Magda Arnaboldi, Eric F. Bell, Leandro Beraldo e Silva, Gurtina Besla, Aparajito Bhattacharya, Souradeep Bhattacharya, Vedant Chandra, Yumi Choi, Michelle L. M. Collins, Emily C. Cunningham, Julianne J. Dalcanton, Ivanna Escala, Hayden R. Foote, Annette M. N. Ferguson, Benjamin J. Gibson, Oleg Y. Gnedin , et al. (28 additional authors not shown)

Abstract: As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($Δt\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the ha… ▽ More As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($Δt\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the halo (10$σ$ detection in F146, F062 of 26.5, 26.1AB mag respectively) and yield proper motions to $\sim$25 microarcsec/year (i.e., $\sim$90 km/s) for all stars brighter than F146 $\approx 23.6$ AB mag (i.e., reaching the red clump stars in the Andromeda halo). This survey will yield (through averaging) high-fidelity proper motions for all satellites and compact substructures in the Andromeda halo and will enable statistical searches for clusters in chemo-dynamical space. Adding a third epoch during the extended mission will improve these proper motions by $\sim t^{-1.5}$, to $\approx 11$ km/s, but this requires obtaining the first epoch in Year 1 of Roman operations. In combination with ongoing and imminent spectroscopic campaigns with ground-based telescopes, this Roman survey has the potential to yield full 3-d space motions of $>$100,000 stars in the Andromeda halo, including (by combining individual measurements) robust space motions of its entire globular cluster and most of its dwarf galaxy satellite populations. It will also identify high-velocity stars in Andromeda, providing unique information on the processes that create this population. These data offer a unique opportunity to study the immigration history, halo formation, and underlying dark matter scaffolding of a galaxy other than our own. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: Submitted in response to the call for Roman Space Telescope Core Community Survey white papers

NANCY: Next-generation All-sky Near-infrared Community surveY

Authors: Jiwon Jesse Han, Arjun Dey, Adrian M. Price-Whelan, Joan Najita, Edward F. Schlafly, Andrew Saydjari, Risa H. Wechsler, Ana Bonaca, David J Schlegel, Charlie Conroy, Anand Raichoor, Alex Drlica-Wagner, Juna A. Kollmeier, Sergey E. Koposov, Gurtina Besla, Hans-Walter Rix, Alyssa Goodman, Douglas Finkbeiner, Abhijeet Anand, Matthew Ashby, Benedict Bahr-Kalus, Rachel Beaton, Jayashree Behera, Eric F. Bell, Eric C Bellm , et al. (184 additional authors not shown)

Abstract: The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GAL… ▽ More The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GALAH, 4MOST, WEAVE, MOONS, PFS, UVEX, NEO Surveyor, etc.). Roman can uniquely provide uniform high-spatial-resolution (~0.1 arcsec) imaging over the entire sky, vastly expanding the science reach and precision of all of these near-term and future surveys. This imaging will not only enhance other surveys, but also facilitate completely new science. By imaging the full sky over two epochs, Roman can measure the proper motions for stars across the entire Milky Way, probing 100 times fainter than Gaia out to the very edge of the Galaxy. Here, we propose NANCY: a completely public, all-sky survey that will create a high-value legacy dataset benefiting innumerable ongoing and forthcoming studies of the universe. NANCY is a pure expression of Roman's potential: it images the entire sky, at high spatial resolution, in a broad infrared bandpass that collects as many photons as possible. The majority of all ongoing astronomical surveys would benefit from incorporating observations of NANCY into their analyses, whether these surveys focus on nearby stars, the Milky Way, near-field cosmology, or the broader universe. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Comments: Submitted to the call for white papers for the Roman Core Community Survey (June 16th, 2023), and to the Bulletin of the AAS

The mass-loss rates of star clusters with stellar-mass black holes: implications for the globular cluster mass function

Authors: Mark Gieles, Oleg Y. Gnedin

Abstract: Stellar-mass black holes (BHs) can be retained in globular clusters (GCs) until the present. Simulations of GC evolution find that the relaxation driven mass-loss rate is elevated if BHs are present, especially near dissolution. We capture this behaviour in a parameterised mass-loss rate, benchmarked by results from $N$-body simulations, and use it to evolve an initial GC mass function (GCMF), sim… ▽ More Stellar-mass black holes (BHs) can be retained in globular clusters (GCs) until the present. Simulations of GC evolution find that the relaxation driven mass-loss rate is elevated if BHs are present, especially near dissolution. We capture this behaviour in a parameterised mass-loss rate, benchmarked by results from $N$-body simulations, and use it to evolve an initial GC mass function (GCMF), similar to that of young massive clusters in the Local Universe, to an age of 12 Gyr. Low-metallicity GCs ([Fe/H]$\lesssim-1.5$) have the highest mass-loss rates, because of their relatively high BH masses, which combined with their more radial orbits and stronger tidal field in the past explains the high turnover mass of the GCMF ($\sim10^5\,{\rm M}_\odot$) at large Galactic radii ($\gtrsim 10\,$kpc). The turnover mass at smaller Galactic radii is similar as the result of the upper mass truncation of the initial GCMF and the lower mass-loss rate because of the higher metallicities. The density profile in the Galaxy of mass lost from massive GCs ($\gtrsim10^{5}\,{\rm M}_\odot$) resembles that of nitrogen-rich stars in the halo, confirming that these stars originated from GCs. We conclude that two-body relaxation is the dominant effect in shaping the GCMF from a universal initial GCMF, because including the effect of BHs reduces the need for additional disruption mechanisms. △ Less

Submitted 26 April, 2023; v1 submitted 7 March, 2023; originally announced March 2023.

Comments: 18 pages, 16 figures, accepted for MNRAS

Formation of globular clusters in dwarf galaxies of the Local Group

Authors: Yingtian Chen, Oleg Y. Gnedin

Abstract: The existence of globular clusters (GCs) in a few satellite galaxies, and their absence in majority of dwarf galaxies, present a challenge for models attempting to understand the origins of GCs. In addition to GC presence appearing stochastic and difficult to describe with average trends, in the smallest satellite galaxies GCs contribute a substantial fraction of total stellar mass. We investigate… ▽ More The existence of globular clusters (GCs) in a few satellite galaxies, and their absence in majority of dwarf galaxies, present a challenge for models attempting to understand the origins of GCs. In addition to GC presence appearing stochastic and difficult to describe with average trends, in the smallest satellite galaxies GCs contribute a substantial fraction of total stellar mass. We investigate the stochasticity and number of GCs in dwarf galaxies using an updated version of our model that links the formation of GCs to the growth of the host galaxy mass. We find that more than 50% of dwarf galaxies with stellar mass $M_\odot\lesssim 2\times10^7 M_\odot$ do not host GCs, whereas dwarfs with $M_\odot\sim10^8 M_\odot$ almost always contain some GCs, with a median number $\sim 10$ at $z=0$. These predictions are in agreement with the observations of the Local Volume dwarfs. We also confirm the near-linear GC system mass--halo mass relation down to $M_{\rm h}\simeq10^8 M_\odot$ under the assumption that GC formation and evolution in galaxies of all mass can be described by the same physical model. A detailed case study of two model dwarfs that resemble the Fornax dwarf spheroidal galaxy shows that observational samples can be notably biased by incompleteness below detection limit and at large radii. △ Less

Submitted 22 May, 2023; v1 submitted 19 January, 2023; originally announced January 2023.

Comments: 16 pages, 9 figures, accepted for publication in MNRAS

Ultrafaint Dwarf Galaxy Candidates in the M81 Group: Signatures of Group Accretion

Authors: Eric F. Bell, Adam Smercina, Paul A. Price, Richard D'Souza, Jeremy Bailin, Roelof S. de Jong, Katya Gozman, In Sung Jang, Antonela Monachesi, Oleg Y. Gnedin, Colin T. Slater

Abstract: The faint and ultrafaint dwarf galaxies in the Local Group form the observational bedrock upon which our understanding of small-scale cosmology rests. In order to understand whether this insight generalizes, it is imperative to use resolved-star techniques to discover similarly faint satellites in nearby galaxy groups. We describe our search for ultrafaint galaxies in the M81 group using deep grou… ▽ More The faint and ultrafaint dwarf galaxies in the Local Group form the observational bedrock upon which our understanding of small-scale cosmology rests. In order to understand whether this insight generalizes, it is imperative to use resolved-star techniques to discover similarly faint satellites in nearby galaxy groups. We describe our search for ultrafaint galaxies in the M81 group using deep ground-based resolved-star data sets from Subaru's Hyper Suprime-Cam. We present one new ultrafaint dwarf galaxy in the M81 group and identify five additional extremely low surface brightness candidate ultrafaint dwarfs that reach deep into the ultrafaint regime to $M_V \sim -6$ (similar to current limits for Andromeda satellites). These candidates' luminosities and sizes are similar to known Local Group dwarf galaxies Tucana B, Canes Venatici I, Hercules, and Boötes I. Most of these candidates are likely to be real, based on tests of our techniques on blank fields. Intriguingly, all of these candidates are spatially clustered around NGC 3077, which is itself an M81 group satellite in an advanced state of tidal disruption. This is somewhat surprising, as M81 itself and its largest satellite M82 are both substantially more massive than NGC 3077 and by virtue of their greater masses, would have been expected to host as many or more ultrafaint candidates. These results lend considerable support to the idea that satellites of satellites are an important contribution to the growth of satellite populations around Milky Way-mass galaxies. △ Less

Submitted 13 September, 2022; originally announced September 2022.

Comments: The Astrophysical Journal Letters; in press. 11 pages, 4 figures, 1 table

Optical and X-ray GRB Fundamental Planes as Cosmological Distance Indicators

Authors: Maria Giovanna Dainotti, Via Nielson, Giuseppe Sarracino, Enrico Rinaldi, Shigehiro Nagataki, Salvatore Capozziello, Oleg Y. Gnedin, Giada Bargiacchi

Abstract: Gamma-Ray Bursts (GRBs), can be employed as standardized candles, extending the distance ladder beyond Supernovae Type Ia (SNe Ia, $z=2.26$). We standardize GRBs using the 3D fundamental plane relation (the Dainotti relation) among the rest-frame end time of the X-ray plateau emission, its corresponding luminosity, and the peak prompt luminosity. Combining SNe Ia and GRBs, we constrain… ▽ More Gamma-Ray Bursts (GRBs), can be employed as standardized candles, extending the distance ladder beyond Supernovae Type Ia (SNe Ia, $z=2.26$). We standardize GRBs using the 3D fundamental plane relation (the Dainotti relation) among the rest-frame end time of the X-ray plateau emission, its corresponding luminosity, and the peak prompt luminosity. Combining SNe Ia and GRBs, we constrain $Ω_{\text{M}}= 0.299 \pm 0.009$ assuming a flat $Λ$CDM cosmology with and without correcting GRBs for selection biases and redshift evolution. Using a 3D optical Dainotti correlation, we find this sample is as efficacious in the determination of $Ω_{\text{M}}$ as the X-ray sample. We trimmed our GRB samples to achieve tighter planes to simulate additional GRBs. We determined how many GRBs are needed as standalone probes to achieve a comparable precision on $Ω_{\text{M}}$ to the one obtained by SNe Ia only. We reach the same error measurements derived using SNe Ia in 2011 and 2014 with 142 and 284 simulated optical GRBs, respectively, considering the errorbars on the variables halved. These error limits will be reached in 2038 and in 2047, respectively. Using a doubled sample (obtained by future machine learning approaches allowing a lightcurve reconstruction and the estimates of GRB redhifts when z is unknown) compared to the current sample, with errorbars halved we will reach the same precision as SNe Ia in 2011 and 2014, now and in 2026, respectively. If we consider the current SNe precision, this will be reached with 390 optical GRBs by 2054. △ Less

Submitted 17 June, 2022; v1 submitted 29 March, 2022; originally announced March 2022.

Comments: 31 pages, 17 figures, 10 tables

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 514, Issue 2, August 2022, Pages 1828-1856

Modeling the kinematics of globular cluster systems

Authors: Yingtian Chen, Oleg Y. Gnedin

Abstract: Globular clusters (GCs) are old massive star clusters that serve as `fossils' of galaxy formation. The advent of Gaia observatory has enabled detailed kinematics studies of the Galactic GCs and revolutionized our understanding of the connections between GC properties and galaxy assembly. However, lack of kinematic measurements of extragalactic GCs limits the sample size of GC systems that we can f… ▽ More Globular clusters (GCs) are old massive star clusters that serve as `fossils' of galaxy formation. The advent of Gaia observatory has enabled detailed kinematics studies of the Galactic GCs and revolutionized our understanding of the connections between GC properties and galaxy assembly. However, lack of kinematic measurements of extragalactic GCs limits the sample size of GC systems that we can fully study. In this work, we present a model for GC formation and evolution, which includes positional and kinematic information of individual GCs by assigning them to particles in the Illustris TNG50-1 simulation based on age and location. We calibrate the three adjustable model parameters using observed properties of the Galactic and extragalactic GC systems, including the distributions of position, systemic velocity, velocity dispersion, anisotropy parameter, orbital actions, and metallicities. We also analyze the properties of GCs from different origins. In outer galaxy, ex-situ clusters are more dominant than the clusters formed in-situ. This leads to the GC metallicities decreasing outwards due to the increasing abundance of accreted, metal-poor clusters. We also find the ex-situ GCs to have greater velocity dispersions and orbital actions, in agreement with their accretion origin. △ Less

Submitted 10 July, 2022; v1 submitted 1 March, 2022; originally announced March 2022.

Comments: 20 pages, 17 figures, accepted for publication in MNRAS

Testing Feedback from Star Clusters in Simulations of the Milky Way Formation

Authors: Gillen Brown, Oleg Y. Gnedin

Abstract: We present a suite of galaxy formation simulations that directly model star cluster formation and disruption. Starting from a model previously developed by our group, here we introduce several improvements to the prescriptions for cluster formation and feedback, then test these updates using a large suite of cosmological simulations of Milky Way mass galaxies. We perform a differential analysis wi… ▽ More We present a suite of galaxy formation simulations that directly model star cluster formation and disruption. Starting from a model previously developed by our group, here we introduce several improvements to the prescriptions for cluster formation and feedback, then test these updates using a large suite of cosmological simulations of Milky Way mass galaxies. We perform a differential analysis with the goal of understanding how each of the updates affects star cluster populations. Two key parameters are the momentum boost of supernova feedback $f_{\mathrm{boost}}$ and star formation efficiency per freefall time $ε_{\mathrm{ff}}$. We find that $f_{\mathrm{boost}}$ has a strong influence on the galactic star formation rate, with higher values leading to less star formation. The efficiency $ε_{\mathrm{ff}}$ does not have a significant impact on the global star formation rate, but dramatically changes cluster properties, with increasing $ε_{\mathrm{ff}}$ leading to a higher maximum cluster mass, shorter age spread of stars within clusters, and higher integrated star formation efficiencies. We also explore the redshift evolution of the observable cluster mass function, finding that most massive clusters have formed at high redshift $z>4$. Extrapolation of cluster disruption to $z=0$ produces good agreement with both the Galactic globular cluster mass function and age-metallicity relation. Our results emphasize the importance of using small-scale properties of galaxies to calibrate subgrid models of star cluster formation and feedback. △ Less

Submitted 25 April, 2022; v1 submitted 1 March, 2022; originally announced March 2022.

Comments: Accepted for publication in MNRAS

Tidal disruption of star clusters in galaxy formation simulations

Authors: Xi Meng, Oleg Y. Gnedin

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

Submitted 12 August, 2022; v1 submitted 24 January, 2022; originally announced January 2022.

Comments: 13 pages, 13 figures, accepted by MNRAS

Faint objects in motion: the new frontier of high precision astrometry

Authors: Fabien Malbet, Céline Boehm, Alberto Krone-Martins, Antonio Amorim, Guillem Anglada-Escudé, Alexis Brandeker, Frédéric Courbin, Torsten Enßlin, Antonio Falcão, Katherine Freese, Berry Holl, Lucas Labadie, Alain Léger, Gary Mamon, Barbara Mcarthur, Alcione Mora, Mike Shao, Alessandro Sozzetti, Douglas Spolyar, Eva Villaver, Ummi Abbas, Conrado Albertus, João Alves, Rory Barnes, Aldo Stefano Bonomo , et al. (61 additional authors not shown)

Abstract: Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the front… ▽ More Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the frontier of precision astrometry from evidence of Earth-mass habitable worlds around the nearest stars, to distant Milky Way objects, and out to the Local Group of galaxies. As we enter the era of the James Webb Space Telescope and the new ground-based, adaptive-optics-enabled giant telescopes, by obtaining these high precision measurements on key objects that Gaia could not reach, a mission that focuses on high precision astrometry science can consolidate our theoretical understanding of the local Universe, enable extrapolation of physical processes to remote redshifts, and derive a much more consistent picture of cosmological evolution and the likely fate of our cosmos. Already several missions have been proposed to address the science case of faint objects in motion using high precision astrometry missions: NEAT proposed for the ESA M3 opportunity, micro-NEAT for the S1 opportunity, and Theia for the M4 and M5 opportunities. Additional new mission configurations adapted with technological innovations could be envisioned to pursue accurate measurements of these extremely small motions. The goal of this White Paper is to address the fundamental science questions that are at stake when we focus on the motions of faint sky objects and to briefly review instrumentation and mission profiles. △ Less

Submitted 16 November, 2021; originally announced November 2021.

Comments: arXiv admin note: substantial text overlap with arXiv:1910.08028, arXiv:1707.01348

Journal ref: Experimental Astronomy, Springer Link, 2021, 51 (3), pp.845-886

Reconstruction of the Density Power Spectrum from Quasar Spectra using Machine Learning

Authors: Maria Han Veiga, Xi Meng, Oleg Y. Gnedin, Nickolay Y. Gnedin, Xun Huan

Abstract: We describe a novel end-to-end approach using Machine Learning to reconstruct the power spectrum of cosmological density perturbations at high redshift from observed quasar spectra. State-of-the-art cosmological simulations of structure formation are used to generate a large synthetic dataset of line-of-sight absorption spectra paired with 1-dimensional fluid quantities along the same line-of-sigh… ▽ More We describe a novel end-to-end approach using Machine Learning to reconstruct the power spectrum of cosmological density perturbations at high redshift from observed quasar spectra. State-of-the-art cosmological simulations of structure formation are used to generate a large synthetic dataset of line-of-sight absorption spectra paired with 1-dimensional fluid quantities along the same line-of-sight, such as the total density of matter and the density of neutral atomic hydrogen. With this dataset, we build a series of data-driven models to predict the power spectrum of total matter density. We are able to produce models which yield reconstruction to accuracy of about 1% for wavelengths $k \leq 2 h Mpc^{-1}$, while the error increases at larger $k$. We show the size of data sample required to reach a particular error rate, giving a sense of how much data is necessary to reach a desired accuracy. This work provides a foundation for developing methods to analyse very large upcoming datasets with the next-generation observational facilities. △ Less

Submitted 19 July, 2021; originally announced July 2021.

Comments: 10 pages, 9 figures

Radii of Young Star Clusters in Nearby Galaxies

Authors: Gillen Brown, Oleg Y. Gnedin

Abstract: We measure the projected half-light radii of young star clusters in 31 galaxies from the Legacy Extragalactic UV Survey (LEGUS). We implement a custom pipeline specifically designed to be robust against contamination, which allows us to measure radii for 6097 clusters. This is the largest sample of young star cluster radii currently available. We find that most (but not all) galaxies share a commo… ▽ More We measure the projected half-light radii of young star clusters in 31 galaxies from the Legacy Extragalactic UV Survey (LEGUS). We implement a custom pipeline specifically designed to be robust against contamination, which allows us to measure radii for 6097 clusters. This is the largest sample of young star cluster radii currently available. We find that most (but not all) galaxies share a common cluster radius distribution, with the peak at around 3 pc. We find a clear mass-radius relation of the form $R_{\mathrm{eff}} \propto M^{0.24}$. This relation is present at all cluster ages younger than 1 Gyr, but with a shallower slope for clusters younger than 10 Myr. We present simple toy models to interpret these age trends, finding that high-mass clusters are more likely to be not tidally limited and expand. We also find that most clusters in LEGUS are gravitationally bound, especially at older ages or higher masses. Lastly, we present the cluster density and surface density distributions, finding a large scatter that appears to decrease with cluster age. The youngest clusters have a typical surface density of 100 $M_\odot$ pc$^{-2}$. △ Less

Submitted 5 October, 2021; v1 submitted 23 June, 2021; originally announced June 2021.

Comments: Accepted by MNRAS. Cluster catalogs available at https://gillenbrown.github.io/LEGUS-sizes/

Improving Performance of Zoom-In Cosmological Simulations using Initial Conditions with Customized Grids

Authors: Gillen Brown, Oleg Y. Gnedin

Abstract: We present a method for customizing the root grid of zoom-in initial conditions used for simulations of galaxy formation. Starting from the white noise used to seed the structures of an existing initial condition, we cut out a smaller region of interest and use this trimmed white noise cube to create a new root grid. This new root grid contains similar structures as the original, but allows for a… ▽ More We present a method for customizing the root grid of zoom-in initial conditions used for simulations of galaxy formation. Starting from the white noise used to seed the structures of an existing initial condition, we cut out a smaller region of interest and use this trimmed white noise cube to create a new root grid. This new root grid contains similar structures as the original, but allows for a smaller box volume and different grid resolution that can be tuned to best suit a given simulation code. To minimally disturb the zoom region, the dark matter particles and gas cells from the original zoom region are placed within the new root grid, with no modification other than a bulk velocity offset to match the systemic velocity of the corresponding region in the new root grid. We validate this method using a zoom-in initial condition containing a Local Group analog. We run collisionless simulations using the original and modified initial conditions, finding good agreement. The dark matter halo masses of the two most massive galaxies at $z=0$ match the original to within 15%. The times and masses of major mergers are reproduced well, as are the full dark matter accretion histories. While we do not reproduce specific satellite galaxies found in the original simulation, we obtain qualitative agreement in the distributions of the maximum circular velocity and the distance from the central galaxy. We also examine the runtime speedup provided by this method for full hydrodynamic simulations with the ART code. We find that reducing the root grid cell size improves performance, but the increased particle and cell numbers can negate some of the gain. We test several realizations, with our best runs achieving a speedup of nearly a factor of two. △ Less

Submitted 18 September, 2020; originally announced September 2020.

Comments: Accepted for publication in New Astronomy

On the response of a star cluster to a tidal perturbation

Authors: Luis A. Martinez-Medina, Mark Gieles, Oleg Y. Gnedin, Hui Li

Abstract: We study the response of star clusters to individual tidal perturbations using controlled $N$-body simulations. We consider perturbations by a moving point mass and by a disc, and vary the duration of the perturbation as well as the cluster density profile. For fast perturbations (i.e. `shocks'), the cluster gains energy in agreement with theoretical predictions in the impulsive limit. For slow di… ▽ More We study the response of star clusters to individual tidal perturbations using controlled $N$-body simulations. We consider perturbations by a moving point mass and by a disc, and vary the duration of the perturbation as well as the cluster density profile. For fast perturbations (i.e. `shocks'), the cluster gains energy in agreement with theoretical predictions in the impulsive limit. For slow disc perturbations, the energy gain is lower, and this has previously been attributed to adiabatic damping. However, the energy gain due to slow perturbations by a point-mass is similar to, or larger than that due to fast shocks, which is not expected because adiabatic damping should be almost independent of the nature of the tides. We show that the geometric distortion of the cluster during slow perturbations is of comparable importance for the energy gain as adiabatic damping, and that the combined effect can qualitatively explain the results. The half-mass radius of the bound stars after a shock increases up to $\sim$7\% for low-concentration clusters, and decreases $\sim$3\% for the most concentrated ones. The fractional mass loss is a non-linear function of the energy gain, and depends on the nature of the tides and most strongly on the cluster density profile, making semi-analytic model predictions for cluster lifetimes extremely sensitive to the adopted density profile. △ Less

Submitted 10 August, 2022; v1 submitted 14 September, 2020; originally announced September 2020.

Comments: 13 pages, 15 figures. Accepted for publication in MNRAS

Evolution of Disc Thickness in Simulated High-Redshift Galaxies

Authors: Xi Meng, Oleg Y. Gnedin

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

Submitted 13 August, 2022; v1 submitted 18 June, 2020; originally announced June 2020.

Comments: 9 pages, 8 figures, accepted for publication in MNRAS

The Strikingly Metal-Rich Halo of the Sombrero Galaxy

Authors: Roger E. Cohen, Paul Goudfrooij, Matteo Correnti, Oleg Y. Gnedin, William E. Harris, Rupali Chandar, Thomas H. Puzia, Ruben Sanchez-Janssen

Abstract: The nature of the Sombrero galaxy (M 104 = NGC 4594) has remained elusive despite many observational studies at a variety of wavelengths. Here we present Hubble Space Telescope imaging of two fields at $\sim$16 and 33 kpc along the minor axis to examine stellar metallicity gradients in the extended spheroid. We use this imaging, extending more than 2 mag below the tip of the red giant branch (TRGB… ▽ More The nature of the Sombrero galaxy (M 104 = NGC 4594) has remained elusive despite many observational studies at a variety of wavelengths. Here we present Hubble Space Telescope imaging of two fields at $\sim$16 and 33 kpc along the minor axis to examine stellar metallicity gradients in the extended spheroid. We use this imaging, extending more than 2 mag below the tip of the red giant branch (TRGB), in combination with artificial star tests to forward model observed color-magnitude diagrams (CMDs), measuring metallicity distribution functions (MDFs) at different radii along the minor axis. An important and unexpected result is that the halo of the Sombrero is strikingly metal-rich: even the outer field, located at $\sim$17 effective radii of the bulge, has a median metallicity [Z/H]$\sim$-0.15 and the fraction of stars with [Z/H]<-1.0 is negligible. This is unprecedented among massive galaxy halos studied to date, even among giant ellipticals. We find significant radial metallicity gradients, characterized by an increase in the fraction of metal-poor stars with radius and a gradient in median metallicity of $\sim$-0.01 dex/kpc. The density profile is well fit by power laws with slopes that exhibit a dependence on metallicity, with flatter slopes for more metal-poor stars. We discuss our results in the context of recent stellar MDF studies of other nearby galaxies and potential formation scenarios for the Sombrero galaxy. △ Less

Submitted 7 January, 2020; v1 submitted 6 January, 2020; originally announced January 2020.

Comments: ApJ accepted

ESA Voyage 2050 white paper -- Faint objects in motion: the new frontier of high precision astrometry

Authors: F. Malbet, U. Abbas, J. Alves, C. Boehm, W. Brown, L. Chemin, A. Correia, F. Courbin, J. Darling, A. Diaferio, M. Fortin, M. Fridlund, O. Gnedin, B. Holl, A. Krone-Martins, A. Léger, L. Labadie, J. Laskar, G. Mamon, B. McArthur, D. Michalik, A. Moitinho, M. Oertel, L. Ostorero, J. Schneider , et al. (6 additional authors not shown)

Abstract: Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the front… ▽ More Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the frontier of precision astrometry from evidence of earth-massed habitable worlds around the nearest starts, and also into distant Milky way objects up to the Local Group of galaxies. As we enter the era of the James Webb Space Telescope and the new ground-based, adaptive-optics-enabled giant telescopes, by obtaining these high precision measurements on key objects that Gaia could not reach, a mission that focuses on high precision astrometry science can consolidate our theoretical understanding of the local universe, enable extrapolation of physical processes to remote redshifts, and derive a much more consistent picture of cosmological evolution and the likely fate of our cosmos. Already several missions have been proposed to address the science case of faint objects in motion using high precision astrometry ESA missions: NEAT for M3, micro-NEAT for S1 mission, and Theia for M4 and M5. Additional new mission configurations adapted with technological innovations could be envisioned to pursue accurate measurements of these extremely small motions. The goal of this white paper is to address the fundamental science questions that are at stake when we focus on the motions of faint sky objects and to briefly review quickly instrumentation and mission profiles. △ Less

Submitted 12 October, 2019; originally announced October 2019.

Comments: White paper for the Voyage 2050 long-term plan in the ESA Science Programme. arXiv admin note: substantial text overlap with arXiv:1707.01348

Origin of Giant Stellar Clumps in High-Redshift Galaxies

Authors: Xi Meng, Oleg Y. Gnedin

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

Submitted 13 August, 2022; v1 submitted 8 October, 2019; originally announced October 2019.

Comments: 13 pages, 5 figures, accepted for publication in MNRAS

Formation and evolution of globular clusters in cosmological simulations

Authors: Hui Li, Oleg Gnedin

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

Submitted 2 August, 2019; originally announced August 2019.

Comments: 6 pages, 4 figures, to appear in Proceedings IAU Symposium No. 351 Star Clusters: From the Milky Way to the Early Universe, 2019

Journal ref: Proc. IAU 14 (2019) 34-39

High Resolution Optical Spectroscopy of Stars in the Sylgr Stellar Stream

Authors: Ian U. Roederer, Oleg Y. Gnedin

Abstract: We observe two metal-poor main sequence stars that are members of the recently-discovered Sylgr stellar stream. We present radial velocities, stellar parameters, and abundances for 13 elements derived from high-resolution optical spectra collected using the Magellan Inamori Kyocera Echelle spectrograph. The two stars have identical compositions (within 0.13 dex or 1.2 sigma) among all elements det… ▽ More We observe two metal-poor main sequence stars that are members of the recently-discovered Sylgr stellar stream. We present radial velocities, stellar parameters, and abundances for 13 elements derived from high-resolution optical spectra collected using the Magellan Inamori Kyocera Echelle spectrograph. The two stars have identical compositions (within 0.13 dex or 1.2 sigma) among all elements detected. Both stars are very metal poor ([Fe/H] = -2.92 +/- 0.06). Neither star is highly enhanced in C ([C/Fe] < +1.0). Both stars are enhanced in the alpha elements Mg, Si, and Ca ([alpha/Fe] = +0.32 +/- 0.06), and ratios among Na, Al, and all Fe-group elements are typical for other stars in the halo and ultra-faint and dwarf spheroidal galaxies at this metallicity. Sr is mildly enhanced ([Sr/Fe] = +0.22 +/- 0.11), but Ba is not enhanced ([Ba/Fe] < -0.4), indicating that these stars do not contain high levels of neutron-capture elements. The Li abundances match those found in metal-poor unevolved field stars and globular clusters (log epsilon (Li) = 2.05 +/- 0.07), which implies that environment is not a dominant factor in determining the Li content of metal-poor stars. The chemical compositions of these two stars cannot distinguish whether the progenitor of the Sylgr stream was a dwarf galaxy or a globular cluster. If the progenitor was a dwarf galaxy, the stream may originate from a dense region such as a nuclear star cluster. If the progenitor was a globular cluster, it would be the most metal-poor globular cluster known. △ Less

Submitted 8 July, 2019; originally announced July 2019.

Comments: Accepted for publication in the Astrophysical Journal (17 pages, 4 figures, 3 tables)

Origins of scaling relations of globular cluster systems

Authors: Nick Choksi, Oleg Y. Gnedin

Abstract: Globular cluster (GC) systems demonstrate tight scaling relations with the properties of their host galaxies. In previous work, we developed an analytic model for GC formation in a cosmological context and showed that it matches nearly all of the observed scaling relations across 4 orders of magnitude in host galaxy mass. Motivated by the success of this model, we investigate in detail the physica… ▽ More Globular cluster (GC) systems demonstrate tight scaling relations with the properties of their host galaxies. In previous work, we developed an analytic model for GC formation in a cosmological context and showed that it matches nearly all of the observed scaling relations across 4 orders of magnitude in host galaxy mass. Motivated by the success of this model, we investigate in detail the physical origins and evolution of these scaling relations. The ratio of the combined mass in GCs $M_{\rm GC}$ to the host dark matter halo mass $M_h$ is nearly constant at all redshifts, but its normalization evolves by a factor of $\sim$10 from birth to $z=0$. The relation is steeper than linear at halo masses $M_h \lesssim 10^{11.5} M_{\odot}$, primarily due to non-linearity in the stellar mass-halo mass relation. The near constancy of the ratio $M_{\rm GC}/M_h$, combined with the shape of the stellar mass-halo mass relation, sets the characteristic $U-$shape of the GC specific frequency as a function of host galaxy mass. The contribution of accreted satellite galaxies to the buildup of GC systems is a strong function of the host galaxy mass, ranging from $\approx$0% at $M_h \approx 10^{11} M_{\odot}$ to 80% at $M_h \approx 10^{15} M_{\odot}$. The metal-poor clusters are significantly more likely to form ex-situ relative to the metal-rich clusters, but a substantial fraction of metal-poor clusters still form in-situ in lower mass galaxies. Similarly, the fraction of red clusters increases from $\approx 10$% at $M_h = 10^{11} M_{\odot}$ to $\approx 60$% at $M_h \approx 10^{13} M_{\odot}$, and flattens at higher $M_h$. Clusters formation occurs essentially continuously at high redshift, while at low redshift galactic mergers become increasingly important for cluster formation. △ Less

Submitted 13 May, 2019; originally announced May 2019.

Comments: 10 pages, 9 figures. Comments welcome

Disruption of giant molecular clouds and formation of bound star clusters under the influence of momentum stellar feedback

Authors: Hui Li, Mark Vogelsberger, Federico Marinacci, Oleg Y. Gnedin

Abstract: Energetic feedback from star clusters plays a pivotal role in shaping the dynamical evolution of giant molecular clouds (GMCs). To study the effects of stellar feedback on the star formation efficiency of the clouds and the dynamical response of embedded star clusters, we perform a suite of isolated GMC simulations with star formation and momentum feedback subgrid models using the moving-mesh hydr… ▽ More Energetic feedback from star clusters plays a pivotal role in shaping the dynamical evolution of giant molecular clouds (GMCs). To study the effects of stellar feedback on the star formation efficiency of the clouds and the dynamical response of embedded star clusters, we perform a suite of isolated GMC simulations with star formation and momentum feedback subgrid models using the moving-mesh hydrodynamics code \textsc{Arepo}. The properties of our simulated GMCs span a wide range of initial mass, radius, and velocity configurations. We find that the ratio of the final stellar mass to the total cloud mass, $ε_{\rm int}$, scales strongly with the initial cloud surface density and momentum feedback strength. This correlation is explained by an analytic model that considers force balancing between gravity and momentum feedback. For all simulated GMCs, the stellar density profiles are systematically steeper than that of the gas at the epochs of the peaks of star formation, suggesting a centrally concentrated stellar distribution. We also find that star clusters are always in a sub-virial state with a virial parameter $\sim0.6$ prior to gas expulsion. Both the sub-virial dynamical state and steeper stellar density profiles prevent clusters from dispersal during the gas removal phase of their evolution. The final cluster bound fraction is a continuously increasing function of $ε_{\rm int}$. GMCs with star formation efficiency smaller than 0.5 are still able to form clusters with large bound fractions. △ Less

Submitted 19 May, 2019; v1 submitted 26 April, 2019; originally announced April 2019.

Comments: 17 pages, 13 figures; Match the published version in MNRAS

All-Sky Near Infrared Space Astrometry

Authors: Barbara McArthur, David Hobbs, Erik Høg, Valeri Makarov, Alessandro Sozzetti, Anthony Brown, Alberto Krone Martins, Jennifer Lynn Bartlett, John Tomsick, Mike Shao, Fritz Benedict, Eduardo Bendek, Celine Boehm, Charlie Conroy, Johan Peter Uldall Fynbo, Oleg Gnedin, Lynne Hillenbrand, Lennart Lindegren, David R. Rodriguez, Rick White, Slava Turyshev, Stephen Unwin, ChengXing Zhai

Abstract: Gaia is currently revolutionizing modern astronomy. However, much of the Galactic plane, center and the spiral arm regions are obscured by interstellar extinction, rendering them inaccessible because Gaia is an optical instrument. An all-sky near infrared (NIR) space observatory operating in the optical NIR, separated in time from the original Gaia would provide microarcsecond NIR astrometry and m… ▽ More Gaia is currently revolutionizing modern astronomy. However, much of the Galactic plane, center and the spiral arm regions are obscured by interstellar extinction, rendering them inaccessible because Gaia is an optical instrument. An all-sky near infrared (NIR) space observatory operating in the optical NIR, separated in time from the original Gaia would provide microarcsecond NIR astrometry and millimag photometry to penetrate obscured regions unraveling the internal dynamics of the Galaxy. △ Less

Submitted 18 April, 2019; originally announced April 2019.

Comments: 7 pages

Journal ref: NASA Astro 2020 White Paper

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

Authors: Hui Li, Oleg Y. Gnedin

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

Submitted 19 May, 2019; v1 submitted 25 October, 2018; originally announced October 2018.

Comments: 15 pages, 10 figures. Match the published version in MNRAS

Structure and stability of high-redshift galaxies in cosmological simulations

Authors: Xi Meng, Oleg Y. Gnedin, Hui Li

Abstract: We investigate the structure of galaxies formed in a suite of high-resolution cosmological simulations. Consistent with observations of high-redshift galaxies, our simulated galaxies show irregular, prolate shapes, which are dominated by turbulent motions instead of rotation. Yet molecular gas and young stars are restricted to a relatively thin plane. We examine the accuracy of applying the Toomre… ▽ More We investigate the structure of galaxies formed in a suite of high-resolution cosmological simulations. Consistent with observations of high-redshift galaxies, our simulated galaxies show irregular, prolate shapes, which are dominated by turbulent motions instead of rotation. Yet molecular gas and young stars are restricted to a relatively thin plane. We examine the accuracy of applying the Toomre linear stability analysis to predict the location and amount of gas available for star formation. We find that the Toomre criterion still works for these irregular galaxies, after correcting for multiple gas and stellar components: the $Q$ parameter in $\rm{H_2}$ rich regions is in the range $0.5-1$, remarkably close to unity. Due to the violent stellar feedback from supernovae and strong turbulent motions, young stars and molecular gas are not always spatially associated. Neither the $Q$ map nor the $\rm{H_2}$ surface density map coincide with recent star formation exactly. We argue that the Toomre criterion is a better indicator of future star formation than a single $\rm{H_2}$ surface density threshold because of the smaller dynamic range of $Q$. The depletion time of molecular gas is below 1~Gyr on kpc scale, but with large scatter. Centering the aperture on density peaks of gas/young stars systematically biases the depletion time to larger/smaller values and increases the scatter. △ Less

Submitted 13 August, 2022; v1 submitted 15 October, 2018; originally announced October 2018.

Comments: 17 pages, 11 figures; accepted for publication in MNRAS

Formation of Globular Cluster Systems II: Impact of the Cutoff of the Cluster Initial Mass Function

Authors: Nick Choksi, Oleg Y. Gnedin

Abstract: Observations of young star clusters reveal that the high-mass end of the cluster initial mass function (CIMF) deviates from a pure power-law and instead truncates exponentially. We investigate the effects of this truncation on the formation of globular cluster (GC) systems by updating our analytic model for cluster formation and evolution, which is based on dark matter halo merger trees coupled to… ▽ More Observations of young star clusters reveal that the high-mass end of the cluster initial mass function (CIMF) deviates from a pure power-law and instead truncates exponentially. We investigate the effects of this truncation on the formation of globular cluster (GC) systems by updating our analytic model for cluster formation and evolution, which is based on dark matter halo merger trees coupled to empirical galactic scaling relations, and has been shown in previous work to match a wide array of observational data. The cutoff masses of $M_c=10^{6.5} M_{\odot}$ or $10^{7}M_{\odot}$ match many scaling relations: between the GC system mass and host halo mass, between the average metallicity of the GC system and host halo mass, and the distribution of cluster masses. This range of $M_c$ agrees with indirect measurements from extragalactic GC systems. Models with $M_c<10^{6.5}M_{\odot}$ cannot reproduce the observed GC metallicity and mass distributions in massive galaxies. The slope of the mass-metallicity relation for metal-poor clusters (blue tilt) for all $M_c$ models is consistent with observations within their errors, when measured using the same method. We introduce an alternative, more robust fitting method, which reveals a trend of increasing tilt slope for lower $M_c$. In our model the blue tilt arises because the metal-poor clusters form in relatively low-mass galaxies which lack sufficient cold gas to sample the CIMF at highest masses. Massive blue clusters form in progressively more massive galaxies and inherit their higher metallicity. The metal-rich clusters do not exhibit such a tilt because they form in significantly more massive galaxies, which have enough cold gas to fully sample the CIMF. △ Less

Submitted 20 March, 2019; v1 submitted 3 October, 2018; originally announced October 2018.

Comments: 14 pages, 11 figures. Updated to match MNRAS accepted version. Minor edits to text and new appendix on the impact of IMF sampling methods

The star clusters that make black hole binaries across cosmic time

Authors: Nick Choksi, Marta Volonteri, Monica Colpi, Oleg Y. Gnedin, Hui Li

Abstract: We explore the properties of dense star clusters that are likely to be the nurseries of stellar black holes pairing in close binaries. We combine a cosmological model of globular cluster formation with analytic prescriptions for the dynamical assembly and evolution of black hole binaries to constrain which types of clusters are most likely to form binaries tight enough to coalesce within a Hubble… ▽ More We explore the properties of dense star clusters that are likely to be the nurseries of stellar black holes pairing in close binaries. We combine a cosmological model of globular cluster formation with analytic prescriptions for the dynamical assembly and evolution of black hole binaries to constrain which types of clusters are most likely to form binaries tight enough to coalesce within a Hubble time. We find that black hole binaries which are ejected and later merge ex-situ form in clusters of a characteristic mass $M_{\rm cl} \sim 10^{5.3}M_{\odot}$, whereas binaries which merge in-situ form in more massive clusters, $M_{\rm cl} \sim 10^{5.7}M_{\odot}$. The clusters which dominate the production of black hole binaries are similar in age and metallicity to the entire population. Finally, we estimate an approximate cosmic black hole merger rate of dynamically assembled binaries using the mean black hole mass for each cluster given its metallicity. We find an intrinsic rate of $\sim 6\,\mathrm{Gpc\,yr^{-1}}$ at $z=0$, a weakly increasing merger rate out to $z=1.5$, and then a decrease out to $z=4$. Our results can be used to provide a cosmological context and choose initial conditions in numerical studies of black hole binaries assembled in star clusters. △ Less

Submitted 4 March, 2019; v1 submitted 4 September, 2018; originally announced September 2018.

Comments: 15 pages, 12 figures. Updated to match ApJ accepted version. Minor edits to text relative to v1

Millisecond pulsars and the gamma-ray excess in Andromeda

Authors: Giacomo Fragione, Fabio Antonini, Oleg Y. Gnedin

Abstract: The Fermi Gamma-Ray Space Telescope has provided evidence for diffuse gamma-ray emission in the central parts of the Milky Way and the Andromeda galaxy. This excess has been interpreted either as dark matter annihilation emission or as emission from thousands of millisecond pulsars (MSPs). We have recently shown that old massive globular clusters may move towards the center of the Galaxy by dynami… ▽ More The Fermi Gamma-Ray Space Telescope has provided evidence for diffuse gamma-ray emission in the central parts of the Milky Way and the Andromeda galaxy. This excess has been interpreted either as dark matter annihilation emission or as emission from thousands of millisecond pulsars (MSPs). We have recently shown that old massive globular clusters may move towards the center of the Galaxy by dynamical friction and carry within them enough MSPs to account for the observed gamma-ray excess. In this paper we revisit the MSP scenario for the Andromeda galaxy, by modeling the formation and disruption of its globular cluster system. We find that our model predicts gamma-ray emission $\sim 2-3$ times larger than for the Milky Way, but still nearly an order of magnitude smaller than the observed Fermi excess in the Andromeda. Our MSP model can reproduce the observed excess only by assuming $\sim 8$ times larger number of old clusters than inferred from galaxy scaling relations. To explain the observations we require either that Andromeda deviates significantly from the scaling relations, or that a large part of its high-energy emission comes from additional sources. △ Less

Submitted 6 January, 2019; v1 submitted 7 August, 2018; originally announced August 2018.

Comments: 5 pages, 2 figures, 1 Table, accepted by ApJ Lett

Nuclear Star Clusters in Cosmological Simulations

Authors: Gillen Brown, Oleg Y. Gnedin, Hui Li

Abstract: We investigate the possible connection between the most massive globular clusters, such as $ω$ Cen and M54, and nuclear star clusters of dwarf galaxies that exhibit similar spreads in age and metallicity. We examine galactic nuclei in cosmological galaxy formation simulations at $z \approx 1.5$ to explore whether their age and metallicity spreads could explain these massive globular clusters. We d… ▽ More We investigate the possible connection between the most massive globular clusters, such as $ω$ Cen and M54, and nuclear star clusters of dwarf galaxies that exhibit similar spreads in age and metallicity. We examine galactic nuclei in cosmological galaxy formation simulations at $z \approx 1.5$ to explore whether their age and metallicity spreads could explain these massive globular clusters. We derive structural properties of these nuclear regions, including mass, size, rotation, and shape. By using theoretical supernova yields to model the supernova enrichment in the simulations, we obtain individual elemental abundances for Fe, O, Na, Mg, and Al. Our nuclei are systematically more metal-rich than their host galaxies, which lie on the expected mass-metallicity relation. Some nuclei have a spread in Fe and age comparable to the massive globular clusters of the Milky Way, lending support to the hypothesis that nuclear star clusters of dwarf galaxies could be the progenitors of these objects. None of our nuclear regions contain the light element abundance spreads that characterize globular clusters, even when a large age spread is present. Our results demonstrate that extended star formation history within clusters, with metal pollution provided solely by supernova ejecta, is capable of replicating the metallicity spreads of massive globular clusters, but still requires another polluter to produce the light element variations. △ Less

Submitted 23 July, 2018; v1 submitted 25 April, 2018; originally announced April 2018.

Comments: 10 pages, 12 figures, accepted to ApJ

Journal ref: Brown et. al. 2018, ApJ, 864, 94

Globular Cluster Formation and Evolution in the Context of Cosmological Galaxy Assembly: Open Questions

Authors: Duncan A. Forbes, Nate Bastian, Mark Gieles, Robert A. Crain, J. M. Diederik Kruijssen, Søren S. Larsen, Sylvia Ploeckinger, Oscar Agertz, Michele Trenti, Annette M. N. Ferguson, Joel Pfeffer, Oleg Y. Gnedin

Abstract: We discuss some of the key open questions regarding the formation and evolution of globular clusters (GCs) during galaxy formation and assembly within a cosmological framework. The current state-of-the-art for both observations and simulations is described, and we briefly mention directions for future research. The oldest GCs have ages $\ge$ 12.5 Gyr and formed around the time of reionisation. Res… ▽ More We discuss some of the key open questions regarding the formation and evolution of globular clusters (GCs) during galaxy formation and assembly within a cosmological framework. The current state-of-the-art for both observations and simulations is described, and we briefly mention directions for future research. The oldest GCs have ages $\ge$ 12.5 Gyr and formed around the time of reionisation. Resolved colour-magnitude diagrams of Milky Way GCs and direct imaging of lensed proto-GCs at z $\sim$ 6 with JWST promise further insight. Globular clusters are known to host multiple populations of stars with variations in their chemical abundances. Recently, such multiple populations have been detected in $\sim$2 Gyr old compact, massive star clusters. This suggests a common, single pathway for the formation of GCs at high and low redshift. The shape of the initial mass function for GCs remains unknown, however for massive galaxies a power-law mass function is favoured. Significant progress has been made recently modelling GC formation in the context of galaxy formation, with success in reproducing many of the observed GC-galaxy scaling relations. △ Less

Submitted 30 January, 2018; v1 submitted 17 January, 2018; originally announced January 2018.

Comments: 28 pages, accepted for publication in the Proceedings of the Royal Society A journal. Version 2 with updated references

Formation of Globular Cluster Systems: From Dwarf Galaxies to Giants

Authors: Nick Choksi, Oleg Gnedin, Hui Li

Abstract: Globular cluster (GC) systems around galaxies of a vast mass range show remarkably simple scaling relations. The combined mass of all GCs is a constant fraction of the total galaxy mass and the mean metallicity and metallicity dispersion of the GC system scale up weakly with galaxy mass. The metallicity of massive, metal-poor ("blue") clusters increases with cluster mass, while that of metal-rich… ▽ More Globular cluster (GC) systems around galaxies of a vast mass range show remarkably simple scaling relations. The combined mass of all GCs is a constant fraction of the total galaxy mass and the mean metallicity and metallicity dispersion of the GC system scale up weakly with galaxy mass. The metallicity of massive, metal-poor ("blue") clusters increases with cluster mass, while that of metal-rich ("red") clusters does not. A significant age-metallicity relation emerges from analysis of resolved stellar populations in Galactic GCs and unresolved populations in nearby galaxies. Remarkably, all these trends can be explained by a simple merger-based model developed in previous work and updated here using recent observations of galaxy scaling relations at high redshift. We show that the increasing dispersion of GC metallicity distributions with galaxy mass is a robust prediction of the model. It arises from more massive galaxies having more mergers that combine satellite GC systems. The average metallicity also increases by 0.6~dex over 3~dex in halo mass. The models show a non-linear trend between the GC system mass and host galaxy mass which is consistent with the data. The model does not consider GC self-enrichment, yet predicts a correlation between cluster mass and metallicity for massive blue clusters. The age-metallicity relation is another robust prediction of the model. Half of all clusters are predicted to form within the redshift range $5<z<2.3$, corresponding to ages of $10.8-12.5$~Gyr, in halos of masses $10^{11}-10^{12.5} M_{\odot}$. △ Less

Submitted 23 August, 2018; v1 submitted 10 January, 2018; originally announced January 2018.

Comments: 15 pages, 10 figures. Current version matches MNRAS-accepted submission. For supplementary data (model catalogs and compilation of obs. data) see http://ugastro.berkeley.edu/~nchoksi/cgl18_supplemental.zip

Journal ref: 2018, MNRAS, 480, 2343-2356

Star Cluster Formation in Cosmological Simulations. II. Effects of Star Formation Efficiency and Stellar Feedback

Authors: Hui Li, Oleg Y. Gnedin, Nickolay Y. Gnedin

Abstract: The implementation of star formation and stellar feedback in cosmological simulations plays a critical role in shaping galaxy properties. In the first paper of the series, we presented a new method to model star formation as a collection of star clusters. In this paper, we improve the algorithm by eliminating accretion gaps, boosting momentum feedback, and introducing a subgrid initial bound fract… ▽ More The implementation of star formation and stellar feedback in cosmological simulations plays a critical role in shaping galaxy properties. In the first paper of the series, we presented a new method to model star formation as a collection of star clusters. In this paper, we improve the algorithm by eliminating accretion gaps, boosting momentum feedback, and introducing a subgrid initial bound fraction, $f_i$, that distinguishes cluster mass from stellar particle mass. We perform a suite of simulations with different star formation efficiency per freefall time $ε_{\rm ff}$ and supernova momentum feedback intensity $f_{\rm boost}$. We find that the star formation history of a Milky Way-sized galaxy is sensitive to $f_{\rm boost}$, which allows us to constrain its value, $f_{\rm boost}\approx5$, in the current simulation setup. Changing $ε_{\rm ff}$ from a few percent to 200\% has little effect on global galaxy properties. However, on smaller scales, the properties of star clusters are very sensitive to $ε_{\rm ff}$. We find that $f_i$ increases with $ε_{\rm ff}$ and cluster mass. Through the dependence on $f_i$, the shape of the cluster initial mass function varies strongly with $ε_{\rm ff}$. The fraction of clustered star formation and maximum cluster mass increase with the star formation rate surface density, with the normalization of both relations dependent on $ε_{\rm ff}$. The cluster formation timescale systematically decreases with increasing $ε_{\rm ff}$. Local variations in the gas accretion history lead to a 0.25~dex scatter for the integral cluster formation efficiency. Joint constraints from all the observables prefer the runs that produce a median integral efficiency of 16\%. △ Less

Submitted 25 June, 2018; v1 submitted 4 December, 2017; originally announced December 2017.

Comments: 18 pages, 14 figure., ApJ in press

Disrupted Globular Clusters and the Gamma-Ray Excess in the Galactic Centre

Authors: Giacomo Fragione, Fabio Antonini, Oleg Y. Gnedin

Abstract: The Fermi Large Area Telescope has provided the most detailed view toward the Galactic Centre (GC) in high-energy gamma rays. Besides the interstellar emission and point-source contributions, the data suggest a residual diffuse gamma-ray excess. The similarity of its spatial distribution with the expected profile of dark matter has led to claims that this may be evidence for dark matter particle a… ▽ More The Fermi Large Area Telescope has provided the most detailed view toward the Galactic Centre (GC) in high-energy gamma rays. Besides the interstellar emission and point-source contributions, the data suggest a residual diffuse gamma-ray excess. The similarity of its spatial distribution with the expected profile of dark matter has led to claims that this may be evidence for dark matter particle annihilation. Here, we investigate an alternative explanation that the signal originates from millisecond pulsars (MSPs) formed in dense globular clusters and deposited at the GC as a consequence of cluster inspiral and tidal disruption. We use a semi-analytical model to calculate the formation, migration, and disruption of globular clusters in the Galaxy. Our model reproduces the mass of the nuclear star cluster and the present-day radial and mass distribution of globular clusters. For the first time, we calculate the evolution of MSPs from disrupted globular clusters throughout the age of the Galaxy and consistently include the effect of the MSP spin-down due to magnetic-dipole breaking. The final gamma-ray amplitude and spatial distribution are in good agreement with the Fermi observations and provide a natural astrophysical explanation for the GC excess. △ Less

Submitted 22 January, 2018; v1 submitted 11 September, 2017; originally announced September 2017.

Comments: 10 pages, 7 figures, 2 tables, submitted to MNRAS

Theia: Faint objects in motion or the new astrometry frontier

Authors: The Theia Collaboration, Celine Boehm, Alberto Krone-Martins, Antonio Amorim, Guillem Anglada-Escude, Alexis Brandeker, Frederic Courbin, Torsten Ensslin, Antonio Falcao, Katherine Freese, Berry Holl, Lucas Labadie, Alain Leger, Fabien Malbet, Gary Mamon, Barbara McArthur, Alcione Mora, Michael Shao, Alessandro Sozzetti, Douglas Spolyar, Eva Villaver, Conrado Albertus, Stefano Bertone, Herve Bouy, Michael Boylan-Kolchin , et al. (74 additional authors not shown)

Abstract: In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15… ▽ More In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision. △ Less

Submitted 2 July, 2017; originally announced July 2017.

Comments: http://theia.phyip3.dur.ac.uk

AGN Activity in Nucleated Galaxies as Measured by Chandra

Authors: Adi Foord, Elena Gallo, Edmund Hodges-Kluck, Brendan P. Miller, Vivienne F. Baldassare, Kayhan Gultekin, Oleg Y. Gnedin

Abstract: Motivated by theoretical expectations that Nuclear Star Clusters (NSCs) in galactic centers may provide a favorable environment for super-massive black holes to form and/or efficiently grow, we set out to measure the fraction of nearby nucleated galaxies that also host an Active Galactic Nucleus (AGN). We targeted a distance-limited sample of 98 objects with the Chandra X-ray Telescope, down to a… ▽ More Motivated by theoretical expectations that Nuclear Star Clusters (NSCs) in galactic centers may provide a favorable environment for super-massive black holes to form and/or efficiently grow, we set out to measure the fraction of nearby nucleated galaxies that also host an Active Galactic Nucleus (AGN). We targeted a distance-limited sample of 98 objects with the Chandra X-ray Telescope, down to a uniform X-ray luminosity threshold of $\sim$10$^{38}$ erg s$^{-1}$. The sample is composed of 47 late-types and 51 early-types, enabling us to further investigate the active fraction as a function of galactic morphology. After correcting for contamination to the nuclear X-ray signal from bright X-ray binaries, we measure an active fraction $f$=11.2$\%^{+7.4}_{-4.9}$ (1$σ$ C.L.) across the whole sample, in agreement with previous estimates based on an heterogeneous combination of optical, X-ray and radio diagnostics, by Seth et al. (2008). After accounting for the different stellar mass distributions in our samples, we find no statistically significant difference in the active fraction of early- vs. late-type nucleated galaxies, with $f$=10.6$\%^{+11.9}_{-4.9}$ and 10.8$\%^{+11.3}_{-6.3}$, respectively. For the early-type nucleated galaxies, we are able to carry out a controlled comparison with a parent sample of non-nucleated galaxies covering the same stellar mass range, finding again no statistically significant difference in the active fraction. Taken at face value, our findings suggest that the presence of a NSC does not facilitate nor enhance accretion-powered emission from a nuclear super-massive black hole. This is true even for late-type nucleated galaxies, home to bluer NSCs and arguably larger gas reservoirs. △ Less

Submitted 12 April, 2017; originally announced April 2017.

Comments: ApJ accepted, 12 pages, 8 Figures, 4 Tables

Globular Cluster Systems in Brightest Cluster Galaxies. III: Beyond Bimodality

Authors: W. E. Harris, S. M. Ciccone, G. M. Eadie, O. Y. Gnedin, D. Geisler, B. Rothberg, J. Bailin

Abstract: We present new deep photometry of the rich globular cluster (GC) systems around the Brightest Cluster Galaxies UGC 9799 (Abell 2052) and UGC 10143 (Abell 2147), obtained with the HST ACS and WFC3 cameras. For comparison, we also present new reductions of similar HST/ACS data for the Coma supergiants NGC 4874 and 4889. All four of these galaxies have huge cluster populations (to the radial limits o… ▽ More We present new deep photometry of the rich globular cluster (GC) systems around the Brightest Cluster Galaxies UGC 9799 (Abell 2052) and UGC 10143 (Abell 2147), obtained with the HST ACS and WFC3 cameras. For comparison, we also present new reductions of similar HST/ACS data for the Coma supergiants NGC 4874 and 4889. All four of these galaxies have huge cluster populations (to the radial limits of our data, comprising from 12000 to 23000 clusters per galaxy). The metallicity distribution functions (MDFs) of the GCs can still be matched by a bimodal-Gaussian form where the metal-rich and metal-poor modes are separated by ~0.8 dex, but the internal dispersions of each mode are so large that the total MDF becomes very broad and nearly continuous from [Fe/H] = -2.4 to Solar. There are, however, significant differences between galaxies in the relative numbers of \emph{metal-rich} clusters, suggesting that they underwent significantly different histories of mergers with massive, gas-rich halos. Lastly, the proportion of metal-poor GCs rises especially rapidly outside projected radii R > 4 R_eff, suggesting the importance of accreted dwarf satellites in the outer halo. Comprehensive models for the formation of GCs as part of the hierarchical formation of their parent galaxies will be needed to trace the systematic change in structure of the MDF with galaxy mass, from the distinctly bimodal form in smaller galaxies up to the broad continuum that we see in the very largest systems. △ Less

Submitted 21 December, 2016; originally announced December 2016.

Comments: In press for Astrophysical Journal

Star cluster formation in cosmological simulations. I. Properties of young clusters

Authors: Hui Li, Oleg Y. Gnedin, Nickolay Y. Gnedin, Xi Meng, Vadim A. Semenov, Andrey V. Kravtsov

Abstract: We present a new implementation of star formation in cosmological simulations, by considering star clusters as a unit of star formation. Cluster particles grow in mass over several million years at the rate determined by local gas properties, with high time resolution. The particle growth is terminated by its own energy and momentum feedback on the interstellar medium. We test this implementation… ▽ More We present a new implementation of star formation in cosmological simulations, by considering star clusters as a unit of star formation. Cluster particles grow in mass over several million years at the rate determined by local gas properties, with high time resolution. The particle growth is terminated by its own energy and momentum feedback on the interstellar medium. We test this implementation for Milky Way-sized galaxies at high redshift, by comparing the properties of model clusters with observations of young star clusters. We find that the cluster initial mass function is best described by a Schechter function rather than a single power law. In agreement with observations, at low masses the logarithmic slope is $α\approx 1.8-2$, while the cutoff at high mass scales with the star formation rate. A related trend is a positive correlation between the surface density of star formation rate and fraction of stars contained in massive clusters. Both trends indicate that the formation of massive star clusters is preferred during bursts of star formation. These bursts are often associated with major merger events. We also find that the median timescale for cluster formation ranges from 0.5 to 4 Myr and decreases systematically with increasing star formation efficiency. Local variations in the gas density and cluster accretion rate naturally lead to the scatter of the overall formation efficiency by an order of magnitude, even when the instantaneous efficiency is kept constant. Comparison of the formation timescale with the observed age spread of young star clusters provides an additional important constraint on the modeling of star formation and feedback schemes. △ Less

Submitted 7 November, 2016; v1 submitted 10 August, 2016; originally announced August 2016.

Comments: 17 pages, 13 figures, ApJ accepted

Globular Cluster Systems in Brightest Cluster Galaxies. II: NGC 6166

Authors: William E. Harris, John P. Blakeslee, Bradley C. Whitmore, Oleg Y. Gnedin, Douglas Geisler, Barry Rothberg

Abstract: We present new deep photometry of the globular cluster system (GCS) around NGC 6166, the central supergiant galaxy in Abell 2199. HST data from the ACS and WFC3 cameras in F475W, F814W are used to determine the spatial distribution of the GCS, its metallicity distribution function (MDF), and the dependence of the MDF on galactocentric radius and on GC luminosity. The MDF is extremely broad, with t… ▽ More We present new deep photometry of the globular cluster system (GCS) around NGC 6166, the central supergiant galaxy in Abell 2199. HST data from the ACS and WFC3 cameras in F475W, F814W are used to determine the spatial distribution of the GCS, its metallicity distribution function (MDF), and the dependence of the MDF on galactocentric radius and on GC luminosity. The MDF is extremely broad, with the classic red and blue subpopulations heavily overlapped, but a double-Gaussian model can still formally match the MDF closely. The spatial distribution follows a Sersic-like profile detectably to a projected radius of at least $R_{gc} = 250$ kpc. To that radius, the total number of clusters in the system is N_{GC} = 39000 +- 2000, the global specific frequency is S_N = 11.2 +- 0.6, and 57\% of the total are blue, metal-poor clusters. The GCS may fade smoothly into the Intra-Cluster Medium of A2199; we see no clear transition from the core of the galaxy to the cD halo or the ICM. The radial distribution, projected ellipticity, and mean metallicity of the red (metal-richer) clusters match the halo light extremely well for R > 15 kpc, both of them varying as σ_{MRGC} ~ σ_{light} ~ R^-1.8. By comparison, the blue (metal-poor) GC component has a much shallower falloff σ_{MPGC} ~ R^-1.0 and a more nearly spherical distribution. This strong difference in their density distributions produces a net metallicity gradient in the GCS as a whole that is primarily generated by the population gradient. With NGC 6166 we appear to be penetrating into a regime of high enough galaxy mass and rich enough environment that the bimodal two-phase description of GC formation is no longer as clear or effective as it has been in smaller galaxies. △ Less

Submitted 26 November, 2015; originally announced November 2015.

Comments: 19 pages, 22 figures. In press for Astrophysical Journal

Proper Motions and Trajectories for 16 Extreme Runaway and Hypervelocity Stars

Authors: Warren R. Brown, Jay Anderson, Oleg Y. Gnedin, Howard E. Bond, Margaret J. Geller, Scott J. Kenyon

Abstract: We measure proper motions with the Hubble Space Telescope for 16 extreme radial velocity stars, mostly unbound B stars in the Milky Way halo. Twelve of these stars have proper motions statistically consistent with zero, and thus have radial trajectories statistically consistent with a Galactic center "hypervelocity star" origin. The trajectory of HE 0437-5439 is consistent with both Milky Way and… ▽ More We measure proper motions with the Hubble Space Telescope for 16 extreme radial velocity stars, mostly unbound B stars in the Milky Way halo. Twelve of these stars have proper motions statistically consistent with zero, and thus have radial trajectories statistically consistent with a Galactic center "hypervelocity star" origin. The trajectory of HE 0437-5439 is consistent with both Milky Way and Large Magellanic Cloud origins. A Galactic center origin is excluded at 3-sigma confidence for two of the lowest radial velocity stars in our sample, however. These two stars are probable disk runaways and provide evidence for ~500 km/s ejections from the disk. We also measure a significant proper motion for the unbound sdO star US 708. Its 1,000 km/s motion is in some tension with proposed supernova ejection models, but can be explained if US 708 was ejected from the stellar halo. In the future, we expect Gaia will better constrain the origin of these remarkable unbound stars. △ Less

Submitted 17 February, 2015; originally announced February 2015.

Comments: 11 pages, 5 figures, submitted to ApJ

Globular Cluster Systems in Brightest Cluster Galaxies: A Near-Universal Luminosity Function?

Authors: William E. Harris, Warren Morningstar, Oleg Y. Gnedin, Heather O'Halloran, John P. Blakeslee, Bradley C. Whitmore, Patrick Cote, Douglas Geisler, Eric W. Peng, Jeremy Bailin, Barry Rothberg, Robert Cockcroft, Regina Barber DeGraaff

Abstract: We present the first results from our HST Brightest Cluster Galaxy (BCG) survey of seven central supergiant cluster galaxies and their globular cluster (GC) systems. We measure a total of 48000 GCs in all seven galaxies, representing the largest single GC database. We find that a log-normal shape accurately matches the observed luminosity function (LF) of the GCs down to the GCLF turnover point, w… ▽ More We present the first results from our HST Brightest Cluster Galaxy (BCG) survey of seven central supergiant cluster galaxies and their globular cluster (GC) systems. We measure a total of 48000 GCs in all seven galaxies, representing the largest single GC database. We find that a log-normal shape accurately matches the observed luminosity function (LF) of the GCs down to the GCLF turnover point, which is near our photometric limit. In addition, the LF has a virtually identical shape in all seven galaxies. Our data underscore the similarity in the formation mechanism of massive star clusters in diverse galactic environments. At the highest luminosities (log L > 10^7 L_Sun) we find small numbers of "superluminous" objects in five of the galaxies; their luminosity and color ranges are at least partly consistent with those of UCDs (Ultra-Compact Dwarfs). Lastly, we find preliminary evidence that in the outer halo (R > 20 kpc), the LF turnover point shows a weak dependence on projected distance, scaling as L_0 ~ R^-0.2, while the LF dispersion remains nearly constant. △ Less

Submitted 23 October, 2014; originally announced October 2014.

Comments: To appear in Astrophysical Journal, December 2014

Modeling the Formation of Globular Cluster Systems in the Virgo Cluster

Authors: Hui Li, Oleg Y. Gnedin

Abstract: The mass distribution and chemical composition of globular cluster (GC) systems preserve fossil record of the early stages of galaxy formation. The observed distribution of GC colors within massive early-type galaxies in the ACS Virgo Cluster Survey (ACSVCS) reveals a multi-modal shape, which likely corresponds to a multi-modal metallicity distribution. We present a simple model for the formation… ▽ More The mass distribution and chemical composition of globular cluster (GC) systems preserve fossil record of the early stages of galaxy formation. The observed distribution of GC colors within massive early-type galaxies in the ACS Virgo Cluster Survey (ACSVCS) reveals a multi-modal shape, which likely corresponds to a multi-modal metallicity distribution. We present a simple model for the formation and disruption of GCs that aims to match the ACSVCS data. This model tests the hypothesis that GCs are formed during major mergers of gas-rich galaxies and inherit the metallicity of their hosts. To trace merger events, we use halo merger trees extracted from a large cosmological N-body simulation. We select 20 halos in the mass range of $2\times 10^{12}$ to $7\times 10^{13}M_\odot$ and match them to 19 Virgo galaxies with K-band luminosity between $3\times 10^{10}$ and $3\times 10^{11}L_\odot$. To set the [Fe/H] abundances, we use an empirical galaxy mass-metallicity relation. We find that a minimal merger ratio of 1:3 best matches the observed cluster metallicity distribution. A characteristic bimodal shape appears because metal-rich GCs are produced by late mergers between massive halos, while metal-poor GCs are produced by collective merger activities of less massive hosts at early times. The model outcome is robust to alternative prescriptions for cluster formation rate throughout cosmic time, but a gradual evolution of the mass-metallicity relation with redshift appears to be necessary to match the observed cluster metallicities. We also affirm the age-metallicity relation, predicted by an earlier model, in which metal-rich clusters are systematically several billion years younger than their metal-poor counterparts. △ Less

Submitted 3 November, 2014; v1 submitted 4 May, 2014; originally announced May 2014.

Comments: 14 pages, 15 figures, Accepted for publication in the Astrophysical Journal

Co-Evolution of Galactic Nuclei and Globular Cluster Systems

Authors: Oleg Y. Gnedin, Jeremiah P. Ostriker, Scott Tremaine

Abstract: We revisit the hypothesis that dense galactic nuclei are formed from inspiraling globular clusters. Recent advances in understanding of the continuous formation of globular clusters over cosmic time and the concurrent evolution of the galaxy stellar distribution allow us to construct a simple model that matches the observed spatial and mass distributions of clusters in the Galaxy and the giant ell… ▽ More We revisit the hypothesis that dense galactic nuclei are formed from inspiraling globular clusters. Recent advances in understanding of the continuous formation of globular clusters over cosmic time and the concurrent evolution of the galaxy stellar distribution allow us to construct a simple model that matches the observed spatial and mass distributions of clusters in the Galaxy and the giant elliptical galaxy M87. In order to compare with observations, we model the effects of dynamical friction and dynamical evolution, including stellar mass loss, tidal stripping of stars, and tidal disruption of clusters by the growing galactic nucleus. We find that inspiraling globular clusters form a dense central structure, with mass and radius comparable to the typical values in observed nuclear star clusters (NSCs) in late-type and low-mass early-type galaxies. The density contrast associated with the NSC is less pronounced in giant elliptical galaxies. Our results indicate that the NSC mass as a fraction of mass of the galaxy stellar spheroid scales as M_NSC/M_* = 0.0025 M_{*,11}^{-0.5}. Thus disrupted globular clusters could contribute most of the mass of NSCs in galaxies with stellar mass below 10^{11} Msun. The inner part of the accumulated cluster may seed the growth of a central black hole via stellar dynamical core collapse, thereby relieving the problem of how to form luminous quasars at high redshift. The seed black hole may reach about 10^5 Msun within 1 Gyr of the beginning of globular cluster formation. △ Less

Submitted 13 January, 2014; v1 submitted 31 July, 2013; originally announced August 2013.

Comments: 15 pages, submitted to ApJ, revised version after referee comments

Revisiting The First Galaxies: The epoch of Population III stars

Authors: Alexander L. Muratov, Oleg Y. Gnedin, Nickolay Y. Gnedin, Marcel Zemp

Abstract: We investigate the transition from primordial Population III (Pop III) star formation to normal Pop II star formation in the first galaxies using new cosmological hydrodynamic simulations. We find that while the first stars seed their host galaxies with metals, they cannot sustain significant outflows to enrich the intergalactic medium, even assuming a top-heavy initial mass function. This means t… ▽ More We investigate the transition from primordial Population III (Pop III) star formation to normal Pop II star formation in the first galaxies using new cosmological hydrodynamic simulations. We find that while the first stars seed their host galaxies with metals, they cannot sustain significant outflows to enrich the intergalactic medium, even assuming a top-heavy initial mass function. This means that Pop III star formation could potentially continue until z~6 in different unenriched regions of the universe, before being ultimately shut off by cosmic reionization. Within an individual galaxy, the metal production and stellar feedback from Pop II stars overtake Pop III stars in 20-200 Myr, depending on galaxy mass. △ Less

Submitted 9 August, 2013; v1 submitted 5 December, 2012; originally announced December 2012.

Comments: 9 pages, 7 figures, published in ApJ

Journal ref: Muratov et al. 2013, ApJ, 773, 19

Revisiting The First Galaxies: The effects of Population III stars on their host galaxies

Authors: Alexander L. Muratov, Oleg Y. Gnedin, Nickolay Y. Gnedin, Marcel Zemp

Abstract: We revisit the formation and evolution of the first galaxies using new hydrodynamic cosmological simulations with the ART code. Our simulations feature a recently developed model for H2 formation and dissociation, and a star formation recipe that is based on molecular rather than atomic gas. Here, we develop and implement a recipe for the formation of metal-free Population III stars in galaxy-scal… ▽ More We revisit the formation and evolution of the first galaxies using new hydrodynamic cosmological simulations with the ART code. Our simulations feature a recently developed model for H2 formation and dissociation, and a star formation recipe that is based on molecular rather than atomic gas. Here, we develop and implement a recipe for the formation of metal-free Population III stars in galaxy-scale simulations that resolve primordial clouds with sufficiently high density. We base our recipe on the results of prior zoom-in simulations that resolved the protostellar collapse in pre-galactic objects. We find the epoch during which Pop III stars dominated the energy and metal budget of the first galaxies to be short-lived. Galaxies which host Pop III stars do not retain dynamical signatures of their thermal and radiative feedback for more than 10^8 yr after the lives of the stars end in pair-instability supernovae, even when we consider the maximum reasonable efficiency of the feedback. Though metals ejected by the supernovae can travel well beyond the virial radius of the host galaxy, they typically begin to fall back quickly, and do not enrich a large fraction of the intergalactic medium. Galaxies with total mass in excess of 3 x 10^6 Msun re-accrete most of their baryons and transition to metal-enriched Pop II star formation. △ Less

Submitted 9 August, 2013; v1 submitted 4 December, 2012; originally announced December 2012.

Comments: 13 pages, 9 figures, published in ApJ

Journal ref: Muratov et al. 2013, ApJ, 772, 106