Source author record

B. W. Holwerda

B. 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

55works
3topics
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

55 published item(s)

preprint2025arXiv

Pushing JWST to the extremes: search and scrutiny of bright galaxy candidates at z$\simeq$15-30

We designed customized Lyman-break color selection techniques to identify galaxy candidates in the redshift ranges $15 \leq z \leq 20$ and $20 \leq z \leq 28$. The selection was performed on the ASTRODEEP-JWST multi-band catalogs of the CEERS, Abell-2744, JADES, NGDEEP, and PRIMER survey fields, covering a total area of $\sim0.2$ sq. deg. We identify five candidates at $15 \leq z \leq 20$, while no objects are found based on the $z\gtrsim20$ color selection criteria. Despite exhibiting a $>$1.5 mag break, all the objects display multimodal redshift probability distributions across different SED-fitting codes and methodologies. The alternative solutions correspond to poorly understood populations of low-mass quiescent or dusty galaxies at z$\sim$3-7. This conclusion is supported by the analysis of five F200W-dropout objects that we find to be interlopers on the basis of NIRSpec PRISM spectra: four dusty star-forming galaxies at z$\sim$2.2-6.6, and a passive galaxy at z=4.91 with log$(M_{\rm star}/{\rm M}_{\odot}) \lesssim$ 9. We measured the UV luminosity function under different assumptions on the contamination level within our sample. We find that if even a fraction of the candidates is indeed at $z\gtrsim15$, the resulting UV LF points to a very mild evolution compared to estimates at $z<15$, implying a significant tension with existing theoretical models. In particular, confirming our bright ($M_{\text{UV}}<-21$) candidates would require substantial revisions to the theoretical framework. In turn, if all these candidates will be confirmed to be interlopers, we conclude that future surveys may need ten times wider areas to select $M_{\text{UV}}\lesssim-20$ galaxies at $z>15$. Observations in the F150W and F200W filters at depths comparable to those in the NIRCam LW bands are also required to mitigate contamination from rare red objects at z$\lesssim$8.

preprint2022arXiv

Galaxy And Mass Assembly (GAMA): Self-Organizing Map Application on Nearby Galaxies

Galaxy populations show bimodality in a variety of properties: stellar mass, colour, specific star-formation rate, size, and Sérsic index. These parameters are our feature space. We use an existing sample of 7556 galaxies from the Galaxy and Mass Assembly (GAMA) survey, represented using five features and the K-means clustering technique, showed that the bimodalities are the manifestation of a more complex population structure, represented by between 2 and 6 clusters. Here we use Self Organizing Maps (SOM), an unsupervised learning technique which can be used to visualize similarity in a higher dimensional space using a 2D representation, to map these five-dimensional clusters in the feature space onto two-dimensional projections. To further analyze these clusters, using the SOM information, we agree with previous results that the sub-populations found in the feature space can be reasonably mapped onto three or five clusters. We explore where the "green valley" galaxies are mapped onto the SOM, indicating multiple interstitial populations within the green valley population. Finally, we use the projection of the SOM to verify whether morphological information provided by GalaxyZoo users, for example, if features are visible, can be mapped onto the SOM-generated map. Voting on whether galaxies are smooth, likely ellipticals, or "featured" can reasonably be separated but smaller morphological features (bar, spiral arms) can not. SOMs promise to be a useful tool to map and identify instructive sub-populations in multidimensional galaxy survey feature space, provided they are large enough.

preprint2022arXiv

Galaxy And Mass Assembly: Galaxy Zoo spiral arms and star formation rates

Understanding the effect spiral structure has on star formation properties of galaxies is important to completing our picture of spiral structure evolution. Previous studies have investigated connections between spiral arm properties with star formation, but the effect that the number of spiral arms has on this process is unclear. Here we use the Galaxy and Mass Assembly (GAMA) survey paired with the citizen science visual classifications from the Galaxy Zoo project to explore galaxies' spiral arm number and how it connects to the star formation process. We use the votes from the GAMA-KiDS GalaxyZoo classification to investigate the link between spiral arm number with stellar mass, star formation rate, and specific star formation rate. We find that galaxies with fewer spiral arms have lower stellar masses and higher sSFRs, while those with more spiral arms tend toward higher stellar masses and lower sSFRs, and conclude that galaxies are less efficient at forming stars if they have more spiral arms. We note how previous studies' findings may indicate a cause for this connection in spiral arm strength or opacity.

preprint2022arXiv

North Ecliptic Pole merging galaxy catalogue

We aim to generate a catalogue of merging galaxies within the 5.4 sq. deg. North Ecliptic Pole over the redshift range $0.0 < z < 0.3$. To do this, imaging data from the Hyper Suprime-Cam are used along with morphological parameters derived from these same data. The catalogue was generated using a hybrid approach. Two neural networks were trained to perform binary merger non-merger classifications: one for galaxies with $z < 0.15$ and another for $0.15 \leq z < 0.30$. Each network used the image and morphological parameters of a galaxy as input. The galaxies that were identified as merger candidates by the network were then visually checked by experts. The resulting mergers will be used to calculate the merger fraction as a function of redshift and compared with literature results. We found that 86.3% of galaxy mergers at $z < 0.15$ and 79.0% of mergers at $0.15 \leq z < 0.30$ are expected to be correctly identified by the networks. Of the 34 264 galaxies classified by the neural networks, 10 195 were found to be merger candidates. Of these, 2109 were visually identified to be merging galaxies. We find that the merger fraction increases with redshift, consistent with literature results from observations and simulations, and that there is a mild star-formation rate enhancement in the merger population of a factor of $1.102 \pm 0.084$.

preprint2022arXiv

Stellar masses, sizes, and radial profiles for 465 nearby early-type galaxies: an extension to the Spitzer Survey of Stellar Structure in Galaxies (S$^{4}$G)

The Spitzer Survey of Stellar Structure in Galaxies (S$^{4}$G) is a detailed study of over 2300 nearby galaxies in the near-infrared (NIR), which has been critical to our understanding of the detailed structures of nearby galaxies. Because the sample galaxies were selected only using radio-derived velocities, however, the survey favored late-type disk galaxies over lenticulars and ellipticals. A follow-up Spitzer survey was conducted to rectify this bias, adding 465 early-type galaxies (ETGs) to the original sample, to be analyzed in a manner consistent with the initial survey. We present the data release of this ETG extension, up to the third data processing pipeline (P3): surface photometry. We produce curves of growth and radial surface brightness profiles (with and without inclination corrections) using reduced and masked Spitzer IRAC 3.6$μ$m and 4.5$μ$m images produced through Pipelines 1 and 2, respectively. From these profiles, we derive the following integrated quantities: total magnitudes, stellar masses, concentration parameters, and galaxy size metrics. We showcase NIR scaling relations for ETGs among these quantities. We examine general trends across the whole S$^{4}$G and ETG extension among our derived parameters, highlighting differences between ETGs and late-type galaxies (LTGs). ETGs are, on average, more massive and more concentrated than LTGs, and also show subtle distinctions among ETG morphological sub-types. We also derive the following scaling relations and compare with previous results in visible light: mass--size (both half-light and isophotal), mass--concentration, mass--surface brightness (central, effective, and within 1 kpc), and mass--color. We find good agreement with previous works, though some relations (e.g., mass--central surface brightness) will require more careful multi-component decompositions to be fully understood.

preprint2022arXiv

The XXL Survey XLV. Linking the ages of optically selected groups to their X-ray emission

We investigate the properties of 232 optical spectroscopically selected groups from the Galaxy And Mass Assembly (GAMA) survey that overlap the XXL X-ray cluster survey. X-ray aperture flux measurements combined with GAMA group data provides the largest available sample of optical groups with detailed galaxy membership information and consistently measured X-ray fluxes and upper limits. 142 of these groups are divided into three subsets based on the relative strength of X-ray and optical emission, and we see a trend in galaxy properties between these subsets: X-ray overluminous groups contain a lower fraction of both blue and star forming galaxies compared with X-ray underluminous systems. X-ray overluminous groups also have a more dominant central galaxy, with a magnitude gap between first and second ranked galaxies on average 0.22 mag larger than in underluminous groups. The central galaxy in overluminous groups also lies closer to the centre of the group. We examine a number of other structural properties of our groups, such as axis ratio, velocity dispersion, and group crossing time and find trends with X-ray emission in some of these properties despite the high stochastic noise from the limited number of group galaxies. We attribute the trends we see to the evolutionary state of groups, with X-ray overluminous systems being more dynamically evolved than underluminous groups. The X-ray overluminous groups have had more time to develop a luminous intragroup medium, quench member galaxies, and build the mass of the central galaxy through mergers compared to underluminous groups. However, a minority of X-ray underluminous groups have properties that suggest them to be dynamically mature. The lack of hot gas in these systems cannot be accounted for by high star formation efficiency, suggesting that high gas entropy resulting from feedback is the likely cause of their weak X-ray emission.

preprint2021arXiv

Galaxy And Mass Assembly: Group and field galaxy morphologies in the star-formation rate - stellar mass plane

We study the environment in which a galaxy lies (i.e. field or group) and its connection with the morphology of the galaxy. This is done by examining the distribution of parametric and non-parametric statistics across the star-formation rate (SFR) - stellar mass (M$_{\star}$) plane and studying how these distributions change with the environment in the local universe ($z<0.15$). We determine the concentration (C), Gini, M$_{20}$, asymmetry, Gini-M$_{20}$ bulge statistic (GMB), 50\% light radius ($r_{50}$), total Sérsic index, and bulge Sérsic index ($n_{Bulge}$) for galaxies from the Galaxy and Mass Assembly (GAMA) survey using optical images from the Kilo Degree Survey. We determine the galaxy environment using the GAMA group catalogue and split the galaxies into field or group galaxies. The group galaxies are further divided by the group halo mass (M$_{h}$) - $11\leq\mathrm{log(M}_{h}/\mathrm{M}_\odot)<12$, $12\leq\mathrm{log(M}_{h}/\mathrm{M}_\odot)<13$, and $13\leq\mathrm{log(M}_{h}/\mathrm{M}_\odot)<14$ - and into central and satellite galaxies. The galaxies in each of these samples are then placed onto the SFR-M$_{\star}$ plane, and each parameter is used as a third dimension. We fit the resulting distributions for each parameter in each sample using two two-dimensional Gaussian distributions: one for star-forming galaxies and one for quiescent galaxies. The coefficients of these Gaussian fits are then compared between environments. Using C and $r_{50}$, we find that galaxies typically become larger as the group mass increases. This change is greater for larger galaxies. There is no indication that galaxies are typically more or less clumpy as the environment changes. Using GMB and $n_{Bulge}$, we see that the star-forming galaxies do not become more bulge or disk dominated as the group mass changes. Asymmetry does not appear to be greatly influenced by environment.

