The influence of ionized gas kinematics on Hii galaxies
The cases of Tol 1004-296 and Tol 0957-278

The dwarf galaxy Tol 1004-296, also known as NGC 3125, ESO 435-G 041, ESO 100419-2941.5, AM 1004-294, MCG 05-24-022, is classified as a BCG with a magnitude of $M_{V}$ = -17.7 (Bergvall & Olofsson, 1986) and a low metallicity of 12 + Log(O/H) = 8.34 (Vacca & Conti, 1992). It is located at a distance of 11.99 Mpc ¹¹1Local Group distance from Nasa Extragalactic Database. Smith et al. (1976) described it as an object with super condensations showing emission lines and a diffuse stellar continuum, from an objective prism survey (Smith, 1975). Spectroscopic survey by Penston et al. (1977) confirmed these objective prism data, exhibiting HeII$\lambda$4686 emission line in one of the condensation. Reif et al. (1982) has detected Hi content and Paturel et al. (2003) converted it to mass according to Pustilnik & Martin (2007), giving an $M_{HI}=3.19\times 10^{8}M_{\odot}$. Kunth et al. (1985) report the presence of Wolf Rayet (WR) stars in the object. Dust content and mass have been estimated to $4.4\times 10^{5}M_{\odot}$ by Engelbrachtet al. (2008), using Spitzer Observatory, also giving a dust temperature of $53.7\pm 0.3$ K.

Tol 1004-296 has been observed in the R, B and H$\alpha$ bands on several occasions. Kunth et al. (1988), Papaderos et al. (2002), Gil de Paz et al. (2003, 2005) present studies of this object based on images in these three bands. Kunth et al. (1988) B image mentioned four intense star formation regions. Gil de Paz et al. (2003) shows both the broad and narrow bands images of the object. It is notable that the extension of the H$\alpha$ emission goes beyond 40″. We can also see that this extension, at lower limits, appears as filaments. A 20″$\times$ 20″zoomed field of view is also presented in Gil de Paz et al. (2003) H$\alpha$ images (slightly different compared to the observed field with SIFS - see below).

The spectrum shown by Terlevich et al. (1991) is dominated by intense emission lines of hydrogen (H$\alpha$ and H$\beta$) and forbidden oxygen lines such as [OIII]$\lambda$5007 and [OIII]$\lambda$4959. Kehrig et al. (2004) also give the fluxes of the main emission lines from the high efficient spectrograph $FEROS$ using one fibre.

Lagos et al. (2007) published an H$\beta$ emission map and equivalent width (EW) image showing EW(H$\beta$) above 90Å. Three intense emission sites are also visible in these maps.

In Schwartz et al. (2006), low-amplitude outflows were shown through spectroscopy conducted with the HST in the two main star-forming regions of the object.

Observations in the near-infrared, using the ESO-VLT SINFONI (Eisenhauer et al., 2003), made it possible to obtain images of the continuum emission in the K-band, where a series of faint compact sources are detected near the bright central region (Vanzi et al., 2011). The SINFONI field was 8 ′$\times$ 8 ′with a resolution of 0.250 $\times$ 0.125 arcsec/pixel. Vanzi et al. (2011) focused specifically on node 1 in the nomenclature of Lagos et al. (2007) and A in the nomenclature of Vacca & Conti (1992) and Stevens et al. (2002).

The galaxy is also composed of starburst regions, as well as being populated by massive star clusters that span a range of ages. Additionally, Vanzi et al. (2011) found that in several areas there is thermal bremsstrahlung emission and molecular hydrogen clouds separate from the Hii regions. These regions of molecular hydrogen are mainly excited in a fluorescent manner, extending beyond the giant Hii regions and outlining the shape of the molecular cloud involved in star formation Vanzi et al. (2011).

Another component of this galaxy is the presence of peaks of [FeII] emission located outside the Hii regions, indicating the presence of supernovae that may have triggered starburst episodes. Vanzi et al. (2011) suggests that these episodes occurred in a chain, with the older one triggering the most recent one. This hypothesis was also raised by Cresci et al. (2010).

Both Cresci et al. (2010) and Vanzi et al. (2011) show regions where starburst episodes are occurring. However, Vanzi et al. (2011) found that the massive starburst episodes occurred in three separate episodes, resulting in the distinct intense emission sources.

