Source author record

Benne W. Holwerda

Benne W. Holwerda appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

48works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

48 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 "mid-z" group at z>1.6 with M_V < -16.0 mag or a "low-z" 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' luminosities fall in the range of the so-called "gap transients" 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 "Ultra-Diffuse Galaxies" 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') 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'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' 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'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' 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 "Nkalakatha"). 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'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'$), 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'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.

preprint2016arXiv

Dependence of GAMA galaxy halo masses on the cosmic web environment from 100 square degrees of KiDS weak lensing data

Galaxies and their dark matter haloes are part of a complex network of mass structures, collectively called the cosmic web. Using the tidal tensor prescription these structures can be classified into four cosmic environments: voids, sheets, filaments and knots. As the cosmic web may influence the formation and evolution of dark matter haloes and the galaxies they host, we aim to study the effect of these cosmic environments on the average mass of galactic haloes. To this end we measure the galaxy-galaxy lensing profile of 91,195 galaxies, within 0.039 < z < 0.263, from the spectroscopic Galaxy And Mass Assembly (GAMA) survey, using ~100 square degrees of overlapping data from the Kilo-Degree Survey (KiDS). In each of the four cosmic environments we model the contributions from group centrals, satellites and neighbouring groups to the stacked galaxy-galaxy lensing profiles. After correcting the lens samples for differences in the stellar mass distribution, we find no dependence of the average halo mass of central galaxies on their cosmic environment. We do find a significant increase in the average contribution of neighbouring groups to the lensing profile in increasingly dense cosmic environments. We show, however, that the observed effect can be entirely attributed to the galaxy density at much smaller scales (within 4 Mpc/h), which is correlated with the density of the cosmic environments. Within our current uncertainties we find no direct dependence of galaxy halo mass on their cosmic environment.

preprint2016arXiv

Galaxy And Mass Assembly (GAMA): $\mathcal{M_\star}-R_{\rm e}$ relations of $z=0$ bulges, discs and spheroids

We perform automated bulge + disc decomposition on a sample of $\sim$7500 galaxies from the Galaxy And Mass Assembly (GAMA) survey in the redshift range of 0.002$<$z$<$0.06 using SIGMA, a wrapper around GALFIT3. To achieve robust profile measurements we use a novel approach of repeatedly fitting the galaxies, varying the input parameters to sample a large fraction of the input parameter space. Using this method we reduce the catastrophic failure rate significantly and verify the confidence in the fit independently of $χ^2$. Additionally, using the median of the final fitting values and the 16$^{th}$ and 84$^{th}$ percentile produces more realistic error estimates than those provided by GALFIT, which are known to be underestimated. We use the results of our decompositions to analyse the stellar mass - half-light radius relations of bulges, discs and spheroids. We further investigate the association of components with a parent disc or elliptical relation to provide definite $z=0$ disc and spheroid $\mathcal{M_\star}-R_{\rm e}$ relations. We conclude by comparing our local disc and spheroid $\mathcal{M_\star}-R_{\rm e}$ to simulated data from EAGLE and high redshift data from CANDELS-UDS. We show the potential of using the mass-size relation to study galaxy evolution in both cases but caution that for a fair comparison all data sets need to be processed and analysed in the same manner.

preprint2016arXiv

Galaxy And Mass Assembly (GAMA): Active Galactic Nuclei in Pairs of Galaxies

There exist conflicting observations on whether or not the environment of broad and narrow line AGN differ and this consequently questions the validity of the AGN unification model. The high spectroscopic completeness of the GAMA survey makes it ideal for a comprehensive analysis of the close environment of galaxies. To exploit this, and conduct a comparative analysis of the environment of broad and narrow line AGN within GAMA, we use a double-Gaussian emission line fitting method to model the more complex line profiles associated with broad line AGN. We select 209 type 1 (i.e., unobscured), 464 type 1.5-1.9 (partially obscured), and 281 type 2 (obscured) AGN within the GAMA II database. Comparing the fractions of these with neighbouring galaxies out to a pair separation of $350\,\text{kpc }h^{-1}$ and $Δz < 0.012$ shows no difference between AGN of different type, except at separations less than $20\,\text{kpc }h^{-1}$ where our observations suggest an excess of type 2 AGN in close pairs. We analyse the properties of the galaxies neighbouring our AGN and find no significant differences in colour or the star formation activity of these galaxies. Further to this we find that $Σ_5$ is also consistent between broad and narrow line AGN. We conclude that the observations presented here are consistent with AGN unification.

preprint2016arXiv

Galaxy And Mass Assembly (GAMA): Stellar mass growth of spiral galaxies in the cosmic web