preprint2021arXiv

GAMA/DEVILS: Constraining the cosmic star-formation history from improved measurements of the 0.3-2.2 micron Extragalactic Background Light

We present a revised measurement of the optical extragalactic background light (EBL), based on the contribution of resolved galaxies to the integrated galaxy light (IGL). The cosmic optical background radiation (COB), encodes the light generated by star-formation, and provides a wealth of information about the cosmic star formation history (CSFH). We combine wide and deep galaxy number counts from the Galaxy And Mass Assembly survey (GAMA) and Deep Extragalactic VIsible Legacy Survey (DEVILS), along with the Hubble Space Telescope (HST) archive and other deep survey datasets, in 9 multi-wavelength filters to measure the COB in the range from 0.35 micron to 2.2 micron. We derive the luminosity density in each band independently and show good agreement with recent and complementary estimates of the optical-EBL from very high-energy (VHE) experiments. Our error analysis suggests that the IGL and Gamma-ray measurements are now fully consistent to within ~10%, suggesting little need for any additional source of diffuse light beyond the known galaxy population. We use our revised IGL measurements to constrain the cosmic star-formation history, and place amplitude constraints on a number of recent estimates. As a consistency check, we can now demonstrate convincingly, that the CSFH, stellar mass growth, and the optical-EBL provide a fully consistent picture of galaxy evolution. We conclude that the peak of star-formation rate lies in the range 0.066-0.076 Msol/yr/Mpc^3 at a lookback time of 9.1 to 10.9 Gyrs.

preprint2021arXiv

The detection of a massive chain of dark HI clouds in the GAMA G23 Field

We report on the detection of a large, extended HI cloud complex in the GAMA G23 field, located at a redshift of $z\,\sim\,0.03$, observed as part of the MeerHOGS campaign (a pilot survey to explore the mosaicing capabilities of MeerKAT). The cloud complex, with a total mass of $10^{10.0}\,M_\odot$, lies in proximity to a large galaxy group with $M_\mathrm{dyn}\sim10^{13.5}\,M_\odot$. We identify seven HI peak concentrations, interconnected as a tenuous 'chain' structure, extending $\sim 400\,\mathrm{kpc}$ from east-to-west, with the largest (central) concentration containing $10{^{9.7}}\,M_\odot$ in HI gas distributed across $50\,\mathrm{kpc}$. The main source is not detected in ultra-violet, optical or infrared imaging. The implied gas mass-to-light ($M_\mathrm{HI}$/$L_\mathrm{r}$) is extreme ($>$1000) even in comparison to other 'dark clouds'. The complex has very little kinematic structure ($110\,\mathrm{km}\,\mathrm{s}^{-1}$), making it difficult to identify cloud rotation. Assuming pressure support, the total mass of the central concentration is $>10^{10.2}\,M_\odot$, while a lower limit to the dynamical mass in the case of full rotational support is $10^{10.4}\,M_\odot$. If the central concentration is a stable structure, it has to contain some amount of unseen matter, but potentially less than is observed for a typical galaxy. It is, however, not clear whether the structure has any gravitationally stable concentrations. We report a faint UV--optical--infrared source in proximity to one of the smaller concentrations in the gas complex, leading to a possible stellar association. The system nature and origins is enigmatic, potentially being the result of an interaction with or within the galaxy group it appears to be associated with.

preprint2020arXiv

Galaxy And Mass Assembly (GAMA): Defining Passive Galaxy Samples and Searching for the UV Upturn

We use data from the GAMA and GALEX surveys to demonstrate that the UV upturn, an unexpected excess of ultraviolet flux from a hot stellar component, seen in the spectra of many early-type galaxies, arises from processes internal to individual galaxies with no measurable influence from the galaxies' larger environment. We first define a clean sample of passive galaxies without a significant contribution to their UV flux from low-level star formation. We confirm that galaxies with the optical colours of red sequence galaxies often have signs of residual star formation, which, without other information, would prevent a convincing demonstration of the presence of UV upturns. However, by including (NUV$-u$) and {\it WISE} (W2-W3) colours, and FUV data where it exists, we can convincingly constrain samples to be composed of non-star-forming objects. Using such a sample, we examine GALEX photometry of low redshift GAMA galaxies in a range of low-density environments, from groups to the general field, searching for UV upturns. We find a wide range of (NUV$-r$) colours, entirely consistent with the range seen -- and attributed to the UV upturn -- in low-redshift red sequence cluster galaxies. The range of colours is independent of group multiplicity or velocity dispersion, with isolated passive galaxies just as likely to have blue UV-to-optical colours, implying significant upturn components, as those in richer groups and in the previous data on clusters. This is supported by equivalent results for (FUV$-r$) colours which are clear indicators of upturn components.

preprint2020arXiv

Galaxy and Mass Assembly (GAMA): Demonstrating the power of WISE in the study of Galaxy Groups to $z<0.1$

Combining high-fidelity group characterisation from the Galaxy and Mass Assembly (GAMA) survey and source-tailored $z<0.1$ photometry from the WISE survey, we present a comprehensive study of the properties of ungrouped galaxies, compared to 497 galaxy groups (4$\leq$ N$_{\rm FoF}$ $\leq$ 20) as a function of stellar and halo mass. Ungrouped galaxies are largely unimodal in WISE color, the result of being dominated by star-forming, late-type galaxies. Grouped galaxies, however, show a clear bimodality in WISE color, which correlates strongly with stellar mass and morphology. We find evidence for an increasing early-type fraction, in stellar mass bins between $10^{10}\lesssim$M$_{\rm stellar} \lesssim10^{11}$ M$_\odot$, with increasing halo mass. Using ungrouped, late-type galaxies with star-forming colors (W2$-$W3$>$3), we define a star-forming main-sequence (SFMS), which we use to delineate systems that have moved below the sequence ("quenched" for the purposes of this work). We find that with increasing halo mass, the relative number of late-type systems on the SFMS decreases, with a corresponding increase in early-type, quenched systems at high stellar mass (M$_{\rm stellar}>{10}^{10.5}$ M$_\odot$), consistent with mass quenching. Group galaxies with masses M$_{\rm stellar}<{10}^{10.5}$ M$_\odot$ show evidence of quenching consistent with environmentally-driven processes. The stellar mass distribution of late-type, quenched galaxies suggests they may be an intermediate population as systems transition from being star-forming and late-type to the "red sequence". Finally, we use the projected area of groups on the sky to extract groups that are (relatively) compact for their halo mass. Although these show a marginal increase in their proportion of high mass and early-type galaxies compared to nominal groups, a clear increase in quenched fraction is not evident.

preprint2020arXiv

Galaxy And Mass Assembly (GAMA): the interplay between galaxy mass, SFR and heavy element abundance in paired galaxy sets

We study the star formation rate (SFR), stellar mass ($M_{\star}$) and the gas metallicity (Z) for 4,636 galaxy pairs using the Galaxy And Mass Assembly (GAMA) survey. Our galaxy pairs lie in a redshift range of 0 $<$ $z$ $<$ 0.35, mass range of 7.5 $<$ log( M$_{\star}$/ M$_\odot$) $<$ 11.5 and $ΔV$ $<$ 1000 km s$^{-1}$. We explore variations in SFR and Z from three point of views: multiplicity, pair separation and dynamics. We define multiplicity as the number of galaxies paired with a single galaxy, and analyzed for the first time variations in SFR and Z for both, single pairs and pairs with higher multiplicity. For the latter, we find SFR enhancements from 0.025-0.15 dex, that would shift the M-SFR relation of single pairs by 27$\%$ to higher SFRs. The effect of Z on the other hand, is of only 4$\%$. We analyze the most and least massive galaxy of major/minor pairs as a function of the pair separation. We define major pairs those with mass ratios of 0.5 $<$ $M_1$/$M_2$ $<$ 2, while pairs with more discrepant mass ratios are classified as minor pairs. We find SFR enhancements of up to 2 and 4 times with respect to their control sample, for major and minor pairs. For the case of Z, we find decrements of up to 0.08 dex for the closest pairs. When we focus on dynamics, Z enhancements are found for minor pairs with high velocity dispersion $(σ_p > 250 \; \mathrm{km\,s ^{-1}})$ and high multiplicity.

preprint2020arXiv

Mergers Do Trigger AGNs out to z $\sim$ 0.6

Aims. The fueling and feedback of Active Galactic Nuclei (AGNs) are important in understanding the co-evolution between black holes and host galaxies. Mergers are thought to have the capability to bring gas inwards and ignite nuclear activity, especially for more powerful AGNs. However, there is still significant ongoing debate on whether mergers can trigger AGNs and, if they do, whether mergers are a significant triggering mechanism. Methods. We select a low-redshift ($0.005<z<0.1$) sample from the Sloan Digital Sky Survey (SDSS) and a high-redshift ($0 < z < 0.6$) sample from the Galaxy And Mass Assembly (GAMA) survey. We take advantage of the convolutional neural network technique to identify mergers. We use mid-infrared (MIR) color cut and optical emission line diagnostics to classify AGNs. We also include Low Excitation Radio Galaxies (LERGs) to investigate the connection between mergers and low accretion rate AGNs. Results. We find that AGNs are more likely to be found in mergers than non-mergers, with an AGN excess up to $1.81\pm{0.16}$, suggesting that mergers can trigger AGNs. We also find the fraction of mergers in AGNs is higher than that in non-AGN controls, for both MIR and optically selected AGNs, as well as LERGs, with values between $16.40\pm{0.5}\%$ and $39.23\pm{2.10}\%$, implying a non-negligible to potentially significant role of mergers in triggering AGNs. This merger fraction in AGNs increases as stellar mass increases which supports the idea that mergers are more important for triggering AGNs in more massive galaxies. In terms of merger fraction as a function of AGN power we find a positive trend for MIR selected AGNs and a complex trend for optically selected AGNs, which we interpret under an evolutionary scenario proposed by previous studies. In addition, obscured MIR selected AGNs are more likely to be hosted in mergers than unobscured MIR selected AGNs.

preprint2020arXiv

WALLABY -- An SKA Pathfinder HI Survey

The Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) is a next-generation survey of neutral hydrogen (HI) in the Local Universe. It uses the widefield, high-resolution capability of the Australian Square Kilometer Array Pathfinder (ASKAP), a radio interferometer consisting of 36 x 12-m dishes equipped with Phased-Array Feeds (PAFs), located in an extremely radio-quiet zone in Western Australia. WALLABY aims to survey three-quarters of the sky (-90 degr < Dec < +30 degr) to a redshift of z < 0.26, and generate spectral line image cubes at ~30 arcsec resolution and ~1.6 mJy/beam per 4 km/s channel sensitivity. ASKAP's instantaneous field of view at 1.4 GHz, delivered by the PAF's 36 beams, is about 30 sq deg. At an integrated signal-to-noise ratio of five, WALLABY is expected to detect over half a million galaxies with a mean redshift of z ~ 0.05 (~200 Mpc). The scientific goals of WALLABY include: (a) a census of gas-rich galaxies in the vicinity of the Local Group; (b) a study of the HI properties of galaxies, groups and clusters, in particular the influence of the environment on galaxy evolution; and (c) the refinement of cosmological parameters using the spatial and redshift distribution of low-bias gas-rich galaxies. For context we provide an overview of previous large-scale HI surveys. Combined with existing and new multi-wavelength sky surveys, WALLABY will enable an exciting new generation of panchromatic studies of the Local Universe. - First results from the WALLABY pilot survey are revealed, with initial data products publicly available in the CSIRO ASKAP Science Data Archive (CASDA).

preprint2019arXiv

Galaxy And Mass Assembly (GAMA): Properties and evolution of red spiral galaxies

We use multi-wavelength data from the Galaxy and Mass Assembly (GAMA) survey to explore the cause of red optical colours in nearby (0.002<z<0.06) spiral galaxies. We show that the colours of red spiral galaxies are a direct consequence of some environment-related mechanism(s) which has removed dust and gas, leading to a lower star formation rate. We conclude that this process acts on long timescales (several Gyr) due to a lack of morphological transformation associated with the transition in optical colour. The sSFR and dust-to-stellar mass ratio of red spiral galaxies is found to be statistically lower than blue spiral galaxies. On the other hand, red spirals are on average $0.9$ dex more massive, and reside in environments 2.6 times denser than their blue counterparts. We find no evidence of excessive nuclear activity, or higher inclination angles to support these as the major causes for the red optical colours seen in >= 47% of all spirals in our sample. Furthermore, for a small subsample of our spiral galaxies which are detected in HI, we find that the SFR of gas-rich red spiral galaxies is lower by ~1 dex than their blue counterparts.

preprint2016arXiv

All NIRspec needs is HST/WFC3 pre-imaging? The use of Milky Way Stars in WFC3 Imaging to Register NIRspec MSA Observations

The James Webb Space Telescope (JWST) will be an exquisite new near-infrared observatory with imaging and multi-object spectroscopy through ESA's NIRspec instrument with its unique Micro-Shutter Array (MSA), allowing for slits to be positioned on astronomical targets by opening specific 0.002"-wide micro shutter doors. To ensure proper target acquisition, the on-sky position of the MSA needs to be verified before spectroscopic observations start. An onboard centroiding program registers the position of pre-identified guide stars in a Target Acquisition (TA) image, a short pre-spectroscopy exposure without dispersion (image mode) through the MSA with all shutters open. The outstanding issue is the availability of Galactic stars in the right luminosity range for TA relative to typical high redshift targets. We explore this here using the stars and $z\sim8$ candidate galaxies identified in the source extractor catalogs of Brightest of Reionizing Galaxies survey (BoRG[z8]), a pure-parallel program with Hubble Space Telescope Wide-Field Camera 3. We find that (a) a single WFC3 field contains enough Galactic stars to satisfy the NIRspec astrometry requirement (20 milli-arcseconds), provided its and the NIRspec TA's are $m_{lim}>24.5$ AB in WFC3 F125W, (b) a single WFC3 image can therefore serve as the pre-image if need be, (c) a WFC3 mosaic and accompanying TA image satisfy the astrometry requirement at $\sim23$ AB mag in WFC3 F125W, (d) no specific Galactic latitude requires deeper TA imaging due to a lack of Galactic stars, and (e) a depth of $\sim24$ AB mag in WFC3 F125W is needed if a guide star in the same MSA quadrant as a target is required. We take the example of a BoRG identified $z\sim8$ candidate galaxy and require a Galactic star within 20" of it. In this case, a depth of 25.5 AB in F125W is required (with $\sim$97% confidence).

preprint2016arXiv

Galaxy And Mass Assembly (GAMA): Accurate Panchromatic Photometry from Optical Priors using LAMBDAR

We present the Lambda Adaptive Multi-Band Deblending Algorithm in R (LAMBDAR), a novel code for calculating matched aperture photometry across images that are neither pixel- nor PSF-matched, using prior aperture definitions derived from high resolution optical imaging. The development of this program is motivated by the desire for consistent photometry and uncertainties across large ranges of photometric imaging, for use in calculating spectral energy distributions. We describe the program, specifically key features required for robust determination of panchromatic photometry: propagation of apertures to images with arbitrary resolution, local background estimation, aperture normalisation, uncertainty determination and propagation, and object deblending. Using simulated images, we demonstrate that the program is able to recover accurate photometric measurements in both high-resolution, low-confusion, and low-resolution, high-confusion, regimes. We apply the program to the 21-band photometric dataset from the Galaxy And Mass Assembly (GAMA) Panchromatic Data Release (PDR; Driver et al. 2016), which contains imaging spanning the far-UV to the far-IR. We compare photometry derived from LAMBDAR with that presented in Driver et al. (2016), finding broad agreement between the datasets. Nonetheless, we demonstrate that the photometry from LAMBDAR is superior to that from the GAMA PDR, as determined by a reduction in the outlier rate and intrinsic scatter of colours in the LAMBDAR dataset. We similarly find a decrease in the outlier rate of stellar masses and star formation rates using LAMBDAR photometry. Finally, we note an exceptional increase in the number of UV and mid-IR sources able to be constrained, which is accompanied by a significant increase in the mid-IR colour-colour parameter-space able to be explored.

preprint2016arXiv

Galaxy And Mass Assembly (GAMA): detection of low-surface-brightness galaxies from SDSS data

We report on a search for new low-surface-brightness galaxies (LSBGs) using Sloan Digital Sky Survey (SDSS) data within the GAMA equatorial fields. The search method consisted of masking objects detected with SDSS photo, combining gri images weighted to maximise the expected signal-to-noise ratio (SNR), and smoothing the images. The processed images were then run through a detection algorithm that finds all pixels above a set threshold and groups them based on their proximity to one another. The list of detections was cleaned of contaminants such as diffraction spikes and the faint wings of masked objects. From these, selecting potentially the brightest in terms of total flux, a list of 343 LSBGs was produced having been confirmed using VISTA Kilo-degree Infrared Galaxy Survey (VIKING) imaging. The photometry of this sample was refined using the deeper VIKING Z band as the aperture-defining band. Measuring their $g-i$ and $J-K$ colours shows that most are consistent with being at redshifts less than 0.2. The photometry is carried out using an AUTO aperture for each detection giving surface brightnesses of $μ_{r} \ge 25 \, \mathrm{mag} \, \mathrm{arcsec}^{-2}$ and magnitudes of $r > 19.8$ mag. None of these galaxies are bright enough to be within the GAMA main survey limit but could be part of future deeper surveys to measure the low-mass end of the galaxy stellar mass function.

preprint2016arXiv

GAMA/H-ATLAS: Common star-formation rate indicators and their dependence on galaxy physical parameters

We compare common star-formation rate (SFR) indicators in the local Universe in the GAMA equatorial fields (around 160 sq. deg.), using ultraviolet (UV) photometry from GALEX, far-infrared (FIR) and sub-millimetre (sub-mm) photometry from H-ATLAS, and Halpha spectroscopy from the GAMA survey. With a high-quality sample of 745 galaxies (median redshift 0.08), we consider three SFR tracers: UV luminosity corrected for dust attenuation using the UV spectral slope beta (SFRUV,corr), Halpha line luminosity corrected for dust using the Balmer decrement (BD) (SFRHalpha,corr), and the combination of UV and IR emission (SFRUV+IR). We demonstrate that SFRUV,corr can be reconciled with the other two tracers after applying attenuation corrections by calibrating IRX (i.e. the IR to UV luminosity ratio) and attenuation in the Halpha (derived from BD) against beta. However, beta on its own is very unlikely to be a reliable attenuation indicator. We find that attenuation correction factors depend on parameters such as stellar mass, z and dust temperature (Tdust), but not on Halpha equivalent width (EW) or Sersic index. Due to the large scatter in the IRX vs beta correlation, when compared to SFRUV+IR, the beta-corrected SFRUV,corr exhibits systematic deviations as a function of IRX, BD and Tdust.

preprint2016arXiv

Milky Way Red Dwarfs in the BoRG Survey; Galactic scale-height and the distribution of dwarfs stars in WFC3 imaging

We present a tally of Milky Way late-type dwarf stars in 68 WFC3 pure-parallel fields (227 arcmin^2) from the Brightest of Reionizing Galaxies (BoRG) survey for high-redshift galaxies. Using spectroscopically identified M-dwarfs in two public surveys, the CANDELS and the ERS mosaics, we identify a morphological selection criterion using the half-light radius (r50), a near-infrared J-H, G-J color region where M-dwarfs are found, and a V-J relation with M-dwarf subtype. We apply this morphological selection of stellar objects, color-color selection of M-dwarfs and optical-near-infrared color subtyping to compile a catalog of 274 M-dwarfs belonging to the disk of the Milky Way with a limiting magnitude of m_F125W < 24. Based on the M-dwarfs statistics, we conclude that (a) the previously identified North/South discrepancy in M-dwarf numbers persists in our sample; there are more M-dwarfs in the Northern fields on average than in Southern ones, (b) the Milky Way's single disk scale-height for M-dwarfs is 0.3-4 kpc, depending on sub-type, (c) {\bf ERRATUM:} we present corrected coordinates (AstroPy) and distances and find a constant $z_0$=600 pc for all types. (d) a second component is visible in the vertical distribution, with a different, much higher scale-height. We report the M-dwarf component of the Sagittarius stream in one of our fields with 11 confirmed M-dwarfs, 7 of which are at the stream's distance. The dwarf scale-height and the relative low incidence in our fields of L- and T-dwarfs in these fields makes it unlikely that these stars will be interlopers in great numbers in color-selected samples of high-redshift galaxies. The relative ubiquity of M-dwarfs however will make them ideal tracers of Galactic Halo substructure with EUCLID and reference stars for JWST observations.

