Researcher profile

Robin Ciardullo

Robin Ciardullo contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2026arXiv

ODIN: Clustering Properties of Ly$α$ Blobs at $z$ $\sim$ 2.4 and 3.1

Spatially extended Ly$α$ nebulae, known as Ly$α$ blobs (LABs), are a rare population at $z > 2$ that are thought to trace proto-groups or the progenitors of massive galaxies in the present-day universe. However, their dark matter halo properties (e.g., halo mass) are still uncertain due to their rarity and strong field-to-field variation. The One-hundred-deg$^2$ DECam Imaging in Narrowbands (ODIN) survey has discovered 103 and 112 LABs in the extended ($\sim$9~\sqdeg) COSMOS field at $z\sim2.4$ and 3.1, respectively, enabling estimation of their bias and host halo masses through clustering analysis. We measure the angular auto-correlation functions (ACFs) of LABs and derive galaxy bias factors of $b$ = $4.0\pm0.8$ and $3.8\pm0.7$, corresponding to minimum halo masses of $2.8^{+3.0}_{-1.8}$ and $0.7^{+0.8}_{-0.5}\times10^{12}~M_\odot$ and median halo masses of $4.2^{+3.8}_{-2.5}$ and $1.1^{+1.1}_{-0.7}\times10^{12}~M_\odot$ at $z\sim2.4$ and 3.1, respectively. LABs occupy $\sim$11$^{+39}_{-8}$\% and $\sim$3$^{+9}_{-2}$\% of all dark matter halos above these minimum halo masses. These findings suggest that LABs inhabit massive dark matter halos, likely tracing proto-group environments that evolve into present-day massive halos ($\sim$10$^{13}~M_\odot$), where massive elliptical galaxies or galaxy groups reside, by $z=0$.

preprint2025arXiv

Towards Precision Cosmology With Improved PNLF Distances Using VLT-MUSE. III. Impact of Stellar Populations in Early-Type Galaxy

Distance measurements using the planetary nebula luminosity function (PNLF) rely on the bright-end power-law cut-off magnitude ($M^*$), which is defined by a number of the [OIII]$\lambda5007$-brightest planetary nebulae (PNe). In early-type galaxies (ETGs), the formation of these PNe is enigmatic; the population is typically too old to form the expected $M^*$ PNe from single star evolution. We aim to give a solution to this problem. We selected five ETGs with known MUSE-PNLF distances. The MUSE instrument allows us to calculate the PNLF and consistently investigate the underlying stellar populations. Using stellar population synthesis, we derive the population age, star formation history, metallicity, and alpha abundance. We compare these parameters to the PNLF variables: $M^*$ and luminosity-specific PN number at the top 0.5 mag of the PNLF ($α_{0.5}$). We also compare our results with PNe In Cosmological Simulations (PICS) model applied to Magneticum Pathfinder analogue galaxies. The average mass-weighted ages and metallicities of our observations are typically old ($9 <\mathrm{Age}< 13.5$ Gyr) and rather metal-rich ($-0.4 <\mathrm{[M/H]}< +0.2$). We find $M^*$ to be independent of age and metallicity in these ages and metallicity intervals. We discover a positive correlation between $α_{0.5}$ values and the mass fraction of stellar population ages of 2--10 Gyr, implying that most of the PNe originate from stars with intermediate ages. Similar trends are also found in the PICS analogue galaxies. We show that the presence of at least $\sim 2\%$ of stellar mass younger than 10 Gyr is, in principle, sufficient to form the $M^*$ PNe in ETGs. We also present observing requirements for an ideal PNLF distance determination in ETGs.

preprint2024arXiv

Absorption Troughs of Lyman Alpha Emitters in HETDEX

The Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) is designed to detect and measure the redshifts of more than one million Ly$α$ emitting galaxies (LAEs) between $1.88 < z < 3.52$. In addition to its cosmological measurements, these data enable studies of Ly$α$ spectral profiles and the underlying radiative transfer. Using the roughly half a million LAEs in the HETDEX Data Release 3, we stack various subsets to obtain the typical Ly$α$ profile for the $z \sim 2-3$ epoch and to understand their physical properties. We find clear absorption wings around Ly$α$ emission, which extend $\sim 2000$ km $\mathrm{s}^{-1}$ both redward and blueward of the central line. Using far-UV spectra of nearby ($0.002 < z < 0.182$) LAEs in the CLASSY treasury and optical/near-IR spectra of $2.8 < z < 6.7$ LAEs in the MUSE-Wide survey, we observe absorption profiles in both redshift regimes. Dividing the sample by volume density shows that the troughs increase in higher density regions. This trend suggests that the depth of the absorption is dependent on the local density of objects near the LAE, a geometry that is similar to damped Lyman-$α$ systems. Simple simulations of Ly$α$ radiative transfer can produce similar troughs due to absorption of light from background sources by HI gas surrounding the LAEs.

preprint2023arXiv

The Spitzer-HETDEX Exploratory Large Area Survey. IV. Model-Based Multi-wavelength Photometric Catalog

We present a 0.3--4.5 $μ$m 16-band photometric catalog for the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers a $\sim 27$ deg$^2$ field within the footprint of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). Here we present new DECam imaging and a $rizK_s$-band-selected catalog of four million sources extracted using a fully model-based approach. We validate our photometry by comparing with the model-based DECam Legacy Survey. We analyze the differences between model-based and aperture photometry by comparing with the previous SHELA catalog, finding that our model-based photometry can measure point sources to fainter fluxes and better capture the full emission of resolved sources. The catalog is $80\%$ ($50\%$) complete at $riz \sim 24.7$ ($25.1$) AB mag, and the optical photometry reaches a $5σ$ depth of $\sim 25.5$ AB mag. We measure photometric redshifts and achieve $1σ$ scatter of $Δz/(1+z)$ of 0.04 with available spectroscopic redshifts at $0 \le z \le 1$. This large area, multi-wavelength photometric catalog, combined with spectroscopic information from HETDEX, will enable a wide range of extragalactic science investigations.

preprint2022arXiv

Lyα Halos around [O III]-Selected Galaxies in HETDEX

We present extended Lyman-α (Lyα) emission out to 800 kpc of 1034 [O III]-selected galaxies at redshifts 1.9<z<2.35 using the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). The locations and redshifts of the galaxies are taken from the 3D-HST survey. The median-stacked surface brightness profile of Lyα emission of the [O III]-selected galaxies agrees well with that of 968 bright Lyα-emitting galaxies (LAEs) at r>40 kpc from the galaxy centers. The surface brightness in the inner parts (r<10 kpc) around the [O III]-selected galaxies, however, is ten times fainter than that of the LAEs. Our results are consistent with the notion that photons dominating the outer regions of the Lyα halos are not produced in the central galaxies but originate outside of them.

preprint2022arXiv

Spectroscopic Confirmation of Two Luminous Post-AGB Stars in the Globular Cluster M19

The visually brightest stars in globular clusters (GCs) are the ones evolving off the asymptotic giant branch (AGB) and passing through spectral types F and A--the &#34;yellow&#34; post-AGB (yPAGB) stars. yPAGB stars are potentially excellent &#34;Population II&#34; standard candles for measuring extragalactic distances. A recent survey of the Galactic GC system, using uBVI photometry to identify stars of low surface gravities with large Balmer discontinuities, discovered a candidate luminous yPAGB star in the GC M19 (NGC 6273), designated ZNG 4. The same survey also identified a bright, hotter candidate blue PAGB star, ZNG 2, lying near the top of the white-dwarf cooling sequence. Both PAGB candidates have proper motions and parallaxes in the recent Gaia Early Data Release 3 consistent with cluster membership, but they still lacked spectroscopic verification. Here we present moderate-resolution spectra of both stars, confirming them as low-gravity objects that are extremely likely to be cluster members. Through comparison with a library of synthetic spectra, we made approximate estimates of the stars&#39; atmospheric parameters. We find that the yPAGB star ZNG 4 has an effective temperature of Teff ~ 6500 K, a surface gravity of log g ~ 1, and a metallicity of [Fe/H] ~ -1.5, similar to that of the host cluster. The blue PAGB star ZNG 2 has Teff ~ 18000 K, log g ~ 3, and an apparently low metallicity in the range of [Fe/H] ~ -2.0 to -2.5. Both stars are bright (V=12.5 and 13.3, respectively). We urge high-dispersion spectroscopic follow-up to determine detailed atmospheric parameters and chemical compositions, and to monitor radial velocities.