We look for correlated changes in stellar mass and star formation rate along filaments in the cosmic web by examining the stellar masses and UV-derived star formation rates (SFR) of 1,799 ungrouped and unpaired spiral galaxies that reside in filaments. We devise multiple distance metrics to characterise the complex geometry of filaments, and find that galaxies closer to the cylindrical centre of a filament have higher stellar masses than their counterparts near the periphery of filaments, on the edges of voids. In addition, these peripheral spiral galaxies have higher specific star formation rates (SSFR) at a given mass. Complementing our sample of filament spiral galaxies with spiral galaxies in tendrils and voids, we find that the average SFR of these objects in different large scale environments are similar to each other with the primary discriminant in SFR being stellar mass, in line with previous works. However, the distributions of SFRs are found to vary with large-scale environment. Our results thus suggest a model in which in addition to stellar mass as the primary discriminant, the large-scale environment is imprinted in the SFR as a second order effect. Furthermore, our detailed results for filament galaxies suggest a model in which gas accretion from voids onto filaments is primarily in an orthogonal direction. Overall, we find our results to be in line with theoretical expectations of the thermodynamic properties of the intergalactic medium in different large-scale environments.

preprint2016arXiv

Galaxy And Mass Assembly (GAMA): the Stellar Mass Budget of Galaxy Spheroids and Disks

We build on a recent photometric decomposition analysis of 7506 Galaxy and Mass Assembly (GAMA) survey galaxies to derive stellar mass function fits to individual spheroid and disk component populations down to a lower mass limit of log(M_*/M_sun)= 8. We find that the spheroid/disk mass distributions for individual galaxy morphological types are well described by single Schechter function forms. We derive estimates of the total stellar mass densities in spheroids (rho_spheroid = 1.24+/-0.49 * 10^8 M_sun Mpc^-3 h_0.7) and disks (rho_disk = 1.20+/-0.45 * 10^8 M_sun Mpc^-3 h_0.7), which translates to approximately 50% of the local stellar mass density in spheroids and 48% in disks. The remaining stellar mass is found in the dwarf "little blue spheroid" class, which is not obviously similar in structure to either classical spheroid or disk populations. We also examine the variation of component mass ratios across galaxy mass and group halo mass regimes, finding the transition from spheroid to disk mass dominance occurs near galaxy stellar mass ~10^11 M_sun and group halo mass ~10^12.5 M_sun/h. We further quantify the variation in spheroid-to-total mass ratio with group halo mass for central and satellite populations as well as the radial variation of this ratio within groups.

preprint2016arXiv

Galaxy And Mass Assembly (GAMA): Understanding the wavelength dependence of galaxy structure with bulge-disc decompositions

With a large sample of bright, low-redshift galaxies with optical$-$near-IR imaging from the GAMA survey we use bulge-disc decompositions to understand the wavelength-dependent behavior of single-Sérsic structural measurements. We denote the variation in single-Sérsic index with wavelength as $\mathcal{N}$, likewise for effective radius we use $\mathcal{R}$. We find that most galaxies with a substantial disc, even those with no discernable bulge, display a high value of $\mathcal{N}$. The increase in Sérsic index to longer wavelengths is therefore intrinsic to discs, apparently resulting from radial variations in stellar population and/or dust reddening. Similarly, low values of $\mathcal{R}$ ($<$ 1) are found to be ubiquitous, implying an element of universality in galaxy colour gradients. We also study how bulge and disc colour distributions vary with galaxy type. We find that, rather than all bulges being red and all discs being blue in absolute terms, both components become redder for galaxies with redder total colours. We even observe that bulges in bluer galaxies are typically bluer than discs in red galaxies, and that bulges and discs are closer in colour for fainter galaxies. Trends in total colour are therefore not solely due to the colour or flux dominance of the bulge or disc.

preprint2015arXiv

Galaxy And Mass Assembly (GAMA): Panchromatic Data Release (far-UV --- far-IR) and the low-z energy budget

We present the GAMA Panchromatic Data Release (PDR) constituting over 230deg$^2$ of imaging with photometry in 21 bands extending from the far-UV to the far-IR. These data complement our spectroscopic campaign of over 300k galaxies, and are compiled from observations with a variety of facilities including: GALEX, SDSS, VISTA, WISE, and Herschel, with the GAMA regions currently being surveyed by VST and scheduled for observations by ASKAP. These data are processed to a common astrometric solution, from which photometry is derived for 221,373 galaxies with r<19.8 mag. Online tools are provided to access and download data cutouts, or the full mosaics of the GAMA regions in each band. We focus, in particular, on the reduction and analysis of the VISTA VIKING data, and compare to earlier datasets (i.e., 2MASS and UKIDSS) before combining the data and examining its integrity. Having derived the 21-band photometric catalogue we proceed to fit the data using the energy balance code MAGPHYS. These measurements are then used to obtain the first fully empirical measurement of the 0.1-500$μ$m energy output of the Universe. Exploring the Cosmic Spectral Energy Distribution (CSED) across three time-intervals (0.3-1.1Gyr, 1.1-1.8~Gyr and 1.8---2.4~Gyr), we find that the Universe is currently generating $(1.5 \pm 0.3) \times 10^{35}$ h$_{70}$ W Mpc$^{-3}$, down from $(2.5 \pm 0.2) \times 10^{35}$ h$_{70}$ W Mpc$^{-3}$ 2.3~Gyr ago. More importantly, we identify significant and smooth evolution in the integrated photon escape fraction at all wavelengths, with the UV escape fraction increasing from 27(18)% at z=0.18 in NUV(FUV) to 34(23)% at z=0.06. The GAMA PDR will allow for detailed studies of the energy production and outputs of individual systems, sub-populations, and representative galaxy samples at $z<0.5$. The GAMA PDR can be found at: http://gama-psi.icrar.org/

