Researcher profile

Benne W. Holwerda

Benne W. Holwerda contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
21works
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

21 published item(s)

preprint2026arXiv

Galaxies in the Epoch of Reionization Are All Bark and No Bite -- Plenty of Ionizing Photons, Low Escape Fractions

Early results from JWST suggest that epoch-of-reionization (EoR) galaxies produce copious ionizing photons, which, if they escape efficiently, could cause reionization to occur too early. We study this problem using \jwst\ imaging and prism spectroscopy for 412 galaxies at 4.5 < z < 9.0. We fit these data simultaneously with stellar-population and nebular-emission models that include a parameter for the fraction of ionizing photons that escape the galaxy, $f_\mathrm{esc}$. We find that the ionization production efficiency, $ξ_\mathrm{ion}$ = Q(H) / L(UV), increases with redshift and decreasing UV luminosity, but shows significant scatter, $σ( \log ξ_\mathrm{ion})$ = 0.3 dex. The inferred escape fractions averaged over the population are low, ranging from $\langle f_\mathrm{esc} \rangle$ = $2.6\pm 1.4$\% at 6 < z < 9 to $6.5\pm 2.2$\% at 4.5 < z < 6 with weak or no indication of evolution with redshift. This implies that in our models most of the ionizing photons need to be absorbed to account for the nebular emission. We compute the impact of our results on reionization, including the distributions for $ξ_\mathrm{ion}$ and $f_\mathrm{esc}$, and the evolution and uncertainty of the UV luminosity function. Considering galaxies brighter than M(UV) < -16 mag, we would produce an IGM hydrogen-ionized fraction of $x_e = 0.5$ at 5.3 < z < 5.8, possibly too late compared to constraints from from QSO sightlines. Including fainter galaxies, M(UV) < -14 mag, we obtain $x_e = 0.5$ at 6.0 < z < 8.1, fully consistent with QSO and CMB data. This implies that EoR galaxies produce plenty of ionizing photons, but these do not efficiently escape. This may be a result of high gas column densities combined with burstier star-formation histories, which limit the time massive stars are able to clear channels through the gas for ionizing photons to escape.

preprint2026arXiv

PEARLS: 21 Transients Found in the Three-Epoch NIRCam Observations in the Continuous Viewing Zone of the James Webb Space Telescope

We present 21 transients from our three-epoch, four-band NIRCam observations covering 14.16 arcmin^2 in the Spitzer IRAC Dark Field (IDF), taken by the JWST Prime Extragalactic Areas for Reionization and Lensing Science program with a time cadence of ~6 months. A separate Hubble Space Telescope program provided Advanced Camera for Surveys optical imaging contemporaneous with the second and third epochs of the NIRCam observations. The NIRSpec spectroscopy on three transients confirmed a Type Ia supernova at z=1.63 and the host galaxies of the other two at z=2.64 and 1.90, respectively. Combining these with the photometric redshifts (z_ph) of the host galaxies in the rest of the sample, we find that the transients are in either a &#34;mid-z&#34; group at z>1.6 with M_V < -16.0 mag or a &#34;low-z&#34; group at z < 0.4 with M_H > -14.0 mag. The mid-z transients are consistent with supernovae. In contrast, the low-z transients&#39; luminosities fall in the range of the so-called &#34;gap transients&#34; between supernovae and novae. However, this latter conclusion is only tentative due to possible catastrophic failures in z_ph that could bias them to low-z. Conversely, if they are indeed at z < 0.4, it would be worth studying similar transients in the future. Our work further demonstrates the power of NIRCam in transient science and also shows that it would be more fruitful to carry out a long-term monitoring program with more passbands, a higher cadence and prompt follw-up spectroscopy. Being in the continuous viewing zone of the JWST, the IDF is an ideal field for this purpose.

preprint2024arXiv

PEARLS: A Potentially Isolated Quiescent Dwarf Galaxy with a TRGB Distance of 30 Mpc

