Quantified Morphology of HI Disks in the Universe

he upcoming new perspective of the high redshift Universe in the 21 cm line of atomic hydrogen opens possibilities to explore topics of spiral disk evolution, hitherto reserved for the optical regime. The growth of spiral gas disks over Cosmic time can be explored with the new generation of radio telescopes, notably the SKA, and its precursors, as accurately as with the Hubble Space Telescope for stellar disks. Since the atomic hydrogen gas is the building block of these disks, it should trace their formation accurately. Morphology of HI disks can now equally be quantified over Cosmic time. In studies of HST deep fields, the optical or UV morphology of high-redshift galaxy disks have been characterized using a few quantities: concentration (C), asymmetry (A), smoothness (S), second-order-moment (M20), the GINI coefficient (G), and Ellipticity (E). We have applied these parameters across wavelengths and compared them to the HI morphology over the THINGS sample. NGC 3184, an unperturbed disk, and NGC 5194, the canonical 3:1 interaction, serve as examples for quantified morphology. We find that morphology parameters determined in HI are as good or better a tracer of interaction compared to those in any other wavelength, notably in Asymmetry, Gini and M20. This opens the possibility of using them in the parameterization pipeline for SKA precursor catalogues to select interacting or harassed galaxies from their HI morphology. Asymmetry, Gini and M20 may be redefined for use on data-cubes rather than HI column density image.Comment: 6 pages, 3 figures, proceeding of the conference "Panoramic Radio Astronomy: Wide-field 1-2 GHz research on galaxy evolution", June 02 - 05 2009, Groningen, update after small edit

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.Comment: 10 pages, 3 figures, accepted by MNRAS, appendixes not include

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.Comment: 26 figures, 4 tables, two appendices. Third appendix (HI maps of all galaxies) omitted. Accepted by MNRA

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.Comment: 12 pages, 6 figures, 4 tables, accepted by MNRAS, appendix not include

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.Comment: 12 pages, 5 figures, 1 table, accepted by MNRAS, appendix not include

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)Comment: 39 pages, 38 figures, Accepted for publication in Astronomy & Astrophysic

High-precision Photometric Redshifts from Spitzer/IRAC: Extreme [3.6]-[4.5] Colors Identify Galaxies in the Redshift Range z~6.6-6.9

One of the most challenging aspects of studying galaxies in the z>~7 universe is the infrequent confirmation of their redshifts through spectroscopy, a phenomenon thought to occur from the increasing opacity of the intergalactic medium to Lya photons at z>6.5. The resulting redshift uncertainties inhibit the efficient search for [C II] in z~7 galaxies with sub-mm instruments such as ALMA, given their limited scan speed for faint lines. One means by which to improve the precision of the inferred redshifts is to exploit the potential impact of strong nebular emission lines on the colors of z~4-8 galaxies as observed by Spitzer/IRAC. At z~6.8, galaxies exhibit IRAC colors as blue as [3.6]-[4.5] ~-1, likely due to the contribution of [O III]+Hb to the 3.6 mum flux combined with the absence of line contamination in the 4.5 mum band. In this paper we explore the use of extremely blue [3.6]-[4.5] colors to identify galaxies in the narrow redshift window z~6.6-6.9. When combined with an I-dropout criterion, we demonstrate that we can plausibly select a relatively clean sample of z~6.8 galaxies. Through a systematic application of this selection technique to our catalogs from all five CANDELS fields, we identify 20 probable z~6.6-6.9 galaxies. We estimate that our criteria select the ~50% strongest line emitters at z~6.8 and from the IRAC colors we estimate a typical [O III]+Hb rest-frame equivalent width of 1085A for this sample. The small redshift uncertainties on our sample make it particularly well suited for follow-up studies with facilities such as ALMA.Comment: In submission to the Astrophysical Journal, updated in response to the referee report, 13 pages, 11 figures, 1 tabl

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 {\mu}m at the 5{\sigma} 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.Comment: 21 pages, 12 figures. Accepted for publication in Astronomy & Astrophysic

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.Comment: 16 pages, 11 figure, table 1, accepted by MNRAS, appendix not include

The Bright End of the z~9 and z~10 UV Luminosity Functions using all five CANDELS Fields

The deep, wide-area (~800-900 arcmin**2) near-infrared/WFC3/IR + Spitzer/IRAC observations over the CANDELS fields have been a remarkable resource for constraining the bright end of high redshift UV luminosity functions (LFs). However, the lack of HST 1.05-micron observations over the CANDELS fields has made it difficult to identify z~9-10 sources robustly, since such data are needed to confirm the presence of an abrupt Lyman break at 1.2 microns. We report here on the successful identification of many such z~9-10 sources from a new HST program (z9-CANDELS) that targets the highest-probability z~9-10 galaxy candidates with observations at 1.05 microns, to search for a robust Lyman-break at 1.2 microns. The potential z~9-10 candidates are preselected from the full HST, Spitzer/IRAC S-CANDELS observations, and the deepest-available ground-based optical+near-infrared observations. We identified 15 credible z~9-10 galaxies over the CANDELS fields. Nine of these galaxies lie at z~9 and 5 are new identifications. Our targeted follow-up strategy has proven to be very efficient in making use of scarce HST time to secure a reliable sample of z~9-10 galaxies. Through extensive simulations, we replicate the selection process for our sample (both the preselection and follow-up) and use it to improve current estimates for the volume density of bright z~9 and z~10 galaxies. The volume densities we find are 5(-2)(+3)x and 8(-3)(+9)x lower, respectively, than found at z~8. When compared with the best-fit evolution (i.e., dlog_{10} rho(UV)/dz=-0.29+/-0.02) in the UV luminosities densities from z~8 to z~4 integrated to 0.3L*(z=3) (-20 mag), these luminosity densities are 2.6(-0.9)(+1.5)x and 2.2(-1.1)(+2.0)x lower, respectively, than the extrapolated trends. Our new results are broadly consistent with the "accelerated evolution" scenario at z>8, as seen in many theoretical models.Comment: 23 pages, 15 figures, 7 tables, updated to match the version in press, including some minor textual corrections identified at the proof stag