preprint2022arXiv

Stellar Populations of Lyman-alpha Emitting Galaxies in the HETDEX Survey I: An Analysis of LAEs in the GOODS-N Field

We present the results of a stellar-population analysis of Lyman-alpha emitting galaxies (LAES) in GOODS-N at 1.9 < z < 3.5 spectroscopically identified by the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). We provide a method for connecting emission-line detections from the blind spectroscopic survey to imaging counterparts, a crucial tool needed as HETDEX builds a massive database of ~1 million Lyman-alpha detections. Using photometric data spanning as many as 11 filters covering 0.4-4.5 microns from the Hubble and Spitzer Space Telescopes, we study the objects&#39; global properties and explore which properties impact the strength of Lyman-alpha emission. We measure a median stellar mass of 0.8 (^+2.9_-0.5) x 10^9 Msol and conclude that the physical properties of HETDEX spectroscopically-selected LAEs are comparable to LAEs selected by previous deep narrow band studies. We find that stellar mass and star formation rate correlate strongly with the Lyman-alpha equivalent width. We then use a known sample of z>7 LAEs to perform a proto-study of predicting Lyman-alpha emission from galaxies in the Epoch of Reionization, finding agreement at the 1-sigma level between prediction and observation for the majority of strong emitters.

preprint2022arXiv

Surface Brightness Profile of Lyman-$α$ Halos out to 320 kpc in HETDEX