A wealth of observations have long suggested that the vast majority of isolated classical dwarf galaxies ($M_*=10^7$-$10^9$ M$_\odot$) are currently star-forming. However, recent observations of the large abundance of &#34;Ultra-Diffuse Galaxies&#34; beyond the reach of previous large spectroscopic surveys suggest that our understanding of the dwarf galaxy population may be incomplete. Here we report the serendipitous discovery of an isolated quiescent dwarf galaxy in the nearby Universe, which was imaged as part of the PEARLS GTO program. Remarkably, individual red-giant branch stars are visible in this near-IR imaging, suggesting a distance of $30\pm4$ Mpc, and a wealth of archival photometry point to an sSFR of $2\times10^{-11}$ yr$^{-1}$ and SFR of $4\times10^{-4}$ M$_\odot$ yr$^{-1}$. Spectra obtained with the Lowell Discovery Telescope find a recessional velocity consistent with the Hubble Flow and ${>}1500$ km/s separated from the nearest massive galaxy in SDSS, suggesting that this galaxy was either quenched from internal mechanisms or had a very high-velocity ($>1000$ km/s) interaction with a nearby massive galaxy in the past. This analysis highlights the possibility that many nearby quiescent dwarf galaxies are waiting to be discovered and that JWST has the potential to resolve them.

preprint2022arXiv

An empirical measurement of the Halo Mass Function from the combination of GAMA DR4, SDSS DR12, and REFLEX II data