preprint2015arXiv

Galaxy And Mass Assembly (GAMA): the Stellar Mass Budget by Galaxy Type

We report an expanded sample of visual morphological classifications from the Galaxy and Mass Assembly (GAMA) survey phase two, which now includes 7,556 objects (previously 3,727 in phase one). We define a local (z <0.06) sample and classify galaxies into E, S0-Sa, SB0-SBa, Sab-Scd, SBab-SBcd, Sd-Irr, and "little blue spheroid" types. Using these updated classifications, we derive stellar mass function fits to individual galaxy populations divided both by morphological class and more general spheroid- or disk-dominated categories with a lower mass limit of log(Mstar/Msun) = 8 (one dex below earlier morphological mass function determinations). We find that all individual morphological classes and the combined spheroid-/bulge-dominated classes are well described by single Schechter stellar mass function forms. We find that the total stellar mass densities for individual galaxy populations and for the entire galaxy population are bounded within our stellar mass limits and derive an estimated total stellar mass density of rho_star = 2.5 x 10^8 Msun Mpc^-3 h_0.7, which corresponds to an approximately 4% fraction of baryons found in stars. The mass contributions to this total stellar mass density by galaxies that are dominated by spheroidal components (E and S0-Sa classes) and by disk components (Sab-Scd and Sd-Irr classes) are approximately 70% and 30%, respectively.

preprint2015arXiv

Galaxy And Mass Assembly (GAMA): the wavelength dependence of galaxy structure versus redshift and luminosity

We study how the sizes and radial profiles of galaxies vary with wavelength, by fitting Sersic functions simultaneously to imaging in nine optical and near-infrared bands. To quantify the wavelength dependence of effective radius we use the ratio, $\mathcal{R}$, of measurements in two restframe bands. The dependence of Sersic index on wavelength, $\mathcal{N}$, is computed correspondingly. Vulcani et al. (2014) have demonstrated that different galaxy populations present sharply contrasting behaviour in terms of $\mathcal{R}$ and $\mathcal{N}$. Here we study the luminosity dependence of this result. We find that at higher luminosities, early-type galaxies display a more substantial decrease in effective radius with wavelength, whereas late-types present a more pronounced increase in Sersic index. The structural contrast between types thus increases with luminosity. By considering samples at different redshifts, we demonstrate that lower data quality reduces the apparent difference between the main galaxy populations. However, our conclusions remain robust to this effect. We show that accounting for different redshift and luminosity selections partly reconciles the size variation measured by Vulcani et al. with the weaker trends found by other recent studies. Dividing galaxies by visual morphology confirms the behaviour inferred using morphological proxies, although the sample size is greatly reduced. Finally, we demonstrate that varying dust opacity and disc inclination can account for features of the joint distribution of $\mathcal{R}$ and $\mathcal{N}$ for late-type galaxies. However, dust does not appear to explain the highest values of $\mathcal{R}$ and $\mathcal{N}$. The bulge-disc nature of galaxies must also contribute to the wavelength-dependence of their structure.

preprint2015arXiv

The Sizes of Candidate $z\sim9-10$ Galaxies: confirmation of the bright CANDELS sample and relation with luminosity and mass