We present the median-stacked Lyman-$α$ surface brightness profile of 968 spectroscopically selected Lyman-$α$ emitting galaxies (LAEs) at redshifts $1.9<z<3.5$ in the early data of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). The selected LAEs are high-confidence Lyman-$α$ detections with large signal-to-noise ratios observed with good seeing conditions (point-spread-function full-width-at-half-maximum $<1.4&#34;$), excluding active galactic nuclei (AGN). The Lyman-$α$ luminosities of the LAEs are $10^{42.4}-10^{43}\, \mathrm{erg}\, \mathrm{s}^{-1}$. We detect faint emission in the median-stacked radial profiles at the level of $(3.6\pm 1.3)\times 10^{-20}\,\mathrm{erg}\,\mathrm{s}^{-1}\,\mathrm{cm}^{-2}\,\mathrm{arcsec}^{-2}$ from the surrounding Lyman-$α$ halos out to $r\simeq 160$ kpc (physical). The shape of the median-stacked radial profile is consistent at $r<80\,\mathrm{kpc}$ with that of much fainter LAEs at $3<z<4$ observed with the Multi Unit Spectroscopic Explorer (MUSE), indicating that the median-stacked Lyman-$α$ profiles have similar shapes at redshifts $2<z<4$ and across a factor of $10$ in Lyman-$α$ luminosity. While we agree with the results from the MUSE sample at $r<80\,\mathrm{kpc}$, we extend the profile over a factor of two in radius. At $r>80\,\mathrm{kpc}$, our profile is flatter than the MUSE model. The measured profile agrees at most radii with that of galaxies in the Byrohl et al. (2021) cosmological radiative transfer simulation at $z=3$. This suggests that the surface brightness of a Lyman-$α$ halo at $r\lesssim 100$ kpc is dominated by resonant scattering of Lyman-$α$ photons from star-forming regions in the central galaxy, whereas at $r > 100$ kpc it is dominated by photons from galaxies in surrounding dark matter halos.

preprint2022arXiv

The Active Galactic Nuclei in the Hobby-Eberly Telescope Dark Energy Experiment Survey (HETDEX) I. Sample selection

We present the first Active Galactic Nuclei (AGN) catalog in the Hobby-Eberly Telescope Dark Energy Experiment Survey (HETDEX) observed between January 2017 and June 2020. HETDEX is an ongoing spectroscopic survey with no pre-selection based on magnitudes, colors or morphologies, enabling us to select AGN based on their spectral features. Both luminous quasars and low-luminosity Seyferts are found in our catalog. AGN candidates are selected with at least two significant AGN emission lines, such as the LyA and CIV line pair, or with single broad emission lines (FWHM > 1000 km/s). Each source is further confirmed by visual inspections. This catalog contains 5,322 AGN, covering an effective sky coverage of 30.61 deg^2. A total of 3,733 of these AGN have secure redshifts, and we provide redshift estimates for the remaining 1,589 single broad-line AGN with no cross matched spectral redshifts from SDSS DR14Q. The redshift range of the AGN catalog is 0.25 < z < 4.32, with a median of z = 2.1. The bolometric luminosity range is 10^9-10^14 Lsun with a median of 10^12 Lsun. The median r-band magnitude of the AGN is 21.6 mag, with 34% of the AGN have r > 22.5, and 2.6% reaching the detection limit at r ~ 26 mag of the deepest imaging surveys we searched. We also provide a composite spectrum of the AGN sample covering 700 AA - 4400 AA.

preprint2022arXiv

The Active Galactic Nuclei in the Hobby-Eberly Telescope Dark Energy Experiment Survey (HETDEX) II. Luminosity Function

We present the LyA emission line luminosity function (LF) of the Active Galactic Nuclei (AGN) in the first release of the Hobby-Eberly Telescope Dark Energy Experiment Survey (HETDEX) AGN catalog (Liu et al. 2022, Paper I). The AGN are selected either by emission-line pairs characteristic of AGN or by single broad emission line, free of any photometric pre-selections (magnitude/color/morphology). The sample consists of 2,346 AGN spanning 1.88<z<3.53, covering an effective area of 30.61 deg^2. Approximately 2.6 of the HETDEX AGN are not detected at $>5σ$ confidence at r~26 in the deepest $r$-band images we have searched. The LyA line luminosity ranges from ~10^42.3 to ~10^45.9 erg s^-1. Our LyA LF shows a turnover luminosity with opposite slopes on the bright end and the faint end: The space density is highest at L_LyA^*=10^43.4 erg s^-1. We explore the evolution of the AGN LF over a broader redshift range (0.8<z<3); constructing the rest-frame ultraviolet (UV) LF with the 1450 AA monochromatic luminosity of the power-law component of the continuum ($\rm M_{1450}$) from M_1450~-18 to ~-27.5. We divide the sample into three redshift bins (z~1.5, 2.1, and 2.6). In all three redshift bins, our UV LFs indicate that the space density of AGN is highest at the turnover luminosity M_1450^* with opposite slopes on the bright end and the faint end. The M_1450 LFs in the three redshift bins can be well-fit with a luminosity-evolution-density-evolution (LEDE) model: the turnover luminosity (M_1450^*) increases and the turnover density (Phi^*) decreases with increasing redshift.

preprint2022arXiv

Yellow Post-Asymptotic-Giant-Branch Stars as Standard Candles. I. Calibration of the Luminosity Function in Galactic Globular Clusters

We use results of a survey for low-surface-gravity stars in Galactic (and LMC) globular clusters to show that &#34;yellow&#34; post-asymptotic-giant-branch (yPAGB) stars are likely to be excellent extragalactic standard candles, capable of producing distances to early-type galaxies that are accurate to a few percent. We show that the mean bolometric magnitude of the 10 known yPAGB stars in globular clusters is <Mbol> = -3.38 +/- 0.03, a value that is ~0.2 mag brighter than that predicted from the latest post-horizontal-branch evolutionary tracks. More importantly, we show that the observed dispersion in the distribution is only 0.10 mag, i.e., smaller than the scatter for individual Cepheids. We describe the physics that can produce such a small dispersion, and show that, if one restricts surveys to the color range 0 < (B-V)0 < 0.5, then samples of non-variable yPAGB stars can be identified quite easily with a minimum of contamination. The bright absolute V magnitudes of these stars (<Mv> = -3.37) make them, by far, the visually brightest objects in old stellar populations and ideal Population II standard candles for measurements out to ~10 Mpc with current instrumentation. A Hubble Space Telescope survey in the halos of galaxies in the M81 and Sculptor groups could therefore serve as an effective cross-check on both the Cepheid and TRGB distance scales.

preprint2021arXiv

A Census of Above-Horizontal-Branch Stars in Galactic Globular Clusters

We have carried out a search for above-horizontal-branch (AHB) stars--objects lying above the horizontal branch (HB) and blueward of the asymptotic giant branch (AGB) in the color-magnitude diagram--in 97 Galactic and seven Magellanic Cloud globular clusters (GCs). We selected AHB candidates based on photometry in the $uBVI$ system, which is optimized for detection of low-gravity stars with large Balmer jumps, in the color range $-0.05\le(B-V)_0\le1.0$. We then used $Gaia$ astrometry and Gaussian-mixture modeling to confirm cluster membership and remove field interlopers. Our final catalog contains 438 AHB stars, classified and interpreted in the context of post-HB evolution as follows: (1) AHB1: 280 stars fainter than $M_V=-0.8$, evolving redward from the blue HB (BHB) toward the base of the AGB. (2) Post-AGB (PAGB): 13 stars brighter than $M_V\simeq-2.75$, departing from the top of the AGB and evolving rapidly blueward. (3) AHB2: 145 stars, with absolute magnitudes between those of the AHB1 and PAGB groups. This last category includes a mixture of objects leaving the extreme BHB and evolving toward the AGB, and brighter ones moving back from the AGB toward higher temperatures. Among the AHB1 stars are 59 RR Lyrae interlopers, observed by chance in our survey near maximum light. PAGB and AHB2 stars (including W Virginis Cepheids) overwhelmingly belong to GCs containing BHB stars, in accordance with predictions of post-HB evolutionary tracks. We suggest that most W Vir variables are evolving toward lower temperatures and are in their first crossings of the instability strip. Non-variable yellow PAGB stars show promise as a Population II standard candle for distance measurement.

preprint2021arXiv

HETDEX [OIII] Emitters I: A spectroscopically selected low-redshift population of low-mass, low-metallicity galaxies

We assemble a sample of 17 low metallicity (7.45 < log(O/H)+12 < 8.12) galaxies with z < 0.1 found spectroscopically, without photometric pre-selection, in early data from the Hobby Eberly Telescope Dark Energy Experiment (HETDEX). Star forming galaxies that occupy the lowest mass and metallicity end of the mass-metallicity relation tend to be under sampled in continuum-based surveys as their spectra are typically dominated by emission from newly forming stars. We search for galaxies with high [OIII]$λ$5007 / [OII]$λ$3727, implying highly ionized nebular emission often indicative of low metallicity systems. With the Second Generation Low Resolution Spectrograph on the Hobby Eberly Telescope we acquired follow-up spectra, with higher resolution and broader wavelength coverage, of each low-metallicity candidate in order to confirm the redshift, measure the H$α$ and [NII] line strengths and, in many cases, obtain deeper spectra of the blue lines. We find our galaxies are consistent with the mass-metallicity relation of typical low mass galaxies. However, galaxies in our sample tend to have similar specific star formation rates (sSFRs) as the incredibly rare &#34;blueberry&#34; galaxies found in (Yang et. al. 2017). We illustrate the power of spectroscopic surveys for finding low mass and metallicity galaxies and reveal that we find a sample of galaxies that are a hybrid between the properties of typical dwarf galaxies and the more extreme blueberry galaxies.

preprint2021arXiv

The Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) Survey Design, Reductions, and Detections