We construct the halo mass function (HMF) from the GAMA galaxy group catalogue over the mass range 10^12.7M_sol to 10^15.5M_sol, and find good agreement with the expectation from LambdaCDM. In comparison to previous studies, this result extends the mass range over which the HMF has now been measured over by an order of magnitude. We combine the GAMA DR4 HMF with similar data from the SDSS DR12 and REFLEX II surveys, and fit a four-parameter Murray-Robotham-Power (MRP) function, valid at z~0.1, yielding: a density normalisation of: log10 (phi Mpc^3)=-3.96[+0.55,-0.82], a high mass turn-over of: log10(M/M_sol)=14.13[+0.43,-0.40], a low mass power law slope of: alpha=-1.68[+0.21,-0.24] , and a high mass softening parameter of: beta= 0.63[+0.25,-0.11]. If we fold in the constraint on Omega_M from Planck 2018 Cosmology, we are able to reduce these uncertainties further, but this relies on the assumption that the power-law trend can be extrapolated from 10^12.7M_sol to zero mass. Throughout, we highlight the effort needed to improve on our HMF measurement: improved halo mass estimates that do not rely on calibration to simulations; reduced halo mass uncertainties needed to mitigate the strong Eddington Bias that arises from the steepness of the HMF low mass slope; and deeper wider area spectroscopic surveys. To our halo mass limit of 10^12.7 M_sol, we are directly resolving (`seeing&#39;) 41+/-5 per cent of the total mass density, i.e. Omega_[M>12.7]=0.128+/-0.016, opening the door for the direct construction of 3D dark matter mass maps at Mpc resolution.

preprint2022arXiv

Beyond the Local Volume. II. Population Scaleheights and Ages of Ultracool Dwarfs in Deep HST/WFC3 Parallel Fields

Ultracool dwarfs represent a significant proportion of stars in the Milky Way,and deep samples of these sources have the potential to constrain the formation history and evolution of low-mass objects in the Galaxy. Until recently, spectral samples have been limited to the local volume (d<100 pc). Here, we analyze a sample of 164 spectroscopically-characterized ultracool dwarfs identified by Aganze et al. (2022) in the Hubble Space Telescope WFC3 Infrared Spectroscopic Parallel (WISP) Survey and 3D-HST. We model the observed luminosity function using population simulations to place constraints on scaleheights, vertical velocity dispersions and population ages as a function of spectral type. Our star counts are consistent with a power-law mass function and constant star formation history for ultracool dwarfs, with vertical scaleheights 249$_{-61}^{+48}$ pc for late M dwarfs, 153$_{-30}^{+56}$ pc for L dwarfs, and 175$_{-56}^{+149}$ pc for T dwarfs. Using spatial and velocity dispersion relations, these scaleheights correspond to disk population ages of 3.6$_{-1.0}^{+0.8}$ for late M dwarfs, 2.1$_{-0.5}^{+0.9}$ Gyr for L dwarfs, and 2.4$_{-0.8}^{+2.4}$ Gyr for T dwarfs, which are consistent with prior simulations that predict that L-type dwarfs are on average a younger and less dispersed population. There is an additional 1-2 Gyr systematic uncertainty on these ages due to variances in age-velocity relations. We use our population simulations to predict the UCD yield in the JWST PASSAGES survey, a similar and deeper survey to WISPS and 3D-HST, and find that it will produce a comparably-sized UCD sample, albeit dominated by thick disk and halo sources.

preprint2022arXiv

Deep Extragalactic VIsible Legacy Survey (DEVILS): The emergence of bulges and decline of disk growth since $z = 1$

We present a complete structural analysis of the ellipticals (E), diffuse bulges (dB), compact bulges (cB), and disks (D) within a redshift range $0 < z < 1$, and stellar mass $\log_{10}(\mathrm{M}_*/\mathrm{M}_\odot) \geq 9.5$ volume-limited sample drawn from the combined DEVILS and HST-COSMOS region. We use the {\sc ProFit} code to profile over $\sim35,000$ galaxies for which visual classification into single or double-component was predefined in Paper-I. Over this redshift range, we see a growth in the total stellar mass density (SMD) of a factor of 1.5. At all epochs we find that the dominant structure, contributing to the total SMD, is the disk, and holds a fairly constant share of $\sim60\%$ of the total SMD from $z = 0.8$ to $z = 0.2$, dropping to $\sim30\%$ at $z = 0.0$ (representing $\sim33\%$ decline in the total disk SMD). Other classes (E, dB, and cB) show steady growth in their numbers and integrated stellar mass densities. By number, the most dramatic change across the full mass range is in the growth of diffuse bulges. In terms of total SMD, the biggest gain is an increase in massive elliptical systems, rising from 20\% at $z = 0.8$ to equal that of disks at $z = 0.0$ (30\%) representing an absolute mass growth of a factor of 2.5. Overall we see a clear picture of the emergence and growth of all three classes of spheroids over the past 8 Gyrs, and infer that in the later half of the Universe&#39;s timeline spheroid forming-processes and pathways (secular evolution, mass-accretion, and mergers) appear to dominate mass transformation over quiescent disk growth.

preprint2022arXiv

Exploring the effect of baryons on the radial distribution of satellite galaxies with GAMA and IllustrisTNG

We explore the radial distribution of satellite galaxies in groups in the Galaxy and Mass Assembly (GAMA) survey and the IllustrisTNG simulations. Considering groups with masses $12.0 \leq \log_{10} (\mathcal{M}_h / h^{-1} \mathrm{M}_{\odot}) < 14.8$ at $z<0.267$, we find a good agreement between GAMA and a sample of TNG300 groups and galaxies designed to match the GAMA selection. Both display a flat profile in the centre of groups, followed by a decline that becomes steeper towards the group edge, and normalised profiles show no dependence on group mass. Using matched satellites from TNG and dark matter-only TNG-Dark runs we investigate the effect of baryons on satellite radial location. At $z=0$, we find that the matched subhaloes from the TNG-Dark runs display a much flatter radial profile: namely, satellites selected above a minimum stellar mass exhibit both smaller halo-centric distances and longer survival times in the full-physics simulations compared to their dark-matter only analogues. We then divide the TNG satellites into those which possess TNG-Dark counterparts and those which do not, and develop models for the radial positions of each. We find the satellites with TNG-Dark counterparts are displaced towards the halo centre in the full-physics simulations, and this difference has a power-law behaviour with radius. For the `orphan&#39; galaxies without TNG-Dark counterparts, we consider the shape of their radial distribution and provide a model for their motion over time, which can be used to improve the treatment of satellite galaxies in semi-analytic and semi-empirical models of galaxy formation.

preprint2022arXiv

Galactic Component Mapping of Galaxy UGC 2885 by Machine Learning Classification

Automating classification of galaxy components is important for understanding the formation and evolution of galaxies. Traditionally, only the larger galaxy structures such as the spiral arms, bulge, and disc are classified. Here we use machine learning (ML) pixel-by-pixel classification to automatically classify all galaxy components within digital imagery of massive spiral galaxy UGC 2885. Galaxy components include young stellar population, old stellar population, dust lanes, galaxy center, outer disc, and celestial background. We test three ML models: maximum likelihood classifier (MLC), random forest (RF), and support vector machine (SVM). We use high-resolution Hubble Space Telescope (HST) digital imagery along with textural features derived from HST imagery, band ratios derived from HST imagery, and distance layers. Textural features are typically used in remote sensing studies and are useful for identifying patterns within digital imagery. We run ML classification models with different combinations of HST digital imagery, textural features, band ratios, and distance layers to determine the most useful information for galaxy component classification. Textural features and distance layers are most useful for galaxy component identification, with the SVM and RF models performing the best. The MLC model performs worse overall but has comparable performance to SVM and RF in some circumstances. Overall, the models are best at classifying the most spectrally unique galaxy components including the galaxy center, outer disc, and celestial background. The most confusion occurs between the young stellar population, old stellar population, and dust lanes. We suggest further experimentation with textural features for astronomical research on small-scale galactic structures.

preprint2022arXiv

Galapagos-2/Galfitm/GAMA -- multi-wavelength measurement of galaxy structure: separating the properties of spheroid and disk components in modern surveys

We present the capabilities of Galapagos--2 and Galfitm in the context of fitting 2-component profiles to galaxies, on the way to providing complete multi-band, multi-component fitting of large samples of galaxies in future surveys. We release both the code and the fit results to 234,239 objects from the DR3 of the Gama survey, a sample significantly deeper than previous works. We use stringent tests on both simulated and real data, as well as comparison to public catalogues to evaluate the advantages of using multi-band over single-band data. We show that multi-band fitting using Galfitm provides significant advantages when trying to decompose galaxies into their individual constituents, as more data are being used, by effectively being able to use the colour information buried in the individual exposures to its advantage. Using simulated data, we find that multi-band fitting significantly reduces the deviations from real parameter values, allows component sizes and Sérsic indices to be recovered more accurately, and, by design, constrains the band-to-band variations of these parameters to more physical values. On both simulated and real data, we confirm that the SEDs of the 2 main components can be recovered to fainter magnitudes compared to using single-band fitting, which tends to recover disks and bulges to - on average - have identical SEDs when the galaxies become too faint, instead of the different SEDs they truly have. By comparing our results to those provided by other fitting codes, we confirm that they agree in general, but measurement errors can be significantly reduced by using the multi-band tools developed by the MegaMorph project. We conclude that the multi-band fitting employed by Galapagos-2 and Galfitm significantly improves the accuracy of structural galaxy parameters and enables much larger samples to be be used in a scientific analysis.

preprint2022arXiv

Galaxy And Mass Assembly (GAMA): Bulge-disk decomposition of KiDS data in the nearby universe

We derive single Sérsic fits and bulge-disk decompositions for 13096 galaxies at redshifts z < 0.08 in the GAMA II equatorial survey regions in the Kilo-Degree Survey (KiDS) g, r and i bands. The surface brightness fitting is performed using the Bayesian two-dimensional profile fitting code ProFit. We fit three models to each galaxy in each band independently with a fully automated Markov-chain Monte Carlo analysis: a single Sérsic model, a Sérsic plus exponential and a point source plus exponential. After fitting the galaxies, we perform model selection and flag galaxies for which none of our models are appropriate (mainly mergers/Irregular galaxies). The fit quality is assessed by visual inspections, comparison to previous works, comparison of independent fits of galaxies in the overlap regions between KiDS tiles and bespoke simulations. The latter two are also used for a detailed investigation of systematic error sources. We find that our fit results are robust across various galaxy types and image qualities with minimal biases. Errors given by the MCMC underestimate the true errors typically by factors 2-3. Automated model selection criteria are accurate to > 90 % as calibrated by visual inspection of a subsample of galaxies. We also present g-r component colours and the corresponding colour-magnitude diagram, consistent with previous works despite our increased fit flexibility. Such reliable structural parameters for the components of a diverse sample of galaxies across multiple bands will be integral to various studies of galaxy properties and evolution. All results are integrated into the GAMA database.

preprint2022arXiv

Galaxy And Mass Assembly (GAMA): Data Release 4 and the z < 0.1 total and z < 0.08 morphological galaxy stellar mass functions

In Galaxy And Mass Assembly Data Release 4 (GAMA DR4), we make available our full spectroscopic redshift sample. This includes 248682 galaxy spectra, and, in combination with earlier surveys, results in 330542 redshifts across five sky regions covering ~250deg^2. The redshift density, is the highest available over such a sustained area, has exceptionally high completeness (95 per cent to r_KIDS=19.65mag), and is well suited for the study of galaxy mergers, galaxy groups, and the low redshift (z<0.25) galaxy population. DR4 includes 32 value-added tables or Data Management Units (DMUs) that provide a number of measured and derived data products including GALEX, ESO KiDS, ESO VIKING, WISE and Herschel Space Observatory imaging. Within this release, we provide visual morphologies for 15330 galaxies to z<0.08, photometric redshift estimates for all 18million objects to r_KIDS~25mag, and stellar velocity dispersions for 111830 galaxies. We conclude by deriving the total galaxy stellar mass function (GSMF) and its sub-division by morphological class (elliptical, compact-bulge and disc, diffuse-bulge and disc, and disc only). This extends our previous measurement of the total GSMF down to 10^6.75 M_sol h^-2_70 and we find a total stellar mass density of rho_*=(2.97+/-0.04)x10^8 M_sol h_70 Mpc^-3 or Omega_*=(2.17+/-0.03)x10^-3 h^-1_70. We conclude that at z<0.1, the Universe has converted 4.9+/-0.1 per cent of the baryonic mass implied by Big Bang Nucleosynthesis into stars that are gravitationally bound within the galaxy population.

preprint2022arXiv

Identifying Galaxy Mergers in Simulated CEERS NIRCam Images using Random Forests

Identifying merging galaxies is an important - but difficult - step in galaxy evolution studies. We present random forest classifications of galaxy mergers from simulated JWST images based on various standard morphological parameters. We describe (a) constructing the simulated images from IllustrisTNG and the Santa Cruz SAM, and modifying them to mimic future CEERS observations as well as nearly noiseless observations, (b) measuring morphological parameters from these images, and (c) constructing and training the random forests using the merger history information for the simulated galaxies available from IllustrisTNG. The random forests correctly classify $\sim60\%$ of non-merging and merging galaxies across $0.5 < z < 4.0$. Rest-frame asymmetry parameters appear more important for lower redshift merger classifications, while rest-frame bulge and clump parameters appear more important for higher redshift classifications. Adjusting the classification probability threshold does not improve the performance of the forests. Finally, the shape and slope of the resulting merger fraction and merger rate derived from the random forest classifications match with theoretical Illustris predictions, but are underestimated by a factor of $\sim 0.5$.

preprint2022arXiv

JWST&#39;s PEARLS: A JWST/NIRCam view of ALMA sources

We report the results of James Webb Space Telescope/NIRCam observations of 19 (sub)millimeter (submm/mm) sources detected by the Atacama Large Millimeter Array (ALMA). The accurate ALMA positions allowed unambiguous identifications of their NIRCam counterparts. Taking gravitational lensing into account, these represent 16 distinct galaxies in three fields and constitute the largest sample of its kind to date. The counterparts&#39; spectral energy distributions from rest-frame ultraviolet to near infrared provide photometric redshifts ($1<z<4.5$) and stellar masses ($M_*>10^{10.5}$ Msol), which are similar to sub-millimeter galaxy (SMG) hosts studied previously. However, our sample is fainter in submm/mm than the classic SMG samples are, and our sources exhibit a wider range of properties. They have dust-embedded star-formation rates as low as 10 Msol yr$^{-1}$, and the sources populate both the star-forming main sequence and the quiescent categories. The deep NIRCam data allow us to study the rest-frame near-IR morphologies. Excluding two multiply imaged systems and one quasar, the majority of the remaining sources are disk-like and show either little or no disturbance. This suggests that secular growth is a potential route for the assembly of high-mass disk galaxies. While a few hosts have large disks, the majority have small disks (median half-mass radius of 1.6 kpc). At this time, it is unclear whether this is due to the prevalence of small disks at these redshifts or some unknown selection effects of deep ALMA observations. A larger sample of ALMA sources with NIRCam observations will be able to address this question.

preprint2022arXiv

LADUMA: Discovery of a luminous OH megamaser at $z > 0.5$

In the local Universe, OH megamasers (OHMs) are detected almost exclusively in infrared-luminous galaxies, with a prevalence that increases with IR luminosity, suggesting that they trace gas-rich galaxy mergers. Given the proximity of the rest frequencies of OH and the hyperfine transition of neutral atomic hydrogen (HI), radio surveys to probe the cosmic evolution of HI in galaxies also offer exciting prospects for exploiting OHMs to probe the cosmic history of gas-rich mergers. Using observations for the Looking At the Distant Universe with the MeerKAT Array (LADUMA) deep HI survey, we report the first untargeted detection of an OHM at $z > 0.5$, LADUMA J033046.20$-$275518.1 (nicknamed &#34;Nkalakatha&#34;). The host system, WISEA J033046.26$-$275518.3, is an infrared-luminous radio galaxy whose optical redshift $z \approx 0.52$ confirms the MeerKAT emission line detection as OH at a redshift $z_{\rm OH} = 0.5225 \pm 0.0001$ rather than HI at lower redshift. The detected spectral line has 18.4$σ$ peak significance, a width of $459 \pm 59\,{\rm km\,s^{-1}}$, and an integrated luminosity of $(6.31 \pm 0.18\,{\rm [statistical]}\,\pm 0.31\,{\rm [systematic]}) \times 10^3\,L_\odot$, placing it among the most luminous OHMs known. The galaxy&#39;s far-infrared luminosity $L_{\rm FIR} = (1.576 \pm 0.013) \times 10^{12}\,L_\odot$ marks it as an ultra-luminous infrared galaxy; its ratio of OH and infrared luminosities is similar to those for lower-redshift OHMs. A comparison between optical and OH redshifts offers a slight indication of an OH outflow. This detection represents the first step towards a systematic exploitation of OHMs as a tracer of galaxy growth at high redshifts.

preprint2022arXiv

The Subaru HSC weak lensing mass-observable scaling relations of spectroscopic galaxy groups from the GAMA survey

We utilize the galaxy shape catalogue from the first-year data release of the Subaru Hyper Suprime-cam Survey (HSC) to study the dark matter content of galaxy groups in the Universe using weak lensing. We use galaxy groups from the Galaxy Mass and Assembly galaxy survey in approximately $100$ sq. degrees of the sky that overlap with the HSC survey as lenses. We restrict our analysis to the $1587$ groups with at least five members. We divide these groups into six bins each of group luminosity and group member velocity dispersion and measure the lensing signal with a signal-to-noise ratio of $55$ and $51$ for these two different selections, respectively. We use a Bayesian halo model framework to infer the halo mass distribution of our groups binned in the two different observable properties and constrain the power-law scaling relation, and the scatter between mean halo masses and the two group observable properties. We obtain a 5 percent constraint on the amplitude of the scaling relation between halo mass and group luminosity with $\avg{M} = (0.81\pm 0.04)\times10^{14}\hinvMsun$ for $L_{\rm grp}=10^{11.5}\hinvsqLsun$, and a power-law index of $α=1.01\pm 0.07$. We constrain the amplitude of the scaling relation between halo mass and velocity dispersion to be $\avg{M}=(0.93\pm 0.05)\times10^{14}\hinvMsun$ for $σ=500 \kms$ and a power-law index to be $α=1.52\pm0.10$. However, these scaling relations are sensitive to the exact cuts applied to the number of group members. Comparisons with similar scaling relations from the literature show that our results are consistent and have significantly reduced errors.

preprint2022arXiv

The variation of the gas content of galaxy groups and pairs compared to isolated galaxies

We measure how the atomic gas (HI) fraction ($f_{HI}={\rm \frac{M_{HI}}{M_{*}}}$) of groups and pairs taken as single units vary with average stellar mass ($\langle {\rm M_*} \rangle$) and average star-formation rate ($\langle {\rm SFR} \rangle$), compared to isolated galaxies. The HI 21 cm emission observation are from (i) archival ALFALFA survey data covering three fields from the GAMA survey (provides environmental and galaxy properties), and (ii) DINGO pilot survey data of one of those fields. The mean $f_{HI}$ for different units (groups/pairs/isolated galaxies) are measured in regions of the log($\langle {\rm M_*} \rangle$) -- log($\langle {\rm SFR} \rangle$) plane, relative to the z $\sim 0$ star-forming main sequence (SFMS) of individual galaxies, by stacking $f_{HI}$ spectra of individual units. For ALFALFA, $f_{HI}$ spectra of units are measured by extracting HI spectra over the full groups/pair areas and dividing by the total stellar mass of member galaxies. For DINGO, $f_{HI}$ spectra of units are measured by co-adding HI spectra of individual member galaxies, followed by division by their total stellar mass. For all units the mean $f_{HI}$ decreases as we move to higher $\langle {\rm M_*} \rangle$ along the SFMS, and as we move from above the SFMS to below it at any $\langle {\rm M_*} \rangle$. From the DINGO-based study, mean $f_{HI}$ in groups appears to be lower compared to isolated galaxies for all $\langle {\rm M_*} \rangle$ along the SFMS. From the ALFALFA-based study we find substantially higher mean $f_{HI}$ in groups compared to isolated galaxies (values for pairs being intermediate) for ${\langle{\rm M_*}\rangle}\lesssim10^{9.5}~{\rm M_{\odot}}$, indicating the presence of substantial amounts of HI not associated with cataloged member galaxies in low mass groups.

preprint2021arXiv

Galaxy and Mass Assembly (GAMA): The Weak Environmental Dependence of Quasar Activity at 0.1<z<0.35

Understanding the connection between nuclear activity and galaxy environment remains critical in constraining models of galaxy evolution. By exploiting extensive catalogued data from the Galaxy and Mass Assembly (GAMA) survey, we identify a representative sample of 205 quasars at 0.1 < z < 0.35 and establish a comparison sample of galaxies, closely matched to the quasar sample in terms of both stellar mass and redshift. On scales <1 Mpc, the galaxy number counts and group membership of quasars appear entirely consistent with those of the matched galaxy sample. Despite this, we find that quasars are ~1.5 times more likely to be classified as the group center, indicating a potential link between quasar activity and cold gas flows or galaxy interactions associated with rich group environments. On scales of ~a few Mpc, the clustering strength of both samples are statistically consistent and beyond 10 Mpc we find no evidence that quasars trace large scale structures any more than the galaxy control sample. Both populations are found to prefer intermediate-density sheets and filaments to either very high- or very low- density environments. This weak dependence of quasar activity on galaxy environment supports a paradigm in which quasars represent a phase in the lifetime of all massive galaxies and in which secular processes and a group-centric location are the dominant trigger of quasars at low redshift.

preprint2021arXiv

Significant Dust-Obscured Star Formation in Luminous Lyman-break Galaxies at $z$$\sim$$7$$-$$8$

We make use of ALMA continuum observations of $15$ luminous Lyman-break galaxies at $z$$\sim$$7$$-$$8$ to probe their dust-obscured star-formation. These observations are sensitive enough to probe to obscured SFRs of $20$ $M_{\odot}$$/$$yr$ ($3σ$). Six of the targeted galaxies show significant ($\geq$$3$$σ$) dust continuum detections, more than doubling the number of known dust-detected galaxies at $z$$>$$6.5$. Their IR luminosities range from $2.7$$\times$$10^{11}$ $L_{\odot}$ to $1.1$$\times$$10^{12}$ $L_{\odot}$, equivalent to obscured SFRs of $20$ to $105$ $M_{\odot}$$/$$yr$. We use our results to quantify the correlation of the infrared excess IRX on the UV-continuum slope $β_{UV}$ and stellar mass. Our results are most consistent with an SMC attenuation curve for intrinsic $UV$-slopes $β_{UV,intr}$ of $-2.63$ and most consistent with an attenuation curve in-between SMC and Calzetti for $β_{UV,intr}$ slopes of $-2.23$, assuming a dust temperature $T_d$ of $50$ K. Our fiducial IRX-stellar mass results at $z$$\sim$$7$$-$$8$ are consistent with marginal evolution from $z$$\sim$$0$. We then show how both results depend on $T_d$. For our six dust-detected sources, we estimate their dust masses and find that they are consistent with dust production from SNe if the dust destruction is low ($<$$90$%). Finally we determine the contribution of dust-obscured star formation to the star formation rate density for $UV$ luminous ($<$$-$$21.5$ mag: $\gtrsim$$1.7$$L_{UV} ^*$) $z$$\sim$$7$$-$$8$ galaxies, finding that the total SFR density at $z$$\sim$$7$ and $z$$\sim$$8$ from bright galaxies is $0.18_{-0.10}^{+0.08}$ dex and $0.20_{-0.09}^{+0.05}$ dex higher, respectively, i.e. $\sim$$\frac{1}{3}$ of the star formation in $\gtrsim$$1.7$$L_{UV} ^*$ galaxies at $z$$\sim$$7$$-$$8$ is obscured by dust.

preprint2020arXiv

Is NGC 300 a pure exponential disk galaxy?

NGC 300 is a low-mass disk galaxy in the Sculptor group. In the literature, it has been identified as a pure exponential disk galaxy, as its luminosity profile could be well fitted with a single exponential law over many disk scale lengths (Type I). We investigate the stellar luminosity distribution of NGC 300 using $Hubble$ $Space$ $Telescope$ (HST) archive data, reaching farther and deeper than any other previous studies. Color magnitude diagrams show a significant population of old red giant branch (RGB) stars in all fields out to $R\sim19$ kpc ($32&#39;$), as well as younger populations in the inner regions. We construct the density profiles of the young, intermediate-aged, and old stellar populations. We find two clear breaks in the density profiles of the old RGB and intermediate-aged stars: one down-bending (Type II) at $R\sim5.9$ kpc, and another up-bending (Type III) at $R\sim8.3$ kpc. Moreover, the old RGB stars exhibit a negative radial color gradient with an up-bending at $R\sim8$~kpc, beyond which the stellar populations are uniformly old ($>$7~Gyr) and metal-poor ($\rm[Fe/H] = -1.6^{+0.2}_{-0.4}$ dex). The outer stellar component at $R\gtrapprox8$ kpc is, therefore, well separated from the inner disk in terms of the stellar density and stellar populations. While our results cast doubt on the currently established wisdom that NGC\,300 is a pure exponential disk galaxy, a more detailed survey should be carried out to identify the outskirts as either a disk or a stellar halo.

preprint2020arXiv

The Sizes of $z\sim9-10$ Galaxies Identified in the BoRG Survey

Redshift $z=9--10$ object selection is the effective limit of Hubble Space Telescope imaging capability, even when confirmed with Spitzer. If only a few photometry data points are available, it becomes attractive to add criteria based on their morphology in these J- and H-band images. One could do so through visual inspection, a size criterion, or alternate morphometrics. We explore a vetted sample of BoRG $z\sim9$ and $z\sim10$ candidate galaxies and the object rejected by Morishita+ (2018) to explore the utility of a size criterion in z=9-10 candidate selection. A stringent, PSF-corrected effective radius criterion ($r_e<0\farcs3$) would result in the rejection of 65-70\% of the interlopers visually rejected by Morishita+. It may also remove up to $\sim20$\% of bona-fide brightest ($L>>L^*$) z=9 or 10 candidates from a BoRG selected sample based on the Mason+ (2015) luminosity functions, assuming the Holwerda+ (2015) $z\sim9$ size-luminosity relation. We argue that including a size constraint in lieu of a visual inspection may serve in wide-field searches for these objects in e.g. EUCLID or HST archival imaging with the understanding that some brightest ($L>>L^*$) candidates may be missed. The sizes of the candidates found by Morishita+ (2018) follow the expected size distribution of $z\sim9$ for bright galaxies, consistent with the lognormal in Shibuya+ (2015) and single objects. Two candidates show high star-formation surface density ($Σ_{SFR} > 25 M_\odot/kpc^2$) and all merit further investigation and follow-up observations.

preprint2020arXiv

Tracing the anemic stellar halo of M101

Models of galaxy formation in a cosmological context predict that massive disk galaxies should have structured extended stellar halos. Recent studies in integrated light, however, report a few galaxies, including the nearby disk galaxy M101, that have no measurable stellar halos to the detection limit. We aim to quantify the stellar content and structure of M101&#39;s outskirts by resolving its stars. We present the photometry of its stars based on deep F606W and F814W images taken with Hubble Space Telescope as part of the GHOSTS survey. The constructed CMDs of stars reach down to two magnitudes below the tip of the red giant branch. We derived radial number density profiles of the bright red giant branch (RGB) stars. The mean color of the RGB stars at $R \sim$ 40 -- 60 kpc is similar to those of metal-poor globular clusters in the Milky Way. We also derived radial surface brightness profiles using the public image data provided by the Dragonfly team. Both the radial number density and surface brightness profiles were converted to radial mass density profiles and combined. We find that the mass density profiles show a weak upturn at the very outer region, where surface brightness is as faint as $μ_g\approx 34$ mag arcsec$^{-1}$. An exponential disk + power-law halo model on the mass density profiles finds the total stellar halo mass of $M_{halo}=8.2_{-2.2}^{+3.5}\times 10^7M_\odot$. The total stellar halo mass does not exceed $M_{halo} = 3.2 \times 10^8$ $M_{\odot}$ when strongly truncated disk models are considered. Combining the halo mass with the total stellar mass of M101, we obtain the stellar halo mass fraction of $M_{halo}/M_{gal} = 0.20_{-0.08}^{+0.10}\%$ with an upper limit of 0.78\%. We compare the halo properties of M101 with those of six GHOSTS survey galaxies as well as the Milky Way and M31 and find that M101 has an anemic stellar halo.