Open Access
Issue
A&A
Volume 635, March 2020
Article Number A69
Number of page(s) 21
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201936556
Published online 10 March 2020

© K. A. Lutz et al. 2020

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

The study of outliers to scaling relations may inform of the physical processes that underlay the scaling relation. In the case of the relations between atomic hydrogen (H I) and stellar content of galaxies, the study of galaxies that are more H I-rich than expected from the stellar or optical properties helps understand how galaxies acquire and maintain their H I content. Possible scenarios are more efficient gas accretion or less efficient star formation than in average galaxies.

Previous surveys of H I-rich galaxies include HIghMass (Huang et al. 2014), BLUEDISK (Wang et al. 2013), H I monsters (Lee et al. 2014), H I excess galaxies (Geréb et al. 2016, 2018), and a sample of H I-rich, massive galaxies (Lemonias et al. 2014). While the HIghMass, H I monsters, and Lemonias et al. (2014) samples investigate generally similar galaxies, BLUEDISK galaxies are systematically less H I-rich than H IX galaxies (see also Lutz et al. 2017). H I-excess systems differ from other H I-rich galaxies as they were selected for being H I-rich for their specific star formation rates. For five HIghMass and all Lemonias et al. (2014) galaxies resolved H I observations were observed, analysed, and published. Out of the five HIghMass galaxies, three galaxies have a higher than average spin parameter, one galaxy is on the brink of a star burst, and one galaxy (a group central) accreted gas from its neighbouring galaxies (Hallenbeck et al. 2014, 2016). The H I column densities in the galaxies of Lemonias et al. (2014) are very low and thus star formation is suppressed. So far the gas-phase metallicity has not been taken into account to investigate possible gas accretion onto H I-rich galaxies.

When comparing the current cold gas content (atomic and molecular hydrogen, i.e. H I and H2) of local star forming galaxies to their star formation activity, Saintonge et al. (2017), Schiminovich et al. (2010), and others have found that these galaxies would use their entire cold gas (H I plus H2) content, on average within 4 Gyr. Hence, for galaxies to remain star formers in the future, they need to replenish their cold gas reservoir. Currently two main avenues of cold gas accretion are commonly considered: gas accretion from the intergalactic medium (Birnboim & Dekel 2003; Dekel & Birnboim 2006; Kereš et al. 2005; van de Voort et al. 2011) or the infall of satellite galaxies, which bring their own cold gas reservoir with them (Di Teodoro & Fraternali 2014; van de Voort et al. 2011). Additionally, the so-called Galactic Fountain could also increase the cold gas content of a galaxy in the following way: through stellar feedback gas that is expelled from the disc into the halo, where hot halo gas condenses onto the expelled parcels of gas and together they rain back onto the galaxy disc (Oosterloo et al. 2007; Fraternali et al. 2007). Furthermore, this mechanism redistributes gas within the galaxy and its halo.

The gas-phase metallicity (12 + log(O/H)) is an important indicator of the mode of accretion. One line of argument that suggests that galaxies should accrete metal-poor gas from the intergalactic medium (IGM) is the comparison of chemical evolution models of galaxies to observationsof the gas-phase metallicity in galaxies. When assuming a closed box model, where the interstellar medium (ISM) of galaxies is not diluted by pristine gas from the IGM, metallicities of modelled galaxies are overestimated with respect to observed metallicities (van den Bergh 1962; Kudritzki et al. 2015). Hence, galaxies need to accrete relatively metal-poor gas from the IGM to dilute their ISM.

An important observational probe into the chemical evolution of galaxies is the mass–metallicity relation. This relation connects thestellar mass of a galaxy to the central gas-phase metallicity (Tremonti et al. 2004). It was shown that the scatter of the mass–metallicity relation is correlated with the star formation rate of galaxies (SFR) (Mannucci et al. 2010) or the H I content of galaxies (Hughes et al. 2013; Bothwell et al. 2013; Lagos et al. 2016; Brown et al. 2018). Growing evidence now points to H I being the primary driver and SFR being a secondary effect due to the dependence of star formation on gas (e.g. Lagos et al. 2016, Brown et al. 2018).

N-body, smoothed particle hydrodynamical simulations (e.g. Davé et al. 2013), cosmological simulations (e.g. Torrey et al. 2019), and simple models (e.g. Forbes et al. 2014) have found that the mass–metallicity relation can be explained if galaxies are seen as systems in equilibrium between gas accretion, star formation, and outflows. Stochastic accretion of pristine gas removes galaxies from the equilibrium: the gas content is increased, the ISM diluted, and star formation triggered. Thus, the galaxy initially moves from the centre of the mass–metallicity relation to a more metal-poor position for its stellar mass. With the triggered star formation, the galaxy now grows in stellar mass and its ISM is gradually enriched again (through feedback). Thus, it moves back to the equilibrium line of the mass–metallicity relation.

The observed mass–metallicity relation was originally defined using only the central metallicities of galaxies. For example the Tremonti et al. (2004) relation is based on SDSS spectra, which were observed with a 3 arcsec fibre. Hence, only the metallicity of the central 3arcsec is included in their relation. However, the metallicity of the disc and outskirts of the galaxy also provide vital information on the chemical evolution. Recently, a number of large samples of galaxies have been observed with long-slit or integral field spectra, examples are the Sydney-AAO Multi-object Integral field spectrograph galaxy survey (SAMI survey; Croom et al. 2012; Fogarty et al. 2012), the Calar Alto Legacy Integral Field spectroscopy Area survey (CALIFA survey; Sánchez et al. 2012), the Mapping Nearby Galaxies at APO survey (MaNGA survey; Bundy et al. 2015), or follow-up of galaxies (Moran et al. 2012) in the GALEX Arecibo SDSS survey (GASS; Catinella et al. 2010, 2013). Using these data, a metallicity gradient (i.e. the gas-phase metallicity change with radius) can be measured. The gas-phase metallicity in regular spiral galaxies decreases towards the edges of the disc. In the CALIFA survey, Sánchez et al. (2014) found a universal metallicity gradient in the sense that the steepness of the gradient is not dependent on the stellar mass1. They suggested that this result indicates a similar chemical evolution in all disc galaxies. Hence, the metallicity at a given radius is set by the evolutionary state at that radius rather than the basic properties of the galaxy. As galaxies form from the inside out, the outskirts of a galaxy would be less evolved and thus more metal-poor. These findings were confirmed by Ho et al. (2015) and Kudritzki et al. (2015), who were able to reproduce the uniform metallicity gradients with a simple chemical evolution model including in- and outflows. These models indicate that the metallicity is dependent on the local stellar-to-gas mass ratio, which has also been observed in MaNGA galaxies (Barrera-Ballesteros et al. 2018) (only using estimates of gas mass from dust attenuation). Other results based on the MaNGA survey, show a metallicity gradient that varies with stellar mass (Belfiore et al. 2017). They suggest that flatter gradients in low mass galaxies show that strong feedback, gas mixing, and wind recycling must also be important in these galaxies (i.e. that a mechanism like the Galactic Fountain flattens their metallicity gradient).