We describe the survey design, calibration, commissioning, and emission-line detection algorithms for the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). The goal of HETDEX is to measure the redshifts of over a million Ly$α$ emitting galaxies between 1.88<z<3.52, in a 540 deg^2 area encompassing a co-moving volume of 10.9 Gpc^3. No pre-selection of targets is involved; instead the HETDEX measurements are accomplished via a spectroscopic survey using a suite of wide-field integral field units distributed over the focal plane of the telescope. This survey measures the Hubble expansion parameter and angular diameter distance, with a final expected accuracy of better than 1%. We detail the project&#39;s observational strategy, reduction pipeline, source detection, and catalog generation, and present initial results for science verification in the COSMOS, Extended Groth Strip, and GOODS-N fields. We demonstrate that our data reach the required specifications in throughput, astrometric accuracy, flux limit, and object detection, with the end products being a catalog of emission-line sources, their object classifications, and flux-calibrated spectra.

preprint2021arXiv

Two Luminous Post-AGB Stars in the Galactic Globular Cluster M19

We report the discovery of a luminous &#34;yellow&#34; post-asymptotic-giant-branch (PAGB) star in the globular cluster (GC) M19 (NGC 6273), identified during our uBVI survey of Galactic GCs. The uBVI photometric system is optimized to detect stars with large Balmer discontinuities, indicating very low surface gravities and high luminosities. The spectral-energy distribution (SED) of the star is consistent with an effective temperature of about 6250 K and a surface gravity of $\log g=0.5$. We use Gaia data to show that the star&#39;s proper motion and radial velocity are consistent with cluster membership. One aim of our program is to test yellow PAGB stars as candidate Population II standard candles for determining extragalactic distances. We derive a visual absolute magnitude of $M_V=-3.39\pm0.09$ for the M19 star. This is in close agreement with the $M_V$ values found for yellow PAGB stars in the GCs omega Cen, NGC 5986, and M79, indicating a very narrow luminosity function. These objects are four magnitudes brighter than RR Lyrae variables, and they can largely avoid the issues of interstellar extinction that are a problem for Population I distance indicators. We also identified a second luminous PAGB object in M19, this one a hotter &#34;UV-bright&#34; star. Its SED is consistent with an effective temperature of about 11,750 K and $\log g=2.0$. The two objects have nearly identical bolometric luminosities, $\log L/L_\odot=3.24$ and 3.22, respectively.