Recently, a small sample of six $z\sim9-10$ candidates was discovered in CANDELS that are $\sim10-20\times$ more luminous than any of the previous $z\sim9-10$ galaxies identified over the HUDF/XDF and CLASH fields. We measure the sizes of these candidates to map out the size evolution of galaxies from the earliest observable times. Their sizes are also used to provide a valuable constraint on whether these unusual galaxy candidates are at high redshift. Using galfit to derive sizes from the CANDELS F160W images of these candidates, we find a mean size of 0.13$\pm$0.02" (or 0.5$\pm$0.1 kpc at $z\sim9-10$). This handsomely matches the 0.6 kpc size expected extrapolating lower redshift measurements to $z\sim9-10$, while being much smaller than the 0.59" mean size for lower-redshift interlopers to $z\sim9-10$ photometric selections lacking the blue IRAC color criterion. This suggests that source size may be an effective constraint on contaminants from $z\sim9-10$ selections lacking IRAC data. Assuming on the basis of the strong photometric evidence that the Oesch et al. 2014 sample is entirely at $z\sim9-10$, we can use this sample to extend current constraints on the size-luminosity, size-mass relation, and size evolution of galaxies to $z\sim10$. We find that the $z\sim9-10$ candidate galaxies have broadly similar sizes and luminosities as $z\sim6$-8 counterparts with star-formation-rate surface densities in the range of $\rm Σ_{SFR}=1-20\, M_\odot~ yr^{-1}\, kpc^{-2}$. The stellar mass-size relation is uncertain, but shallower than those inferred for lower-redshift galaxies. In combination with previous size measurements at z=4-7, we find a size evolution of $(1+z)^{-m}$ with $m=1.0\pm0.1$ for $>0.3L^*_{z=3}$ galaxies, consistent with the evolution previously derived from $2 < z < 8$ galaxies.

preprint2014arXiv

Evidence for the concurrent growth of thick discs and central mass concentrations from S$^4$G imaging

We have produced $3.6μ{\rm m}+4.5μ{\rm m}$ vertically integrated radial luminosity profiles of 69 edge-on galaxies from the Spitzer Survey of Stellar Structure in Galaxies (S$^4$G). We decomposed the luminosity profiles into a disc and a central mass concentration (CMC). These fits, combined with thin/thick disc decompositions from our previous studies, allow us to estimate the masses of the CMCs, the thick discs, and the thin discs ($\mathcal{M}_{\rm CMC}$, $\mathcal{M}_{\rm T}$, and $\mathcal{M}_{\rm T}$). We obtained atomic disc masses ($\mathcal{M}_{\rm g}$) from the literature. We then consider the CMC and the thick disc to be dynamically hot components and the thin disc and the gas disc to be dynamically cold components. We find that the ratio between the mass of the hot components and that of the cold components, $(\mathcal{M}_{\rm CMC}+\mathcal{M}_{\rm T})/(\mathcal{M}_{\rm t}+\mathcal{M}_{\rm g})$, does not depend on the total galaxy mass as described by circular velocities ($v_{\rm c}$). We also find that the higher the $v_{\rm c}$, the more concentrated the hot component of a galaxy. We suggest that our results are compatible with having CMCs and thick discs built in a short and early high star forming intensity phase. These components were born thick because of the large scale height of the turbulent gas disc in which they originated. Our results indicate that the ratio between the star forming rate in the former phase and that of the formation of the thin disc is of the order of 10, but the value depends on the duration of the high star forming intensity phase.

preprint2014arXiv

Galaxy And Mass Assembly (GAMA): Stellar mass functions by Hubble type

We present an estimate of the galaxy stellar mass function and its division by morphological type in the local (0.025 < z < 0.06) Universe. Adopting robust morphological classifications as previously presented (Kelvin et al.) for a sample of 3,727 galaxies taken from the Galaxy And Mass Assembly survey, we define a local volume and stellar mass limited sub-sample of 2,711 galaxies to a lower stellar mass limit of M = 10^9.0 M_sun. We confirm that the galaxy stellar mass function is well described by a double Schechter function given by M* = 10^10.64 M_sun, α1 = -0.43, ϕ*1 = 4.18 dex^-1 Mpc^-3, α2 = -1.50 and ϕ*2 = 0.74 dex^-1 Mpc^-3. The constituent morphological-type stellar mass functions are well sampled above our lower stellar mass limit, excepting the faint little blue spheroid population of galaxies. We find approximately 71+3-4% of the stellar mass in the local Universe is found within spheroid dominated galaxies; ellipticals and S0-Sas. The remaining 29+4-3% falls predominantly within late type disk dominated systems, Sab-Scds and Sd-Irrs. Adopting reasonable bulge-to-total ratios implies that approximately half the stellar mass today resides in spheroidal structures, and half in disk structures. Within this local sample, we find approximate stellar mass proportions for E : S0-Sa : Sab-Scd : Sd-Irr of 34 : 37 : 24 : 5.

preprint2014arXiv

IRAC Near-Infrared Features in the Outer Parts of S4G Galaxies