Conversely, Moran et al. (2012) measured the radial metallicity gradient in massive galaxies from the GASS survey and found the steepest declining metallicity gradients in galaxies at their lower stellar mass limit of log M [M] > 10. They furthermore found that about 10% of their sample have a sharp downturn in metallicity at large radii. This decline is correlated with the H I content of these galaxies and was interpreted as a sign for a phase of active gas inflow and disc-building. Similarly, Sánchez Almeida et al. (2014a) reported some star forming regions in dwarf galaxies to be more metal-poor than their surroundings. They and Ceverino et al. (2016) subsequently suggested that these could be star forming regions that formed in parcels of recently accreted metal-poor gas.

To date not many samples in addition to the GASS sample and the sample in this work have actually measured gas properties. In particular, the overlap between large integral field unit (IFU) surveys (SAMI, MaNGA, and CALIFA) and H I surveys is still limited for single-dish surveys and close to non-existent for interferometric (i.e. spatially resolved) H I data.

We previously compiled a sample of galaxies that contain at least 2.5 times more H I than expected from their optical R-band luminosity. For these H I eXtreme galaxies (H IX galaxies) and a CONTROL sample, we analysed the star formation activity (Lutz et al. 2017) and the H I kinematic properties (Lutz et al. 2018). We found that H IX galaxies maintain their H I reservoir by retaining large amounts of H I outside the stellar disc. There the gas cannot move towards central, denser regions and is stabilised against star formation due to its high specific angular momentum, which is likely inherited from a high spin halo. To investigate possible accretion of metal-poor gas or recent mergers with the H IX galaxies, we acquired optical integral field spectra within the stellar discs of H IX galaxies. The data presented in this paper allow us to compare spatially resolved gas-phase oxygen abundance measurements to spatially resolved H I observations as well as the kinematics of different galaxy components.

In addition to the high angular momentum scenario, another scenario that might explain the H I-richness of H IX galaxies is very effective or active gas accretion, which increases the H I content of H IX galaxies. This might either be a temporarily high accretion rate of gas from the IGM or a recent gas-rich merger. To probe these scenarios, this paper compares Hα and stellar kinematics to the H I kinematics. Any misalignments in the rotation of these components can point towards recent mergers (Corsini 2014). We note, however, that the lack of misalignments cannot exclude recent merger or accretion events. Furthermore, the central gas-phase metallicity of H IX galaxies is examined using the mass–metallicity relation and the steepness of the gas-phase metallicity gradient in H IX galaxies is compared to the literature. Both measures can also inform of recent inflow of pristine gas.

This paper is structured as follows: in Sect. 2, details of observations and the data analysis are described. In Sect. 3, results from the analysis of the optical integral field spectra (IFS) data are presented. These results are discussed and summarised in Sect. 4.

2 Galaxy samples and data analysis

2.1 The H IX galaxy survey: sample selection

The sample selection was described in detail in Lutz et al. (2017); here we give a brief summary: H IX galaxies are selected from the H I Parkes All Sky Survey (HIPASS, Barnes et al. 2001). In particular, we used a high quality, high fidelity catalogue of galaxies with both H I mass measurement and optical photometry (Dénes et al. 2014), which was based on the HIPASS catalogues (Barnes et al. 2001; Meyer et al. 2004; Zwaan et al. 2004; Koribalski et al. 2004) and the optical counterparts in the HOPCAT catalogue (Doyle et al. 2005).

Based on the Dénes et al. (2014) scaling relation between H I mass and R-band luminosity, H IX galaxies were selected to fulfil the following criteria:

  • they host at least 2.5 times more H I than expected from their R-band luminosity, i.e. they lie at least 1.4 σ above the scaling relation;

  • they are located south of Dec <−30 deg for good observability with the Australian Telescope Compact Array;

  • they are brighter in absolute Ks-band magnitude than −22 mag to restrict the sample to massive spiral galaxies.

For comparison, a CONTROL sample was compiled from the same catalogue with the same selection criteria except that galaxies should have between 1.6 times less and more H I than expected from the Dénes et al. (2014) R-band scaling relation (i.e. within ± 0.7σ).

An overview of basic properties of H IX and CONTROL galaxies examined in this paper is given in Table 1.

2.2 Observations of H IX galaxies with WiFeS

We used the WIde FiEld Spectrograph (WIFES, Dopita et al. 2007) on the ANU 2.3 m telescope in Siding Spring, Australia to obtain optical integral field spectra of 10 out of 12 H IX galaxies and one CONTROL galaxy.

WIFES is an image slicing integral field unit, meaning that it consists of 25 slitlets, each 1 × 36 arcsec in size. Combining the length of the slitlets and their number, WIFES has a 36 × 25 arcsec2 field of view. This is smaller than the optical disc of the H IX sample galaxies, which have an average 25 mag arcsec−2 isophote radius of 47 arcsec. Therefore, the aim was to observe every galaxy with multiple pointings to cover the centre and the outer regions of the stellar disc. Figure 1 shows an example of multiple pointings towards ESO111-G014. In the left panel the WIFES pointings are shown, while the right panel details the final binning of the pointings.

The aimwas to to obtain 60 min (90 min) of on-source time per pointing for galaxies with average surface brightness brighter (fainter) than 22 mag arcsec−2. The actual total on-source exposure times of each pointing are between 15 and 90 min depending on the signal strength. The total on-source time per pointing is split into single exposures of 15 min. Stacking multiple short exposures helps to remove cosmic rays and small image errors. To subtract the night sky foreground, all science spectra are taken in the so-called nod-and-shuffle mode, in which the telescope nods between the science target (i.e. the galaxy) and an “empty” part of the sky. This way a spectrum of the sky and a sum of the sky spectrum and the galaxy is obtained. The separate sky spectrum is then subtracted from the intermingled sky and galaxy spectrum to obtain a pure galaxy spectrum.

Incoming spectra were split by a dichroic at 560 nm into a red and a blue half. Thus, a wide wavelength range can be covered in one shot. WIFES wasset up with R = 3000 gratings in the red and the blue arm for all observations. At a wavelength of 660 nm this is equivalent to a wavelength step of 0.22 nm per pixel or a velocity stepof 100 km s−1.

Every night of observations was completed with standard calibration images including bias, sky flat field, dome flat field, wire imaging for centring the slitlets, and NeAr and CuAr arc lamp spectra for wavelength calibration. Two to three times a night a standard star was observed for flux calibration. Dark current information was obtained from overscan regions (i.e. parts of the CCD that were not exposed). The observations were conducted between August 2014 and April 2016 under photometric conditions. Details about the number of pointings and on-source exposure times are given in Table 2.

2.3 Data reduction of WiFeS observations

Childress et al. (2014) provide the fully automated PYWIFES pipeline for WIFES data. This pipeline includes bad pixel repair, bias and dark current subtraction, flat fielding, wavelength calibration, sky subtraction, flux calibration, and data cube creation. The resulting data cubes are 70 × 25 pixels (35 × 25 arcsec) in size and their wavelength range is set from 650.0 to 685.0 nm and 460.0 to 525.0 nm for the red and blue cubes, respectively. We note that these measures imply a pixel size of 1 arcsec in one spatial direction and 0.5 arcsec in the other spatial direction.