preprint2020arXiv

Exploring AGN and Star Formation Activity of Massive Galaxies at Cosmic Noon

We investigate the relation between AGN and star formation (SF) activity at $0.5 < z < 3$ by analyzing 898 galaxies with X-ray luminous AGN ($L_X > 10^{44}$ erg s$^{-1}$) and a large comparison sample of $\sim 320,000$ galaxies without X-ray luminous AGN. Our samples are selected from a large (11.8 deg$^2$) area in Stripe 82 that has multi-wavelength (X-ray to far-IR) data. The enormous comoving volume ($\sim 0.3$ Gpc$^3$) at $0.5 < z < 3$ minimizes the effects of cosmic variance and captures a large number of massive galaxies ($\sim 30,000$ galaxies with $M_* > 10^{11} \ M_{\odot}$) and X-ray luminous AGN. While many galaxy studies discard AGN hosts, we fit the SED of galaxies with and without X-ray luminous AGN with Code Investigating GALaxy Emission (CIGALE) and include AGN emission templates. We find that without this inclusion, stellar masses and star formation rates (SFRs) in AGN host galaxies can be overestimated, on average, by factors of up to $\sim 5$ and $\sim 10$, respectively. The average SFR of galaxies with X-ray luminous AGN is higher by a factor of $\sim 3$ to $10$ compared to galaxies without X-ray luminous AGN at fixed stellar mass and redshift, suggesting that high SFRs and high AGN X-ray luminosities may be fueled by common mechanisms. The vast majority ($> 95 \%$) of galaxies with X-ray luminous AGN at $z=0.5-3$ do not show quenched SF: this suggests that if AGN feedback quenches SF, the associated quenching process takes a significant time to act and the quenched phase sets in after the highly luminous phases of AGN activity.