We present a catalogue and images of visually detected features, such as asymmetries, extensions, warps, shells, tidal tails, polar rings, and obvious signs of mergers or interactions, in the faint outer regions (at and outside of R_25) of nearby galaxies. This catalogue can be used in future quantitative studies that examine galaxy evolution due to internal and external factors. We are able to reliably detect outer region features down to a brightness level of 0.03 MJy/sr per pixel at 3.6 microns in the Spitzer Survey of Stellar Structure in Galaxies (S4G). We also tabulate companion galaxies. We find asymmetries in the outer isophotes in 22+/-1 per cent of the sample. The asymmetry fraction does not correlate with galaxy classification as an interacting galaxy or merger remnant, or with the presence of companions. We also compare the detected features to similar features in galaxies taken from cosmological zoom re-simulations. The simulated images have a higher fraction (33 per cent) of outer disc asymmetries, which may be due to selection effects and an uncertain star formation threshold in the models. The asymmetries may have either an internal (e.g., lopsidedness due to dark halo asymmetry) or external origin.

preprint2014arXiv

The Luminosity Function at z~8 from 97 Y-band dropouts: Inferences About Reionization

[Abbreviated] We present the largest search to date for $z\sim8$ Lyman break galaxies (LBGs) based on 350 arcmin$^2$ of HST observations in the V-, Y-, J- and H-bands from the Brightest of Reionizing Galaxies (BoRG) survey. The BoRG dataset includes $\sim$50 arcmin$^2$ of new data and deeper observations of two previous BoRG pointings, from which we present 9 new $z\sim8$ LBG candidates, bringing the total number of BoRG LBGs to 38 with $25.5\leqslant m_{J} \leqslant 27.6$ (AB system). We introduce a new Bayesian formalism for estimating the galaxy luminosity function (LF), which does not require binning (and thus smearing) of the data and includes a likelihood based on the formally correct binomial distribution as opposed to the often used approximate Poisson distribution. We demonstrate the utility of the new method on a sample of $97$ LBGs that combines the bright BoRG galaxies with the fainter sources published in Bouwens et al. (2012) from the HUDF and ERS programs. We show that the $z\sim8$ LF is well described by a Schechter function with a characteristic magnitude $M^\star = -20.15^{+0.29}_{-0.38}$, a faint-end slope of $α= -1.87^{+0.26}_{-0.26}$, and a number density of $\log_{10} ϕ^\star [\textrm{Mpc}^{-3}] = -3.24^{+0.25}_{-0.24}$. Integrated down to $M=-17.7$ this LF yields a luminosity density, $\log_{10} ε[\textrm{erg}/\textrm{s/Hz/Mpc}^{3}] = 25.52^{+0.05}_{-0.05}$. Our LF analysis is consistent with previously published determinations within 1$σ$. We discuss the implication of our study for the physics of reionization. By assuming theoretically motivated priors on the clumping factor and the photon escape fraction we show that the UV LF from galaxy samples down to $M=-17.7$ can ionize only 10-50% of the neutral hydrogen at $z\sim8$. Full reionization would require extending the LF down to $M=-15$.

preprint2014arXiv

The Ultraviolet Attenuation Law in Backlit Spiral Galaxies

(Abridged) The effective extinction law (attenuation behavior) in galaxies in the emitted ultraviolet is well known only for actively star-forming objects and combines effects of the grain properties, fine structure in the dust distribution, and relative distributions of stars and dust. We use GALEX, XMM Optical Monitor, and HST data to explore the UV attenuation in the outer parts of spiral disks which are backlit by other UV-bright galaxies, starting with candidates provided by Galaxy Zoo participants. Our analysis incorporates galaxy symmetry, using non-overlapping regions of each galaxy to derive error estimates on the attenuation measurements. The entire sample has an attenuation law close to the Calzetti et al. (1994) form; the UV slope for the overall sample is substantially shallower than found by Wild et al. (2011), a reasonable match to the more distant galaxies in our sample but not to the weighted combination including NGC 2207. The nearby, bright spiral NGC 2207 alone gives accuracy almost equal to the rest of our sample, and its outer arms have a very low level of foreground starlight. This "grey" law can be produced from the distribution of dust alone, without a necessary contribution from differential escape of stars from dense clouds. The extrapolation needed to compare attenution between backlit galaxies at moderate redshifts, and local systems from SDSS data, is mild enough to allow use of galaxy overlaps to trace the cosmic history of dust. For NGC 2207, the covering factor of clouds with small optical attenuation becomes a dominant factor farther into the ultraviolet, which opens the possibility that widespread diffuse dust dominates over dust in star-forming regions deep into the ultraviolet. Comparison with published radiative-transfer models indicates that the role of dust clumping dominates over differences in grain populations, at this spatial resolution.

preprint2013arXiv

Constraining the age of the NGC 4565 HI Disk Warp: Determining the Origin of Gas Warps

