1. Introduction
Neutral hydrogen (Hi) exists in multiple phases in a galaxy’s interstellar medium (ISM). Hi is observed in two long-lived phases, the warm neutral medium (WNM, 5000–10,000 K) and cold neutral medium (CNM, 20–200K) (McKee & Ostriker, Reference McKee and Ostriker1977). In a pressure equilibrium these two phases will coexist, but that equilibrium is dependent on metallicity (Wolfire et al., Reference Wolfire, Hollenbach, McKee, Tielens and Bakes1995, Reference Wolfire, McKee, Hollenbach and Tielens2003; Bialy & Sternberg, 2019). Turbulence and colliding flows can drive the formation of CNM (Hennebelle & Audit, Reference Hennebelle and Audit2007; Kim & Ostriker, Reference Kim and Ostriker2017). The CNM is a precursor to the formation of the dense cores of molecular hydrogen (
$\mathrm{H}_2$
) from which stars form. The fraction of Hi in the CNM state is an important metric for the efficiency of star formation and galaxy evolution, along with the molecular gas fraction (Krumholz et al., Reference Krumholz, McKee and Tumlinson2009; Kennicutt & Evans, Reference Kennicutt and Evans2012).
We can use the spin temperature (
$T_{\mathrm{S}}$
) of Hi to assess the fraction of cold gas. The spin temperature of Hi is the excitation temperature of the Hi 21-cm spin-flip transition. As the gas is thermalised by collisions in the dense CNM, the spin temperature will be equal to the kinetic temperature of the gas in this environment (Field, Reference Field1958). In the WNM, it will be a lower limit for the kinetic temperature (Liszt, Reference Liszt2001).
By combining absorption observations with adjacent emission observations we are able to measure the spin temperature of the gas (Murray et al., Reference Murray, Stanimirović and Goss2015; Jameson et al., Reference Jameson, McClure-Griffiths and Liu2019). Hi absorption against continuum sources allows us to directly detect the presence of cold gas clouds. These clouds are otherwise difficult to detect as, due to their low spin temperature, they produce only weak 21-cm emission.
The Small Magellanic Cloud (SMC) is the perfect laboratory to study cold gas formation at high resolution in a low metallicity environment. The SMC is a nearby (
$61\ \mathrm{kpc}$
; Graczyk et al., Reference Graczyk, Pietrzyński and Thompson2014) low mass galaxy, part of the interacting Magellanic System, along with the Large Magellanic Cloud, the Leading Arm and the trailing Magellanic Stream. The SMC has metallicity of 0.2 solar (Russell & Dopita, Reference Russell and Dopita1992) and thus will have a typical temperature range for the CNM of
$50 \le {\mathrm{T_S}} \le 100\ \mathrm{K}$
(Bialy & Sternberg, 2019).
Dickey et al. (Reference Dickey, Mebold, Stanimirović and Staveley-Smith2000) presented the first survey of Hi absorption across the SMC, examining emission and absorption towards 32 continuum sources. They achieved optical depth noise levels of
$ 0.05 \leq \sigma_{\tau} \leq 0.203$
at
$0.825\ \mathrm{km\ s}^{-1}$
spectral resolution and detected significant absorption in 13 of these spectra. The sources were selected for their strong continuum flux. Jameson et al. (Reference Jameson, McClure-Griffiths and Liu2019) built on this by examining emission and absorption towards 55 continuum sources in 21 fields. These fields were selected for the greatest likelihood of finding absorption, with strong continuum flux (
$S_{\mathrm{cont}} > 50\ \mathrm{mJy}$
) and towards high column density regions (
$N_{\mathrm{HI}} > 4 \times 10^{20}\ \mathrm{cm}^{-2}$
). In their
$\ge 10$
hr dwell time per target they reached noise levels of
$ 0.01 \leq \sigma_{\tau} \leq 1.28$
and detected absorption in 37 of the spectra. Their spectra were imaged at
$0.2\ \mathrm{km\ s}^{-1}$
but smoothed to
$0.6\ \mathrm{km\ s}^{-1}$
for analysis. Both surveys used the Australia Telescope Compact Array with baselines up to 6 km (
$\mathrm{FWHM}_{\mathrm{beam}} \approx 5\ \mathrm{arcsec}$
) for absorption and compared them to emission from the Stanimirovic et al. (Reference Stanimirovic, Staveley-Smith, Dickey, Sault and Snowden1999) survey of the SMC (
$\mathrm{FWHM}_{\mathrm{beam}} \approx 98\ \mathrm{arcsec}$
). Table 1 summarises the observational parameters of these surveys and the survey presented in this work.
Comparison of SMC absorption survey parameters
a Dickey et al. (Reference Dickey, Mebold, Stanimirović and Staveley-Smith2000) and Jameson et al. (Reference Jameson, McClure-Griffiths and Liu2019) were targeted surveys which affects the detection rate.
Using the Australian Square Kilometre Array Pathfinder (ASKAP) telescope, the Galactic ASKAP (GASKAP; Dickey et al., Reference Dickey, McClure-Griffiths and Gibson2013) survey will observe Hi and OH in the Galactic Plane and the Magellanic System with unprecedented detail. The observations use high angular resolution (
$\mathrm{FWHM}_{\mathrm{beam}} \approx 16\ \mathrm{arcsec}$
) and high spectral resolution (
$0.24 \ \mathrm{km\ s}^{-1}$
per channel). Planned observations include long dwell times on the Magellanic Clouds and the low latitude Galactic Plane and shorter dwell times on the Magellanic Stream and Bridge. The repeated observations of key fields will provide Hi absorption spectra with a flux density sensitivity of
$\sigma_{\mathrm{S}} = 0.5 \ \mathrm{mJy}$
(Dickey et al., Reference Dickey, McClure-Griffiths and Gibson2013). The survey is planned to cover
$13,020 \ \mathrm{deg}^2$
in total. With a rate of
$\approx 10$
sources per square degree, up to 130,200 Hi absorption spectra are expected in the GASKAP-HI survey. With such high volumes of spectra, a repeatable process is essential.
In this work, we present the GASKAP Hi absorption pipeline and use it with the pilot phase I SMC Hi observations (see Pingel et al., Reference Pingel, Dempsey and McClure-Griffiths2022) to explore the distribution of cold gas in the SMC in an unbiased way. In section 2, we describe the GASKAP pilot observations of the SMC. In section 3, we present our Hi absorption pipeline along with the processing parameters used. We describe the observed absorption in the SMC and surrounds in section 4. In section 5, we discuss the results and their implications. Finally, in section 6, we summarise our findings.
2. Observations
The SMC was the first of three adjacent Magellanic fields targeted during the GASKAP Pilot Phase I observations. Two 12-hour observations of the field (ASKAP scheduling blocks 10941 and 10944) were taken in December 2019 using the standard GASKAP-HI observing configuration. The closepack-36 phased-array feed (PAF) footprint was used with a pitch of
$0.9$
deg and 3 interleaves for even coverage across the 25 square degree field. The field was centred on J2000
$ \mathrm{RA} = 00^{\mathrm{h}}58^{\mathrm{m}}43.280^{\mathrm{s}}$
,
$\mathrm{Dec} = -72^{\mathrm{d}}31^{\mathrm{m}}49.03^{\mathrm{s}}$
. The zoom-16 mode, with
$\Delta\nu = 1.15$
kHz, was used to provide a spectral resolution of
$\Delta v \sim 0.24\,\mathrm{km\ s}^{-1}$
. The observed band covered 18.5 MHz centred on 1419.85 MHz with 15558 channels, however only the 2048 channels covering the Milky Way and SMC velocity ranges were processed. The flagging was even across the field, providing a change in RMS across the field of
$\leq 7.5 \%$
.
As noted in Dickey et al. (Reference Dickey, McClure-Griffiths and Gibson2013), ASKAP’s bimodal baseline distribution provides excellent capabilities for measuring both Hi emission and absorption. In these observations, there was good coverage of baselines longer than 2000m, which provide fine spatial resolution. Baselines up to 6000m were present. The expected sensitivity from the combined observation is 3.3 mJy/beam after accounting for flagging and excluded baselines (see Section 3.2.1).
The data were flagged and calibrated using the standard ASKAPSoft pipeline (Hotan et al., Reference Hotan, Bunton and Chippendale2021) with configuration suitable for wide-field emission. A continuum image and a continuum source catalogue were also produced by the ASKAPSoft pipeline for each observation. The observations and the initial processing with ASKAPSoft are described in further detail in Pingel et al. (Reference Pingel, Dempsey and McClure-Griffiths2022).
3. Hi Absorption Pipeline
To measure the absorption against continuum sources we need a spectral line cube for the region surrounding each source. We have two potential strategies to produce these source cubes:
-
a) producing a large, continuum included cube covering the entire field and then extracting sub-cubes for each source, or
-
b) producing sub-cubes for each source position directly from the measurement sets.
Note, we cannot use the GASKAP emission cube as the continuum subtraction will limit our ability to accurately measure deep absorption, and the cube contains large-scale emission which will add noise to the absorption spectrum. We have chosen to take approach (b) as we know in advance the positions we wish to image. This approach also has several advantages for our use case. This saves us from imaging the unused regions (
$> 99$
%) of the cube not sampled by the sparse and small continuum sources. Moreover, we are able to set the phase-centre to the source position for each sub-cube, thus avoiding w-term effects which can impact wide fields (Cornwell & Perley, Reference Cornwell and Perley1992). Finally we are able to optimise the cube production to obtain the most accurate absorption spectra possible.
We have developed a pipeline (Dempsey, Reference Dempsey2022) to produce Hi absorption spectra for a large number of sources from the calibrated ASKAP measurement sets. As GASKAP is a wide-field survey, and not targeted at specific sources, we have the opportunity to take an unbiased sample of the cold Hi absorption in every field. We obtain absorption spectra towards all continuum sources above a given brightness (e.g.,
$S_{\mathrm{cont}} \geq 15\ \mathrm{mJy}$
). For each source, we produce a small cube around the source, sufficient to include all continuum emission from the source. From this we extract an integrated absorption spectrum. We also measure the emission surrounding the source from the GASKAP emission cube for the field (Pingel et al., Reference Pingel, Dempsey and McClure-Griffiths2022). Finally, we then process all spectra to find any significant absorption detections and produce catalogues of spectra and absorption features.
The pipeline takes as input:
-
1. A continuum source catalogue, the ASKAP Selavy (Whiting & Humphreys, Reference Whiting and Humphreys2012) component catalogue
-
2. A collection of calibrated ASKAP beam measurement sets.
The pipeline is split into three phases: preparation, imaging, and spectra extraction, each described below.
3.1. Preparation Phase
In the preparation phase we identify the target sources and the input data for the pipeline run. We take the continuum component catalogue produced as part of the ASKAPSoft processing of the field (Hotan et al., Reference Hotan, Bunton and Chippendale2021). We designate all sources with
$S_{\mathrm{cont}} \geq 15\ \mathrm{mJy}$
as target sources. A typical GASKAP-HI observation for 20 hrs has a
$\sigma_{\mathrm{S}} = 1.89\ \mathrm{mJy/beam}$
(Pingel et al., Reference Pingel, Dempsey and McClure-Griffiths2022), and after the baseline cutoff described in Section 3.2
$\sigma_{\mathrm{S}} = 2.80\ \mathrm{mJy/beam}$
. For sources with
$S_{\mathrm{cont}} = 15\ \mathrm{mJy}$
we have a theoretical optical depth noise level of
$\sigma_{\mathrm{cont}} = 0.19$
. Examining the 47 sources near this threshold (
$15 \leq {\mathrm{S_{cont}}} \le 17\ \mathrm{mJy}$
) in the SMC pilot data, we see a median optical depth noise level of
$\sigma_{\mathrm{cont,median}} = 0.23$
with a range of
$0.07 \leq \sigma_{\mathrm{cont}} \leq 1.34$
, showing that we are at the limit of useful spectra in these data. For the 200 hr integrations on the Magellanic Clouds planned in the full survey, we could expect to push this limit down to as low as
$S_{\mathrm{cont}} \sim 5\ \mathrm{mJy}$
with similar optical depth noise levels.
For each target source we use all beams within 0.55 deg of the source (or 0.8 deg for edge sources where no beams are closer). From these we produce a list of all target sources and the beams which will be used to produce the cutout for each source.
3.2. Imaging Phase
The imaging phase is the most computationally intensive. As discussed earlier, we produce a sub-cube around each source. We process this phase in parallel with a job per sub-cube. In the pipeline, we dynamically schedule the sub-cube extraction jobs to minimise competing use of measurement sets and thus improve overall I/O throughput while providing high parallelism within the resources of the compute infrastructure.
For each cube we use the CASA (McMullin et al., Reference McMullin, Waters, Schiebel, Young and Golap2007) task tclean to image a 50 arcsec x 50 arcsec region around the source at 1 arcsec per pixel and with
$0.24\ \mathrm{km\,s}^{-1}$
LSRK velocity resolution (the native spectral resolution). This images the beam measurement sets selected for each source in the preparation phase. We use a 1.5 k
$\lambda$
baseline length cutoff, natural weighting, primary beam correction and light cleaning (1000 total iterations), as discussed in Section 3.2.1. This gives a typical synthesised beam size of
$16 \times 14\ \mathrm{arcsec}^2$
. We use a primary beam model of a 12
$\,\textrm{m}$
dish with a 0.75
$\,\textrm{m}$
blockage.
3.2.1. Imaging Parameters for GASKAP Absorption
For absorption studies we wish to exclude emission from the spectra while maintaining high signal-to-noise ratios. In this pipeline we have achieved that primarily by setting a minimum baseline length when producing the cutout cube. Figure 1 shows the results of our trials of a range of baseline length cutoffs on a sample of seven spectra selected to represent the diverse range of spectra seen in the full dataset. Based on these trials, we found that the optimum balance of reduced emission and reduced noise came from a 1.5 k
$\lambda$
(315 m) baseline length cutoff. This retains 574 of the 630 ASKAP baselines but provides sensitivity only to features of
$\lesssim3$
arcmin in size. This is well suited for analysis of compact extra-Galactic sources, but means these data products would not be suitable for analysis of larger Galactic structures such as supernova remnants and Hii regions.
Box-plot comparison of maximum levels of noise from emission (where
$e^{-\tau} > 1$
; top) and optical depth noise (bottom) for each of the sample spectra for a range of baseline length cutoffs. In these plots the central horizontal line is the median, the ends of the box are the 25th and 75th percentiles and the top and bottom lines show the maximum and minimum values respectively. Diamonds show outliers based on their distance from the interquartile range.
Similarly, we tested a variety of weighting parameters, as shown in Figure 2. Based on these results we found ‘natural’ weighting produced cutout cubes with the lowest optical depth noise and emission noise. Additionally, we used primary beam correction and light cleaning (1000 total iterations) when imaging the data.
Box-plot comparison of maximum levels of noise from emission (where
$e^{-\tau} > 1$
; top) and optical depth noise (bottom) for each of the sample spectra for a range of weighting parameters. See Fig. 1 for details of the ranges.
3.3. Spectra Extraction Phase
From the cutout cubes for each source, we extract absorption spectra for every target source. We use the technique developed in Dempsey et al. (Reference Dempsey, McClure-Griffiths, Jameson and Buckland-Willis2020), but we describe the process here for completeness. The continuum source catalogue defines a source ellipse for each component. We combine all pixels within this source ellipse to produce the spectrum. We define a line-free region of the cube (with a velocity range
$-100 \leq \mathrm{v}_{\mathrm{LSRK}} \leq -60\ \mathrm{km\,s}^{-1}$
for the SMC) and measure the mean brightness within this range for each pixel. We then weight each pixel’s contribution to the spectrum by the square of the mean brightness, as described in Dickey et al. (Reference Dickey, Brinks and Puche1992). Lastly, we divide the combined spectrum by its mean brightness within the line-free region to produce the absorption spectrum,
$e^{-\tau}$
.
We estimate the noise in the spectrum using a combination of the noise in the off-line region and emission in the primary beam of the dish, as described in Jameson et al. (Reference Jameson, McClure-Griffiths and Liu2019). The standard deviation of the spectrum in the line-free region is taken as the base noise for the spectrum. In order to model the increase in system temperature due to emission received by the antenna at different frequencies, we then measure the emission level in the ASKAP primary beam of 62 arcmin using data from the Parkes Galactic All-Sky Survey (GASS; McClure-Griffiths et al., Reference McClure-Griffiths, Pisano and Calabretta2009; Kalberla & Haud, Reference Kalberla and Haud2015). We average the GASS emission across a 7 pixel (33 arcmin) radius annulus centred on the source position, with a 1 pixel exclusion at the centre. The 1
$\sigma$
noise envelope for the GASKAP absorption spectrum is then calculated as:
\begin{equation}\sigma_{\tau}(v) = \sigma_{\mathrm{cont}} \frac{T_{\mathrm{sys}} + \eta_{\mathrm{ant}}T_{\mathrm{em}}(v)}{T_{sys}},\end{equation}
where the system temperature
$T_{\mathrm{sys}} = 50\ \mathrm{K}$
and antenna efficiency
$\eta_{\mathrm{ant}} = 0.67$
(Hotan et al., Reference Hotan, Bunton and Chippendale2021),
$\sigma_{\mathrm{cont}}$
is the standard deviation of the line-free region of the spectrum, and the mean emission
$T_{\mathrm{em}}(v)$
at each velocity step is as measured in GASS.
We identify any absorption features in the spectra during this phase. We classify absorption features as those having one channel of
$3 \sigma$
absorption and an adjacent channel of
$\geq 2.8 \sigma$
. We expand the feature to include any adjacent channels of at least
$2.8 \sigma$
significance. These criteria were chosen empirically to detect even shallower absorption whilst still avoiding noise spikes being detected as features. Statistically, these criteria give an average of less than one false positive absorption feature per
$\approx 370$
spectra.
With a large number of spectra expected in each field, it is important to consistently identify which spectra are of sufficient quality to be analysed. We use a subset of the Brown et al. (Reference Brown, Dickey, Dawson and McClure-Griffiths2014) rating system to classify the spectra quality level from A to D. Rating A spectra pass all tests, while each spectrum’s rating is reduced by one step for each failed test until rating D spectra fail all tests. The tests described in Brown et al. aim to flag unphysical or noisy spectra.
The three tests we use are:
-
Continuum noise—The
$1\sigma$
noise level (
$\sigma_{\mathrm{cont}}$
) in optical depth (
$e^{-\tau}$
) must be less than
$1/3$
-
Max signal to max noise—The ratio of the deepest absorption to the highest emission noise (i.e.,
$e^{-\tau} > 1$
) in the spectrum:
$(1-{\mathrm{min}}( e^{-\tau}))/({\mathrm{max}}( e^{-\tau}) -1) \ge 3$
-
Optical depth range—The range from the deepest absorption to highest emission noise in the spectrum:
$ {\mathrm{max}}(e^{-\tau}) - {\mathrm{min}}(e^{-\tau}) < 1.5$
.
We do not use the other two Brown et al. tests as they are not applicable to our spectral extraction method.
To make measurements of spin temperature it is also necessary to have an emission spectrum for each source. The emission spectra are calculated from the full field continuum-subtracted GASKAP emission cube produced for the same observations (Pingel et al., Reference Pingel, Dempsey and McClure-Griffiths2022). We calculate the emission spectrum of a source as the mean value of an annulus of radius 56 arcsec (8 pixels,
$\approx 2$
beam widths) around each source with a central exclusion of 28 arcsec (4 pixels,
$\approx 1$
beam width). This allows us to measure the mean emission in the close vicinity of the source without measuring the source itself, in all but the most extended cases. The high resolution gives us as close an approximation for the emission in the line of sight of the source as possible. We use the standard deviation of the annulus as the emission noise level.
The pipeline produces a series of outputs describing the spectra towards each source. For each source we produce a spectrum votable file with the absorption and, where available, the emission spectra, along with the noise in both the absorption and emission spectra. We output catalogues of the spectra, (see Table 2), as well as the absorption features detected (see Table 3).
Sample of the GASKAP spectrum catalogue. This is a sample of the key fields from the GASKAP-HI absorption spectrum catalogue for 14 sources. The full catalogue of all 229 sources is available in the dataset (Dempsey et al., Reference Dempsey, Mcclure-Griffiths and Murray2022).
Sample of the GASKAP absorption feature catalogue. This is a sample of the key fields for those Hi absorption features detected in the spectra listed in Table 1. Note that multiple features are detected in some of the spectra, while other spectra have no detectable features. The full catalogue of all 130 features is available in the dataset (Dempsey et al., Reference Dempsey, Mcclure-Griffiths and Murray2022).
4. SMC Hi Absorption
For each of the 373 continuum sources which satisfy our criteria of peak flux density
$\ge 15\ \mathrm{mJy}$
(see Sec. 3.1), we produced an absorption spectrum using the process described in Section 3.3. However, we averaged the data to
$1\ \mathrm{km\ s}^{-1}$
spectral resolution in the imaging phase to reduce noise. We selected the sources from the Selavy (Whiting & Humphreys, Reference Whiting and Humphreys2012) continuum source catalogue from the SB 10941 observation. We then excluded any sources with
$\sigma_{\mathrm{cont}} > 0.3$
or at the edges of the field where beam power
$< 80\%$
, leaving 229 sources. The source density is 10 sources per square degree, a rate which we expect to be replicated in most GASKAP fields of similar noise. The distribution of all sources by noise is shown in Figure 3. Notably, all of the excluded noisy sources are towards the edges of the field.
Distribution of continuum sources showing their optical depth noise against the SMC Hi column density map from GASKAP (Pingel et al., Reference Pingel, Dempsey and McClure-Griffiths2022). Triangles are sources excluded due to either high noise or being on the edges of the cube, squares are sources against which absorption was detected, and circles are other sources. Darker colours indicate lower optical depth noise.
An example spectrum for continuum source
$\mathrm{J}005556-722605$
is shown in Figure 4. In the top panel, we can see that the spectral bandpass has a slope of
$\approx0$
and is consistent with random noise around a constant continuum level. The emission has been successfully excluded from the spectrum, with the spectrum above the continuum level consistent with the continuum noise envelope. In the bottom panel, we show the brightness temperature spectrum from the GASKAP emission cube. Note that the velocity range of the emission data is restricted to the range of the SMC emission cube. Overall, the approach of excluding short baselines from the imaging has been successful in excluding emission while maintaining a high signal-to-noise ratio. The full set of spectra are available in the dataset (Dempsey et al., Reference Dempsey, Mcclure-Griffiths and Murray2022).
Absorption (top) and emission (bottom) spectra for source
$\mathrm{J}005556-722605$
in the SMC velocity range. The measured absorption is shown as a black line, the continuum level is shown as a red line, the 1
$\sigma$
noise envelope is shaded grey, and the dotted orange line is the 3
$\sigma$
absorption level. Regions of detected absorption features are shaded in both the absorption and emission spectra. In the emission plot, the black line shows the brightness temperature, and the 1
$\sigma$
uncertainty in the brightness temperature is shown as a grey envelope. Note that the SMC emission cube does not cover the velocity range
$\mathrm{v}_{\mathrm{LSRK}} > 250\ \mathrm{km\,s}^{-1}$
.
A subset of sources are described in Table 2, with the full set of sources provided in the dataset. For each source we provide
-
a) ID—the id of the source within the field
-
b) Source—name based on the detected location of the source
-
c) RA—right ascension of the source
-
d) Dec—declination of the source
-
e) Peak Flux—peak flux density of the source from the Selavy continuum catalogue
-
f) Rating—quality rating of the spectrum (see Sec 3.3)
-
g)
$\sigma_{\mathrm{cont}}$
- 1
$\sigma$
noise level of the spectrum in absorption units. Does not include emission noise -
h) Peak
$\tau$
—the maximum
$\tau$
value within the SMC velocity range. Where the spectrum is saturated (
$e^{-\tau} \le 0$
) a minimum limit is specified based on the
$1 \sigma$
noise level of the channel in the spectrum -
i)
$N_{\mathrm{HI,uncorr}}$
—column density towards the continuum source, excluding Milky Way velocities, from the GASS survey under the assumption that the Hi is optically thin -
j)
$\langle T_{\mathrm{S}} \rangle$
—the density-weighted mean spin temperature of the sight-line (see Section 5.3) -
k)
$\mathcal{R}_{\mathrm{HI}}$
—column density correction ratio for the sight-line (see Section 5.2).
Using the criteria described in Sec. 3.3, we detect absorption features at SMC velocities (
$\mathrm{v}_{\mathrm{LSRK}} \geq 50\ \mathrm{km s}^{-1}$
) in 65 (28%) of the spectra. Thirty-six of these spectra have multiple features, giving a total of 122 features. Note that a full decomposition of these features into components is beyond the scope of this work. A subset of absorption features is shown in Table 3, with the full set of features provided in the dataset. For each feature we provide
-
a) Source—name of the source based on the location of the source
-
b) Feature—name of the feature based on the source name and the minimum velocity of the feature
-
c) Min Velocity—the minimum velocity bound of the feature
-
d) Max Velocity—the maximum velocity bound of the feature
-
e) Width—the number of
$1\ \mathrm{km\ s}^{-1}$
channels the feature spans -
f) Peak Absorption - the maximum measured absorption (
$1-e^{-\tau}$
) of the feature with the uncertainty in absorption spectrum at that velocity, values greater than 1 are saturated -
g) Peak
$\tau$
—the maximum
$\tau$
value of the feature. Where the spectrum is saturated (
$e^{-\tau} \le 0$
) a minimum limit is specified based on the noise level in the peak channel -
h) Significance—the highest single-channel significance of the feature as measured in absorption (
$e^{-\tau}$
) -
i) Equivalent Width—the integral of
$\tau$
for the feature. This will be a lower limit for saturated spectra as we use a minimum limit based on noise for each saturated velocity channel in the same way as peak
$\tau$
above.
The distribution of optical depth of these features is shown in Figure 5. Sixteen of the features are saturated (
$e^{-\tau} < 0$
) and 6 features have sufficient noise that they could be saturated (above the noise limit in Fig. 5). The
$\tau$
values of the 100 remaining features range from 0.04 to as high as 3.3, with a median
$\tau = 0.5$
. The trend to have higher
$\tau$
values for higher noise spectra is a side-effect of our minimum significance requirement.
Distribution of the peak optical depth of detected absorption features against the noise in optical depth. Blue dots are non-saturated features in the body of the SMC, orange triangles are non-saturated features outside the body of the SMC. Blue plus signs show saturated (
$e^{-\tau} < 0$
) features in the body of the SMC and the green line shows our
$\tau$
sensitivity limit due to noise.
Following Jameson et al. (Reference Jameson, McClure-Griffiths and Liu2019), we define an arbitrary column density limit for the body of the SMCFootnote a (
$N_{\mathrm{HI,uncorr}} > 2 \times 10^{21} \ \mathrm{cm}^{-1}$
) and show this as the contour in Figure 6. Within the body of the SMC, we see a much higher detection rate, as shown in Table 4, with absorption features in 62% of spectra. Thus 75% of all spectra with absorption are found in the body of the SMC. Notably, of the 23 bright (
$S_{\mathrm{cont}} \ge 50\ \mathrm{mJy}$
) continuum sources within the body of the SMC only one,
$\mathrm{J}003754{-}725157$
, does not have any detectable absorption features in its spectrum. This spectrum has
$\tau_{\mathrm{peak}} = 0.11\pm0.06$
. So in these higher density regions we see absorption in almost all low-noise spectra. The detection rate also drops with source strength, with only 41% of the 51 SMC spectra against faint (
$S_{\mathrm{cont}} < 30\ \mathrm{mJy}$
) continuum sources having detectable absorption.
Detection statistics for different regions
Locations of the 122 SMC absorption features described in Sec. 4, plotted against
$10 \ \mathrm{km\,s}^{-1}$
slices of emission. Each absorption feature is shown as a circle centred on its position with the colour reflecting the depth of absorption. The scale for the absorption is shown in the top right of each panel, with darker colours indicating deeper absorption and lighter colours shallower absorption. Features are plotted in the velocity slice containing the centre of the absorption feature’s velocity range. In each slice, the column density for that velocity range, as measured in GASKAP emission (Pingel et al., Reference Pingel, Dempsey and McClure-Griffiths2022), is shown as a linear grayscale, with the scale shown on the right of each row. The green contour shows the 2
$\times 10^{21}$
column density limit for all SMC velocity ranges, as detected in GASKAP emission (Pingel et al., Reference Pingel, Dempsey and McClure-Griffiths2022), representing the outline of the SMC. Locations of the SMC absorption features plotted against
$10\ \mathrm{kms}^{-1}$
slices of emission continued.
The 16 features outside the body of the SMC have generally low
$\tau$
values within a smaller range (min=0.09, med=0.5 max=1.3) than within the body, and none are saturated. In contrast to the body of the SMC, only one of the 16 spectra outside the SMC has multiple features. Of the 71 bright continuum sources outside the SMC, only 7 show detectable absorption features. As shown in Figure 3, some detections such as
$\mathrm{J}013218{-}715348$
and
$\mathrm{J}005715{-}704046$
are well away from the SMC. Also, as seen in Figure 6, particularly at the extremes of velocity, we detect some shallow absorption features in regions with little emission. These are discussed in further detail in Section 4.1.
4.1. Absorption in Low Density Regions
Six of the detected absorption features are in velocity regions with little emission (
$T_{\mathrm{B}} < 5\ \mathrm{K}$
), indicating low column density (see Fig. 6). This is surprising as cold Hi is normally associated with a higher density of warm Hi (e.g., Kanekar et al., Reference Kanekar, Braun and Roy2011). Three of these features are reliable detections. The first feature,
$\mathrm{J}005732{-}741243\_137$
is a very clear detection (significance
$\approx$
61
$\sigma$
) with deep absorption through low column density (see Fig. 7), and is also reported in Jameson et al. (Reference Jameson, McClure-Griffiths and Liu2019). The feature is also apparent as absorption in the GASKAP emission cube. Its location, seen in Fig. 6h, is in the lower centre of the image below the SMC body. The second feature,
$\mathrm{J}005715{-}704046\_255$
is outside the velocity range for which we have emission data and thus is not shown in Fig. 6. It is located well to the North of the SMC body. The feature is only two channels wide, but has adjacent channels which also show less significant absorption, making this likely to be a real absorption feature. A further absorption feature,
$\mathrm{J}011134{-}711414\_114$
, was detected within a velocity region with
$T_{\mathrm{B}} = 9\pm4\ \mathrm{K}$
. This source is shown in Fig. 6f and is located in the upper left of the image slightly to the North-East of the SMC body. With two adjacent channels over 3
$\sigma$
significance this is a reliable detection. Of these three features, the latter two are shallow, similar to the low-N(Hi) clouds found by Stanimirović et al. (Reference Stanimirović, Heiles, Kanekar and Goss2007) in the Milky Way. They found these type of clouds to be the lower column density population of the CNM but likely to be short-lived due to their small size and lack of shielding.
Absorption (top) and emission (bottom) spectra for source
$\mathrm{J}005732{-}741243$
in the SMC velocity range. See Figure 4 for details.
The other three features are less firm detections. In two of these cases (
$\mathrm{J}013134{-}700042\_240$
which is not plotted, and
$\mathrm{J}005116{-}734000\_65$
in Fig. 6a), the features are marginal detections and consistent with noise. Notably, the spectrum
$\mathrm{J}013134{-}700042$
was flagged as having a potentially underestimated noise level. Finally, the feature
$\mathrm{J}010524{-}722524\_87$
(see Fig. 6c) has a large emission noise spike in the next channel (
$v=89\ \mathrm{km s}^{-1}$
) to the feature. This may indicate that the absorption is a noise artefact also.
5. Discussion
5.1. Column density thresholds and detection rate
Using the GASKAP emission data we can calculate the Hi column density for each sight-line, under the assumption that the Hi gas is optically thin (
$\tau \ll 1$
) (Dickey & Lockman, Reference Dickey and Lockman1990, eq. 3):
\begin{equation} N_{\mathrm{HI,uncorr}} = 1.823 \times 10^{18} \int T_{\mathrm{B}}(\text{v})\,d\text{v}\,\text{cm}^{-2}\end{equation}
We use a Monte-Carlo method to establish the uncertainties in
$N_{\mathrm{HI,uncorr}}$
. For each velocity channel of each spectrum we draw 1000 samples of both the optical depth and the brightness temperature. The samples are randomly drawn from a normal distribution utilising the actual spectrum value as the mean and the 1-
$\sigma$
noise envelope as the standard deviation. We then perform all calculations across all samples, and use the median of the result for each sample as the final value, and the 15th and 85th percentiles as the 1-
$\sigma$
uncertainty ranges.
With an unbiased sample of absorption in the field we have the opportunity to examine where the cold gas is present. Figure 8 shows a comparison of the Hi column density (uncorrected for optical depth, see Sec. 5.2) for both spectra with and without absorption detections. We can see that we have some detections in regions with column densities as low as
$1.25 \times 10^{20}\ \mathrm{cm}^{-2}$
. However, it is not until
$2.5 \times 10^{21}\ \mathrm{cm}^{-2}$
that we have a high rate of detection. Most absorption detections are in sight-lines with a column density of
$3 \times 10^{21}\ \mathrm{cm}^{-2}$
or greater. From
$6 \times 10^{21}\ \mathrm{cm}^{-2}$
almost all sight-lines, even for faint continuum sources, show absorption.
Comparison of the distribution of spectra with and without absorption detections compared with the column density of the lines of sight as measured by GASKAP under the assumption that the Hi is optically thin. The blue bars show the count of spectra without detected absorption features in each column density bin, while the orange bars show the count of spectra with absorption features in those bins. The vertical dashed green line marks the
$2\times10^{21}\ \mathrm{cm}^{-2}$
column density limit of the SMC body.
Kanekar et al. (Reference Kanekar, Braun and Roy2011) presented a relationship between the observables of uncorrected column density and the integral of
$\tau$
(the equivalent width) in Milky Way gas clouds. The following equation (their Equation 3) models the presence of cold gas within warm Hi envelopes, which shield the cold Hi from surrounding hot gas and radiation.
\begin{equation} N_{\mathrm{HI,uncorr}} = N_0 e^{-\tau^{\prime}_{c}} + N_{\infty} (1-e^{-\tau^{\prime}_{c}})\end{equation}
Here,
$N_0$
is the threshold column density for the formation of cold gas,
$N_{\infty}$
is the column density at which the cold gas saturates (
$\tau \to \infty$
) and
$\tau^{\prime}_{c}$
is the effective opacity, or the mean observed optical depth over the velocity range (
$\Delta V$
). They further define a lower and upper limit using Eq. 3 on the expected range in which absorption forms in the Milky Way. The lower limit has
$N_0 = 10^{20}\ \mathrm{cm}^{-2}$
,
$N_{\infty} = 5.0 \times 10^{21}\ \mathrm{cm}^{-2}$
and
$\Delta V = 20\ \mathrm{km\,s}^{-1}$
. The upper limit has
$N_0 = 2 \times 10^{20}\ \mathrm{cm}^{-2}$
,
$N_{\infty} = 10^{22}\ \mathrm{cm}^{-2}$
and
$\Delta V = 10 \ \mathrm{km\,s}^{-1}$
.
Comparing our detections against the Kanekar et al. (Reference Kanekar, Braun and Roy2011) limits (see Figure 9), we see outliers on both sides of the envelope, despite most sources lying within the limits. The two sources to the left of the lower limit hint at cold gas formation below the column density limit seen in the Milky Way. In the current data we do not have the sensitivity to probe this lower column density regime in detail. To the right of the upper limit we find eight sources. Four of these sources have low equivalent width noise levels and thus support a higher limit for the condensation of dense Hi into molecular
$\mathrm{H}_2$
saturation in the SMC. This higher limit reflects the findings of both Bialy & Sternberg (Reference Bialy and Sternberg2016) and Krumholz et al. (Reference Krumholz, McKee and Tumlinson2009) that higher column density is required at low metallicities for the formation of molecular gas, and of Bolatto et al. (Reference Bolatto, Leroy and Jameson2011) that the SMC has a very low fraction of molecular gas as compared to its atomic gas content. With the addition of future GASKAP observations of the SMC we should be able to test these bounds more rigorously.
Comparison of uncorrected column density with integrated optical depth (Equivalent Width). The curved lines are the lower (blue) and upper (red) limits defined in Kanekar et al. (Reference Kanekar, Braun and Roy2011) section 3.
5.2. Corrected column density for the SMC field
The column density measurements from emission data are made under the assumption that the gas is optically thin (
$\tau \ll 1$
). However, if the Hi is not optically thin we need to account for the absorption to get an accurate estimation of the column density. For the sight-lines we have observed, we can correct these measurements for the actual optical depth of the gas. Using the assumption that the gas is isothermal, that is each velocity channel of gas only has a single temperature, we can calculate the corrected column density using the formula (Dickey & Benson, Reference Dickey and Benson1982, eq 5 and Chengalur et al., Reference Chengalur, Kanekar and Roy2013, eq 8):
\begin{equation}N_{\mathrm{HI,corr,iso}} = 1.823 \times 10^{18} \int \frac{T_{\mathrm{B}}({\mathrm{v}}) \times \tau({\mathrm{v}})}{1-e^{-\tau({\mathrm{v}})}}\,d{\mathrm{v}}\,{\mathrm{cm}}^{-2}\end{equation}
where
$T_{\mathrm{B}}({\mathrm{v}})$
and
$\tau({\mathrm{v}})$
are the brightness temperature and absorption respectively for a velocity channel. The column density correction ratio is then
\begin{equation} \mathcal{R}_{\mathrm{HI}} = \frac{N_{\mathrm{HI,corr,iso}}}{N_{\mathrm{HI,uncorr,iso}}},\end{equation}
We use a Monte-Carlo method to establish the uncertainties in
$N_{\mathrm{H,corr,iso}}$
and
$\mathcal{R}_{\mathrm{HI}}$
, as described in Sec. 5.1.
There is a large spread in the correction factors as compared to uncorrected column density, as shown in Fig. 10. Strikingly, the range of correction factors is very different between the wing and the bar of the SMC. The bar shows generally low correction factors with an increase in correction factor with increasing column density. The exception is the sight-line towards
$\mathrm{J}010029{-}713826$
, which shows both a deep and a wide absorption component in a region with relatively low column density for the bar region. In contrast, the wing shows a great variety of correction factors unrelated to column density. This likely reflects that the bar has a longer line of sight than the wing, and so multiple features would be blurred. In the shorter line of sight through the wing, we are more likely to pick out individual features such as diffuse tidal structures and regions of intense star forming. The cluster of four sight-lines with the highest column density and highest correction factor in the wing are all close to NGC 460 and NGC 465. The other two sight-lines with
$\mathcal{R}_{\mathrm{HI}} > 1.3$
are towards the SE end of the wing, another region of active star formation. The detections outside the wing and the bar, all at lower column densities, show small correction factors as expected. The one exception is the sight-line towards
$\mathrm{J}012924{-}733153$
, in which the correction factor is boosted by noise (
$\sigma_{cont} = 0.211$
) as well as the single narrow detection and a potential shallow, wide component. For comparison we also show the correction factors of low-noise (
$\sigma_{cont} < 0.1$
) sight-lines where no absorption was detected. These are mostly at lower column densities and have minimal correction, as expected.
Comparison of uncorrected column density with the correction factor due to optical depth measured for sight-lines with detected absorption.
Overall, the correction factors range from
$\mathcal{R}_{\mathrm{HI}} \approx 1$
at the low column densities, to
$\mathcal{R}_{\mathrm{HI}} = 1.49$
in the higher density regions of the wing, close to known star formation regions. Dickey et al. (Reference Dickey, Mebold, Stanimirović and Staveley-Smith2000) found a correction factor relationship with column density for the SMC of
$\mathcal{R}_{\mathrm{HI}} \simeq 1 + 0.667(\text{log}_{10}N - 21.4)$
. With our broader view of the SMC optical depth, we find a slightly shallower linear rise in correction factor with column density for just the SMC bar of
$\mathcal{R}_{\mathrm{HI}} \simeq 1 + 0.51(\text{log}_{10}N - 21.43)$
. The wide range of correction factors in the SMC wing means that a single line cannot be fit to that region. The wide range for the SMC wing is similar to the results of Nguyen et al. (Reference Nguyen, Dawson and Lee2019) who found large regional variation in
$\mathcal{R}_{\mathrm{HI}}$
around five giant molecular clouds.
The shallower relationship of column density and correction factor than found in Dickey et al. (Reference Dickey, Mebold, Stanimirović and Staveley-Smith2000) has two drivers. Firstly, as noted in Sec. 5.4, we detect a higher proportion of shallower absorption features due to the unbiased selection of continuum sources. Secondly, the smaller emission beam size will reduce beam smearing of small-scale dense structure, so in these cases we will measure a higher uncorrected column density. One of the great advantages GASKAP has over past absorption studies is that we have a very closely matched beam size between the emission (
$30 \times 30\ \mathrm{arcsec}^2$
) and absorption (
$15 \times 15 \ \mathrm{arcsec}^2$
). This means both that our emission sample is taken from as close as possible to the target continuum source and that the scale of structures detected by the emission and absorption are similar.
Comparing the cumulative distribution of correction factors between the field outside the SMC and the body of the SMC, we find that the body has far higher correction factors (see Fig. 13, top). The higher correction factors reflect the column density threshold apparent at
$N_{\mathrm{HI,uncorr}} < 2 \times 10^{21}\ \mathrm{cm}^{-2}$
in Fig. 10, below which there are much lower correction factors. The uncertainties in this top panel are generated from the
$\mathcal{R}_{\mathrm{HI}}$
uncertainties for each source.
5.3. Temperature distribution
To measure the density-weighted mean spin temperature we can compare the integrals of brightness temperature and optical depth using the following formula (Dickey et al., Reference Dickey, Mebold, Stanimirović and Staveley-Smith2000, Eq. 4):
\begin{equation} \langle T_{\mathrm{S}} \rangle = \frac{\int T_{\mathrm{B}}(\text{v})\,d\text{v}}{\int 1-e^{-\tau(\text{v})}\,d\text{v} } = \int \frac{n(\text{s})}{[n(\text{s})/T_{\mathrm{S}}(\text{s})]} d\text{s} ,\end{equation}
where
$n(\text{s})$
is the line-of-sight volume density. We use a Monte-Carlo method to establish the uncertainties in
$ \langle T_{\mathrm{S}} \rangle$
, as described in Sec. 5.1.
Following Dickey et al. (Reference Dickey, Mebold, Stanimirović and Staveley-Smith2000) we measure the integrals across the velocity range where
$T_B \geq 3\ \mathrm{K}$
so as to reduce the noise in the integrals. Note that in some cases the integral range will contain regions where the emission level dips below this threshold (e.g., where the sight-line has a bimodal emission distribution). In those cases we still include the emission from this inner velocity region in the integral.
Within the SMC (
$N_{\mathrm{HI,GASKAP,uncorr}} \geq 2\times10^{21}\ \mathrm{cm}^{-2}$
) we find an inverse noise-weighted mean spin temperature of
$\langle T_{\mathrm{S}} \rangle = 245\pm2\ \mathrm{K}$
. This is lower than the 350 K found by Dickey et al. (Reference Dickey, Mebold, Stanimirović and Staveley-Smith2000). Jameson et al. (Reference Jameson, McClure-Griffiths and Liu2019) found an arithmetic mean spin temperature of
$\langle T_{\mathrm{S}} \rangle = 117.2 \pm 101.7\ \mathrm{K}$
, reflecting their targeting of cold gas. Outside the SMC we find
$\langle T_{\mathrm{S}} \rangle = 156^{+9}_{-7}\ \mathrm{K}$
. The spread of mean spin temperatures against corrected column density for all sight-lines with detected absorption is shown in Fig. 11.
(Top) Comparison of column density with density-weighted mean spin temperatures for sight-lines with detected absorption. The column density is calculated using GASKAP emission data and corrected for absorption, see Sec. 5.2. (Bottom) Distribution of the mean spin temperatures against the SMC Hi column density map from GASKAP (Pingel et al., Reference Pingel, Dempsey and McClure-Griffiths2022).
The distribution both by temperature and by corrected column density also varies by region, as shown in Fig. 11. As was done in McClure-Griffiths et al. (Reference McClure-Griffiths, Dénes and Dickey2018), we split the SMC into two regions, the bar running from N to SW and the wing from N to SE (see Fig. 11 bottom). The bar has a large spread in both column density and in mean spin temperature. It contains the highest densities and the highest temperatures of sight-lines in the SMC. However, there is no significant trend in temperature as column density increases in the bar. The wing has a lower column density than much of the bar and both column density and mean spin temperatures are more tightly clustered than the bar. Despite the greatly varying environment, the distribution of spin temperatures in the wing is remarkably flat (
$163 < \langle T_{\mathrm{S,wing}} \rangle < 546\ \mathrm{K}$
) with no strong trend with column density. In the lower column density regions outside the SMC, there is a much wider spread of spin temperatures and again there is no significant trend of temperature with column density.
In two cases (
$\mathrm{J}013134{-}700042$
and
$\mathrm{J}013701{-}730415$
) our spin temperatures are unconstrained. If there is sufficient emission pollution of the absorption spectrum to drive the integral of
$1-e^{-\tau}$
negative, then the calculated mean spin temperature will be negative. As noted in section 4.1, the absorption in
$\mathrm{J}013134{-}700042$
is likely a noise artefact and this noise also results in narrow spikes of emission in the absorption spectrum. For
$\mathrm{J}013701{-}730415$
, while the absorption is a clear detection, there are also more random emission spikes than expected in the spectrum.
5.4. Fraction of Cold Gas
By assuming a temperature of the cool gas (
$T_{\mathrm{c}}$
) we can estimate the fraction of cool gas in each sight-line as follows (Dickey et al., Reference Dickey, Mebold, Stanimirović and Staveley-Smith2000, Eq. 7):
\begin{equation} f_{\mathrm{c}} \equiv \frac{N_{\mathrm{c}}}{N_{w} + N_{c}} \simeq \frac{T_{\mathrm{c}}}{\langle T_{\mathrm{s}} \rangle} ,\end{equation}
where
$N_{\mathrm{c}}$
is the column density of the CNM and
$N_{\mathrm{w}}$
is the WNM column density.
Following Jameson et al. (Reference Jameson, McClure-Griffiths and Liu2019), we assume
$T_{\mathrm{c}} = 30\ \mathrm{K}$
. This gives a median cold gas fraction across the field of
$f_{\mathrm{c}} \approx 11\%$
with the distribution as shown in Fig. 12. The distribution within and outside the body of the SMC is not significantly different. However, this is significantly lower than the
$f_{\mathrm{c}} \approx 20\%$
found by Jameson et al. (Reference Jameson, McClure-Griffiths and Liu2019). Dickey et al. (Reference Dickey, Mebold, Stanimirović and Staveley-Smith2000) reported a median cold gas fraction of 13% for
$T_{\mathrm{c}} = 55\ \mathrm{K}$
, which is equivalent to
$f_{\mathrm{c}} \approx 7\%$
for
$T_{\mathrm{c}} = 30\ \mathrm{K}$
. Our result sits between these two earlier SMC results, similarly to the mean spin temperature results in Sec. 5.3, as would be expected.
The higher mean spin temperatures and lower cold gas fractions than those found in the targeted Jameson et al. (Reference Jameson, McClure-Griffiths and Liu2019) survey are driven by our detection of shallower absorption in the unbiased sight-lines of the present study. We can expect that with the addition of future GASKAP observations of the SMC we will see further shallow absorption detections and thus an even lower median cold gas fraction. The key results of this survey and the earlier SMC surveys are shown in Table 5.
Comparison of the results of SMC absorption surveys
a See Sec. 5.3 for details of uncertainties in
$\langle{T_{s}}\rangle$
Histogram of the distribution of cold gas fractions for sight-lines with detected absorption.
5.5. Comparison with the Milky Way
We compare the results for the field and the body to recent absorption studies of the Milky Way (see Fig. 13, bottom). The field shows a similar, but slightly denser, CDF to the MACH survey (Murray et al., Reference Murray, Stanimirović and Heiles2021). The MACH survey probed a window of very low density Galactic gas at high latitudes. In contrast, the SMC CDF is closest to the Perseus cloud (Lee et al., Reference Lee, Stanimirović, Murray, Heiles and Miller2015) at lower
$\mathcal{R}_{\mathrm{HI}}$
but lacks the higher
$\mathcal{R}_{\mathrm{HI}}$
, likely as our sensitivity is not currently sufficient to probe the highest correction factors that they find. Overall this suggests correction factors similar to the envelopes of Giant Molecular Clouds.
(top) The cumulative distribution function (CDF) with uncertainties of column density correction factor for sight-lines within the body of the SMC and outside the SMC. (bottom) Comparison of column density correction factor CDF for the two samples against (Murray et al., Reference Murray, Stanimirović and Goss2018, “21-SPONGE”), (Stanimirović et al., Reference Stanimirović, Murray, Lee, Heiles and Miller2014, “Perseus”), (Heiles & Troland, Reference Heiles and Troland2003, “HT03”), and (Murray et al., Reference Murray, Stanimirović and Heiles2021, “MACH”).
We also compare the distribution of mean spin temperatures with the same Galactic Hi absorption studies. Fig. 14 shows the mean spin temperatures of all sight-lines with detected Hi absorption, but not saturated, for the SMC body and the field surrounding the SMC, along with comparisons against the recent Galactic Hi absorption samples. Over 70% of temperatures in the body of the SMC are in the
$200 < \langle T_{\mathrm{S}} \rangle < 500\ \mathrm{K}$
range, with a slightly lower proportion in the field. Note that our sensitivity in this study limits our ability to detect high spin temperatures (
$T_{\mathrm{S}} \ge 1150\ \mathrm{K}$
). The two distributions are still remarkably similar, likely indicating that the cold gas we detect in the field is associated with outflows from the SMC. Both the SMC body and the field show similar temperatures to the mid-high latitude (
$|b| > 10^{\circ}$
) HT03 (Heiles & Troland, Reference Heiles and Troland2003) and 21-SPONGE (Murray et al., Reference Murray, Stanimirović and Goss2018) samples, but lacking the higher spin temperatures. These two studies surveyed a variety of environments across the Milky Way, rather than focusing on a single region.
(top) The cumulative distribution function (CDF) with uncertainties of mean spin temperature (
$\langle T_{\mathrm{S}} \rangle$
) for sight-lines within the body of the SMC and outside the SMC where absorption was detected but were not saturated. (bottom) Comparison of mean spin temperature CDF for the two samples against Murray et al. (Reference Murray, Stanimirović and Goss2018, “21-SPONGE”), Stanimirović et al. (Reference Stanimirović, Murray, Lee, Heiles and Miller2014, “Perseus”), Heiles & Troland, (Reference Heiles and Troland2003, “HT03”), and Murray et al. (Reference Murray, Stanimirović and Heiles2021, “MACH”).
5.6. Mean Spectrum of the SMC
The large number of sight-lines sampled through the SMC allow us to build a mean spectrum across the galaxy. This provides a view of the where the cold gas is, and is not, across the SMC. In Figure 15 we show the unweighted mean of spectra for the 43 sight-lines with
$\sigma_{\mathrm{cont}} \le 0.1$
within the body of the SMC. The spectrum has not been corrected for the rotation of the SMC. Notably, there is still fine velocity structure in the mean SMC spectrum.
Weighted mean spectrum for all sight-lines through the SMC. The absorption spectrum is weighted by the inverse square of the continuum noise level for each spectrum. The emission spectrum is unweighted. No correction for the SMC rotation has been applied. See Figure 4 for details.
At the primary emission peak (
$155 \le v_{\mathrm{LSR}} \le 165 \ \mathrm{km\,s}^{-1}$
) 14% of the gas is cold. While at the secondary emission peak (
$120 \le v_{\mathrm{LSR}} \le 130 \ \mathrm{km\,s}^{-1}$
) only 7% of the gas is cold. The absorption falls away much faster from the peak at
$160 \ \mathrm{km\,s}^{-1}$
towards higher velocities, predominantly associated with the wing, than towards lower velocities, which are associated with the bar. Overall, this mean spectrum has a total cold gas fraction of
$f_{\mathrm{c}} \approx 9\%$
and a column density correction factor of
$\mathcal{R}_{\mathrm{HI}} \approx 1.06$
.
5.7. Comparing velocity distributions of absorption and emission
$21\,{\mathrm{cm}}$
emission in the SMC is known to exhibit a complex, multi-peaked velocity structure (Stanimirovic et al., Reference Stanimirovic, Staveley-Smith, Dickey, Sault and Snowden1999). Here we investigate whether Hi absorption is found preferentially in one velocity component or another, or if it follows the same velocity distribution of Hi emission.
To identify the velocities of detected absorption components, we take the first numerical derivative of all spectra (here in the form of
$1-e^{-\tau(v)}$
). As the noise in the spectrum is amplified by the derivative operation, we first smooth each spectrum by a Gaussian kernel of standard deviation of two channels. Detected components are identified by their central channels (
$v_{\mathrm{c}}$
) if they pass the following criteria: (1) the first derivative crosses the zero line; (2) the feature is a maximum (the derivative decreases between adjacent channels); (3) the feature is detected above
$3\sigma$
at
$v_c$
and above
$2.8\sigma$
for adjacent channels; and (4) the feature falls within the SMC velocity range (
$70<v_c<250\,\mathrm{km\,s}^{-1}$
). We repeat the same process for the emission spectra.
Figure 16 displays maps of the positions of detected absorption (a) and emission (b) components, coloured by
$v_c$
. Multiple components per line of sight are displayed as concentric circles.
Comparing the positions of detected absorption (panel a) and emission (panel b) components, coloured by their central velocities (
$v_c$
), identified by their smoothed spectral derivatives, overlaid on an Hi column density map. Lines of sight featuring multiple components are shown as concentric circles.
In Figure 17, we compare the cumulative distribution functions of
$v_c$
for absorption and emission. Uncertainties on the distributions are computed by bootstrapping each sample with replacement over
$10^4$
trials, and represent the
$1^{\mathrm{st}}$
through
$99^{\mathrm{th}}$
percentiles. Although the distributions are indistinguishable between
$100<v_c<150\mathrm{km\,s}^{-1}$
, absorption components are not found in these data beyond
$v_c \sim 175\,\mathrm{km\,s}^{-1}$
whereas emission components extend to
$v_c>200\,\mathrm{km\,s}^{-1}$
. Using a two-sided Kolmogorov-Smirnov test across
$10^4$
bootstrapped trials we find a p-value
$\sim0.17$
. This preliminary result may indicate that the absorption components are more likely to be found in low-velocity gas in the SMC relative to emission components. Future, deeper GASKAP observations of the SMC will provide the opportunity to examine this result further.
Cumulative distribution functions (CDFs) of the central velocities (
$v_c$
) of absorption (blue) and emission (orange) components. Uncertainties on the CDFs are computed by bootstrapping each sample with replacement over
$10^4$
trials, and represent the
$1^{\mathrm{st}}$
through
$99^{\mathrm{th}}$
percentiles.
6. Conclusion
In this paper we have presented the first untargeted survey of Hi absorption in the SMC, alongside the pipeline used to process these and future GASKAP Hi absorption data. We have produced Hi absorption spectra against 337 continuum sources, matched those to emission spectra from the Pingel et al. (Reference Pingel, Dempsey and McClure-Griffiths2022) SMC emission cube, and analysed 229 of those which meet our quality requirements. This represents a 275% increase on sight-lines over previous studies of the SMC region. We have also described the major choices for the pipeline and demonstrated the effectiveness of the pipeline using the SMC data as an example.
We have found 122 absorption features across 65 of the spectra, with no features detected in the other 164 spectra. Within the body of the SMC (
$N_{\mathrm{HI,uncorr}} > 2 \times 10^{21}\ \mathrm{cm}^{-1}$
) we found absorption features against 49 of 79 (62%) spectra, including in all sight-lines with
$\mathrm{S}_{\mathrm{cont}} \ge 100\ \mathrm{mJy}$
. Outside the SMC body we found absorption features against only 16 of 150 (11%) of spectra.
From this first work exploring GASKAP-HI’s unbiased view of absorption across the SMC and its surrounds, we have made two primary findings:
-
The median fraction of Hi in the SMC in the CNM is
$f_{\mathrm{c}} \approx 11\%$
, which is lower than found in Jameson et al. (Reference Jameson, McClure-Griffiths and Liu2019) but higher than the fraction reported in Dickey et al. (Reference Dickey, Mebold, Stanimirović and Staveley-Smith2000). This value is driven by the observation of more shallow absorption in less dense regions. This new fraction better represents the cold gas across the whole of the SMC, rather than just the denser regions observed in Jameson et al. (Reference Jameson, McClure-Griffiths and Liu2019). However, it is potentially still an over estimate due to our sensitivity limits on detection of shallow absorption. Future more sensitive GASKAP-HI observations of the field will refine this figure. -
The range of column density correction factors for optically thick gas varies greatly between the wing and the bar of the SMC. The bar shows a linear increase in correction factor with column density (
$\mathcal{R}_{\mathrm{HI}} \simeq 1 + 0.51[\text{log}_{10}N - 21.43]$
) whilst the wing shows a large variety of correction values with no relationship to uncorrected column density. Overall, the SMC is closer in optical thickness to Galactic molecular clouds, such as the Perseus cloud, whilst the field surrounding the SMC is similar to low density fields at high Galactic latitudes.
In addition we have found that:
-
a) The SMC has an inverse noise-weighted mean spin temperature
$\langle T_{\mathrm{S}} \rangle = 245\ \mathrm{K}$
; which is lower than that found by Dickey et al. (Reference Dickey, Mebold, Stanimirović and Staveley-Smith2000) -
b) There is no significant trend of mean spin temperature with column density
-
c) The equivalent widths and uncorrected column densities we see in the SMC are predominantly within the Milky Way cold gas formation thresholds proposed by Kanekar et al. (Reference Kanekar, Braun and Roy2011), however we see a higher limit for the condensation of dense Hi into molecular
$\mathrm{H}_2$
saturation in the SMC, due to the lower metallicity of the SMC, and we see hints of cold gas formation below their threshold -
d) Fine velocity structure is present even in the mean SMC spectrum, which implies that this structure may also be seen in extra-Galactic experiments
-
e) There are indications that absorption components may be preferentially found in the lower velocity gas of the SMC as compared to the emission components.
Further exploration of the pilot SMC absorption data is planned in future GASKAP-HI papers. In addition, a much longer observation (up to 200 hrs) of the SMC in the full GASKAP-HI survey will obtain Hi absorption spectra up to 3 times as sensitive as the spectra presented here. We will have the opportunity to explore the column density threshold for the formation of cold Hi proposed by Kanekar et al. (Reference Kanekar, Braun and Roy2011) in more detail. We also expect that we will see further clear detections of shallow absorption features, and thus slightly refine the cold gas fraction for the SMC.
Acknowledgement
The Parkes radio telescope is part of the Australia Telescope National Facility which is funded by the Australian Government for operation as a National Facility managed by CSIRO. The Australian SKA Pathfinder is part of the Australia Telescope National Facility which is managed by CSIRO. Operation of ASKAP is funded by the Australian Government with support from the National Collaborative Research Infrastructure Strategy. ASKAP uses the resources of the Pawsey Supercomputing Centre. Establishment of ASKAP, the Murchison Radio-astronomy Observatory and the Pawsey Supercomputing Centre are initiatives of the Australian Government, with support from the Government of Western Australia and the Science and Industry Endowment Fund. We acknowledge the Wajarri Yamatji people as the traditional owners of the Observatory site. This paper includes archived data obtained through the CSIRO ASKAP Science Data Archive, CASDA (http://data.csiro.au).
This research was supported by the Australian Research Council (ARC) through grant DP190101571. N.M.-G. acknowledges the support of the ARC through Future Fellowship FT150100024. S.S. gratefully acknowledges support by the NSF grant AST-2108370. We thank the anonymous reviewer for their thoughtful comments that have improved this paper.
Open access