preprint2020arXiv

MCSED: A flexible spectral energy distribution fitting code and its application to $z \sim 2$ emission-line galaxies

We present MCSED, a new spectral energy distribution (SED)-fitting code, which mates flexible stellar evolution calculations with the Markov Chain Monte Carlo algorithms of the software package emcee. MCSED takes broad, intermediate, and narrow-band photometry, emission-line fluxes, and/or absorption line spectral indices, and returns probability distributions and co-variance plots for all model parameters. MCSED includes a variety of dust attenuation curves with parameters for varying the UV slopes and bump strengths, a prescription for continuum and PAH emission from dust, models for continuum and line emission from ionized gas, options for fixed and variable stellar metallicity, and a selection of star formation rate (SFR) histories. The code is well-suited for exploring parameter inter-dependencies in sets of galaxies with known redshifts, for which there is multi-band photometry and/or spectroscopy. We apply MCSED to a sample of $\sim2000$ $1.90<z<2.35$ galaxies in the five CANDELS fields, which were selected via their strong [O III] $\lambda5007$ emission, and explore the systematic behavior of their SEDs. We find the galaxies become redder with stellar mass, due to both increasing internal attenuation and a greater population of older stars. The slope of the UV extinction curve also changes with stellar mass, and at least some galaxies exhibit an extinction excess at 2175 Angstroms. Finally, we demonstrate that below $M\lesssim10^9\,M_{\odot}$), the shape of the star-forming galaxy main sequence is highly dependent on the galaxies&#39; assumed SFR history, as calculations which assume a constant SFR produce stellar masses that are $\sim1$ dex smaller than those found using more realistic SFR histories.

preprint2020arXiv

On the (Lack of) Evolution of the Stellar Mass Function of Massive Galaxies from $z$=1.5 to 0.4

