Search

The stellar masses of galaxies are measured using integrated light via several methods -- however, few of these methods were designed for low-mass ($M_{\star}\lesssim10^{8}\rm{M_{\odot}}$) "dwarf" galaxies, whose properties (e.g., stochastic star formation, low metallicity) pose unique challenges for estimating stellar masses. In this work, we quantify the precision and accuracy at which stellar masses of low-mass galaxies can be recovered using UV/optical/IR photometry. We use mock observations of 469 low-mass galaxies from a variety of models, including both semi-empirical models (GRUMPY, UniverseMachine-SAGA) and cosmological baryonic zoom-in simulations (MARVELous Dwarfs and FIRE-2), to test literature color-$M_\star/L$ relations and multi-wavelength spectral energy distribution (SED) mass estimators. We identify a list of "best practices" for measuring stellar masses of low-mass galaxies from integrated photometry. These include updated prescriptions for stellar mass based on $g-r$ color and WISE 3.4 $\mu$m luminosity, which are less systematically biased than literature calibrations and can recover true stellar masses of low-mass galaxies with $\sim0.1$ dex precision. When using SED fitting to estimate stellar mass, we find that the form of the assumed star formation history can induce significant biases: parametric SFHs can underestimate stellar mass by as much as $\sim0.4$ dex, while non-parametric SFHs recover true stellar masses with insignificant offset ($-0.03\pm0.11$ dex). However, we also caution that non-informative dust attenuation priors may introduce $M_\star$ uncertainties of up to $\sim0.6$ dex., Comment: 28 pages including references, 9 figures; submitted to ApJ

We compute parametric measurements of the Einstein-radius-enclosed total mass for 177 cluster-scale strong gravitational lenses identified by the ChicagO Optically-selected Lenses Located At the Margins of Public Surveys (COOL-LAMPS) collaboration with lens redshifts ranging from $0.2 \lessapprox z \lessapprox 1.0$ using only two measured parameters in each lensing system: the Einstein radius, and the brightest-cluster-galaxy (BCG) redshift. We then constrain the Einstein-radius-enclosed luminosity and stellar mass by fitting parametric spectral energy distributions (SEDs) with aperture photometry from the Dark Energy Camera Legacy Survey (DECaLS) in the $g$, $r$, and $z$-band Dark Energy Camera (DECam) filters. We find that the BCG redshift, enclosed total mass, and enclosed luminosity are strongly correlated and well described by a planar relationship in 3D space. We also find that the enclosed total mass and stellar mass are correlated with a logarithmic slope of $0.443\pm0.035$, and the enclosed total mass and stellar-to-total mass fraction are correlated with a logarithmic slope of $-0.563\pm0.035$. The correlations described here can be used to validate strong lensing candidates in upcoming imaging surveys -- such as Rubin/Legacy Survey of Space and Time (LSST) -- in which an algorithmic treatment of lensing systems will be needed due to the sheer volume of data these surveys will produce., Comment: 17 pages, 5 figures, 2 tables. Submitted to The Astrophysical Journal. v3: updated authors, formatting, grammar, and references

The flux-based statistical theory of the non-hierarchical three-body system predicts that the chaotic outcome distribution reduces to the chaotic emissivity function times a known function, the asymptotic flux. Here, we measure the chaotic emissivity function (or equivalently, the absorptivity) through simulations. More precisely, we follow millions of scattering events only up to the point when it can be decided whether the scattering is regular or chaotic. In this way, we measure a tri-variate absorptivity function. Using it, we determine the flux-based prediction for the chaotic outcome distribution over both binary binding energy and angular momentum, and we find good agreement with the measured distribution. This constitutes a detailed confirmation of the flux-based theory, and demonstrates a considerable reduction in computation to determine the chaotic outcome distribution., Comment: 10 pages, 8 figures

