WIYN Open Cluster Study XVI: Optical/Infrared Photometry and Comparisons With Theoretical Isochrones

We present combined optical/near-IR photometry (BVIK) for six open clusters - M35, M37, NGC 1817, NGC 2477, NGC 2420, and M67. The open clusters span an age range from 150 Myr to 4 Gyr and have metal abundances from [Fe/H] = -0.3 to +0.09 dex. We have utilized these data to test the robustness of theoretical main sequences constructed by several groups as denoted by the following designations - Padova, Baraffe, Y^2, Geneva, and Siess. The comparisons of the models with the observations have been performed in the [Mv, (B-V)o], [Mv, (V-I)o], and [Mv, (V-K)o] colour-magnitude diagrams as well as the distance-independent [(V-K)o, (B-V)o] and [(V-K)o, (V-I)o] two-colour diagrams. We conclude that none of the theoretical models reproduce the observational data in a consistent manner over the magnitude and colour range of the unevolved main sequence. In particular, there are significant zeropoint and shape differences between the models and the observations. We speculate that the crux of the problem lies in the precise mismatch between theoretical and observational colour-temperature relations. These results underscore the importance of pursuing the study of stellar structure and stellar modelling with even greater intensity.Comment: Accepted for publication in MNRAS. 12 pages, 37 figures, 4 tables. High resolution figures available from http://www.astro.ufl.edu/~aaron/opt_ir_figs

Comparison of Alpha-Element Enhanced Simple Stellar Population Models with Milky Way Globular Clusters

We present simple stellar population (SSP) models with scaled-solar and alpha-element enhanced abundances. The SSP models are based on the Dartmouth Stellar Evolution Database, our library of synthetic stellar spectra, and a detailed systematic variation of horizontal-branch (HB) morphology with age and metallicity. In order to test the relative importance of a variety of SSP model ingredients, we compare our SSP models with integrated spectra of 41 Milky Way Globular Clusters (MWGCs) from Schiavon et al. (2005). Using the Mg b and Ca4227 indices, we confirm that Mg and Ca are enhanced by about +0.4 and +0.2 dex, respectively, in agreement with results from high resolution spectra of individual stars in MWGCs. Balmer lines, particularly Hgamma and Hdelta, of MWGCs are reproduced by our alpha-enhanced SSP models not only because of the combination of isochrone and spectral effects but also because of our reasonable HB treatment. Moreover, it is shown that the Mg abundance significantly influences Balmer and iron line indices. Finally, the investigation of power-law initial mass function (IMF) variations suggests that an IMF much shallower than Salpeter is unrealistic because the Balmer lines are too strong on the metal-poor side to be compatible with observations.Comment: 38 pages, 18 figures, AJ accepted. Models are available from http://astro.wsu.edu/hclee/sp_LWD09.html and http://astro.wsu.edu/hclee/sp.htm

The DESI One-Percent Survey: Evidence for Assembly Bias from Low-Redshift Counts-in-Cylinders Measurements

We explore the galaxy-halo connection information that is available in low-redshift samples from the early data release of the Dark Energy Spectroscopic Instrument (DESI). We model the halo occupation distribution (HOD) from z=0.1-0.3 using Survey Validation 3 (SV3; a.k.a., the One-Percent Survey) data of the DESI Bright Galaxy Survey (BGS). In addition to more commonly used metrics, we incorporate counts-in-cylinders (CiC) measurements, which drastically tighten HOD constraints. Our analysis is aided by the Python package, galtab, which enables the rapid, precise prediction of CiC for any HOD model available in halotools. This methodology allows our Markov chains to converge with much fewer trial points, and enables even more drastic speedups due to its GPU portability. Our HOD fits constrain characteristic halo masses tightly and provide statistical evidence for assembly bias, especially at lower luminosity thresholds: the HOD of central galaxies in $z\sim0.15$ samples with limiting absolute magnitude $M_r < -20.0$ and $M_r < -20.5$ samples is positively correlated with halo concentration with a significance of 99.9% and 99.5%, respectively. Our models also favor positive central assembly bias for the brighter $M_r < -21.0$ sample at $z\sim0.25$ (94.8% significance), but there is no significant evidence for assembly bias with the same luminosity threshold at $z\sim0.15$. We provide our constraints for each threshold sample's characteristic halo masses, assembly bias, and other HOD parameters. These constraints are expected to be significantly tightened with future DESI data, which will span an area 100 times larger than that of SV3

The DESI One-percent Survey: Evidence for Assembly Bias from Low-redshift Counts-in-cylinders Measurements