The massive starburst episodes populate the galaxy with hot and young stars. However, this doesn’t imply the presence of older stars in the galaxy: one characteristic of Hii galaxies is the spatial dominance of young stars, making it challenging to detect older stellar populations (Westera et al., 2004).

Tol 1004-296 exhibits the strongest known emission in [HeII]$\lambda 1640$ of stellar origin, providing important information about the most massive stars in the galaxy. In Chanda et al. (2004), a study of the massive stellar content (Wolf-Rayet stars and O stars) in Tol 1004-296 was conducted, revealing a significant predominance of WN late-type Wolf-Rayet stars. Westera et al. (2004) and Cuisinier et al. (2006) have estimated characteristics of young, intermediate, and old stellar populations. The ratio between young + intermediate to old stellar population mass is 1:100 and the age is between 10 to 20 Myrs. And the young stellar mass to the total stellar mass ratio is 0.008. The young population age is between 1 to 2 Myrs and the old population to the total stellar mass ratio is 0.99.

The galaxy Tol 0957-278, also known as Tol 02, ESO 435-IG 020, AM 0957-275, PGC 028863, was initially observed by Smith et al. (1976) and classified as a medium-strong emission line emitter. Kunth et al. (1985); Schaerer et al. (1999); Méndez & Esteban et al. (2000) also report the presence of Wolf-Rayet (WR) stars, similar to Tol 1004-296. Broadband CCD imaging from Kunth et al. (1988) and Doublier et al. (1999) confirms the clumpy shape of the object. Campbell et al. (1988)’s modeling confirms a low metallicity of 12 + log(O/H) = 8.1, consistent with what has been reported in more recent studies (Vacca & Conti, 1992; Masegosa et al., 1994; Buckalew et al., 2005; Engelbrachtet al., 2008).

Terlevich et al. (1991) and Masegosa et al. (1994) published emission line fluxes and ion abundances for different samples of Hii galaxies, including both objects in our study. A more recent study by Kehrig et al. (2004) also provides emission line fluxes, which we used for the selection of our two targets.

Tol 0957-278 is neighbor to several other Hii galaxies as mentioned by Telles & Terlevich (1995) and seems to form a pair system with ESO 435-G 016. Sung et al. (2002) saw this pair as an interacting system. Kim et al. (2017) studied the Hi content of this system, estimating Tol 0957-278’s extended Hi gas mass to be $1.11\times 10^{8}M_{\odot}$. They also concluded that the system has been tidally disturbed, and the interaction is responsible for the recent star formation activities.

Tol 0957-278 has a distance of 9.93 Mpc²²2Local Group distance from NASA Extragalactic Database and an angular extension of $\approx~{}60\arcsec\times 40\arcsec$ (Torres-Campos et al., 2017). It is also classified as a BCD galaxy due to its low B-band magnitude $M_{B}$ = -16.26 mag (Kunth et al., 1988; Doublier et al., 1999; Gil de Paz et al., 2003). Gil de Paz et al. (2003) observed it in the H$\alpha$ band, and, similar to the previous object, Tol 0957-278 shows extended H$\alpha$ emission (25″) with a filamentary aspect on the outskirts.

Lagos et al. (2007) also published H$\beta$ emission and equivalent width (EW) images showing EW(H$\beta$) above 150Å. Four intense emission sites are also visible in these maps.

An underlying older stellar population mixed with a much younger stellar population has been known for quite some time, but it is still difficult to separate (Westera et al., 2004). In their study, Westera et al. (2004) and Cuisinier et al. (2006) found the following characteristics of this older population. The mass ratio between the young + intermediate population/old population is 1 to 100, the mass ratio between the young population and the total stellar mass is between 0.0081 and 0.0161. Finally, they estimated that the ratio between the old population and the total stellar mass is 0.96. The age of the young population is between 1 and 2 Myrs, and the age of the intermediate population is between 50 and 200 Myrs.