We study the evolution in the number density of the highest mass galaxies over $0.4<z<1.5$ (covering 9 Gyr). We use the Spitzer/HETDEX Exploratory Large-Area (SHELA) Survey, which covers 17.5 $\mathrm{deg}^2$ with eight photometric bands spanning 0.3-4.5 $μ$m within the SDSS Stripe 82 field. This size produces the lowest counting uncertainties and cosmic variance yet for massive galaxies at $z\sim1.0$. We study the stellar mass function (SMF) for galaxies with $\log(M_\ast/M_\odot)>10.3$ using a forward-modeling method that fully accounts for statistical and systematic uncertainties on stellar mass. From $z$=0.4 to 1.5 the massive end of the SMF shows minimal evolution in its shape: the characteristic mass ($M^\ast$) evolves by less than 0.1 dex ($\pm$0.05 dex); the number density of galaxies with $\log (M_\ast/M_\odot) >11$ stays roughly constant at $\log (n/\mathrm{Mpc}^{-3})$ $\simeq$ $-$3.4 ($\pm$0.05), then declines to $\log (n/\mathrm{Mpc}^{-3})$=$-$3.7 ($\pm$0.05) at $z$=1.5. We discuss the uncertainties in the SMF, which are dominated by assumptions in the star formation history and details of stellar population synthesis models for stellar mass estimations. For quiescent galaxies, the data are consistent with no (or slight) evolution ($\lesssim0.1$ dex) in the characteristic mass nor number density from $z\sim 1.5$ to the present. This implies that any mass growth (presumably through &#34;dry&#34; mergers) of the quiescent massive galaxy population must balance the rate of mass losses from late-stage stellar evolution and the formation of quenching galaxies from the star-forming population. We provide a limit on this mass growth from $z=1.0$ to 0.4 of $ΔM_\ast/M_\ast\leq$ 45% (i.e., $\simeq0.16$ dex) for quiescent galaxies more massive than $10^{11}$ $M_\odot$.

preprint2020arXiv

The Atacama Cosmology Telescope: SZ-based masses and dust emission from IR-selected cluster candidates in the SHELA survey

We examine the stacked thermal Sunyaev-Zel\text{&#39;}dovich (SZ) signals for a sample of galaxy cluster candidates from the Spitzer-HETDEX Exploratory Large Area (SHELA) Survey, which are identified in combined optical and infrared SHELA data using the redMaPPer algorithm. We separate the clusters into three richness bins, with average photometric redshifts ranging from 0.70 to 0.80. The richest bin shows a clear temperature decrement at 148 GHz in the Atacama Cosmology Telescope data, which we attribute to the SZ effect. All richness bins show an increment at 220 GHz, which we attribute to dust emission from cluster galaxies. We correct for dust emission using stacked profiles from Herschel Stripe 82 data, and allow for synchrotron emission using stacked profiles created by binning source fluxes from NVSS data. We see dust emission in all three richness bins, but can only confidently detect the SZ decrement in the highest richness bin, finding $M_{500}$ = $8.7^{+1.7}_{-1.3} \times 10^{13} M_\odot$. Neglecting the correction for dust depresses the inferred mass by 26 percent, indicating a partial fill-in of the SZ decrement from thermal dust and synchrotron emission by the cluster member galaxies. We compare our corrected SZ masses to two redMaPPer mass--richness scaling relations and find that the SZ mass is lower than predicted by the richness. We discuss possible explanations for this discrepancy, and note that the SHELA richnesses may differ from previous richness measurements due to the inclusion of IR data in redMaPPer.

preprint2019arXiv

Exploring the High-Mass End of the Stellar Mass Function of Star Forming Galaxies at Cosmic Noon

We present the high-mass end of the galaxy stellar mass function using the largest sample to date (5,352) of star-forming galaxies with $M_{\star} > 10^{11} M_{\odot}$ at cosmic noon, $1.5 < z < 3.5$. This sample is uniformly selected across 17.2 deg$^2$ ($\sim$0.44 Gpc$^3$ comoving volume from $1.5 < z < 3.5$), mitigating the effects of cosmic variance and encompassing a wide range of environments. This area, a factor of 10 larger than previous studies, provides robust statistics at the high-mass end. Using multi-wavelength data in the Spitzer/HETDEX Exploratory Large Area (SHELA) footprint we find that the SHELA footprint star-forming galaxy stellar mass function is steeply declining at the high-mass end probing values as high as $\sim$$10^{-4}$ Mpc$^3$/dex and as low as $\sim$5$\times$$10^{-8}$ Mpc$^3$/dex across a stellar mass range of log($M_\star$/$M_\odot$) $\sim$ 11 - 12. We compare our empirical star-forming galaxy stellar mass function at the high mass end to three types of numerical models: hydrodynamical models from IllustrisTNG, abundance matching from the UniverseMachine, and three different semi-analytic models (SAMs; SAG, SAGE, GALACTICUS). At redshifts $1.5 < z < 3.5$ we find that results from IllustrisTNG and abundance matching models agree within a factor of $\sim$2 to 10, however the three SAMs strongly underestimate (up to a factor of 1,000) the number density of massive galaxies. We discuss the implications of these results for our understanding of galaxy evolution.