We report the discovery of five bright strong gravitationally lensed galaxies at $3 < z < 4$: COOLJ0101$+$2055 ($z = 3.459$), COOLJ0104$-$0757 ($z = 3.480$), COOLJ0145$+$1018 ($z = 3.310$), COOLJ0516$-$2208 ($z = 3.549$), and COOLJ1356$+$0339 ($z = 3.753$). These galaxies have magnitudes of $r_{\rm AB}, z_{\rm AB} < 21.81$ mag and are lensed by galaxy clusters at $0.26 < z < 1$. This sample nearly doubles the number of known bright lensed galaxies with extended arcs at $3 < z < 4$. We characterize the lensed galaxies using ground-based grz/giy imaging and optical spectroscopy. We report model-based magnitudes and derive stellar masses, dust content, and star-formation rates via stellar population synthesis modeling. Building lens models based on ground-based imaging, we estimate source magnifications in the range $\sim$29 to $\sim$180. Combining these analyses, we derive demagnified stellar masses in the range $\rm log_{10}(M_{*}/M_{\odot}) \sim 9.69 - 10.75$ and star formation rates in the youngest age bin ranging from $\rm log_{10}(SFR/(M_{\odot}\cdot yr^{-1})) \sim 0.39 - 1.46$, placing the sample galaxies on the massive end of the star-forming main sequence in this redshift interval. In addition, three of the five galaxies have strong Ly$\alpha$ emissions, offering unique opportunities to study Ly$\alpha$ emitters at high redshift in future work., Comment: 20 pages, 10 figures

We examine the spatial distribution and orbital pole correlations of satellites in a suite of zoom-in high-resolution dissipationless simulations of Miky Way (MW) sized haloes. We use the measured distribution to estimate the incidence of satellite configurations as flattened and as correlated in their orbital pole distribution as satellite system of the Milky Way. We confirm that this incidence is sensitive to the radial distribution of subhaloes and thereby to the processes that affect it, such as artificial disruption due to numerical effects and disruption due to the central disk. Controlling for the resolution effects and bracketing the effects of the disk, we find that the MW satellite system is somewhat unusual (at the $\approx 2-3\sigma$ level) but is statistically consistent with the $\Lambda$CDM model, in general agreement with results and conclusions of other recent studies., Comment: 10 pages, 7 figures, submitted to MNRAS

We use the Auriga simulations to probe different satellite quenching mechanisms operating at different mass scales ($10^5 M_\odot \lesssim M_\star \lesssim 10^{11} M_\odot$) in Milky Way-like hosts. Our goal is to understand the origin of the satellite colour distribution and star-forming properties in both observations and simulations. We find that the satellite populations in the Auriga simulations, which was originally designed to model Milky Way-like host galaxies, resemble the populations in the Exploration of Local VolumE Satellites (ELVES) Survey and the Satellites Around Galactic Analogs (SAGA) survey in their luminosity function in the luminosity range $-12 \lesssim M_V \lesssim -15$ and resemble ELVES in their quenched fraction and colour--magnitude distribution in the luminosity range $-12 \lesssim M_g \lesssim -15$. We find that satellites transition from blue colours to red colours at the luminosity range $-15 \lesssim M_g \lesssim -12$ in both the simulations and observations and we show that this shift is driven by environmental effects in the simulations. We demonstrate also that the colour distribution in both simulations and observations can be decomposed into two statistically distinct populations based on their morphological type or star-forming status that are statistically distinct. In the simulations, these two populations also have statistically distinct infall time distributions. The comparison presented here seems to indicate that the tension between the quenched fraction in SAGA and simulations is resolved by the improved target selection of ELVES, but there are still tensions in understanding the colours of faint galaxies, of which ELVES appears to have a significant population of faint blue satellites not recovered in Auriga., Comment: 15 pages, 10 figures. Published in MNRAS on December 15, 2022

We use \texttt{GRUMPY}, a simple regulator-type model for dwarf galaxy formation and evolution, to forward model the dwarf galaxy satellite population of the Milky Way (MW) using the Caterpillar zoom-in simulation suite. We show that luminosity and distance distributions of the model satellites are consistent with the distributions measured in the DES, PS1 and SDSS surveys, even without including a model for the orphan galaxies. We also show that our model for dwarf galaxy sizes can simultaneously reproduce the observed {\it distribution} of stellar half-mass radii, $r_{1/2}$, of the MW satellites and the overall $r_{1/2}-M_\star$ relation exhibited by observed dwarf galaxies. The model predicts that some of the observed faint stellar systems with $r_{1/2}<10$ pc are ultra-faint dwarf galaxies. Scaling of the stellar mass $M_\star$ and peak halo mass $M_{\rm peak}$ for the model satellites is not described by a power law, but has a clear flattening of $M_\star-M_{\rm peak}$ scaling at $M_{\rm peak}<10^8\,M_\odot$ imprinted by reionization. As a result, the fraction of low mass haloes ($M_{\rm peak} < 10^8\, M_\odot$) hosting galaxies with $M_V<0$ is predicted to be 50% at $M_{\rm peak} \sim 3.6 \times 10^7\,M_\odot$. We find that such high fraction at that halo mass helps to reproduce the number of dwarf galaxies discovered recently in the HSC-SSP survey. Using the model we forecast that there should be the total of $440^{+201}_{-147}$ (68\% confidence interval) MW satellites with $M_V < 0$ and $r_{1/2} > 10$ pc within 300 kpc and make specific predictions for the HSC-SSP, DELVE-WIDE and LSST surveys., Comment: MNRAS in press, 28 pages, 22 figures. Code to reproduce figures is available at https://github.com/kibokov/grumpy_mw_forward_model