Torres-Campos et al. (2017) conducted a photometric study of the object in optical, NIR broadbands filters, and narrow H$\alpha$, [OIII]$\lambda$5007 filters. They identified eight star cluster complexes (SCCs). These SCCs are delineated by the H$\alpha$ emission, which is located in the central zone of the galaxy and extends towards the north and south regions. They have focused on decomposing the galaxy into a disk - or host galaxy formed by old stellar populations - and the Star Formation region (Low Surface Brightness Component - LSBC and Star Formation Component). They also estimated the ages and masses of these components. The age and masses of the LSBCs span between 3.7 Myrs to 2.5 Gyrs and $2.4\times 10^{6}M_{\odot}$ to $93.3\times 10^{6}M_{\odot}$. The age and masses of the SCCs span between 3.1 Myrs to 11 Myrs and $35.5\times 10^{3}M_{\odot}$ to $269.1\times 10^{3}M_{\odot}$.

These two objects were observed at the SOAR Observatory ³³3Southern Observatory for Astrophysical Research using the SIFS (Lepine et al., 2003; Fraga, 2018) instrument ⁴⁴4SOAR Integral Field Spectrometer. SIFS is an Integral Field Unit (IFU) spectrograph. The instrument has a spatial scale of 0.3″  per fiber, and the field of view is 15″$\times$ 7.8″. The fibers feed a pseudo slit to a bench spectrograph. The detector is a 4096 x 4112 pixel CCD231-84 from E2V with a pixel physical size of 15 $\mu$m². The Read-Out Noise (RON) of the detector is 2.0 e-/ADU.

Several observations were necessary to cover the targets. We used four different fields for each target resulting to a final field of view of 30″$\times$ 15 ″. Each target was observed with two different gratings: the 700R (with a resolution of 4200@5500) was used for obtaining monochromatic maps for different emission lines, and the 1500M3 (with a resolution of 9500@5500) for the velocity and velocity dispersion maps. Table 1 shows the observation log. Figure 1 illustrates the different fields using an H$\alpha$ image from the Irénée du Pont telescope at La Campanas Observatory ⁵⁵5These images were downloaded from the NED website for both galaxies. We also show the main emission knots.

We estimated a precision equivalent to a spaxel (0.3″) in the reconstruction of the mosaics. Observations couldn’t be completed in a single semester and had to be spread out over two semesters, a year apart.

During the night, two spectrophotometric standard stars were also observed: LTT 3218 and EG 274 (Hamuy et al., 1993), for flux calibration purposes.

We are also presenting data from the MUSE (Bacon et al., 2010) instrument attached to the ESO-VLT for Tol 1004-296. The observation was part of proposal ID 094.B-0745 (PI - Garcia-Benito, R.). The MUSE instrument has the great advantage of having a large Field of View (FoV) equal to $\approx$ 1.0 ′$\times$ 1.0 ′, allowing the observation of the entire object in a single 2880s observation. The spatial sampling is 0.2″  per pixel. On the other hand, the available grating only provides a spectral resolving power of 1770@4800Å.

The subsequent step involved extracting sub-cubes for each of the targeted emission lines. Depending on the lines, the data cube depth had between 30 or 40 spectral pixels, sufficient to determine the continuum level. To account for seeing effects, we applied a light spatial Gaussian smoothing of 3 spaxels (equivalent to 0.9″) Full Width at Half Maximum (FWHM).

Next, we generated various moment maps using two different techniques: an in-house macro and a Gaussian profile fitting technique.

• •

  Monochromatic and Continuum Maps: For each sub-cube, we identified the peak emission and isolated the profile. Then, we performed a linear fit of the remaining continuum and removed it from the sub-cube. Subsequently, we collapsed it to obtain monochromatic maps in different emission lines. We also generated continuum maps by removing the prominent emission lines and collapsing the cube, particularly for the MUSE data where the quality was better.

• •

  Radial Velocity Maps: To produce the velocity map, we isolated the profile as described before and estimated the barycenter ⁷⁷7It is estimated from histogram of line width, with the weight being the area of each rectangle forming the profile. More details in Amram (1991) of the profile to account for possible asymmetries.

• •

  Velocity Dispersion Maps: The velocity dispersion map was derived from the FWHM of the isolated profile. The FWHM was then transformed into dispersion using the formula $\sigma=\frac{FWHM}{2\sqrt{2ln2}}$.
  The velocity dispersion was corrected for thermal⁸⁸8See Rybicki & Lightman (2004) and instrumental broadening following the formula.

  $$\sigma=\sqrt{\sigma_{obs}^{2}-\sigma_{th}^{2}-\sigma_{inst}^{2}}$$

  The thermal dispersion, $\sigma_{th}$, was estimated to be 9.32 km s^-1 for hydrogen at a temperature of $10^{4}$ K. We estimated the instrumental broadening using the calibration lamp observation, giving $\sigma_{inst}$ = 13.0 $\pm$ 3 km s^-1.