preprint2016arXiv

The Size and Shape of the Milky Way Disk and Halo from M-type Brown Dwarfs in the BoRG Survey

We have identified 274 M-type Brown Dwarfs in the Hubble Space Telescope's Wide Field Camera 3 (WFC3) pure parallel fields from the Brightest of Reionizing Galaxies (BoRG) survey for high redshift galaxies. These are near-infrared observations with multiple lines-of-sight out of our Milky Way. Using these observed M-type Brown Dwarfs we fitted a Galactic disk and halo model with a Markov chain Monte Carlo (MCMC) analysis. This model worked best with the scale length of the disk fixed at $h$ = 2.6 kpc. For the scale height of the disk, we found $z_0 = 0.29^{+0.02}_{-0.019}$ kpc and for the central number density $ρ_0 = 0.29^{+0.20}_{-0.13}$ \#/pc$^3$. For the halo we derived a flattening parameter $κ$ = 0.45$\pm{0.04}$ and a power-law index $p$ = 2.4$\pm{0.07}$. We found the fraction of M-type brown dwarfs in the local density that belong to the halo to be $f_{h}$ = 0.0075$^{+0.0025}_{-0.0019}$. We found no correlation between subtype of M-dwarf and any model parameters. The total number of M-type Brown Dwarfs in the disk and halo was determined to be $58.2^{+9.81}_{-6.70} \times10^{9}$. We found an upper limit for the fraction of M-type Brown Dwarfs in the halo of 7$^{+5}_{-4}$\%. The upper limit for the total Galactic Disk mass in M-dwarfs is $4.34^{+0.73}_{-0.5}\times10^{9}$ $M_{\odot}$, assuming all M-type Brown Dwarfs have a mass of $80 M_J$.

preprint2016arXiv

WISE x SuperCOSMOS photometric redshift catalog: 20 million galaxies over 3pi steradians