We explore the galaxy-halo connection information that is available in low-redshift samples from the early data release of the Dark Energy Spectroscopic Instrument (DESI). We model the halo occupation distribution (HOD) from z = 0.1 to 0.3 using Survey Validation 3 (SV3; a.k.a., the One-Percent Survey) data of the DESI Bright Galaxy Survey. In addition to more commonly used metrics, we incorporate counts-in-cylinders (CiC) measurements, which drastically tighten HOD constraints. Our analysis is aided by the Python package, galtab, which enables the rapid, precise prediction of CiC for any HOD model available in halotools. This methodology allows our Markov chains to converge with much fewer trial points, and enables even more drastic speedups due to its GPU portability. Our HOD fits constrain characteristic halo masses tightly and provide statistical evidence for assembly bias, especially at lower luminosity thresholds: the HOD of central galaxies in z ∼ 0.15 samples with limiting absolute magnitude M r < −20.0 and M r < −20.5 samples is positively correlated with halo concentration with a significance of 99.9% and 99.5%, respectively. Our models also favor positive central assembly bias for the brighter M r < −21.0 sample at z ∼ 0.25 (94.8% significance), but there is no significant evidence for assembly bias with the same luminosity threshold at z ∼ 0.15. We provide our constraints for each threshold sample’s characteristic halo masses, assembly bias, and other HOD parameters. These constraints are expected to be significantly tightened with future DESI data, which will span an area 100 times larger than that of SV3

The DESI One-Percent survey: constructing galaxy-halo connections for ELGs and LRGs using auto and cross correlations

In the current Dark Energy Spectroscopic Instrument (DESI) survey, emission line galaxies (ELGs) and luminous red galaxies (LRGs) are essential for mapping the dark matter distribution at $z \sim 1$. We measure the auto and cross correlation functions of ELGs and LRGs at $0.8<z\leq 1.0$ from the DESI One-Percent survey. Following Gao et al. (2022), we construct the galaxy-halo connections for ELGs and LRGs simultaneously. With the stellar-halo mass relation (SHMR) for the whole galaxy population (i.e. normal galaxies), LRGs can be selected directly by stellar mass, while ELGs can also be selected randomly based on the observed number density of each stellar mass, once the probability $P_{\mathrm{sat}}$ of a satellite galaxy becoming an ELG is determined. We demonstrate that the observed small scale clustering prefers a halo mass-dependent $P_{\mathrm{sat}}$ model rather than a constant. With this model, we can well reproduce the auto correlations of LRGs and the cross correlations between LRGs and ELGs at $r_{\mathrm{p}}>0.1$ $\mathrm{Mpc}\,h^{-1}$. We can also reproduce the auto correlations of ELGs at $r_{\mathrm{p}}>0.3$ $\mathrm{Mpc}\,h^{-1}$ ($s>1$ $\mathrm{Mpc}\,h^{-1}$) in real (redshift) space. Although our model has only seven parameters, we show that it can be extended to higher redshifts and reproduces the observed auto correlations of ELGs in the whole range of $0.8<z<1.6$, which enables us to generate a lightcone ELG mock for DESI. With the above model, we further derive halo occupation distributions (HODs) for ELGs which can be used to produce ELG mocks in coarse simulations without resolving subhalos.Comment: 27 pages, 16 figures, accepted by Ap

Detecting and Characterizing Mg II absorption in DESI Survey Validation Quasar Spectra

In this paper we will present findings on the detection of Magnesium II (MgII, lambda = 2796 {\AA}, 2803 {\AA}) absorption systems observed in data from the Early Data Release (EDR) of the Dark Energy Spectroscopic Instrument (DESI). DESI is projected to obtain spectroscopy of approximately 3 million quasars (QSOs), of which over 99% are anticipated to be found at redshifts greater than z < 0.3, such that DESI would be able to observe an associated or intervening Mg II absorber illuminated by the background QSO. We have developed an autonomous supplementary spectral pipeline that detects such systems through an initial line-fitting process and then confirms line properties using a Markov Chain Monte Carlo (MCMC) sampler. Based upon both a visual inspection and the reanalysis of coadded observations, we estimate this sample of absorption systems to have a completeness of 82.56% and purity of 99.08%. As the spectra in which Mg II systems are detected are the result of coadding multiple observations, we can determine the sensitivity, and therefore completeness, of the sample by searching for known Mg II systems in coadded data with fewer observations (and therefore lower signal-to-noise). From a parent catalog containing 83,207 quasars, we detect a total of 23,921 Mg II absorption systems following a series of quality cuts. Extrapolating from this occurrence rate of 28.75% implies a catalog at the completion of the five-year DESI survey that contains over eight hundred thousand Mg II absorbers. The cataloging of these systems will enable significant further research as they carry information regarding circumgalactic medium (CGM) environments, the distribution of intervening galaxies, and the growth of metallicity across the redshift range 0.3 < z < 2.5.Comment: 12 pages, 7 figure

The DESI One-Percent Survey: Exploring the Halo Occupation Distribution of Luminous Red Galaxies and Quasi-Stellar Objects with AbacusSummit