Additionally, we used a Gaussian fit macro to cross-verify our results. However, we opted to retain the technique described above for estimating the velocity maps and velocity dispersion, as it allows for the consideration of possible asymmetries in the profiles.

We were able to produce monochromatic maps of the following emission lines: H$\beta\lambda$4861, [OIII]$\lambda\lambda$4959, 5007, HeI$\lambda$5876, [NII]$\lambda$6548, H$\alpha\lambda$6563, [NII]$\lambda$6584, [SII]$\lambda\lambda$6716, 6732, HeI$\lambda$6678, [ArIII]$\lambda$7136. Figure 2 to Figure 4 present these different monochromatic maps for both objects (except that for Tol 0957-278, where [HeI]$\lambda$5876, [NII]$\lambda$6548, [HeI]$\lambda$6678 and [ArIII]$\lambda$7136 lines were too weak to be extracted). Table 2 and Table 3 summarize the fluxes of the different emission lines in the different emission knots (they also include several others quantities such as diagnostic line ratios, oxygen abundances, interstellar extinction, radial velocity and velocity dispersion). These knots are symbolically shown in Figure 2 (H$\alpha$ map) and Figure 4 (H$\alpha$ map). We estimated these regions by limiting the extension where the SNR is above 5 in the $H\alpha$ line. In Figure 19, we have added several histograms representing the distribution of several parameters ([SII]$\lambda$6716/[SII]$\lambda$6731, EW(H$\beta$), E(B-V) and diagnostics diagrams line ratios).

In Tol 1004-296 case, we recognised the two large knots (Knot 1 and 2) separated by 10.3″ (661 pc), visible in Figure 1.

It is noticeable that only the two largest bright sites appear weaker in others emission wavelengths, such as [NII]$\lambda$6584, [SII]$\lambda$6716, [SII]$\lambda$6716, [HeI]$\lambda$5876, [NII]$\lambda$6548, [HeI]$\lambda$6678 and [ArIII]$\lambda$7136. We can also notice that the Knot 7 has a more elongated shape or formed by two sub nuclei in [OIII]$\lambda$5007 or H$\alpha$.

The monochromatic maps from the MUSE data are available in Figures 20. Extension of the ionized gas is, as expected, much larger compared to our SIFS data. For the brightest emission lines (H$\alpha$, H$\beta$, [OIII]$\lambda\lambda$ 4959, 5007), the extension is 50″$\times$ 50″  (we limit the emission to a Signal to Noise Ratio (SNR) superior to 8, in order to keep the filamentary character of the emission, in the outskirt, without the sky background noise). It is more or less the same extension showed by Gil de Paz et al. (2003) with the narrow band imaging of the H$\alpha$ + [NII] emission. The MUSE data show the following emission lines: H$\beta~{}\lambda$4861, [OIII]$\lambda\lambda$4959, 5007, HeI$\lambda$5876, [OI]$\lambda$6300, [NII]$\lambda$6548, H$\alpha\lambda$6563, [NII]$\lambda$6584, [SII]$\lambda\lambda$6716, 6732, HeI$\lambda$6678, [ArIII]$\lambda$7136. We easily recognise the different Knots we have presented previously.

Tol 0957-278 maps show a more complex morphology (Figure 4), nine emission sites can be seen on the $[OIII]\lambda 4959$ and $[OIII]\lambda 5007$ lines. H$\alpha$ map in Figure  4 shows the location of emission sites in the galaxy. Four knots (1, 2, 3 and 4) are brighter than the other five knots (5, 6, 7, 8 and 9).

Dust extinction has to be estimated and then applied to the observed fluxes in order to have correct fluxes. Interstellar extinction in nebula is computed using the observed Balmer decrement. The intrinsic luminosity, $L_{int}$ in a given wavelength can be written as:

with $L_{obs}(\lambda)$ the observed luminosity, $k(\lambda)$ the extinction law and $E(B-V)$ the color excess (Calzetti et al., 1994; Domínguez et al., 2013). We used the Calzetti et al. (2000) extinction law for $k(\lambda)$ and the color excess $E(B-V)$ has been derived from the Balmer decrement using $H\alpha$ and $H\beta$ lines by:

with $k(\lambda_{H\beta})$ and $k(\lambda_{H\alpha})$ the extinction coefficients for $H\beta$ and $H\alpha$.

Adopting the theoretical ratio of 2.86 (case B recombination at a temperature of $10^{4}$ K; see Osterbrock & Ferland (2006)), the color excess is:

Figure 6 shows $H\alpha$/$H\beta$ ratio maps for both objects. For Tol 1004-296, the ratio is more or less uniform around 3.8 $\pm$ 0.4 with some fluctuation. Some hot spots are visible where the emission is higher (where the strong emission sites are), as for the Knot 1 (where the ratio is 4.3 $\pm$ 0.2), the Knot 2 (where the ratio is 4.2 $\pm$ 0.1) and the Knot 3 (where the ratio is 4.0 $\pm$ 0.3).

For Tol 0957-278, the average ratio is around 4.3 $\pm$ 0.9. It is notable that the Knot 1 shows a value of 5.7 $\pm$ 0.5 and 3.4″East, the ratio goes down to 2 $\pm$ 0.4. It is possible that this difference is coming from a dust lane, present in this area.

Figure 21 shows the $H\alpha/H\beta$ line ratio for Tol 1004-296 observed with MUSE. The line ratio inside the SIFS FoV is consistent with the map from the SIFS data.

It wasn’t possible to use the direct T_e method (e.g. Dinerstein et al., 1990; McCall et al., 1985), to determine the oxygen abundance because 1] the auroral lines such as [OIII]$\lambda$4363 and [NII]$\lambda$5755 were too weak to be detected and 2] the [OII]$\lambda\lambda$3727, 3729 doublet and [OII]$\lambda\lambda$7320, 7331 lines are outside of the SIFS spectral range. However, several studies (e.g. Pilyugin et al., 2012; Pilyugin & Grebel, 2016; Marino et al., 2013; Denicoló et al., 2002) offer a series of calibrators in order to estimate the oxygen abundance using a variety of available emission lines and based on statistical studies.

Denicoló et al. (2002) and Marino et al. (2013) give the so-called $N2$ calibrator based on the [NII]$\lambda$6584 / H$\alpha$ line ratio, with slightly different regression coefficients. First introduced by Alloin et al. (1979), Marino et al. (2013) detailed the so-called $O3N2$ calibrator based on the [OIII]$\lambda$5007 / H$\beta$ $\times$ H$\alpha$ / [NII]$\lambda$6584 using the CALIFA survey. Introduced by Pilyugin & Grebel (2016), the so-called $S$ calibrator used the ratios $([NII]6548+6583)/H\beta$, $([SII]6717+6732)/H\beta$ and $([OIII]4959+5007)/H\beta$ to derive the oxygen abundance. It’s important to note from its definition that the $S$ calibrator depends on the reddening correction.

For Tol 1004-296, we were able to implement the three methods. Figures 7a, b, and c are respectively using the $S$ calibrator, the $O3N2$ calibrator, and $N2$ calibrator. Even if the overall aspect of the maps are similar, showing gradient between the different knots, we can notice some differences between the different methods. The $S$ calibrator (Fig. 7a) shows the greater contrast in the map between different parts of the galaxy. Even if it’s not the minimum of metallicity, the Knots 1, 2, and 3 show respectively a metallicity of 8.14 $\pm$ 0.04 dex, 8.14 $\pm$ 0.02 dex, and 8.10 $\pm$ 0.06 dex. On the outskirts (at 2.8 ″from the Knots centers) the metallicity is higher, from 8.20 $\pm$ 0.05 dex, to 8.55 $\pm$ 0.05 dex. On the other hand, the $O3N2$ calibrator (Fig. 7b) gives a lower amplitude between the metallicity in the emission knots and around them (the highest metallicity values do not exceed 8.22 $\pm$ 0.04 dex,). Finally, the $N2$ calibrator shows, in a clearer way, that higher metallicity is surrounding the emission knots.