We cross-match the two currently largest all-sky photometric catalogs, mid-infrared WISE and SuperCOSMOS scans of UKST/POSS-II photographic plates, to obtain a new galaxy sample that covers 3pi steradians. In order to characterize and purify the extragalactic dataset, we use external GAMA and SDSS spectroscopic information to define quasar and star loci in multicolor space, aiding the removal of contamination from our extended-source catalog. After appropriate data cleaning we obtain a deep wide-angle galaxy sample that is approximately 95% pure and 90% complete at high Galactic latitudes. The catalog contains close to 20 million galaxies over almost 70% of the sky, outside the Zone of Avoidance and other confused regions, with a mean surface density of over 650 sources per square degree. Using multiwavelength information from two optical and two mid-IR photometric bands, we derive photometric redshifts for all the galaxies in the catalog, using the ANNz framework trained on the final GAMA-II spectroscopic data. Our sample has a median redshift of z_{med} = 0.2 but with a broad dN/dz reaching up to z>0.4. The photometric redshifts have a mean bias of |delta_z|~10^{-3}, normalized scatter of sigma_z = 0.033 and less than 3% outliers beyond 3sigma_z. Comparison with external datasets shows no significant variation of photo-z quality with sky position. Together with the overall statistics, we also provide a more detailed analysis of photometric redshift accuracy as a function of magnitudes and colors. The final catalog is appropriate for `all-sky' 3D cosmology to unprecedented depths, in particular through cross-correlations with other large-area surveys. It should also be useful for source pre-selection and identification in forthcoming surveys such as TAIPAN or WALLABY.

preprint2015arXiv

ARRAKIS: Atlas of Resonance Rings As Known In the S4G

Resonance rings are the consequence of secular evolution processes that redistribute material and angular momentum in discs. We produced a Catalogue and an Atlas of the rings detected in the Spitzer Survey of Stellar Structure in Galaxies (S4G) and to conduct a statistical study of the data in the Catalogue. We traced the contours of rings previously identified by Buta et al. (in preparation) and fitted them with ellipses. We found the orientation of bars by studying the galaxy ellipse fits from S4G's Pipeline 4. We used the galaxy orientation data obtained by S4G's Pipeline 4 to obtain intrinsic ellipticities and orientations of rings and the bars. ARRAKIS contains data on 724 ringed galaxies in the S4G. The frequency of resonance rings in the S4G is of 16+-1% and 35+-1% for outer and inner features, respectively. Outer rings are mostly found in Hubble stages -1<=T<=4. Inner rings are found in a distribution that covers the range -1<=T<=7. We confirm that outer rings have two preferred orientations, parallel and perpendicular to the bar. We confirm a tendency for inner rings to be oriented parallel to the bar, but we find that a significant fraction (~50%) of them have random orientations with respect to the bar. These misaligned inner rings are mostly found in late-type galaxies (T>=4). This may be due to spiral modes decoupled from the bar dominating the Fourier amplitude spectrum at the radius of the inner ring. We find that the fraction of barred galaxies hosting outer (inner) rings is ~1.7 times (~1.3 times) that in unbarred galaxies. The fact that rings are only mildly favoured by bars suggests that those in unbarred galaxies either formed due to weak departures from the axisymmetry of the galactic potential or that they are born because of bars that have been destroyed after the ring formation.

preprint2015arXiv

Bright galaxies at Hubble's redshift detection frontier: Preliminary results and design from the redshift z~9-10 BoRG pure-parallel HST survey

We present the first results and design from the redshift z~9-10 Brightest of the Reionizing Galaxies {\it Hubble Space Telescope} survey BoRG[z9-10], aimed at searching for intrinsically luminous unlensed galaxies during the first 700 Myr after the Big Bang. BoRG[z9-10] is the continuation of a multi-year pure-parallel near-IR and optical imaging campaign with the Wide Field Camera 3. The ongoing survey uses five filters, optimized for detecting the most distant objects and offering continuous wavelength coverage from λ=0.35μm to λ=1.7μm. We analyze the initial ~130 arcmin$^2$ of area over 28 independent lines of sight (~25% of the total planned) to search for z>7 galaxies using a combination of Lyman break and photometric redshift selections. From an effective comoving volume of (5-25) $times 10^5$ Mpc$^3$ for magnitudes brighter than $m_{AB}=26.5-24.0$ in the $H_{160}$-band respectively, we find five galaxy candidates at z~8.3-10 detected at high confidence (S/N>8), including a source at z~8.4 with mAB=24.5 (S/N~22), which, if confirmed, would be the brightest galaxy identified at such early times (z>8). In addition, BoRG[z9-10] data yield four galaxies with $7.3 \lesssim z \lesssim 8$. These new Lyman break galaxies with m$\lesssim26.5$ are ideal targets for follow-up observations from ground and space based observatories to help investigate the complex interplay between dark matter growth, galaxy assembly, and reionization.

preprint2015arXiv

Galaxy And Mass Assembly (GAMA) Blended Spectra Catalog: Strong Galaxy-Galaxy Lens and Occulting Galaxy Pair Candidates

We present the catalogue of blended galaxy spectra from the Galaxy And Mass Assembly (GAMA) survey. These are cases where light from two galaxies are significantly detected in a single GAMA fibre. Galaxy pairs identified from their blended spectrum fall into two principal classes: they are either strong lenses, a passive galaxy lensing an emission-line galaxy; or occulting galaxies, serendipitous overlaps of two galaxies, of any type. Blended spectra can thus be used to reliably identify strong lenses for follow-up observations (high resolution imaging) and occulting pairs, especially those that are a late-type partly obscuring an early-type galaxy which are of interest for the study of dust content of spiral and irregular galaxies. The GAMA survey setup and its autoz automated redshift determination were used to identify candidate blended galaxy spectra from the cross-correlation peaks. We identify 280 blended spectra with a minimum velocity separation of 600 km/s, of which 104 are lens pair candidates, 71 emission-line-passive pairs, 78 are pairs of emission-line galaxies and and 27 are pairs of galaxies with passive spectra. We have visually inspected the candidates in the Sloan Digital Sky Survey (SDSS) and Kilo Degree Survey (KiDS) images. Many blended objects are ellipticals with blue fuzz (Ef in our classification). These latter "Ef" classifications are candidates for possible strong lenses, massive ellipticals with an emission-line galaxy in one or more lensed images. The GAMA lens and occulting galaxy candidate samples are similar in size to those identified in the entire SDSS. This blended spectrum sample stands as a testament of the power of this highly complete, second-largest spectroscopic survey in existence and offers the possibility to expand e.g., strong gravitational lens surveys.

preprint2015arXiv

Galaxy And Mass Assembly (GAMA): end of survey report and data release 2

The Galaxy And Mass Assembly (GAMA) survey is one of the largest contemporary spectroscopic surveys of low-redshift galaxies. Covering an area of ~286 deg^2 (split among five survey regions) down to a limiting magnitude of r < 19.8 mag, we have collected spectra and reliable redshifts for 238,000 objects using the AAOmega spectrograph on the Anglo-Australian Telescope. In addition, we have assembled imaging data from a number of independent surveys in order to generate photometry spanning the wavelength range 1 nm - 1 m. Here we report on the recently completed spectroscopic survey and present a series of diagnostics to assess its final state and the quality of the redshift data. We also describe a number of survey aspects and procedures, or updates thereof, including changes to the input catalogue, redshifting and re-redshifting, and the derivation of ultraviolet, optical and near-infrared photometry. Finally, we present the second public release of GAMA data. In this release we provide input catalogue and targeting information, spectra, redshifts, ultraviolet, optical and near-infrared photometry, single-component Sérsic fits, stellar masses, H$α$-derived star formation rates, environment information, and group properties for all galaxies with r < 19.0 mag in two of our survey regions, and for all galaxies with r < 19.4 mag in a third region (72,225 objects in total). The database serving these data is available at http://www.gama-survey.org/.

preprint2015arXiv

Galaxy And Mass Assembly (GAMA): growing up in a bad neighbourhood - how do low-mass galaxies become passive?

Both theoretical predictions and observations of the very nearby Universe suggest that low-mass galaxies (log$_{10}$[M$_{*}$/M$_{\odot}$]<9.5) are likely to remain star-forming unless they are affected by their local environment. To test this premise, we compare and contrast the local environment of both passive and star-forming galaxies as a function of stellar mass, using the Galaxy and Mass Assembly survey. We find that passive fractions are higher in both interacting pair and group galaxies than the field at all stellar masses, and that this effect is most apparent in the lowest mass galaxies. We also find that essentially all passive log$_{10}$[M$_{*}$/M$_{\odot}$]<8.5 galaxies are found in pair/group environments, suggesting that local interactions with a more massive neighbour cause them to cease forming new stars. We find that the effects of immediate environment (local galaxy-galaxy interactions) in forming passive systems increases with decreasing stellar mass, and highlight that this is potentially due to increasing interaction timescales giving sufficient time for the galaxy to become passive via starvation. We then present a simplistic model to test this premise, and show that given our speculative assumptions, it is consistent with our observed results.

preprint2015arXiv

Galaxy And Mass Assembly (GAMA): The Bright Void Galaxy Population in the Optical and Mid-IR

We examine the properties of galaxies in the Galaxies and Mass Assembly (GAMA) survey located in voids with radii $>10~h^{-1}$ Mpc. Utilising the GAMA equatorial survey, 592 void galaxies are identified out to z~0.1 brighter than $M_{r} = -18.4$, our magnitude completeness limit. Using the $W_{\rm{Hα}}$ vs. [NII]/H$α$ (WHAN) line strength diagnostic diagram, we classify their spectra as star forming, AGN, or dominated by old stellar populations. For objects more massive than $5\times10^{9}$ M$_{\odot}$, we identify a sample of 26 void galaxies with old stellar populations classed as passive and retired galaxies in the WHAN diagnostic diagram, else they lack any emission lines in their spectra. When matched to WISE mid-IR photometry, these passive and retired galaxies exhibit a range of mid-IR colour, with a number of void galaxies exhibiting [4.6]-[12] colours inconsistent with completely quenched stellar populations, with a similar spread in colour seen for a randomly drawn non-void comparison sample. We hypothesise that a number of these galaxies host obscured star formation, else they are star forming outside of their central regions targeted for single fibre spectroscopy. When matched to a randomly drawn sample of non-void galaxies, the void and non-void galaxies exhibit similar properties in terms of optical and mid-IR colour, morphology, and star formation activity, suggesting comparable mass assembly and quenching histories. A trend in mid-IR [4.6]-[12] colour is seen, such that both void and non-void galaxies with quenched/passive colours <1.5 typically have masses higher than $10^{10}$ M$_{\odot}$, where internally driven processes play an increasingly important role in galaxy evolution.

preprint2015arXiv

HERschel Observations of Edge-on Spirals (HEROES). II: Tilted-ring modelling of the atomic gas disks

Context. Edge-on galaxies can offer important insights in galaxy evolution as they are the only systems where the distribution of the different components can be studied both radially and vertically. The HEROES project was designed to investigate the interplay between the gas, dust, stars and dark matter (DM) in a sample of 7 massive edge-on spiral galaxies. Aims. In this second HEROES paper we present an analysis of the atomic gas content of 6 out of 7 galaxies in our sample. The remaining galaxy was recently analysed according to the same strategy. The primary aim of this work is to constrain the surface density distribution, the rotation curve and the geometry of the gas disks in a homogeneous way. In addition we identify peculiar features and signs of recent interactions. Methods. We construct detailed tilted-ring models of the atomic gas disks based on new GMRT 21-cm observations of NGC 973 and UGC 4277 and re-reduced archival HI data of NGC 5907, NGC 5529, IC 2531 and NGC 4217. Potential degeneracies between different models are resolved by requiring a good agreement with the data in various representations of the data cubes. Results. From our modelling we find that all but one galaxy are warped along the major axis. In addition, we identify warps along the line of sight in three galaxies. A flaring gas layer is required to reproduce the data only for one galaxy, but (moderate) flares cannot be ruled for the other galaxies either. A coplanar ring-like structure is detected outside the main disk of NGC 4217, which we suggest could be the remnant of a recent minor merger event. We also find evidence for a radial inflow of 15 +- 5 km/s in the disk of NGC 5529, which might be related to the ongoing interaction with two nearby companions. (Abridged)

preprint2015arXiv

The Dependence of the $A_V$ Prior for SN\,Ia on Host Mass and Disk Inclination

Supernovae type Ia (SNIa) are used as "standard candles" for cosmological distance scales. To fit their light curve shape -- absolute luminosity relation, one needs to assume an intrinsic color and a likelihood of host galaxy extinction or a convolution of these, a color distribution prior. The host galaxy extinction prior is typically assumed to be an exponential drop-off for the current supernova programs ($P(A_V) \propto e^{-A_V/τ_0}$). We explore the validity of this prior using the distribution of extinction values inferred when two galaxies accidentally overlap (an occulting galaxy pair). We correct the supernova luminosity distances from the SDSS-III Supernova projects (SDSS-SN) by matching the host galaxies to one of three templates from occulting galaxy pairs based on the host galaxy mass and the $A_V$-bias - prior-scale ($τ_0$) relation from Jha et al. (2007). We find that introducing an $A_V$ prior that depends on host mass results in lowered luminosity distances for the SDSS-SN on average but it does not reduce the scatter in individual measurements. This points, in our view, to the need for many more occulting galaxy templates to match to SNIa host galaxies to rule out this possible source of scatter in the SNIa distance measurements. We match occulting galaxy templates based on both mass and projected radius and we find that one should match by stellar mass first with radius as a secondary consideration. We discuss the caveats of the current approach and our aim is to convince the reader that a library of occulting galaxy pairs observed with HST will provide sufficient priors to improve (optical) SNIa measurements to the next required accuracy in Cosmology.

preprint2015arXiv

The Spitzer Survey of Stellar Structure in Galaxies (S$^4$G): Precise Stellar Mass Distributions from Automated Dust Correction at 3.6 microns

The mid-infrared is an optimal window to trace stellar mass in nearby galaxies and the 3.6$μm$ IRAC band has been exploited to this effect, but such mass estimates can be biased by dust emission. We present our pipeline to reveal the old stellar flux at 3.6$μm$ and obtain stellar mass maps for more than 1600 galaxies available from the Spitzer Survey of Stellar Structure in Galaxies (S$^{4}$G). This survey consists of images in two infrared bands (3.6 and 4.5$μm$), and we use the Independent Component Analysis (ICA) method presented in Meidt et al. (2012) to separate the dominant light from old stars and the dust emission that can significantly contribute to the observed 3.6$μm$ flux. We exclude from our ICA analysis galaxies with low signal-to-noise ratio (S/N < 10) and those with original [3.6]-[4.5] colors compatible with an old stellar population, indicative of little dust emission (mostly early Hubble types, which can directly provide good mass maps). For the remaining 1251 galaxies to which ICA was successfully applied, we find that as much as 10-30% of the total light at 3.6$μm$ typically originates from dust, and locally it can reach even higher values. This contamination fraction shows a correlation with specific star formation rates, confirming that the dust emission that we detect is related to star formation. Additionally, we have used our large sample of mass estimates to calibrate a relationship of effective mass-to-light ratio ($M/L$) as a function of observed [3.6]-[4.5] color: $\log(M/L)=-0.339 (\pm 0.057) \times ([3.6]-[4.5]) -0.336 (\pm 0.002)$. Our final pipeline products have been made public through IRSA, providing the astronomical community with an unprecedentedly large set of stellar mass maps ready to use for scientific applications.

preprint2015arXiv

What can the Occult do for you? STarlight Attenuation & Reddening Survey of Multiple Occulting Galaxies (STARSMOG)

Interstellar dust is still the dominant uncertainty in Astronomy, limiting precision in e.g., cosmological distance estimates and models of how light is re-processed within a galaxy. When a foreground galaxy serendipitously overlaps a more distant one, the latter backlights the dusty structures in the nearer foreground galaxy. Such an overlapping or occulting galaxy pair can be used to measure the distribution of dust in the closest galaxy with great accuracy. The STARSMOG program uses HST observation of occulting galaxy pairs to accurately map the distribution of dust in foreground galaxies in fine ($<$100 pc) detail. Furthermore, Integral Field Unit observations of such pairs will map the effective extinction curve in these occulting galaxies, disentangling the role of fine-scale geometry and grain composition on the path of light through a galaxy. The overlapping galaxy technique promises to deliver a clear understanding of the dust in galaxies: the dust geometry, a probability function of the amount of dimming as a function of galaxy type, its dependence on wavelength, and evolution of all these properties with cosmic time using distant, high-redshift pairs.

preprint2014arXiv

Galaxy And Mass Assembly (GAMA): autoz spectral redshift measurements, confidence and errors

The Galaxy And Mass Assembly (GAMA) survey has obtained spectra of over 230000 targets using the Anglo-Australian Telescope. To homogenise the redshift measurements and improve the reliability, a fully automatic redshift code was developed (autoz). The measurements were made using a cross-correlation method for both absorption-line and emission-line spectra. Large deviations in the high-pass filtered spectra are partially clipped in order to be robust against uncorrected artefacts and to reduce the weight given to single-line matches. A single figure of merit (FOM) was developed that puts all template matches onto a similar confidence scale. The redshift confidence as a function of the FOM was fitted with a tanh function using a maximum likelihood method applied to repeat observations of targets. The method could be adapted to provide robust automatic redshifts for other large galaxy redshift surveys. For the GAMA survey, there was a substantial improvement in the reliability of assigned redshifts and in the lowering of redshift uncertainties with a median velocity uncertainty of 33 km/s.

preprint2014arXiv

Galaxy And Mass Assembly (GAMA): Galaxy close-pairs, mergers, and the future fate of stellar mass

We use a highly complete subset of the GAMA-II redshift sample to fully describe the stellar mass dependence of close-pairs and mergers between 10^8 Msun and 10^12 Msun. Using the analytic form of this fit we investigate the total stellar mass accreting onto more massive galaxies across all mass ratios. Depending on how conservatively we select our robust merging systems, the fraction of mass merging onto more massive companions is 2.0%-5.6%. Using the GAMA-II data we see no significant evidence for a change in the close-pair fraction between redshift $z = 0.05-0.2$. However, we find a systematically higher fraction of galaxies in similar mass close-pairs compared to published results over a similar redshift baseline. Using a compendium of data and the function $γ_M =A(1+z)m$ to predict the major close-pair fraction, we find fitting parameters of $A = 0.021 \pm 0.001$ and $m = 1.53 \pm 0.08$, which represents a higher low-redshift normalisation and shallower power-law slope than recent literature values. We find that the relative importance of in-situ star-formation versus galaxy merging is inversely correlated, with star-formation dominating the addition of stellar material below Mstar and merger accretion events dominating beyond Mstar. We find mergers have a measurable impact on the whole extent of the GSMF, manifest as a deepening of the dip in the GSMF over the next Gyr and an increase in Mstar by as much as 0.01-0.05 dex.

preprint2014arXiv

Galaxy And Mass Assembly (GAMA): Mid-Infrared Properties and Empirical Relations from $WISE$

The Galaxy And Mass Assembly (GAMA) survey furnishes a deep redshift catalog that, when combined with the Wide-field Infrared Explorer ($WISE$), allows us to explore for the first time the mid-infrared properties of $> 110, 000$ galaxies over 120 deg$^2$ to $z\simeq 0.5$. In this paper we detail the procedure for producing the matched GAMA-$WISE$ catalog for the G12 and G15 fields, in particular characterising and measuring resolved sources; the complete catalogs for all three GAMA equatorial fields will be made available through the GAMA public releases. The wealth of multiwavelength photometry and optical spectroscopy allows us to explore empirical relations between optically determined stellar mass (derived from synthetic stellar population models) and 3.4micron and 4.6micron WISE measurements. Similarly dust-corrected Halpha-derived star formation rates can be compared to 12micron and 22micron luminosities to quantify correlations that can be applied to large samples to $z<0.5$. To illustrate the applications of these relations, we use the 12micron star formation prescription to investigate the behavior of specific star formation within the GAMA-WISE sample and underscore the ability of WISE to detect star-forming systems at $z\sim0.5$. Within galaxy groups (determined by a sophisticated friends-of-friends scheme), results suggest that galaxies with a neighbor within 100$\,h^{-1} $kpc have, on average, lower specific star formation rates than typical GAMA galaxies with the same stellar mass.

preprint2014arXiv

SNIa Host Galaxy Properties and the Dust Extinction Distribution

Supernovae Type Ia display a complex relation with their host galaxies. An important prior to the fit of the supernovae's lightcurve is the distribution of host galaxy extinction values that can be encountered. The SDSS-SN project has published light curve fits using both MLCS2k2 and SALT2. We use the former fits extinction parameter ($A_V$) to map this distribution of extinction values. We explore the dependence of this distribution on four observables; the inclination of the host galaxy disk, radial position of the supernova, redshift of the supernova and host, and the level of star-formation in the host galaxy. The distribution of $A_V$ values encountered by supernovae is typically characterised by: $\rm N_0 ~ e^{-A_V/τ}$, with $τ$= 0.4 or 0.33. We find that the inclination correction using an infinitely thin disk for the SNIa is sufficient, resulting in similar exponential $A_V$ distributions for high- and low-inclination disks. The $A_V$ distribution also depends on the radial position in the disk, consistent with previous results on the transparency of spiral disks. The distribution of $A_V$ values narrows with increased star-formation, possibly due to the destruction or dispersion of the dusty ISM by stellar winds prior to the ignition of the supernova. In future supernova searches, certainly the inclination of the host galaxy disk, should be considered in the construction of the \av \ prior with $τ=0.4/cos(i)$ as the most likely prior in each individual host galaxy's case.

preprint2013arXiv

HERschel Observations of Edge-on Spirals (HEROES). I: Far-infrared morphology and dust mass determination

Context. Edge-on spiral galaxies with prominent dust lanes provide us with an excellent opportunity to study the distribution and properties of the dust within them. The HEROES project was set up to observe a sample of seven large edge-on galaxies across various wavelengths for this investigation. Aims. Within this first paper, we present the Herschel observations and perform a qualitative and quantitative analysis on them, and we derive some global properties of the far infrared and submillimetre emission. Methods. We determine horizontal and vertical profiles from the Herschel observations of the galaxies in the sample and describe the morphology. Modified black-body fits to the global fluxes, measured using aperture photometry, result in dust temperatures and dust masses. The latter values are compared to those that are derived from radiative transfer models taken from the literature. Results. On the whole, our Herschel flux measurements agree well with archival values. We find that the exponential horizontal dust distribution model often used in the literature generally provides a good description of the observed horizontal profiles. Three out of the seven galaxies show signatures of extended vertical emission at 100 and 160 μm at the 5σ level, but in two of these it is probably due to deviations from an exactly edge-on orientation. Only for NGC 4013, a galaxy in which vertically extended dust has already been detected in optical images, we can detect vertically extended dust, and the derived scaleheight agrees with the value estimated through radiative transfer modelling. Our analysis hints at a correlation between the dust scaleheight and its degree of clumpiness, which we infer from the difference between the dust masses as calculated from modelling of optical data and from fitting the spectral energy distribution of Herschel datapoints.

preprint2013arXiv

Morphological Parameters of Spitzer Survey of Stellar Structure in Galaxies (S4G)

The morphology of galaxies can be quantified to some degree using a set of scale-invariant parameters. Concentration (C), Asymmetry (A), Smoothness (S), the Gini index (G), relative contribution of the brightest pixels to the second order moment of the flux (M20), ellipticity (E), and the Gini index of the second order moment (GM) have all been applied to morphologically classify galaxies at various wavelengths. Here we present a catalog of these parameters for the Spitzer Survey of Stellar Structure in Galaxies (S4G), a volume-limited near-infrared imaging survey of nearby galaxies using the 3.6 and 4.5 micron channels of the IRAC camera. Our goal is to provide a reference catalog of near-infrared quantified morphology for high-redshift studies and galaxy evolution models with enough detail to resolve stellar mass morphology. We explore where normal, non-interacting galaxies -those typically found on the Hubble tuning fork- lie in this parameter space and show that there is a tight relation between Concentration and M20 for normal galaxies. M20 can be used to classify galaxies into earlier and later types (e.g., to separate spirals from irregulars). Several criteria using these parameters exist to select systems with a disturbed morphology, i.e., those that appear to be undergoing a tidal interaction. We examine the applicability of these criteria to Spitzer near-infrared imaging. We find that four relations, based on the parameters A&S, G&M20, GM, and C&M20, respectively, select outliers in the morphological parameter space, but each selects different subsets of galaxies. Two criteria (GM > 0.6, G > -0.115 x M20 + 0.384) seem most appropriate to identify possible mergers and the merger fraction in near-infrared surveys. We find no strong relation between lopsidedness and most of these morphological parameters, except for a weak dependence of lopsidedness on Concentration and M20.

preprint2013arXiv

Quantified HI Morphology VII: star-formation and tidal influence on local dwarf HI morphology

Scale-invariant morphology parameters applied to atomic hydrogen maps (HI) of galaxies can be used to quantify the effects of tidal interaction or star-formation on the ISM. Here we apply these parameters, Concentration, Asymmetry, Smoothness, Gini, M20, and the GM parameter, to two public surveys of nearby dwarf galaxies, the VLA-ANGST and LITTLE-THINGS survey, to explore whether tidal interaction or the ongoing or past star-formation is a dominant force shaping the HI disk of these dwarfs. Previously, HI morphological criteria were identified for ongoing spiral-spiral interactions. When we apply these to the Irregular dwarf population, they either select almost all or none of the population. We find that only the Asymmetry-based criteria can be used to identify very isolated dwarfs (i.e., these have a low tidal indication). Otherwise, there is little or no relation between the level of tidal interaction and the HI morphology. We compare the HI morphology to three star-formation rates based on either Halpha, FUV or the resolved stellar population, probing different star-formation time-scales. The HI morphology parameters that trace the inequality of the distribution, the Gini, GM, and M20 parameters, correlate weakly with all these star-formation rates. This is in line with the picture that local physics dominates the ISM appearance and not tidal effects. Finally, we compare the SDSS measures of star-formation and stellar mass to the HI morphological parameters for all four HI surveys. In the two lower-resolution HI surveys (12"), there is no relation between star-formation measures and HI morphology. The morphology of the two high-resolution HI surveys (6"), the Asymmetry, Smoothness, Gini, M20, and GM, do show a link to the total star-formation, but a weak one.

preprint2013arXiv

The Occulting Galaxy Pair UGC 3995; Dust properties from HST and CALIFA data

UGC 3995 is an interacting and occulting galaxy pair. UGC 3995B is a foreground face-on spiral and UGC 3995A a bright background spiral with an AGN. We present analysis of the dust in the disc of UGC 3995B based on archival Hubble Space Telescope WFPC2 and PPAK IFU data from the CALIFA survey's first data release. From the HST F606W image, we construct an extinction map by modeling the isophotes of the background galaxy UGC 3995A and the resulting transmission through UGC 3995B. This extinction map of UGC 3995B shows several distinct spiral extinction features. The radial distribution of $A_V$ values declines slowly with peaks corresponding to the spiral structures. The distribution of $A_V$ values in the HST extinction map peaks near $A_V = 0.3-0.4$. Beyond this point, the distribution of $A_V$ values drops like an exponential: $N(A_V) = N_0 \times e^{(-A_V/0.5)}$. The 0.5 value is higher than typical for a spiral galaxy. The outer arms may be tidally distended; the extinction in the corresponding interarm regions is small to an unusually small radius. To analyse the PPAK IFU data, we take the ratio of a fibre spectrum in the overlap region and the corresponding background fiber spectrum to construct an extinction curve. We fit the Cardelli, Clayton and Mathis (CCM) curve to the extinction curve of each fiber element in the overlap region. A map of the extinction constructed from PPEX IFU data-cubes shows the same spiral structure of the HST extinction map but the some differences in the distribution of the normalization of the CCM fits ($A_V$). The inferred extinction slopes ($R_V$) maps do not display any structure and a range of values partly due to the sampling effects of the disc by fibers, sometimes due to bad fits, and possibly partly due to some reprocessing of dust grains in the interacting disc. (abridged)

preprint2013arXiv

VLT/VIMOS Observations of an Occulting Galaxy Pair: Redshifts and Effective Extinction Curve

We present VLT/VIMOS IFU observations of an occulting galaxy pair previously discovered in HST observations. The foreground galaxy is a low-inclination spiral disk, which causes clear attenuation features seen against the bright bulge and disk of the background galaxy. We find redshifts of $z=0.064 \pm0.003$ and z=0.065 for the foreground and background galaxy respectively. This relatively small difference does not rule out gravitational interaction between the two galaxies. Emission line ratios point to a star-forming, not AGN-dominated foreground galaxy. We fit the Cardelli, Clayton & Mathis (CCM) extinction law to the spectra of individual fibres to derive slope ($R_V$) and normalization ($A_V$). The normalization agrees with the HST attenuation map and the slope is lower than the Milky Way relation ($R_V<3.1$), which is likely linked to the spatial sampling of the disk. We speculate that the values of $R_V$ point to either coherent ISM structures in the disk larger than usual ($\sim9$ kpc) or higher starting values of $R_V$, indicative of recent processing of the dust. The foreground galaxy is a low stellar mass spiral ($M_* \sim 3 \times 10^9 M_\odot$) with a high dust content ($M_{\rm dust} \sim 0.5 \times 10^6 M_\odot$). The dust disk geometry visible in the HST image would explain the observed SED properties of smaller galaxies: a lower mean dust temperature, a high dust-to-stellar mass ratio but relatively little optical attenuation. Ongoing efforts to find occulting pairs with a small foreground galaxies will show how common this geometry is.

preprint2012arXiv

Evolution in the Dust Lane Fraction of Edge-on L* Spiral Galaxies since z=0.8

The presence of a well-defined and narrow dust lane in an edge-on spiral galaxy is the observational signature of a thin and dense molecular disk, in which gravitational collapse has overcome turbulence. Using a sample of galaxies out to z~1 extracted from the COSMOS survey, we identify the fraction of massive disks that display a dust lane. Our goal is to explore the evolution in the stability of the molecular ISM disks in spiral galaxies over a cosmic timescale. We check the reliability of our morphological classifications against changes in restframe wavelength, resolution, and cosmic dimming with (artificially redshifted) images of local galaxies from SDSS. We find that the fraction of L* disks with dust lanes in COSMOS is consistent with the local fraction (~80%) out to z~0.7. At z=0.8, the dust lane fraction is only slightly lower. A somewhat lower dust lane fraction in starbursting galaxies tentatively supports the notion that a high specific star formation rate can efficiently destroy or inhibit a dense molecular disk. A small subsample of higher redshift COSMOS galaxies display low internal reddening (E[B-V]), as well as a low incidence of dust lanes. These may be disks in which the growth of the dusty ISM disk lags behind that of the stellar disk. We note that at z=0.8, the most massive galaxies display a lower dust lane fraction than lower mass galaxies. A small contribution of recent mergers or starbursts to this most massive population may be responsible. The fact that the fraction of galaxies with dust lanes in COSMOS is consistent with little or no evolution implies that models to explain the Spectral Energy Distribution or the host galaxy dust extinction of supernovae based on local galaxies are still applicable to higher redshift spirals. It also suggests that dust lanes are long lived phenomena or can be reformed over very short time-scales.

preprint2012arXiv

Herschel/SPIRE Observations of the Dusty Disk of NGC 4244

We present Herschel/SPIRE images at 250, 350, and 500 μm of NGC 4244, a typical low-mass, disk-only and edge-on spiral galaxy. The dust disk is clumpy and shows signs of truncation at the break radius of the stellar disk. This disk coincides with the densest part of the Hi disk. We compare the Spectral Energy Distribution, including the new SPIRE fluxes, to 3D radiative transfer models; a smooth model disk and a clumpy model with embedded heating. Each model requires a very high value for the dust scale-length (h(dust) = 2 - 5 h(stars)), higher dust masses than previous models of NGC 4244 (Md = 0.47 - 1.39 \times 10e7 Msun) and a face-on optical depth of τ(V) = 0.4 - 1.12, in agreement with previous disk opacity studies. The vertical scales of stars and dust are similar. The clumpy model much better mimics the general morphology in the submm images and the general SED. The inferred gas-to-dust mass ratio is compatible with those of similar low-mass disks. The relatively large radial scale-length of the dust disk points to radial mixing of the dusty ISM within the stellar disk. The large vertical dust scale and the clumpy dust distribution of our SED model are both consistent with a scenario in which the vertical structure of the ISM is dictated by the balance of turbulence and self-gravity.

preprint2012arXiv

Quantified H i Morphology VII: The Morphology of Extended Disks in UV and H i

Extended UltraViolet (xuv) disks have been found in a substantial fraction of late-type --S0, spiral and irregular-- galaxies. Similarly, most late-type spirals have an extended gas disk, observable in the 21cm radio line (HI). The morphology of galaxies can be quantified well using a series of scale-invariant parameters; Concentration- Asymmetry-Smoothness (CAS), Gini, M20, and GM parameters. In this paper, we compare the quantified morphology and effective radius (R50) of the Westerbork observations of neutral Hydrogen in Irregular and Spiral galaxies Project (WHISP) HI maps to those of far-and near-ultraviolet images obtained with galex, to explore how close the morphology and scales of HI and UV in these disks correlate. We find that xuv disks do not stand out by their effective radii in UV or HI. However, the concentration index in FUV appears to select some xuv disks. And known xuv disks can be identified via a criterion using Asymmetry and M20; 80% of xuv disks are included but with 55% contamination. This translates into 61 candidate xuv disk out of our 266 galaxies, (23%) consistent with previous findings. We consider three scenarios; tidal features from major mergers, the typical extended Hi disk is a photo- dissociation product of the xuv regions and both Hi and UV features originate in cold flows fueling the main galaxy. We define extended HI and UV disks based on their concentration (CHI > 5 and CFUV > 4 respectively), but note that these two subsamples never overlap in the WHISP sample. This appears to discount a simple photo-dissociation origin of the outer HI disk. Previously, we identified the morphology space occupied by ongoing major mergers. Known xuv disks rarely reside in the merger dominated part of HI morphology space but those that do are Type 1. This suggests cold flows as the origin for the xuv complexes and their surrounding HI structures.

preprint2012arXiv

The Brightest of Reionizing Galaxies Survey: Constraints on the Bright End of the z~8 Luminosity Function

We report the discovery of 33 Lyman break galaxy (LBG) candidates at z~8 detected in Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) imaging as part of the Brightest of Reionizing Galaxies (BoRG) pure-parallel survey. The ongoing BoRG survey currently has the largest area (274 arcmin^2) with Y_098 (or Y_105), J_125, and H_160 band coverage needed to search for z~8 galaxies, about three times the current CANDELS area, and slightly larger than what will be the final CANDELS wide component with Y_105 data. Our sample of 33 relatively bright Y_098-dropout galaxies have J_125 band magnitudes between 25.5 and 27.4 mag. This is the largest sample of bright (J_125 <~ 27.4) z~8 galaxy candidates presented to date. Combining our dataset with the Hubble Ultra-Deep Field (HUDF09) dataset, we constrain the rest-frame ultraviolet galaxy luminosity function at z~8 over the widest dynamic range currently available. The combined datasets are well fitted by a Schechter function, i.e. ϕ(L) = ϕ_* (L/L_*)^α\ e^{-(L/L_*)}, without evidence for an excess of sources at the bright end. At 68% confidence, we derive ϕ_* = (4.3^{+3.5}_{-2.1}) \times 10^{-4} Mpc^{-3}, M_* = -20.26^{+0.29}_{-0.34}, and a very steep faint-end slope α= -1.98^{+0.23}_{-0.22}. While the best-fit parameters still have a strong degeneracy, especially between ϕ_* and M_*, our improved coverage at the bright end has reduced the uncertainty of the faint-end power-law slope at z~8 compared to the best previous determination at +/-0.4. With a future expansion of the BoRG survey, combined with planned ultradeep WFC3/IR observations, it will be possible to further reduce this uncertainty and clearly demonstrate the steepening of the faint-end slope compared to measurements at lower redshift, thereby confirming the key role played by small galaxies in the reionization of the universe.

preprint2012arXiv

The Opacity of Spiral Galaxy Disks IX; Dust and Gas Surface Densities

Our aim is to explore the relation between gas, atomic and molecular, and dust in spiral galaxies. Gas surface densities are from atomic hydrogen and CO line emission maps. To estimate the dust content, we use the disk opacity as inferred from the number of distant galaxies identified in twelve HST/WFPC2 fields of ten nearby spiral galaxies. The observed number of distant galaxies is calibrated for source confusion and crowding with artificial galaxy counts and here we verify our results with sub-mm surface brightnesses from archival Herschel-SPIRE data. We find that the opacity of the spiral disk does not correlate well with the surface density of atomic (Hi) or molecular hydrogen (H2) alone implying that dust is not only associated with the molecular clouds but also the diffuse atomic disk in these galaxies. Our result is a typical dust-to-gas ratio of 0.04, with some evidence that this ratio declines with galactocentric radius, consistent with recent Herschel results. We discuss the possible causes of this high dust-to-gas ratio; an over-estimate of the dust surface-density, an under-estimate of the molecular hydrogen density from CO maps or a combination of both. We note that while our value of the mean dust-to-gas ratio is high, it is consistent with the metallicity at the measured radii if one assumes the Pilyugin & Thuan calibration of gas metallicity.

preprint2011arXiv

Looking At the Distant Universe with the MeerKAT Array (LADUMA)

The MeerKAT (64 x 13.5m dish radio interferometer) is South Africa's precursor instrument for the Square Kilometre Array (SKA), exploring dish design, instrumentation, and the characteristics of a Karoo desert site and is projected to be on sky in 2016. One of two top-priority, Key Projects is a single deep field, integrating for 5000 hours total with the aim to detect neutral atomic hydrogen through its 21 cm line emission out to redshift unity and beyond. This first truly deep HI survey will help constrain fueling models for galaxy assembly and evolution. It will measure the evolution of the cosmic neutral gas density and its distribution over galaxies over cosmic time, explore evolution of the gas in galaxies, measure the Tully-Fisher relation, measure OH maser counts, and address many more topics. Here we present the observing strategy and envisaged science case for this unique deep field, which encompasses the Chandra Deep Field-South (and the footprints of GOODS, GEMS and several other surveys) to produce a singular legacy multi-wavelength data-set.

preprint2011arXiv

New HErschel Multi-wavelength Extragalactic Survey of Edge-on Spirals (NHEMESES)

Edge-on spiral galaxies offer a unique perspective on the vertical structure of spiral disks, both stars and the iconic dark dustlanes. The thickness of these dustlanes can now be resolved for the first time with Herschel in far-infrared and sub-mm emission. We present NHEMESES, an ongoing project that targets 12 edge-on spiral galaxies with the PACS and SPIRE instruments on Herschel. These vertically resolved observations of edge-on spirals will impact on several current topics. First and foremost, these Herschel observations will settle whether or not there is a phase change in the vertical structure of the ISM with disk mass. Previously, a dramatic change in dustlane morphology was observed as in massive disks the dust collapses into a thin lane. If this is the case, the vertical balance between turbulence and gravity dictates the ISM structure and consequently star-formation and related phenomena (spiral arms, bars etc.). We specifically target lower mass nearby edge-ons to complement existing Herschel observations of high-mass edge-on spirals (the HEROES project). Secondly, the combined data-set, together with existing Spitzer observations, will drive a new generation of spiral disk Spectral Energy Distribution models. These model how dust reprocesses starlight to thermal emission but the dust geometry remains the critical unknown. And thirdly, the observations will provide an accurate and unbiased census of the cold dusty structures occasionally seen extending out of the plane of the disk, when backlit by the stellar disk. To illustrate the NHEMESES project, we present early results on NGC 4244 and NGC 891, two well studies examples of a low and high-mass edge-on spiral.

preprint2011arXiv

Quantified HI Morphology I: Multi-Wavelengths Analysis of the THINGS Galaxies

Galaxy evolution is driven to a large extent by interactions and mergers with other galaxies and the gas in galaxies is extremely sensitive to the interactions. One method to measure such interactions uses the quantified morphology of galaxy images. Well-established parameters are Concentration, Asymmetry, Smoothness, Gini, and M20 of a galaxy image. Thus far, the application of this technique has mostly been restricted to restframe ultra-violet and optical images. However, with the new radio observatories being commissioned (MeerKAT, ASKAP, EVLA, WSRT/APERTIF, and ultimately SKA), a new window on the neutral atomic hydrogen gas (HI) morphology of a large numbers of galaxies will open up. The quantified morphology of gas disks of spirals can be an alternative indicator of the level and frequency of interaction. The HI in galaxies is typically spatially more extended and more sensitive to low-mass or weak interactions. In this paper, we explore six morphological parameters calculated over the extent of the stellar (optical) disk and the extent of the gas disk for a range of wavelengths spanning UV, Optical, Near- and Far-Infrared and 21 cm (HI) of 28 galaxies from The HI Nearby Galaxy Survey (THINGS). Though the THINGS sample is small and contains only a single ongoing interaction, it spans both non-interacting and post-interacting galaxies with a wealth of multi-wavelength data. We find that the choice of area for the computation of the morphological parameters is less of an issue than the wavelength at which they are measured. The signal of interaction is as good in the HI as in any of the other wavelengths in which morphology has been used to trace the interaction rate to date, mostly star-formation dominated ones (near- and far-ultraviolet). The Asymmetry and M20 parameters are the ones which show the most promise as tracers of interaction in 21 cm line observations.

preprint2011arXiv

Quantified HI Morphology II : Lopsidedness and Interaction in WHISP Column Density Maps

Lopsidedness of the gaseous disk of spiral galaxies is a common phenomenon in disk morphology, profile and kinematics. Simultaneously, the asymmetry of a galaxy's stellar disk, in combination with other morphological parameters, has seen extensive use as an indication of recent merger or interaction in galaxy samples. Quantified morphology of stellar spiral disks is one avenue to determine the merger rate over much of the age of the Universe. In this paper, we measure the quantitative morphology parameters for the HI column density maps from the Westerbork observations of neutral Hydrogen in Irregular and SPiral galaxies (WHISP). These are Concentration, Asymmetry, Smoothness, Gini, M20, and one addition of our own, the Gini parameter of the second order moment (GM). Our aim is to determine if lopsided or interacting disks can be identified with these parameters. Our sample of 141 HI maps have all previous classifications on their lopsidedness and interaction. We find that the Asymmetry, M20, and our new GM parameter correlate only weakly with the previous morphological lopsidedness quantification. These three parameters may be used to compute a probability that an HI disk is morphologically lopsided but not unequivocally to determine it. However, we do find that that the question whether or not an HI disk is interacting can be settled well using morphological parameters. Parameter cuts from the literature do not translate from ultraviolet to HI directly but new selection criteria using combinations of Asymmetry and M20 or Concentration and M20, work very well. We suggest that future all-sky HI surveys may use these parameters of the column density maps to determine the merger fraction and hence rate in the local Universe with a high degree of accuracy.

preprint2011arXiv

Quantified HI Morphology III: Merger Visibility Times from HI in Galaxy Simulations

Major mergers of disk galaxies are thought to be a substantial driver in galaxy evolution. To trace the fraction and the rate galaxies are in mergers over cosmic times, several observational techniques, including morphological selection criteria, have been developed over the last decade. We apply this morphological selection of mergers to 21 cm radio emission line (HI) column density images of spiral galaxies in nearby surveys. In this paper, we investigate how long a 1:1 merger is visible in HI from N-body simulations. We evaluate the merger visibility times for selection criteria based on four parameters: Concentration, Asymmetry, M20, and the Gini parameter of second order moment of the flux distribution (GM). Of three selection criteria used in the literature, one based on Concentration and M20 works well for the HI perspective with a merger time scale of 0.4 Gyr. Of the three selection criteria defined in our previous paper, the GM performs well and cleanly selects mergers for 0.69 Gyr. The other two criteria (A-M20 and C-M20), select isolated disks as well, but perform best for face-on, gas-rich disks (T(merger) ~ 1 Gyr). The different visibility scales can be combined with the selected fractions of galaxies in any large HI survey to obtain merger rates in the nearby Universe. All-sky surveys such as WALLABY with ASKAP and the Medium Deep Survey with the APETIF instrument on Westerbork are set to revolutionize our perspective on neutral hydrogen and will provide an accurate measure of the merger fraction and rate of the present epoch.

preprint2011arXiv

Quantified HI Morphology IV: The Merger Fraction and Rate in WHISP

The morphology of the atomic hydrogen (HI) disk of a spiral galaxy is the first component to be disturbed by a gravitational interaction such as a merger between two galaxies. We use a simple parametrisation of the morphology of HI column density maps of Westerbork HI Spiral Project (WHISP) to select those galaxies that are likely undergoing a significant interaction. Merging galaxies occupy a particular part of parameter space defined by Asymmetry (A), the relative contribution of the 20% brightest pixels to the second order moment of the column density map (M20) and the distribution of the second order moment over all the pixels (GM). Based on their HI morphology, we find that 13% of the WHISP galaxies are in an interaction (Concentration-M20) and only 7% based on close companions in the data-cube. This apparent discrepancy can be attributed to the difference in visibility time scales: mergers are identifiable as close pairs for 0.5 Gyr but ~1 Gyr by their disturbed HI morphology. Expressed as volume merger rates, the two estimates agree very well: 7 and 6.8 x 10^-3 mergers Gyr^-1 Mpc^-3 for paired and morphologically disturbed HI disks respectively. The consistency of our merger fractions to those published for bigger surveys such as the Sloan Digital Sky Survey, shows that HI morphology can be a very viable way to identify mergers in a large HI survey. The relatively high value for the volume merger rate may be a bias in the selection or WHISP volume. The expected boon in high-resolution HI data by the planned MeerKAT, ASKAP and WSRT/APERTIF radio observatories will reveal the importance of mergers in the local Universe and, with the advent of SKA, over cosmic times.

preprint2011arXiv

Quantified HI Morphology V: HI Disks in the Virgo Cluster

We explore the quantified morphology of atomic hydrogen (HI) disks in the Virgo cluster. These galaxies display a wealth of phenomena in their Hi morphology, e.g., tails, truncation and warps. These morphological disturbances are related to the ram-pressure stripping and tidal interaction that galaxies undergo in this dense cluster environment. To quantify the morphological transformation of the HI disks, we compute the morphological parameters of CAS, Gini, and M20 and our own GM for 51 galaxies in 48 HI column density maps from the VIVA project. Some morphological phenomena can be identified in this space of relatively low resolution HI data. Truncation of the HI disk can be cleanly identified via the Concentration parameter (C<1) and Concentration can also be used to identify HI deficient disks (1<C<5). Tidal interaction is typically identified using combinations of these morphological parameters, applied to (optical) images of galaxies. We find that some selection criteria (Gini-M20, Asymmetry, and a modified Concentration-M20) are still applicable for the coarse (~15" FWHM) VIVA HI data. The phenomena of tidal tails can be reasonably well identified using the Gini-M20 criterion (60% of galaxies with tails identified but with as many contaminants). Ram-pressure does move HI disks into and out of most of our interaction criteria: the ram-pressure sequence identified by Vollmer et al. (2009) tracks into and out of some of these criteria (Asymmetry based and the Gini-M20 selections, but not the Concentration-M20 or the GM based ones). Therefore, future searches for interaction using HI morphologies should take ram-pressure into account as a mechanism to disturb HI disks enough to make them appear as gravitationally interacting. One mechanism would be to remove all the HI deficient (C<5) disks from the sample, as these have undergone more than one HI removal mechanism.

preprint2010arXiv

Trumpeting the Vuvuzela: UltraDeep HI observations with MeerKAT

The MeerKAT UltraDeep HI Survey aims to observe the 21 cm emission line of neutral hydrogen gas out to a redshift of z=1 and beyond. From both direct detections and stacked signal, we will address the HI mass function, the cosmic neutral gas density of the Universe (Omega_HI) and their evolution over cosmic times, as well as galaxy evolution via e.g., the Tully-Fisher relation, the relation between HI mass and Hubble Type or stellar mass, and the Schmidt-Kennicutt star-formation law. We propose to observe two fields, the COSMOS and Chandra Deep Field South (CDF-S) for 1000 hours each, adding an additional 4000 hours to one of these fields in 2015 when the full instantaneous bandwidth of MeerKAT (0.58-2.5 GHz) will be realised.

preprint2009arXiv

Spectrophotometrically Identified stars in the PEARS-N and PEARS-S fields

Deep ACS slitless grism observations and identification of stellar sources are presented within the Great Observatories Origins Deep Survey (GOODS) North and South fields which were obtained in the Probing Evolution And Reionization Spectroscopically (PEARS) program. It is demonstrated that even low resolution spectra can be a very powerful means to identify stars in the field, especially low mass stars with stellar types M0 and later. The PEARS fields lay within the larger GOODS fields, and we used new, deeper images to further refine the selection of stars in the PEARS field, down to a magnitude of mz = 25 using a newly developed stellarity parameter. The total number of stars with reliable spectroscopic and morphological identification was 95 and 108 in the north and south fields respectively. The sample of spectroscopically identified stars allows constraints to be set on the thickness of the Galactic thin disk as well as contributions from a thick disk and a halo component. We derive a thin disk scale height, as traced by the population of M4 to M9 dwarfs along two independent lines of sight, of h_thin = 370 +60/-65 pc. When including the more massive M0 to M4 dwarf population, we derive h_thin = 300 +/- 70pc. In both cases, we observe that we must include a combination of thick and halo components in our models in order to account for the observed numbers of faint dwarfs. The required thick disk scale height is typically h_thick=1000 pc and the acceptable relative stellar densities of the thin disk to thick disk and the thin disk to halo components are in the range of 0.00025<f_halo<0.0005 and 0.05<f_thick<0.08 and are somewhat dependent on whether the more massive M0 to M4 dwarfs are included in our sample.