We have mapped the distribution of young and old stars in the gaseous HI warp of NGC 4565. We find a clear correlation of young stars (<600 Myr) with the warp, but no coincident old stars (>1 Gyr), which places an upper limit on the age of the structure. The formation rate of the young stars, which increased ~300 Myr ago relative to the surrounding regions, is (6.3 +2.5/-1.5) x 10^-5 M_sol/yr/kpc^2. This implies a ~60+/-20 Gyr depletion time of the HI warp, similar to the timescales calculated for the outer HI disks of nearby spiral galaxies. While some stars associated with the warp fall into the asymptotic giant branch (AGB) region of the color magnitude diagram, where stars could be as old as 1 Gyr, further investigation suggests that they may be interlopers rather than real AGB stars. We discuss the implications of these age constraints for the formation of HI warps, and the gas fueling of disk galaxies.

preprint2013arXiv

Emission-Line Galaxies from the Hubble Space Telescope Probing Evolution and Reionization Spectroscopically (PEARS) Grism Survey. II: The Complete Sample

We present a full analysis of the Probing Evolution And Reionization Spectroscopically (PEARS) slitess grism spectroscopic data obtained with the Advanced Camera for Surveys on HST. PEARS covers fields within both the Great Observatories Origins Deep Survey (GOODS) North and South fields, making it ideal as a random survey of galaxies, as well as the availability of a wide variety of ancillary observations to support the spectroscopic results. Using the PEARS data we are able to identify star forming galaxies within the redshift volume 0< z<1.5. Star forming regions in the PEARS survey are pinpointed independently of the host galaxy. This method allows us to detect the presence of multiple emission line regions (ELRs) within a single galaxy. 1162 Ha, [OIII] and/or [OII] emission lines have been identified in the PEARS sample of ~906 galaxies down to a limiting flux of ~1e-18 erg/s/cm^2. The ELRs have also been compared to the properties of the host galaxy, including morphology, luminosity, and mass. From this analysis we find three key results: 1) The computed line luminosities show evidence of a flattening in the luminosity function with increasing redshift; 2) The star forming systems show evidence of disturbed morphologies, with star formation occurring predominantly within one effective (half-light) radius. However, the morphologies show no correlation with host stellar mass; and 3) The number density of star forming galaxies with M_* > 1e9} M_sun decreases by an order of magnitude at z<0.5 relative to the number at 0.5<z<0.9 in support of the argument for galaxy downsizing.

preprint2013arXiv

GHOSTS I: A New Faint very Isolated Dwarf Galaxy at D = 12 +/- 2 Mpc

We report the discovery of a new faint dwarf galaxy, GHOSTS I, using HST/ACS data from one of our GHOSTS (Galaxy Halos, Outer disks, Substructure, Thick disk, and Star clusters) fields. Its detected individual stars populate an approximately one magnitude range of its luminosity function (LF). Using synthetic color-magnitude diagrams (CMDs) to compare with the galaxy's CMD, we find that the colors and magnitudes of GHOSTS I's individual stars are most consistent with being young helium-burning and asymptotic giant branch stars at a distance of 12 +/- 2 Mpc. Morphologically, GHOSTS I appears to be actively forming stars, so we tentatively classify it as a dwarf irregular (dIrr) galaxy, although future HST observations deep enough to resolve a larger magnitude range in its LF are required to make a more secure classification. GHOSTS I's absolute magnitude is $M_V = -9.85^{+ 0.40}_{- 0.33}$, making it one of the least luminous dIrr galaxies known, and its metallicity is lower than [Fe/H] =-1.5 dex. The half-light radius of GHOSTS I is 226 +/- 38 pc and its ellipticity is 0.47 +/- 0.07, similar to Milky Way and M31 dwarf satellites at comparable luminosity. There are no luminous massive galaxies or galaxy clusters within ~ 4 Mpc from GHOSTS I that could be considered as its host, making it a very isolated dwarf galaxy in the Local Universe.

preprint2013arXiv

On the Origin of Lopsidedness in Galaxies as Determined from the Spitzer Survey of Stellar Structure in Galaxies (S4G)

We study the m=1 distortions (lopsidedness) in the stellar components of 167 nearby galaxies that span a wide range of morphologies and luminosities. We confirm the previous findings of 1) a high incidence of lopsidedness in the stellar distributions, 2) increasing lopsidedness as a function of radius out to at least 3.5 exponential scale lengths, and 3) greater lopsidedness, over these radii, for galaxies of later type and lower surface brightness. Additionally, the magnitude of the lopsidedness 1) correlates with the character of the spiral arms (stronger arm patterns occur in galaxies with less lopsidedness), 2) is not correlated with the presence or absence of a bar, or the strength of the bar when one is present, 3) is inversely correlated to the stellar mass fraction, f_*, within one radial scale length, and 4) correlates directly with f_* measured within the radial range over which we measure lopsidedness. We interpret these findings to mean that lopsidedness is a generic feature of galaxies and does not, generally, depend on a rare event, such as a direct accretion of a satellite galaxy onto the disk of the parent galaxy. While lopsidedness may be caused by several phenomena, moderate lopsidedness (<A_1>_i + <A_1>_o)/2 < 0.3) is likely to reflect halo asymmetries to which the disk responds or a gravitationally self-generated mode . We hypothesize that the magnitude of the stellar response depends both on how centrally concentrated the stars are with respect to the dark matter and whether there are enough stars in the region of the lopsidedness that self-gravity is dynamically important.