We present the first comprehensive Halo Occupation Distribution (HOD) analysis of the DESI One-Percent survey Luminous Red Galaxy (LRG) and Quasi-Stellar Object (QSO) samples. We constrain the HOD of each sample and test possible HOD extensions by fitting the redshift-space galaxy 2-point correlation functions in 0.15 < r < 32 Mpc/h in a set of fiducial redshift bins. We use AbacusSummit cubic boxes at Planck 2018 cosmology as model templates and forward model galaxy clustering with the AbacusHOD package. We achieve good fits with a standard HOD model with velocity bias, and we find no evidence for galaxy assembly bias or satellite profile modulation at the current level of statistical uncertainty. For LRGs in 0.4 < z < 0.6, we infer a satellite fraction of fsat = 11+-1%, a mean halo mass of log10 Mh = 13.40+0.02-0.02, and a linear bias of blin = 1.93+0.06-0.04. For LRGs in 0.6 < z < 0.8, we find fsat = 14+-1%, log10 Mh = 13.24+0.02-0.02, and blin = 2.08+0.03-0.03. For QSOs, we infer fsat = 3+8-2%, log10 Mh = 12.65+0.09-0.04, and blin = 2.63+0.37-0.26 in redshift range 0.8 < z < 2.1. Using these fits, we generate a large suite of high-fidelity galaxy mocks. We also study the redshift-evolution of the DESI LRG sample from z = 0.4 up to z = 1.1, revealing significant and interesting trends in mean halo mass, linear bias, and satellite fraction.Comment: Submitted to MNRAS, comments welcom

Local primordial non-Gaussianity from the large-scale clustering of photometric DESI luminous red galaxies

We use angular clustering of luminous red galaxies from the Dark Energy Spectroscopic Instrument (DESI) imaging surveys to constrain the local primordial non-Gaussianity parameter fNL. Our sample comprises over 12 million targets, covering 14,000 square degrees of the sky, with redshifts in the range 0.2< z < 1.35. We identify Galactic extinction, survey depth, and astronomical seeing as the primary sources of systematic error, and employ linear regression and artificial neural networks to alleviate non-cosmological excess clustering on large scales. Our methods are tested against log-normal simulations with and without fNL and systematics, showing superior performance of the neural network treatment in reducing remaining systematics. Assuming the universality relation, we find fNL $= 47^{+14(+29)}_{-11(-22)}$ at 68\%(95\%) confidence. With a more aggressive treatment, including regression against the full set of imaging maps, our maximum likelihood value shifts slightly to fNL$\sim 50$ and the uncertainty on fNL increases due to the removal of large-scale clustering information. We apply a series of robustness tests (e.g., cuts on imaging, declination, or scales used) that show consistency in the obtained constraints. Despite extensive efforts to mitigate systematics, our measurements indicate fNL > 0 with a 99.9 percent confidence level. This outcome raises concerns as it could be attributed to unforeseen systematics, including calibration errors or uncertainties associated with low-\ell systematics in the extinction template. Alternatively, it could suggest a scale-dependent fNL model--causing significant non-Gaussianity around large-scale structure while leaving cosmic microwave background scales unaffected. Our results encourage further studies of fNL with DESI spectroscopic samples, where the inclusion of 3D clustering modes should help separate imaging systematics.Comment: 19 pages, 15 figures, 6 tables (Appendix excluded). Submitted to MNRA

Target Selection and Validation of DESI Emission Line Galaxies

The Dark Energy Spectroscopic Instrument (DESI) will precisely constrain cosmic expansion and the growth of structure by collecting ∼40 million extragalactic redshifts across ∼80% of cosmic history and one-third of the sky. The Emission Line galaxy (ELG) sample, which will comprise about one-third of all DESI tracers, will be used to probe the universe over the 0.6 &lt; z &lt; 1.6 range, including the 1.1 &lt; z &lt; 1.6 range, which is expected to provide the tightest constraints. We present the target selection for the DESI Survey Validation (SV) and Main Survey ELG samples, which relies on the imaging of the Legacy Surveys. The Main ELG selection consists of a g-band magnitude cut and a (g − r) versus (r − z) color box, while the SV selection explores extensions of the Main selection boundaries. The Main ELG sample is composed of two disjoint subsamples, which have target densities of about 1940 deg−2 and 460 deg−2, respectively. We first characterize their photometric properties and density variations across the footprint. We then analyze the DESI spectroscopic data that have been obtained from 2020 December to 2021 December in the SV and Main Survey. We establish a preliminary criterion for selecting reliable redshifts, based on the [O ii] flux measurement, and assess its performance. Using this criterion, we are able to present the spectroscopic efficiency of the Main ELG selection, along with its redshift distribution. We thus demonstrate that the Main selection 1940 deg−2 subsample alone should provide 400 deg−2 and 460 deg−2 reliable redshifts in the 0.6 &lt; z &lt; 1.1 and the 1.1 &lt; z &lt; 1.6 ranges, respectively