Single galaxy exposures are spatially offset from each other. To account for this, PYWIFES was run on each observed galaxy exposure individually. For each pointing, data cubes were then median stacked after the full data reduction by PYWIFES. To align single data cubes of one pointing, the distributions of Hα emission in the data cubes were compared. For an example see Fig. 2. The locations of bright centres of star forming regions were matched up between the different exposures of one pointing. Once the pixel offset between the single data cubes per pointing was determined, the cubes were median stacked.

After stacking, cubes were reshaped such that both spatial dimensions have pixel sizes of 1 arcsec. For this task the SCIPY (Jones et al. 2001) function scipy.interpolate.griddata was used. Data were linearly interpolated to the new pixel size of 1 × 1 arcsec.

Initially, observations of this WIFES project do not include any world coordinate system (WCS) information in their file headers. To add this information in the file headers the following procedure was applied: right before the WIFES observations were taken, an image with the acquisition and guiding (A&G) camera was taken. In the centre of this image the shadow of the WIFES IFU is visible. We applied an astrometric solution to this A&G camera image using bright stars and the Aladin software (Bonnarel et al. 2000). Since the pixel size of the IFU cube (1 arcsec) and the location of the IFU within the A&G camera image (due to its shadow) are both known, we then simply transferred the astrometric solution from the A&G camera image to the IFU cubes. After this procedure, we measured the position of prominent features (e.g. peaks of star formation, stars) in collapsed maps of the WiFeS cubes and compared them to the positions of these features measured on the SuperCOSMOS Bj -band images. We find an average distance between the two position measurements of (3 ± 6) arcsec.

This procedure resulted in two data cubes for each pointing: one containing the red half of the spectrum and the other the blue half. To increase the signal-to-noise ratio (S/N) of spectra, the spaxels in both cubes were binned such that one bin encompasses a star forming region and the Hα and N[II] λ 658.3 nm emission lines are visible in the resulting spectrum. These star forming regions were defined by hand on the Hα emission planes of the WIFES data cubes to be a connected region of Hα emission (see Fig. 1).

In pointings towards the edges of stellar discs, the one-dimensional spectra of individual star forming regions do not pick up emission from the stellar continuum. In those cases the background was estimated as a constant from emission and sky line free parts of the spectrum (Wisnioski et al. 2015). For our data set we chose the following wavelength ranges: indicating with λα and λN the redshifted wavelengths of Hα and N[II] λ 658.3 nm lines, respectively, we chose the red window as [λN + 2 nm, λN + 7 nm] if λα < 665 nm, and [λα − 8 nm, λα − 3 nm] otherwise. The blue window is always [λβ − 8 nm, λβ − 3 nm], where λβ is the redshifted wavelength of the Hβ line. We then took all the available data within these windows and calculated the 40th and 60th percentiles. The mean of the data in the 40th to 60th percentile range was then set as the constant background value for those spectra that do not show any stellar continuum.

If stellar continuum is detected in a spectrum of a star forming region, then the full spectrum (i.e. the blue and the red half) was used to fit a stellar population synthesis model and gas emission lines with the penalized pixel-fitting (pPXF) method and the PYTHON script by Cappellari & Emsellem (2004) and Cappellari (2017). pPXF used the stellar population synthesis models from the Vazdekis et al. (2010) library plus Gaussian emission lines. From the location of stellar absorption features, the stellar recession velocity in this star forming region was measured. pPXF provides the line strength for each modelled emission line (Hα, Hβ, O[III] λ 500.7 nm, and N[II] λ 658.3 nm). For the recession velocity of the gas it was assumed that there is only one kinematic component, and the recession velocity was modelled from all four emission lines simultaneously (i.e. one recession velocity measurement from the four emission lines).

The emission line fluxes were corrected for internal dust attenuation assuming the Calzetti et al. (2000) attenuation curve and RV=4.05$R&#x0027;_{\textrm{V}}=4.05$ (see also Moran et al. 2012). From the corrected Hα, Hβ, O[III] λ 500.7 nm, and N[II] λ 658.3 nm emission line fluxes, the gas-phase oxygen abundances were determined with the O 3N 2 and N2 method as described by Pettini & Pagel (2004).

We tested several options and found for spectra, in which the stellar continuum could be modelled with pPXF, that metallicity measurements based on the subtraction of a constant background and the N2 method agree with metallicities based on subtraction of a pPXF modelled stellar continuum and the O 3N 2 method. Hence, wherever possible we use O3N2 metallicities calculated from pPXF modelled emission line fluxes, otherwise the N2 method with the emission lines measured after subtraction of a constant background.