preprint2013arXiv

Testing galaxy formation models with the GHOSTS survey: The color profile of M81's stellar halo

We study the properties of the stellar populations in M81's outermost part, which hereafter we will term the stellar halo, using HST ACS/WFC observations of 19 fields from the GHOSTS survey. The observed fields probe the stellar halo out to a projected distance of ~ 50 kpc from the galactic center. Each field was observed in both F606W and F814W filters. The 50% completeness levels of the color magnitude diagrams (CMDs) are typically at 2 mag below the tip of the red giant branch (TRGB). Fields at distances closer than 15 kpc show evidence of disk-dominated populations whereas fields at larger distances are mostly populated by halo stars. The RGB of the M81's halo CMDs is well matched with isochrones of ~ 10 Gyr and metallicities [Fe/H] ~ -1.2 dex, suggesting that the dominant stellar population of M81's halo has a similar age and metallicity. The halo of M81 is characterized by a color distribution of width ~ 0.4 mag and an approximately constant median value of (F606W - F814W) ~ 1 mag measured using stars within the magnitude range 23.7 < F814W < 25.5. When considering only fields located at galactocentric radius R > 15 kpc, we detect no color gradient in the stellar halo of M81. We place a limit of 0.03+/-0.11 mag difference between the median color of RGB M81 halo stars at ~ 15 and at 50 kpc, corresponding to a metallicity difference of 0.08+/-0.35 dex over that radial range for an assumed constant age of 10 Gyr. We compare these results with model predictions for the colors of stellar halos formed purely via accretion of satellite galaxies. When we analyze the cosmologically motivated models in the same way as the HST data, we find that they predict no color gradient for the stellar halos, in good agreement with the observations.

preprint2013arXiv

X-ray nuclear activity in S4G barred galaxies: No link between bar strength and co-occurrent supermassive black hole fueling

Stellar bars can lead to gas inflow toward the center of a galaxy and stimulate nuclear star formation. However, there is no compelling evidence on whether they also feed a central supermassive black hole: by measuring the fractions of barred active and inactive galaxies, previous studies have yielded conflicting results. In this paper, we aim to understand the lack of observational evidence for bar-driven active galactic nucleus (AGN) activity by studying a sample of 41 nearby (d < 35 Mpc) barred galaxies from the Spitzer Survey for Stellar Structure in Galaxies. We use Chandra observations to measure nuclear 2--10 keV X-ray luminosities and estimate Eddington ratios, together with Spitzer 3.6um imaging to quantify the strength of the stellar bar in two independent ways: (1) from its structure, as traced by its ellipticity and boxiness, and (2) from its gravitational torque Q_b, taken as the maximum ratio of the tangential force to the mean background radial force. In this way, rather than discretizing the presence of both stellar bars and nuclear activity, we are able to account for the continuum of bar strengths and degrees of AGN activity. We find nuclear X-ray sources in 31 out of 41 galaxies with median X-ray luminosity and Eddington ratio of L_X=4.3x10^{38} erg/s and L_bol/L_Edd=6.9x10^{-6} respectively, consistent with low-luminosity AGN activity. Including upper limits for those galaxies without nuclear detections, we find no significant correlation between any of the bar strength indicators and the degree of nuclear activity, irrespective of galaxy luminosity, stellar mass, Hubble type, or bulge size. Strong bars do not favor brighter or more efficient nuclear activity, implying that at least for the low-luminosity regime, supermassive black hole fueling is not closely connected to large scale features.

preprint2012arXiv

A unified picture of breaks and truncations in spiral galaxies from SDSS and S^{4}G imaging

The mechanism causing breaks in the radial surface brightness distribution of spiral galaxies is not yet well known. Despite theoretical efforts, there is not a unique explanation for these features and the observational results are not conclusive. In an attempt to address this problem, we have selected a sample of 34 highly inclined spiral galaxies present both in the Sloan Digital Sky Survey and in the Spitzer Survey of Stellar Structure in Galaxies. We have measured the surface brightness profiles in the five Sloan optical bands and in the 3.6$μm$ Spitzer band. We have also calculated the color and stellar surface mass density profiles using the available photometric information, finding two differentiated features: an innermost break radius at distances of $\sim 8 \pm 1$ kpc [$0.77 \pm 0.06$ $R_{25}$] and a second characteristic radius, or truncation radius, close to the outermost optical extent ($\sim 14 \pm 2$ kpc [$1.09 \pm 0.05$ $R_{25}$]) of the galaxy. We propose in this work that the breaks might be a phenomena related to a threshold in the star formation, while truncations are more likely a real drop in the stellar mass density of the disk associated with the maximum angular momentum of the stars.