We present a simple regulator-type framework designed specifically for modelling formation of dwarf galaxies. Despite its simplicity, when coupled with realistic mass accretion histories of haloes from simulations and reasonable choices for model parameter values, the framework can reproduce a remarkably broad range of observed properties of dwarf galaxies over seven orders of magnitude in stellar mass. In particular, we show that the model can simultaneously match observational constraints on the stellar mass--halo mass relation, as well as observed relations between stellar mass and gas phase and stellar metallicities, gas mass, size, and star formation rate, as well as general form and diversity of star formation histories (SFHs) of observed dwarf galaxies. The model can thus be used to predict photometric properties of dwarf galaxies hosted by dark matter haloes in $N$-body simulations, such as colors, surface brightnesses, and mass-to-light ratios and to forward model observations of dwarf galaxies. We present examples of such modelling and show that colors and surface brightness distributions of model galaxies are in good agreement with observed distributions for dwarfs in recent observational surveys. We also show that in contrast with the common assumption, the absolute magnitude-halo mass relation is generally predicted to have a non-power law form in the dwarf regime, and that the fraction of haloes that host detectable ultrafaint galaxies is sensitive to reionization redshift ($z_{\rm rei}$) and is predicted to be consistent with observations for $z_{\rm rei}\lesssim 9$., Comment: 26 pages, 20 figures. The code implementing GRUMPY galaxy formation model is available at https://github.com/kibokov/GRUMPY

We present an extensive comparison between the statistical properties of non-hierarchical three-body systems and the corresponding three-body theoretical predictions. We perform and analyze 1 million realizations for each different initial condition considering equal and unequal mass three-body systems to provide high accuracy statistics. We measure 4 quantities characterizing the statistical distribution of ergodic disintegrations: escape probability of each body, the characteristic exponent for escapes by a narrow margin, predicted absorptivity as a function of binary energy and binary angular momentum, and, finally, the lifetime distribution. The escape probabilities are shown to be in agreement down to the 1% level with the emissivity-blind, flux-based theoretical prediction. This represents a leap in accuracy compared to previous three-body statistical theories. The characteristic exponent at the threshold for marginally unbound escapes is an emissivity-independent flux-based prediction, and the measured values are found to agree well with the prediction. We interpret both tests as strong evidence for the flux-based three-body statistical formalism. The predicted absorptivity and lifetime distributions are measured to enable future tests of statistical theories., Comment: 17 pages, 11 figures

We study chaos and L\'evy flights in the general gravitational three-body problem. We introduce new metrics to characterize the time evolution and final lifetime distributions, namely Scramble Density $\mathcal{S}$ and the LF index $\mathcal{L}$, that are derived from the Agekyan-Anosova maps and homology radius $R_{\mathcal{H}}$. Based on these metrics, we develop detailed procedures to isolate the ergodic interactions and L\'evy flight interactions. This enables us to study the three-body lifetime distribution in more detail by decomposing it into the individual distributions from the different kinds of interactions. We observe that ergodic interactions follow an exponential decay distribution similar to that of radioactive decay. Meanwhile, L\'evy flight interactions follow a power-law distribution. L\'evy flights in fact dominate the tail of the general three-body lifetime distribution, providing conclusive evidence for the speculated connection between power-law tails and L\'evy flight interactions. We propose a new physically-motivated model for the lifetime distribution of three-body systems and discuss how it can be used to extract information about the underlying ergodic and L\'evy flight interactions. We discuss mass ejection probabilities in three-body systems in the ergodic limit and compare it to previous ergodic formalisms. We introduce a novel mechanism for a three-body relaxation process and discuss its relevance in general three-body systems., Comment: 20 pages, 14 figures