To de-project the on-sky galactocentric radius r, first the angle θ between the galaxy semi-major axis and the line between galaxy centre and star forming region was determined. The de-projected radius rdpj was then set to rdpj=r×(bcosθ)2+(asinθ)2b,\begin{equation*} r_{\textrm{dpj}} = \frac{r \times \sqrt{(b \cos \theta)^2 &#x002B; (a \sin \theta)^2}} {b}, \end{equation*}(1)

where a and b are the semi-major and semi-minor axis, respectively. Metallicity gradients were measured by fitting a line to the metallicities as a function of de-projected galactocentric radius. The slope of this line is the metallicity gradient.

Table 1

Basic properties of the H IX and CONTROL (below the line) galaxies.

thumbnail Fig. 1

Example of  WIFES data for ESO111-G014. Left: red-scale image in the background shows the Bj -band optical image. The overlaid blue squares are of the size of the WIFES aperture and are located where the telescope was pointed. Right: colour map in the background shows the near-UV image, the overlaid blue squares indicate where “single” star forming regions were defined (similar figures for all galaxies can be found in Appendix A).

Table 2

WIFES observations of the H IX sample and CONTROL galaxy IC 4857.

thumbnail Fig. 2

Colour-scale images in the four left panels show the plane with Hα emission in the four single data cubes of the southern pointing of ESO111-G014. The black circles indicate the feature that was used to align the four data cubes. Right panel: plane with Hα emission in the final, stacked data cube.

Table 3

Summary of the samples and data used in this paper.

2.4 Comparison samples

As only one of the original CONTROL galaxies was observed with WIFES, data from the Sloan Digital Sky Survey (SDSS, York et al. 2000) data release 7 (Abazajian et al. 2009), and the (extended) GASS (Catinella et al. 2018, 2013, 2010; Moran et al. 2012)are used for comparison. An overview of all comparison samples is given in Table 3.

For the SDSS data, line emission and stellar mass measurements (Kauffmann et al. 2003; Salim et al. 2007) from the MPA-JHU SDSS DR72 catalogue are used. For galaxies within a redshift range of 0.002 < z < 0.06 central metallicities were calculated from 5 σ detected emission lines using again the O3N2 method. These data are used as a benchmark mass–metallicity relation in Sect. 3.2.

We also include data from the optical follow-up of GASS galaxies. These data are based on long-slit spectra rather than IFU spectra, andthey cover radii to 1 Petrosian 90% radius. The GASS survey measured the integrated atomic gas content of ~ 1200, stellar-mass selected galaxies. In addition, for about 200 galaxies optical long-slit spectra were obtained. From these spectra metallicity gradients were measured (Moran et al. 2012). Below, metallicity gradients, stellar masses, and H I mass fractions of star forming galaxies (i.e. near-UV-r colour bluer than 4.3) are used for comparison.To mimic the previous control sample, we restricted the sample to galaxies which are within ± 0.7 σ of the HIPASS running average in the H I mass fraction versus stellar mass plane (see Fig. 3). These GASS galaxies are shown as yellow dots in Fig. 3.

thumbnail Fig. 3

H I gas mass fraction as a function of stellar mass for H IX galaxies(blue stars),  CONTROL galaxy IC 4857 (turquois circle), and those galaxies from GASS with an optical long-slit spectrum (yellow dots). The grey line is the running average of the  HIPASS parent sample and the grey shaded area the corresponding 1 σ scatter.

3 Results

3.1 Comparison of stellar, Hα and H I velocities

To compare the kinematics of stars, and ionised and atomic gas, we consider their recession velocities minus the systemic velocity measured from the titled ring models in Lutz et al. (2018). A pairwise comparison of the three components is shown in Fig. 4. This figure shows for every spaxel how the local velocity of H I compares to that of the stars (left panel), to that of the Hα line (middle panel), and the velocities of stars to the Hα line (right panel). The horizontal and vertical lines indicate ± 100km s−1 of the systemicvelocity, which is the velocity resolution of the used setup of the WIFES spectrograph. Spaxels that are located in the upper left or lower right corner of this figure (red-shaded area) would indicate approaching movement in one component and receding movement in the other and thus counter-rotation. As can been seen, hardly any spaxels are located in the red regions. If the sign of the velocities of two components is different, this is usually within 100 km s−1, which is likely due to the low velocity resolution of the optical observations. Detailed kinematic maps for every galaxy can be found in Appendix A.

3.2 H IX galaxies on the mass–metallicity relation

Figure 5 shows the mass–metallicity relation for the H IX and CONTROL galaxies that were observed with WIFES. As mentioned above, the scatter of the mass–metallicity relation is driven by the H I content of the galaxies, where more H I-rich galaxies are more metal-poor at a given stellar mass (Bothwell et al. 2013; Brown et al. 2018). In Fig. 5, the black contours enclose 64 and 90% of the full SDSS DR7 sample (as described in Sect. 2.4). The data points showing the H IX galaxies (stars) and the CONTROL galaxy (large circle) are colour-coded according to H I mass fraction. For these galaxies, we only consider the most central metallicity measurement, which does not show line-ratios that are consistent with active galactic nuclei (AGN) emission.

There are no central WIFES pointing (closer than ≈ 11 kpc to the centre) available for ESO417-G018. Due to a metallicity gradient, ESO417-G018 is therefore located below the SDSS 90% contour. The H IX galaxies with the lowest central metallicity is, however, ESO290-G035. Apart from these two galaxies andESO208-G026 (H IX galaxy with lowest stellar mass), H IX galaxies are generally located in the lower half of the SDSS scatter (black dotted lines). To sum it up, H IX galaxies are generally expected to have lower metallicities in the centres, and indeed we find them to be predominantly located in the low metallicity half of the SDSS scatter.

3.3 Gas-phase metallicity distribution

Exploiting the spatial information from the IFU spectral data, we extracted radial metallicity profiles and measured metallicity gradients. For a consistent comparison between our samples and the literature, the metallicity gradients are measured in units of dex kpc−1. Metallicity gradients as a function of stellar mass are shown in Fig. 6 and data points are colour-coded by H I mass fraction. The grey line is the average Sánchez et al. (2014) gradient based on CALIFA data and the grey shaded area their 1 σ scatter, within which all H IX galaxies are located. Hence, H IX galaxies do not show steeper metallicity profiles than average local spiral galaxies. In fact, they tend to be located at the flat end of the Sánchez et al. (2014) scatter and their average metallicity gradient is (− 0.016 ± 0.009) dex kpc−1. This means that the average Sánchez et al. (2014) gradient does not agree with the average H IX gradient and its 1 σ scatter, but they do agree within both their errors.

More recent work on metallicity gradients with SAMI and MaNGA galaxies generally agrees with Sánchez et al. (2014), despite detecting anti-correlations between metallicity gradients and stellar masses at low stellar masses (log M [M ] < 9.6) (Poetrodjojo et al. 2018; Belfiore et al. 2017). We note, however, that diffuse emission significantly affects metallicity measurements in IFS data with poor spectral resolution, such as SAMI and MaNGA (Poetrodjojo et al. 2019). Furthermore, these surveys preferentially (only) report quantitative gradients in units of dex per effective radius or per isophotal radius. We thus refrain from more detailed comparison.

The CONTROL galaxy IC 4857 and the galaxies from the GASS (Moran et al. 2012) sample are also mostly located within the 1 σ scatter of CALIFA (circles in Fig. 6). We do not see the variation with stellar mass that was reported by Moran et al. (2012), perhaps because we only look at a small subsample. For the GASS galaxies we find an average metallicity gradient of (− 0.01 ± 0.01) dex kpc−1, which is consistent with the average H IX gradient, but again smaller than the average CALIFA gradient.

Moran et al. (2012) furthermore suggested that large, abrupt drops (of about 0.25 dex) in metallicity at the outskirts of the optical disc in massive galaxies can be caused by the inflow and distribution of pristine gas from the outskirts throughout the entire galaxy disc. Similarly, Sánchez Almeida et al. (2014a) found in a sample of dwarf galaxies, star forming regions with very low metallicities compared to other star forming regions within the same dwarf galaxy. They and Ceverino et al. (2016) argue that these metal-poor star forming regions are induced by pristine gas accretion. There are neither metallicity drops at the outskirts nor particularly metal-poor star forming regions observed in the H IX and galaxies or in IC 4857 (for the detailed data see Appendix A).

thumbnail Fig. 4

Comparison of the recession velocities of different galaxy components: stars, H I, and ionised gas. Every panel compares two components. To account for the different redshifts of the different galaxies, the systemic velocity of the respective galaxy (see Table 1) was subtracted from the recession velocity measurement. No clear sign of counter rotation is visible.

thumbnail Fig. 5

Mass–metallicity relation for the H IX galaxies: we show the gas-phase metallicity (determined with the O 3N 2 parameter following Pettini & Pagel 2004) as a function of stellar mass. The contours enclose 64 and 90% of SDSS galaxies. The circle indicates the CONTROL galaxy IC 4857 and H IX galaxies are marked with stars. All symbols are also colour-coded by the H I mass fraction.

thumbnail Fig. 6

Metallicity gradient in dex kpc−1 as a function of stellar mass. The grey line indicates the average gradient by Sánchez et al. (2014) and the grey shaded area the 1 σ scatter. Stars indicate the metallicity gradient measured in H IX galaxies, and the circles the gradients of the control galaxies (IC 4857 and GASS). All symbols are colour-coded by the H I mass fraction of the galaxy. Overall, H IX and control galaxies behave similarly to the galaxies in the Sánchez et al. (2014) sample.

4 Discussion and conclusion: optical spectra of H IX galaxies

We have presented the analysis of optical spectra of ten H IX galaxies as well as a set of control galaxies. In a first step, measurements of the Hα, stellar (where available) and H I recession velocities in the H IX galaxies and IC 4857 have been compared. Generally, the kinematics of all three components are similar. While the detection of counter-rotation, in particular between the stellar and gaseous components, would be strong evidence of recent accretion of large amounts of gas (Corsini 2014), the opposite (no counter-rotation means no accretion) is not necessarily true. Hence, from the kinematics point of view, the evidence for recent accretion of large amounts of gas in H IX galaxies is still inconclusive.

Using strong emission lines, gas-phase metallicities have been measured throughout the detected stellar discs of H IX, CONTROL, and GASS galaxies. These measurements have been used to place these galaxies on the mass–metallicity relation and to compare their metallicity gradients, amongst each other and to the average of Sánchez et al. (2014).

On the mass–metallicity relation, most H IX galaxies are located below the ridge of the relation. This is expected from and consistent with previous studies that found that H I-rich galaxies are generally metal poor (e.g. Bothwell et al. 2013; Hughes et al. 2013; Lagos et al. 2016; Brown et al. 2018). These studies, based on observations and on simulations, suggest that galaxies accrete gas stochastically, and that galaxies below the ridge of the mass–metallicity relation have recently accreted gas, which increased their gas content and decreased their gas-phase metallicity.

The radial distribution of the gas-phase metallicity (within the stellar disc) of H IX galaxies appears similar to the one of galaxies with average H I contents. All measured gradients in H IX galaxies are negative, i.e. they are more metal poor in the outskirts than in the centres, which is in agreement with the scenario where no major merger occurred recently. We have furthermore compared the H IX metallicity gradients to the metallicity gradients of galaxies with less massive H I discs (drawn from the CONTROL sample and GASS). Regardless of their overall H I content, there is no significant difference in the metallicity gradients of these galaxies. We note again that all H IX galaxies except ESO055-G013 have flatter gradients than the average Sánchez et al. (2014) gradient, yet all H IX galaxies are within their 1σ scatter. However, in the future we expect larger surveys, such as MaNGA or SAMI paired with dedicated H I observations to find an answer to this question.

The chemical evolution model of Ho et al. (2015) and Kudritzki et al. (2015) is able to reproduce measured metallicity gradients. This model suggests that the local metallicity is mostly dependent on the local stellar to total (H I + H2) gas mass ratio. This has been observationally confirmed by Barrera-Ballesteros et al. (2018). At the galactocentric radii where gas-phase metallicities can be measured, both the H I and the stellar mass column densities in H IX galaxies are similar to those of control galaxies. It is only at larger radii, that the H I column densities in control galaxies (or generally less H I-rich galaxies) are smaller than in H IX galaxies, simply because the H IX H I discs extend much further out. The stellar surface densities on the other hand are similar at all radii in H IX and the comparison galaxies (similar stellar effective radii; Lutz et al. 2018). This means that the pace of radial metallicity variation in H IX and comparison galaxies should be the same, and this is indeed what we observe. What is not observed in H IX galaxies are large drops in metallicity or star forming regions with a very low metallicity compared to their surroundings. Thus, major accretion events are not likely (Moran et al. 2012; Sánchez Almeida et al. 2014b).

The analysis presented in previous papers (Lutz et al. 2018, 2017) suggests that H IX galaxies are H I-rich because they are able to stabilise their H I disc against star formation due to a higher baryonic specific angular momentum. Examining the dark matter halo spin and H I content of galaxies simulated with a semi-analytic model, further showed that galaxies with a large H I disc tend to reside in dark matter halos with a large spin (i.e. angular momentum).

Based on the general interpretation of the scatter of the mass–metallicity relation, the new results presented here would suggest that despite the high baryonic specific angular momentum, metal-poor (accreted) gas would make its way from the outskirts to the centres of H IX galaxies, where it dilutes the local gas. This would imply that H IX galaxies are indeed those galaxies that recently accreted gas in a stochastic accretion process. However, this would also imply the following two points. First, their extreme H I mass would require them to accrete more gas than they contained before because they were selected to host at least 2.5 times more H I than expected from their optical properties. Second, the accreted gas must have a very high angular momentum and must increase the spin of the dark matter halo as well to comply with the findings from the semi-analytic models.

While Stewart et al. (2011, 2017) find in simulations that gas accreted in the cold mode would indeed come in with a high angular momentum and form a large co-rotating cold disc (at least in galaxies at redshift >1), the two points madeabove appear very extreme. Furthermore, Lara-López et al. (2013) find based on observations of H I content,specific star formation rate and metallicity that there are hardly any metal-poor, H I-rich galaxies at intermediate stellar mass and specific star formation rate. These are, however, exactly the characteristics of the H IX galaxies.

One way to reconcile these findings is provided by the semi-analytic models of Boissier & Prantzos (2000). These models investigate the chemical evolution of galaxies with respect to the circular rotation velocity and the halo spin parameter. They find that galaxies residing in higher spin halos evolve more slowly than galaxies in lower spin halos. With a slower evolution comes less enrichment of the ISM and thus a lower gas-phase metallicity, much like the effect that makes metallicity gradients universal (Sánchez et al. 2014). These findings have been verified by SDSS observations (Cervantes-Sodi & Hernández 2009). Catinella & Cortese (2015) propose a similar scenario for the most H I-rich galaxies detected at redshift 0.2.

Thus, the overall picture for the H IX galaxies is the following: they reside in halos that have higher spins than the halos of average galaxies in terms of HI content, but not as high as necessary to host low surface brightness galaxies (Boissier et al. 2016, 2003; Kim & Lee 2013). The high spin allows them to host larger H I discs than average. At the same time the high spin slows their evolution such that their ISM is (not yet) as enriched as it is in galaxies, which reside in average spin halos.

Acknowledgements

We would like to thank the anonymous referee for helpful comments that improved the paper. K.L. would like to thank Nicola Pastorello for helpful advice on the WIFES data reduction, Natasha Maddox for fruitful discussions, the technical staff at Siding Springs Observatory for their support during the observations, and the ANU service desk team for help with the WIFES data archive. Parts of this research were supported by the Australian Research Council Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), through project number CE170100013. Besides services and tools already mentioned, data used in this paper were provided by and analysed with TOPCAT (Taylor 2005), VizieR3 (Ochsenbein et al. 2000), AAO Data Central4, astropy (Astropy Collaboration 2013), astroquery (Ginsburg et al. 2019), and matplotlib (Hunter 2007).

Appendix A Results for individual H IX galaxies

Figures A.1A.10 show detailed results for each of the H IX galaxies individually. The maps of stellar and ionised gas kinematics, as well as the maps of metallicity (panels d, e, and g), were obtained with the following procedure: for each star forming region (their borders are shown in panel b) we derived one measurement of recession velocity or metallicity. We thus assign this measurement to the entire star forming region.

thumbnail Fig. A.1

ESO111-G014: (a): pointings are overlaid on SuperCOSMOS Bj-band images. (b): star forming regions marked on a GALEX FUV images (again Bj-band image for ESO208-G026). (c): H I moment 0 (intensity) map. Contours are at 0.5, 1.0, 2.5, 5.0, 7.5 M pc−2. (d): stellar kinematics of single star forming regions overlaid on Bj-band image. (e): Hα recession velocities of single star forming regions overlaid on Bj-band image. (f): H I moment 1 (velocity) map. (g): metallicity measurements in single star forming regions overlaid on Bj-band image. (h): radial profileof metallicities with linear fit to the profile, whose slope is the metallicity gradient. The grey dashed line marks the 25 mag arcsec−2 isophotal radius. (i): radial profile of H I column density. Dark blue circles and dashed line show the measured values corrected for inclination. Light blue points and dotted line show the profile rescaled to a 1 M pc−2 isophotal radius if the galaxy had average H I content (for more details see text). The grey dashed line again indicates the 25 mag arcsec−2 isophotal radius. (a)–(g): all seven panels have the same spatial size; north is up and east is left. (d)–(e): colour scale of the velocities is the same in all three panels.

thumbnail Fig. A.2

NGC 289. Panels as in Fig. A.1.

thumbnail Fig. A.3

ESO245-G010. Panels as in Fig. A.1.

thumbnail Fig. A.4

ESO417-G018. Panels as in Fig. A.1.

thumbnail Fig. A.5

ESO055-G013. Panels as in Fig. A.1.

thumbnail Fig. A.6

ESO208-G026. Panels as in Fig. A.1.

thumbnail Fig. A.7

ESO378-G003. Panels as in Fig. A.1.

thumbnail Fig. A.8

ESO381-G005. Panels as in Fig. A.1.

thumbnail Fig. A.9

ESO075-G006. Panels as in Fig. A.1.

thumbnail Fig. A.10

ESO290-G035. Panels as in Fig. A.1.

The rescaled profile in panel (i) was obtained in the following way. First the H I mass expected from the Dénes et al. (2014) relation was calculated. Then the expected H I radius (RHI, exp, a 1 M pc−2 isophotal radius) was calculated from the relation by Wang et al. (2016). Then the radii of the profile were rescaled with a factor RHI, exp/RHI, mea, where RHI, mea is the measured H I radius from Lutz et al. (2018).

Appendix B Results for IC 4857

Figure B.1 shows the detailed results of IC 4857, which is the only H IX-control galaxy that was observed with WIFES. Figure B.1 shows the same data as Figs. A.1A.10.

thumbnail Fig. B.1

IC 4857. Panels as in Fig. A.1.

References

  1. Abazajian, K. N., Adelman-McCarthy, J. K., Agüeros, M. A., et al. 2009, ApJS, 182, 543 [NASA ADS] [CrossRef] [Google Scholar]
  2. Astropy Collaboration (Robitaille, T., et al.) 2013, A&A, 558, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Barnes, D. G., Staveley-Smith, L., de Blok, W. J. G., et al. 2001, MNRAS, 322, 486 [NASA ADS] [CrossRef] [Google Scholar]
  4. Barrera-Ballesteros, J. K., Heckman, T., Sanchez, S. F., et al. 2018, ApJ, 852, 74 [NASA ADS] [CrossRef] [Google Scholar]
  5. Belfiore, F., Maiolino, R., Tremonti, C., et al. 2017, MNRAS, 469, 151 [NASA ADS] [CrossRef] [Google Scholar]
  6. Birnboim, Y., & Dekel, A. 2003, MNRAS, 345, 349 [NASA ADS] [CrossRef] [Google Scholar]
  7. Boissier, S., & Prantzos, N. 2000, MNRAS, 312, 398 [NASA ADS] [CrossRef] [Google Scholar]
  8. Boissier, S., Ragaigne, D. M., Prantzos, N., et al. 2003, MNRAS, 343, 653 [NASA ADS] [CrossRef] [Google Scholar]
  9. Boissier, S., Boselli, A., Ferrarese, L., et al. 2016, A&A, 593, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Bonnarel, F., Fernique, P., Bienaymé, O., et al. 2000, A&AS, 143, 33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Bothwell, M. S., Maiolino, R., Kennicutt, R., et al. 2013, MNRAS, 433, 1425 [NASA ADS] [CrossRef] [Google Scholar]
  12. Brown, T., Cortese, L., Catinella, B., & Kilborn, V. 2018, MNRAS, 473, 1868 [Google Scholar]
  13. Bundy, K., Bershady, M. A., Law, D. R., et al. 2015, ApJ, 798, 7 [NASA ADS] [CrossRef] [Google Scholar]
  14. Calzetti, D., Armus, L., Bohlin, R. C., et al. 2000, ApJ, 533, 682 [NASA ADS] [CrossRef] [Google Scholar]
  15. Cappellari, M. 2017, MNRAS, 466, 798 [NASA ADS] [CrossRef] [Google Scholar]
  16. Cappellari, M., & Emsellem, E. 2004, PASP, 116, 138 [NASA ADS] [CrossRef] [Google Scholar]
  17. Catinella, B., & Cortese, L. 2015, MNRAS, 446, 3526 [NASA ADS] [CrossRef] [Google Scholar]
  18. Catinella, B., Schiminovich, D., Kauffmann, G., et al. 2010, MNRAS, 403, 683 [NASA ADS] [CrossRef] [Google Scholar]
  19. Catinella, B., Schiminovich, D., Cortese, L., et al. 2013, MNRAS, 436, 34 [NASA ADS] [CrossRef] [Google Scholar]
  20. Catinella, B., Saintonge, A., Janowiecki, S., et al. 2018, MNRAS, 476, 875 [NASA ADS] [CrossRef] [Google Scholar]
  21. Cervantes-Sodi, B., & Hernández, X. 2009, Rev. Mex. Astron. Astrofis., 45, 75 [NASA ADS] [Google Scholar]
  22. Ceverino, D., Sánchez Almeida, J., Muñoz Tuñón, C., et al. 2016, MNRAS, 457, 2605 [NASA ADS] [CrossRef] [Google Scholar]
  23. Childress, M., Vogt, F., Nielsen, J., & Sharp, R. 2014, Ap&SS, 349, 617 [NASA ADS] [CrossRef] [Google Scholar]
  24. Corsini, E. 2014, ASP Conf. Ser., 486, 51 [NASA ADS] [Google Scholar]
  25. Croom, S. M., Lawrence, J. S., Bland-Hawthorn, J., et al. 2012, MNRAS, 421, 872 [NASA ADS] [Google Scholar]
  26. Davé, R., Katz, N., Oppenheimer, B., Kollmeier, J., & Weinberg, D. 2013, MNRAS, 434, 2645 [NASA ADS] [CrossRef] [Google Scholar]
  27. Dekel, A., & Birnboim, Y. 2006, MNRAS, 368, 2 [NASA ADS] [CrossRef] [Google Scholar]
  28. Dénes, H., Kilborn, V., & Koribalski, B. 2014, MNRAS, 444, 667 [NASA ADS] [CrossRef] [Google Scholar]
  29. Di Teodoro, E. M., & Fraternali, F. 2014, A&A, 567, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Dopita, M., Hart, J., McGregor, P., et al. 2007, Ap&SS, 310, 255 [NASA ADS] [CrossRef] [Google Scholar]
  31. Doyle, M., Drinkwater, M., Rohde, D., et al. 2005, MNRAS, 361, 34 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  32. Fogarty, L. M. R., Bland-Hawthorn, J., Croom, S. M., et al. 2012, ApJ, 761, 169 [NASA ADS] [CrossRef] [Google Scholar]
  33. Forbes, J., Krumholz, M., Burkert, A., & Dekel, A. 2014, MNRAS, 443, 168 [NASA ADS] [CrossRef] [Google Scholar]
  34. Fraternali, F., Binney, J., Oosterloo, T., & Sancisi, R. 2007, New Astron. Rev., 51, 95 [NASA ADS] [CrossRef] [Google Scholar]
  35. Geréb, K., Catinella, B., Cortese, L., et al. 2016, MNRAS, 462, 382 [NASA ADS] [CrossRef] [Google Scholar]
  36. Geréb, K., Janowiecki, S., Catinella, B., Cortese, L., & Kilborn, V. 2018, MNRAS, 476, 896 [NASA ADS] [Google Scholar]
  37. Ginsburg, A., Sipőcz, B. M., Brasseur, C. E., et al. 2019, AJ, 157, 98 [NASA ADS] [CrossRef] [Google Scholar]
  38. Hallenbeck, G., Huang, S., Spekkens, K., et al. 2014, AJ, 148, 69 [NASA ADS] [CrossRef] [Google Scholar]
  39. Hallenbeck, G., Huang, S., Spekkens, K., et al. 2016, AJ, 152, 225 [NASA ADS] [CrossRef] [Google Scholar]
  40. Ho, I.-T., Kudritzki, R.-P., Kewley, L. J., et al. 2015, MNRAS, 448, 2030 [NASA ADS] [CrossRef] [Google Scholar]
  41. Huang, S., Haynes, M., Giovanelli, R., et al. 2014, ApJ, 793, 40 [NASA ADS] [CrossRef] [Google Scholar]
  42. Hughes, T. M., Cortese, L., Boselli, A., Gavazzi, G., & Davies, J. I. 2013, A&A, 550, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Hunter, J. 2007, Comput. Sci. Eng., 9, 90 [Google Scholar]
  44. Jones, E., Oliphant, T., Peterson, T., et al. 2001, SciPy: Open source scientific tools for Python [Google Scholar]
  45. Kauffmann, G., Heckman, T. M., White, S. D. M., et al. 2003, MNRAS, 341, 33 [NASA ADS] [CrossRef] [Google Scholar]
  46. Kereš, D., Katz, N., Weinberg, D., & Davé, R. 2005, MNRAS, 363, 2 [NASA ADS] [CrossRef] [Google Scholar]
  47. Kim, J.-h., & Lee, J. 2013, MNRAS, 432, 1701 [NASA ADS] [CrossRef] [Google Scholar]
  48. Koribalski, B. S., Staveley-Smith, L., Kilborn, V. A., et al. 2004, AJ, 128, 16 [NASA ADS] [CrossRef] [Google Scholar]
  49. Kudritzki, R.-P., Ho, I.-T., Schruba, A., et al. 2015, MNRAS, 450, 342 [Google Scholar]
  50. Lagos, C., Theuns, T., Schaye, J., et al. 2016, MNRAS, 459, 2632 [NASA ADS] [CrossRef] [Google Scholar]
  51. Lara-López, M. A., Hopkins, A. M., López-Sánchez, A. R., et al. 2013, MNRAS, 433, L35 [NASA ADS] [CrossRef] [Google Scholar]
  52. Lauberts, A.,& Valentijn, E. 1989, The Surface Photometry Catalogue of the ESO-Uppsala Galaxies (Garching: European Southern Observatory) [Google Scholar]
  53. Lee, C., Chung, A., Yun, M., et al. 2014, MNRAS, 441, 1363 [CrossRef] [Google Scholar]
  54. Lemonias, J., Schiminovich, D., Catinella, B., Heckman, T., & Moran, S. 2014, ApJ, 790, 27 [NASA ADS] [CrossRef] [Google Scholar]
  55. Lutz, K., Kilborn, V., Catinella, B., et al. 2017, MNRAS, 467, 1083 [NASA ADS] [Google Scholar]
  56. Lutz, K. A., Kilborn, V. A., Koribalski, B. S., et al. 2018, MNRAS, 476, 3744 [NASA ADS] [CrossRef] [Google Scholar]
  57. Mannucci, F., Cresci, G., Maiolino, R., Marconi, A., & Gnerucci, A. 2010, MNRAS, 408, 2115 [NASA ADS] [CrossRef] [Google Scholar]
  58. Meyer, M., Zwaan, M., Webster, R., et al. 2004, MNRAS, 350, 1195 [NASA ADS] [CrossRef] [Google Scholar]
  59. Moran, S., Heckman, T., Kauffmann, G., et al. 2012, ApJ, 745, 66 [NASA ADS] [CrossRef] [Google Scholar]
  60. Ochsenbein, F., Bauer, P., & Marcout, J. 2000, A&AS, 143, 23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Oosterloo, T., Fraternali, F., & Sancisi, R. 2007, AJ, 134, 1019 [NASA ADS] [CrossRef] [Google Scholar]
  62. Pettini, M., & Pagel, B. 2004, MNRAS, 348, L59 [NASA ADS] [CrossRef] [Google Scholar]
  63. Poetrodjojo, H., Groves, B., Kewley, L. J., et al. 2018, MNRAS, 479, 5235 [NASA ADS] [CrossRef] [Google Scholar]
  64. Poetrodjojo, H., D’Agostino, J. J., Groves, B., et al. 2019, MNRAS, 487, 79 [NASA ADS] [CrossRef] [Google Scholar]
  65. Saintonge, A., Catinella, B., Tacconi, L. J., et al. 2017, ApJS, 233, 22 [Google Scholar]
  66. Salim, S., Rich, R. M., Charlot, S., et al. 2007, ApJS, 173, 267 [NASA ADS] [CrossRef] [Google Scholar]
  67. Sánchez, S. F., Kennicutt, R. C., Gil de Paz, A., et al. 2012, A&A, 538, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Sánchez, S. F., Rosales-Ortega, F. F., Jungwiert, B., et al. 2013, A&A, 554, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Sánchez, S. F., Rosales-Ortega, F. F., Iglesias-Páramo, J., et al. 2014, A&A, 563, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Sánchez Almeida, J., Morales-Luis, A., Muñoz-Tuñón, C., et al. 2014a, ApJ, 783, 45 [NASA ADS] [CrossRef] [Google Scholar]
  71. Sánchez Almeida, J., Elmegreen, B. G., Muñoz-Tuñón, C., & Elmegreen, D. M. 2014b, A&ARv, 22, 71 [NASA ADS] [CrossRef] [Google Scholar]
  72. Schiminovich, D., Catinella, B., Kauffmann, G., et al. 2010, MNRAS, 408, 919 [NASA ADS] [CrossRef] [Google Scholar]
  73. Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [NASA ADS] [CrossRef] [Google Scholar]
  74. Stewart, K., Kaufmann, T., Bullock, J., et al. 2011, ApJ, 738, 39 [NASA ADS] [CrossRef] [Google Scholar]
  75. Stewart, K., Maller, A., Oñorbe, J., et al. 2017, ApJ, 843, 47 [NASA ADS] [CrossRef] [Google Scholar]
  76. Taylor, M. B. 2005, in Astronomical Data Analysis Software and Systems XIV, eds. P. Shopbell, M. Britton, & R. Ebert (San Francisco: ASP), 347 29 [Google Scholar]
  77. Torrey, P., Vogelsberger, M., Marinacci, F., et al. 2019, MNRAS, 484, 5587 [NASA ADS] [Google Scholar]
  78. Tremonti, C. A., Heckman, T. M., Kauffmann, G., et al. 2004, ApJ, 613, 898 [NASA ADS] [CrossRef] [Google Scholar]
  79. van de Voort, F., Schaye, J., Booth, C., Haas, M., & Dalla Vecchia C. 2011, MNRAS, 414, 2458 [NASA ADS] [CrossRef] [Google Scholar]
  80. van den Bergh, S. 1962, AJ, 67, 486 [NASA ADS] [CrossRef] [Google Scholar]
  81. Vazdekis, A., Sánchez-Blázquez, P., Falcón-Barroso, J., et al. 2010, MNRAS, 404, 1639 [NASA ADS] [Google Scholar]
  82. Wang, J., Kauffmann, G., Józsa, G., et al. 2013, MNRAS, 433, 270 [NASA ADS] [CrossRef] [Google Scholar]
  83. Wang, J., Koribalski, B. S., Serra, P., et al. 2016, MNRAS, 460, 2143 [Google Scholar]
  84. Wen, X.-Q., Wu, H., Zhu, Y.-N., et al. 2013, MNRAS, 433, 2946 [NASA ADS] [CrossRef] [Google Scholar]
  85. Wisnioski, E., Schreiber, N. M. F., Wuyts, S., et al. 2015, ApJ, 799, 209 [NASA ADS] [CrossRef] [Google Scholar]
  86. York, D. G., Adelman, J., Anderson, John E., J., et al. 2000, AJ, 120, 1579 [CrossRef] [Google Scholar]
  87. Zwaan, M., Meyer, M., Webster, R., et al. 2004, MNRAS, 350, 1210 [NASA ADS] [CrossRef] [Google Scholar]

1

The galaxies in this sample cover a stellar mass range of 9.0 ≤ log M[M] ≤ 11.2, but the majority are more massive than 1010 M (Sánchez et al. 2013).

All Tables

Table 1

Basic properties of the H IX and CONTROL (below the line) galaxies.

Table 2

WIFES observations of the H IX sample and CONTROL galaxy IC 4857.

Table 3

Summary of the samples and data used in this paper.

All Figures

thumbnail Fig. 1

Example of  WIFES data for ESO111-G014. Left: red-scale image in the background shows the Bj -band optical image. The overlaid blue squares are of the size of the WIFES aperture and are located where the telescope was pointed. Right: colour map in the background shows the near-UV image, the overlaid blue squares indicate where “single” star forming regions were defined (similar figures for all galaxies can be found in Appendix A).

In the text
thumbnail Fig. 2

Colour-scale images in the four left panels show the plane with Hα emission in the four single data cubes of the southern pointing of ESO111-G014. The black circles indicate the feature that was used to align the four data cubes. Right panel: plane with Hα emission in the final, stacked data cube.

In the text
thumbnail Fig. 3

H I gas mass fraction as a function of stellar mass for H IX galaxies(blue stars),  CONTROL galaxy IC 4857 (turquois circle), and those galaxies from GASS with an optical long-slit spectrum (yellow dots). The grey line is the running average of the  HIPASS parent sample and the grey shaded area the corresponding 1 σ scatter.

In the text
thumbnail Fig. 4

Comparison of the recession velocities of different galaxy components: stars, H I, and ionised gas. Every panel compares two components. To account for the different redshifts of the different galaxies, the systemic velocity of the respective galaxy (see Table 1) was subtracted from the recession velocity measurement. No clear sign of counter rotation is visible.

In the text
thumbnail Fig. 5

Mass–metallicity relation for the H IX galaxies: we show the gas-phase metallicity (determined with the O 3N 2 parameter following Pettini & Pagel 2004) as a function of stellar mass. The contours enclose 64 and 90% of SDSS galaxies. The circle indicates the CONTROL galaxy IC 4857 and H IX galaxies are marked with stars. All symbols are also colour-coded by the H I mass fraction.

In the text
thumbnail Fig. 6

Metallicity gradient in dex kpc−1 as a function of stellar mass. The grey line indicates the average gradient by Sánchez et al. (2014) and the grey shaded area the 1 σ scatter. Stars indicate the metallicity gradient measured in H IX galaxies, and the circles the gradients of the control galaxies (IC 4857 and GASS). All symbols are colour-coded by the H I mass fraction of the galaxy. Overall, H IX and control galaxies behave similarly to the galaxies in the Sánchez et al. (2014) sample.

In the text
thumbnail Fig. A.1

ESO111-G014: (a): pointings are overlaid on SuperCOSMOS Bj-band images. (b): star forming regions marked on a GALEX FUV images (again Bj-band image for ESO208-G026). (c): H I moment 0 (intensity) map. Contours are at 0.5, 1.0, 2.5, 5.0, 7.5 M pc−2. (d): stellar kinematics of single star forming regions overlaid on Bj-band image. (e): Hα recession velocities of single star forming regions overlaid on Bj-band image. (f): H I moment 1 (velocity) map. (g): metallicity measurements in single star forming regions overlaid on Bj-band image. (h): radial profileof metallicities with linear fit to the profile, whose slope is the metallicity gradient. The grey dashed line marks the 25 mag arcsec−2 isophotal radius. (i): radial profile of H I column density. Dark blue circles and dashed line show the measured values corrected for inclination. Light blue points and dotted line show the profile rescaled to a 1 M pc−2 isophotal radius if the galaxy had average H I content (for more details see text). The grey dashed line again indicates the 25 mag arcsec−2 isophotal radius. (a)–(g): all seven panels have the same spatial size; north is up and east is left. (d)–(e): colour scale of the velocities is the same in all three panels.

In the text
thumbnail Fig. A.2

NGC 289. Panels as in Fig. A.1.

In the text
thumbnail Fig. A.3

ESO245-G010. Panels as in Fig. A.1.

In the text
thumbnail Fig. A.4

ESO417-G018. Panels as in Fig. A.1.

In the text
thumbnail Fig. A.5

ESO055-G013. Panels as in Fig. A.1.

In the text
thumbnail Fig. A.6

ESO208-G026. Panels as in Fig. A.1.

In the text
thumbnail Fig. A.7

ESO378-G003. Panels as in Fig. A.1.

In the text
thumbnail Fig. A.8

ESO381-G005. Panels as in Fig. A.1.

In the text
thumbnail Fig. A.9

ESO075-G006. Panels as in Fig. A.1.

In the text
thumbnail Fig. A.10

ESO290-G035. Panels as in Fig. A.1.

In the text
thumbnail Fig. B.1

IC 4857. Panels as in Fig. A.1.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.