We can notice that the metallicity of the different emission knots is between 8.06 $\pm$ 0.03 dex, and 8.16 $\pm$ 0.04 dex, except in Knot 7 (where it is 8.26 $\pm$ 0.3 dex). It is worth notice that a lower metallicity seems to be present where the emission is higher. This is more noticeable for the Knot 1 and 2 and when looking at the calibrators O3N2 and N2.

Figure 21 shows the metallicity map (using the S calibrator) for the MUSE data of Tol 1004-296. It can be noticed that the metallicity is slightly higher on the outskirts of the galaxy. In the center of the galaxy (essentially the field of view of the SIFS data), the metallicity ranges between 8.07 $\pm$ 0.04 dex and 8.30 $\pm$ 0.04 dex. As mentioned before, we also observed that lower metallicity corresponds to high emission knots. All seven emission knots (see Fig.  2) correspond to a minimum metallicity value.

In Tol 0957-278 case, we have implemented the N2 and O3N2 calibrators from Marino et al. (2013) (because the [NII]$\lambda$ 6548 line was too weak to be detected). Figure 7d shows the metallicity deduced from the $N2$ calibrator. It is clear here that the metallicity is lower where the emission is higher. The lowest metallicity is 7.84 $\pm$ 0.16 dex in Knot 2. The brightest site (Knot 1) has a metallicity of 8.09 $\pm$ 0.02 dex. Around Knot 1 (at a radius of 1.8 ″from the brightest spot), the metallicity ranges between 8.18 $\pm$ 0.02 dex and 8.25 $\pm$ 0.03 dex.

Photoionization from UV photons coming from hot stars, photoionization from an AGN source, or shocks from collisions (from stellar winds and supernovae (Dopita et al., 2003)), are the three ionization mechanisms could be present in an emission-line system. Baldwin et al. (1981) and Veilleux & Osterbrock (1987) came up with a series of excitation-dependent line ratios, the so-called ionization diagnostic diagrams. The three diagnostic plots usually are [OIII]$\lambda$5007 / H$\beta$ versus [NII]$\lambda$6584 / H$\alpha$, [OIII]$\lambda$5007 / H$\beta$ versus [S II]$\lambda\lambda$6717, 6731 / H$\alpha$, and [OIII]$\lambda$5007 / H$\beta$ versus [O I]$\lambda$6300 / H$\alpha$. The way it works is to mix in the plot hard ionizing emission lines ratio ([OIII]$\lambda$5007 / H$\beta$) with low ionization lines ([NII]$\lambda$6584, [S II]$\lambda\lambda$6717, 6731) known to be present where hydrogen is partially ionized. High [NII]$\lambda$6584 / H$\alpha$, or [S II]$\lambda\lambda 6717,6731$ / H$\alpha$ ratios are usually indicators that the ionization mechanism is more susceptible to be shock and/or AGN than photoionization from hot star UV photons. Then the three diagnostic diagrams are divided into three zones: photoionization from hot star UV photons (commonly noted Hii regions), AGN, and Shocks (noted as LINER). The diagnostic diagrams, commonly called BPT diagrams, from the name of the authors of Baldwin et al. (1981), are presented in Figure 8, Figure 22, and Figure 9. Except for the MUSE data, we were not able to detect the $[OI]\lambda 6300$ line, so for data coming from SIFS - SOAR, we only present the first two BPT diagrams.

The other interest of this study is the kinematical analysis aspect. In that regard, we observed both objects using the most resolving grating, M1500M3 available for SIFS. The M1500M3 grating gives a resolution of 9500@5500Å⁹⁹9https://noirlab.edu/science/programs/ctio/instruments/sifs/characteristics. The spectral sampling at H$\alpha$ is 0.2877Å (13 km s^-1). Even with this grating, we also used the data observed with the R700M grating to produce velocity and velocity dispersion maps (not presented here) in order to check our findings. In this section, we will describe the velocity and velocity dispersion maps for both galaxies. We also look at the kinematical diagnostic diagrams (Bordalo et al., 2009), applying statistical methods in order to determine possible systematic motions. We will also explore in more detail the line emission profiles in order to describe different degrees of symmetry. As a note, observations done using the M1500M3 grating do not match exactly the fields location showed in Figure 1. Small differences can be noted.