preprint2012arXiv

Galaxy Zoo: A Catalog of Overlapping Galaxy Pairs for Dust Studies

Analysis of galaxies with overlapping images offers a direct way to probe the distribution of dust extinction and its effects on the background light. We present a catalog of 1990 such galaxy pairs selected from the Sloan Digital Sky Survey (SDSS) by volunteers of the Galaxy Zoo project. We highlight subsamples which are particularly useful for retrieving such properties of the dust distribution as UV extinction, the extent perpendicular to the disk plane, and extinction in the inner parts of disks. The sample spans wide ranges of morphology and surface brightness, opening up the possibility of using this technique to address systematic changes in dust extinction or distribution with galaxy type. This sample will form the basis for forthcoming work on the ranges of dust distributions in local disk galaxies, both for their astrophysical implications and as the low-redshift part of a study of the evolution of dust properties. Separate lists and figures show deep overlaps, where the inner regions of the foreground galaxy are backlit, and the relatively small number of previously-known overlapping pairs outside the SDSS DR7 sky coverage.

preprint2012arXiv

Outer-Disk Populations in NGC 7793: Evidence for Stellar Radial Migration

We analyzed the radial surface brightness profile of the spiral galaxy NGC 7793 using HST/ACS images from the GHOSTS survey and a new HST/WFC3 image across the disk break. We used the photometry of resolved stars to select distinct populations covering a wide range of stellar ages. We found breaks in the radial profiles of all stellar populations at 280" (~5.1 kpc). Beyond this disk break, the profiles become steeper for younger populations. This same trend is seen in numerical simulations where the outer disk is formed almost entirely by radial migration. We also found that the older stars of NGC 7793 extend significantly farther than the underlying HI disk. They are thus unlikely to have formed entirely at their current radii, unless the gas disk was substantially larger in the past. These observations thus provide evidence for substantial stellar radial migration in late-type disks.

preprint2011arXiv

Thick disks of edge-on galaxies seen through the Spitzer Survey of Stellar Structure in Galaxies (S4G): Lair of missing baryons?

Most, if not all, disk galaxies have a thin (classical) disk and a thick disk. In most models thick disks are thought to be a necessary consequence of the disk formation and/or evolution of the galaxy. We present the results of a study of the thick disk properties in a sample of carefully selected edge-on galaxies with types ranging from T=3 to T=8. We fitted one-dimensional luminosity profiles with physically motivated functions - the solutions of two stellar and one gaseous isothermal coupled disks in equilibrium - which are likely to yield more accurate results than other functions used in previous studies. The images used for the fits come from the Spitzer Survey of Stellar Structure in Galaxies (S4G). We found that thick disks are on average more massive than previously reported, mostly due to the selected fitting function. Typically, the thin and the thick disk have similar masses. We also found that thick disks do not flare significantly within the observed range in galactocentric radii and that the ratio of thick to thin disk scaleheights is higher for galaxies of earlier types. Our results tend to favor an in situ origin for most of the stars in the thick disk. In addition the thick disk may contain a significant amount of stars coming from satellites accreted after the initial build-up of the galaxy and an extra fraction of stars coming from the secular heating of the thin disk by its own overdensities. Assigning thick disk light to the thin disk component may lead to an underestimate of the overall stellar mass in galaxies, because of different mass to light ratios in the two disk components. On the basis of our new results, we estimate that disk stellar masses are between 10% and 50% higher than previously thought and we suggest that thick disks are a reservoir of "local missing baryons".

preprint2010arXiv

The Herschel-SPIRE Legacy Survey (HSLS): the scientific goals of a shallow and wide submillimeter imaging survey with SPIRE

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

preprint2010arXiv

The thick disk in the galaxy NGC 4244 from S4G imaging

If thick disks are ubiquitous and a natural product of disk galaxy formation and/or evolution processes, all undisturbed galaxies which have evolved during a significant fraction of a Hubble time should have a thick disk. The late-type spiral galaxy NGC 4244 has been reported as the only nearby edge-on galaxy without a confirmed thick disk. Using data from the Spitzer Survey of Stellar Structure in Galaxies (S4G) we have identified signs of two disk components in this galaxy. The asymmetries between the light profiles on both sides of the mid-plane of NGC 4244 can be explained by a combination of the galaxy not being perfectly edge-on and a certain degree of opacity of the thin disk. We argue that the subtlety of the thick disk is a consequence of either a limited secular evolution in NGC 4244, a small fraction of stellar material in the fragments which built the galaxy, or a high amount of gaseous accretion after the formation of the galaxy.