Abstract
We investigate the kinematics of neutral gas in the Small Magellanic Cloud (SMC) and test the hypothesis that it is rotating in a disk. To trace the 3D motions of the neutral gas distribution, we identify a sample of young, massive stars embedded within it. These are stars with radial velocity measurements from spectroscopic surveys and proper motion measurements from Gaia, whose radial velocities match with dominant H i components. We compare the observed radial and tangential velocities of these stars with predictions from the state-of-the-art rotating disk model based on high-resolution 21 cm observations of the SMC from the Australian Square Kilometer Array Pathfinder telescope. We find that the observed kinematics of gas-tracing stars are inconsistent with disk rotation. We conclude that the kinematics of gas in the SMC are more complex than can be inferred from the integrated radial velocity field. As a result of violent tidal interactions with the Large Magellanic Cloud, nonrotational motions are prevalent throughout the SMC, and it is likely composed of distinct substructures overlapping along the line of sight.
Export citation and abstract BibTeX RIS
1. Introduction
The Small Magellanic Cloud (SMC) is a gas-rich, low-metallicity (Z ∼ 0.2 Z⊙; Russell & Dopita 1992) dwarf irregular satellite of the Milky Way (MW). At a distance of ∼62 kpc (e.g., Scowcroft et al. 2016), the SMC provides an ideal laboratory for studying the structure and kinematics of dwarf galaxies at high sensitivity and resolution.
As a result of its proximity to the MW and the nearby Large Magellanic Cloud (LMC), the SMC has been strongly influenced by ongoing dynamical interactions. The SMC and LMC likely experienced a direct collision ∼150 Myr ago (e.g., Zivick et al. 2018) and their interactions have produced the broader, gaseous Magellanic System, including the Bridge (Hindman et al. 1963), Stream (Mathewson et al. 1974), and Leading Arm (Putman et al. 1998). In fact, low metallicities observed in the Stream and Leading Arm suggest that these structures originated, in part, from the SMC (Fox et al. 2013; D'Onghia & Fox 2016), and this premise is consistent with most models of the Magellanic System (Gardiner & Noguchi 1996; Connors et al. 2006; Besla et al. 2012; Diaz & Bekki 2012).
The complex dynamical history of the SMC has made it challenging to understand its present-day structure and kinematics. The neutral hydrogen (H i) in the system is highly disturbed (e.g., Stanimirović et al. 1999), featuring multi-peaked velocity profiles indicative of supergiant shells (Hindman 1967), two separate velocity subsystems (e.g., Mathewson et al. 1988), and/or hundreds of expanding shells (Staveley-Smith et al. 1997). However, despite this elaborate structure, the integrated velocity field of the SMC exhibits a strong gradient suggestive of a nicely rotating disk (Kerr et al. 1954; Hindman et al. 1963; Stanimirović et al. 2004; Di Teodoro et al. 2019).
To complicate the picture further, observations of distinct stellar populations provide additional interpretations for the structure of the SMC. The oldest stars are spherically distributed within a radius of 7–12 kpc (Subramanian & Subramaniam 2012), with little evidence for rotation from radial velocities or proper motions (PMs; Harris & Zaritsky 2006; Kallivayalil et al. 2013; Gaia Collaboration et al. 2018; Niederhofer et al. 2018; Zivick et al. 2018). Red clump stars and Cepheid variables are highly elongated along the line of sight (LOS; ∼20 kpc depth; Mathewson et al. 1988; Nidever et al. 2013; Jacyszyn-Dobrzeniecka et al. 2016; Scowcroft et al. 2016; Ripepi et al. 2017). Furthermore, although younger stellar populations and red giant branch (RGB) stars display a velocity gradient, it appears perpendicular to the gradient seen in H i (Evans & Howarth 2008; Dobbie et al. 2014). Similarly, the distance gradient observed in Cepheids from the northeast to southwest is ∼90° from the minor axis, as implied by a rotating gas disk (Scowcroft et al. 2016).
Ultimately, although the observed distribution and kinematics of distinct mass tracers in the SMC have resulted in a diversity of interpretations for its structure, the rotating disk model derived from the H i velocity field remains a fundamental benchmark for our theoretical understanding of the SMC and the broader Magellanic System. Basic properties including its center of mass, halo composition, and total dynamical mass—fundamental for numerical models—are derived from the rotation curve analysis (Stanimirović et al. 2004; Di Teodoro et al. 2019). Furthermore, the SMC LOS depth inferred from the H i velocity field is routinely invoked to estimate the surface densities of molecular gas and star formation in the system. For example, although the star formation efficiency of molecular gas is the same in the SMC as within large spiral galaxies, the ratio of molecular to atomic gas is an order of magnitude lower (Jameson et al. 2016) implying that the SMC is "strikingly bad" at converting its atomic gas reservoir into stars (Bolatto et al. 2011). Although the interaction history between the SMC and the LMC and MW may have affected this capability by removing gas from the system, this result is invoked to explain similarly suppressed star formation rate efficiencies of H i-dominated galaxies at z ∼ 1–3 (e.g., Wolfe & Chen 2006; Rafelski et al. 2016). Although the result can be alleviated by varying uncertain disk parameters, such as its inclination, a preferable model remains elusive.
Ultimately, reconciling the disparate kinematics of gas and stars in the SMC and testing the rotating disk hypothesis requires constraints for the 3D kinematics of gas in the system. Previous models were based on the the radial velocity alone, as this is the only observable velocity component for the gas. In this paper, we constrain the 3D motions of gas within the SMC using the measured 3D motions of embedded young stars (O and B type). In Section 2, we present the pertinent observations and models; in Section 3, we discuss the methods used to analyze the observations; in Section 4, we present the results of comparing the observed kinematics with predictions from the rotating disk model; and in Section 5, we discuss the implications for the SMC system.
2. Data
Although the radial motions of the SMC gas distribution are well measured by radio observations of neutral gas tracers (e.g., H i), the transverse motions are unconstrained. We assume that massive stars whose radial velocities match with the gaseous components formed recently enough (∼1–10 s Myr lifetimes) to trace the 3D motions of the same gas that they are embedded within, and we use the measured transverse velocities of these stars to trace the transverse motions of the underlying gas. To avoid confusion with stars that have migrated away from their original gas clouds (Oey et al. 2018), we make a very conservative selection to generate a small sample of 143 stars whose velocities match closely with clearly defined peaks in the 21 cm spectra in their directions. This velocity matching step rejects some 88% of the original sample of stars, but the remaining stars are the most likely to show the same disk rotation that the gas shows.
2.1. Neutral Gas
As the SMC is rich in neutral atomic gas (e.g., Leroy et al. 2007), we use H i to trace the gas distribution. New observations of the SMC at 21 cm were recently obtained by the Australian Square Kilometer Array Pathfinder (ASKAP; DeBoer et al. 2009) during Commissioning and Early Science observations in 2017 November. These data were combined with single-dish observations of the SMC from the Parkes Telescope as part of the HI4PI survey (HI4PI Collaboration et al. 2016). For details concerning the observations and data reduction for the ASKAP+Parkes cube, we refer the reader to McClure-Griffiths et al. (2018). The final data cube we use has an rms noise in 21 cm brightness temperature (TB) of , a pixel size of 7'' (35'' × 27'' in angular resolution), and velocity channel width of
(McClure-Griffiths et al. 2018; Di Teodoro et al. 2019). For comparison with stellar velocities, we convert the H i velocities from the local standard of rest to the heliocentric reference frame. Hereafter, all quoted velocities will be in the heliocentric reference frame.
In Figure 1, we display a map of the H i column density (N(H i)) of the SMC from the ASKAP+Parkes data cube (panel (a)), the first moment map (panel (b)), and a position–velocity slice along the major axis inferred from kinematic modeling of the H i data cube by Di Teodoro et al. (2019); panel (c), see Section 2.4). The moment map is masked as in Di Teodoro et al. (2019), wherein a flood-fill algorithm applied to a smoothed version of the data cube (smoothed with a Gaussian kernel of FWHM 70'') identified pixels with emission brighter than ∼10
Figure 1. (a) H i column density map of the SMC from the ASKAP+Parkes data cube. The dotted line delineates the main "Bar" from the "Wing" of the SMC, and arrows indicate the directions of the Magellanic Stream and the LMC (following McClure-Griffiths et al. 2018). (b) The first moment map, masked as described in Section 2.1. The contours (gray) display N(H i) at levels of . The dotted line (black) indicates the major axis inferred from kinematic modeling of the ASKAP+Parkes data cube (center:
237,
273, position angle: 66°; Di Teodoro et al. 2019). (c) The position–velocity slice through the center of the SMC and along the major axis (shown in panel (b)). Contours are drawn at TB = 3 × 2n K.
Download figure:
Standard image High-resolution image2.2. Stars
We select a subsample of SMC stars with measured radial velocities from the spectroscopic survey of the SMC by the 2-degree-field instrument at the Anglo-Australian Telescope (2045 stars; Evans & Howarth 2008) and from the Runaways and Isolated O-Type Star Spectroscopic Survey of the SMC (RIOTS4), a uniformly selected survey of young field stars in the SMC using the Inamori-Magellan Areal Camera and Spectrograph (IMACS) and Magellan Inamori Kyocera Echelle (MIKE) on the Magellan telescopes (374 stars; Lamb et al. 2016). Of this sample, 1484 stars are of spectral type O or B and have significant radial velocity (vr) measurements (defined as , where
is the uncertainty on the radial velocity measurement).
To determine the transverse velocities of the selected stars, we cross-match them with SMC stars from the second data release (DR2) of the Gaia mission (Gaia Collaboration et al. 2018). The targets observed by Evans & Howarth (2008) and Lamb et al. (2016) were selected from the original catalog by Massey (2002), who reported a positional uncertainty of 1'' per target. We begin with a conservative search radius of 3 and
where
and
are the uncertainties in R.A. and decl. proper motions, respectively). The matched stars represent the youngest, brightest end of the population observed by Gaia; for example, the median Gaia G magnitude for the matched sample is G = 15.3, compared with the median magnitude for the full SMC sample of G = 19. We further refine the sample by comparing the observed flux densities and remove 249 objects with discrepant G − B colors (defined as being further from the mean G − B color of 0.16 by more than ±1
2.3. Systemic Proper Motion
To estimate the PMs due to internal motions of the SMC, we subtract the systemic motion of the galaxy from the measured Gaia PMs. For consistent comparison with the latest dynamical modeling of the system (see Section 2.4), we adopt the measured center of mass (COM) PM derived from Hubble Space Telescope photometry by Kallivayalil et al. (2013), who estimated and
. We note that many studies have estimated the COM PM of the SMC (e.g., Piatek et al. 2008; Kallivayalil et al. 2013; van der Marel & Sahlmann 2016; Niederhofer et al. 2018; Zivick et al. 2018), and all estimates are largely consistent within uncertainties.
2.4. Dynamical Model
For comparison with the observed kinematics, we use the results of a recent dynamical model in the SMC by (Di Teodoro et al. 2019, hereafter DT19) based on the ASKAP+Parkes H i data cube. DT19 modeled the SMC as a disk undergoing circular rotation. Using a Markov Chain Monte Carlo (MCMC) approach, they fitted the observed velocity field (e.g., Figure 1(b)) to derive global properties of the system, including kinematic center, position and inclination angles, transverse and systemic velocities, and then applied a tilted-ring analysis to decompose the velocity field and extract the rotation velocity as a function of radius.
Using the DT19 results, we predict the velocities for stars within a rotating SMC based on their coordinates. To compute each velocity component, we follow van der Marel & Cioni (2001, hereafter vdM2001) and van der Marel et al. (2002, hereafter vdM2002). We assume the same global SMC parameters as derived in the DT19 model, including: the SMC COM (), distance (
), COM PM (
; Kallivayalil et al. 2013), systemic LOS velocity (Heliocentric,
), inclination (i = 55° ± 9°), position angle of the line of nodes (PA = 66° ± 8°; measured counterclockwise from north) and precession and nutation rate (
). From the assumed COM PM, we also compute the transverse velocity of the COM and its angle on the sky.
Given these global parameters, for each star, we compute its angular coordinates on the sky (Equations (1)–(3), vdM2001), its radius in the galaxy frame (Equation (9), vdM2002) and its ensuing predicted rotation velocity from the DT19 model. We then follow vdM2002 to compute the individual velocity contributions from the SMC COM (Equation (13), vdM2002), its precession and nutation (Equation (16), vdM2002) and internal motions (Equation (21), vdM2002). The sum of these components comprise the predicted velocities in the frame of the SMC. Finally, we extract the predicted radial velocity () and the PMs to the west and north (
; Equation (9), vdM2002) for comparison with the observed sample. As an estimate of the uncertainties in these predictions, we propagate the quoted uncertainties in the nine global SMC parameters from DT19 through the analysis.
3. Analysis
3.1. Line-of-sight Complexity
Although the radial velocity structure of H i in the SMC is complex, it has been shown that a significant fraction of OB stars have coincident radial velocities with H i structures along the same LOS and are, therefore, embedded within the H i distribution (Lamb et al. 2016). Although the majority of young stars are associated with some H i emission, our goal is to construct a pure sample of targets that are most likely to trace the motions of dominant H i components.
We begin by identifying significant H i components along each LOS in the ASKAP+Parkes data cube via the derivative of . First, we resample each
spectrum to
resolution and smooth with a Gaussian kernel eight channels in width to suppress the influence of noise on selecting spurious peaks. We then compute the first numerical derivative of each spectrum and identify components as local maxima with significant emission (
). Figure 2 displays the total number of significant H i components along each LOS. Although the outskirts of the SMC typically feature a single strong H i component, within the main body of the SMC the radial velocity structure of the gas is more complex and there are at least two H i components. This is evident in Figure 1(c); in addition to the "main disk" of the SMC along the major axis inferred by rotation models, there are several significant velocity components evident in the position–velocity diagram.
Figure 2. Map of the total number of significant H i components (Ncomp), masked as described in Section 2.1. White contours denote N(H i) at levels of .
Download figure:
Standard image High-resolution image3.2. Matching Radial Velocities
To identify which stars are moving with the identified H i components, we compare the observed stellar radial velocities () with the radial velocities of the H i components (
). To select the most likely sources to trace the dominant H i along each LOS, we design a method for matching stars with H i components which accounts for the uncertainties in
, which can be large (
; Evans & Howarth 2008; Lamb et al. 2016).
For each star, we select N = 1000 velocities distributed normally with mean and
. For each H i component, we construct similar normally distributed velocity distributions with mean
and an assumed
. We then use a two-sample Anderson–Darling test (A–D; Anderson & Darling 1954) to test for the null hypothesis that the stellar and H i velocity distributions are drawn from the same parent population. A match is achieved if the null hypothesis cannot be rejected with high confidence (>99%), and if the brightness temperature of the H i peak is above a certain fraction of the max TB along the LOS (rpeak, where
), to ensure that matches are with dominant H i peaks.
To select the optimal values for the matching parameters () and select the values that maximize the number of matches while minimizing the contamination fraction. To quantify the contamination of the sample due to random interlopers with coincidental radial velocities to the H i peaks, we repeat the matching process using simulated samples of 1202 stars with random velocities and uncertainties drawn from the observed distributions over 1000 trials. The estimated number of spurious matches is computed as the median over all trials. We find that the choice of
and
.
We note that this approach produces statistically indistinguishable results from other methods of matching and vH i, including: (1) by-eye matching, (2) applying a range of absolute cuts to
, and (3) computing the area of overlap between the
and
distributions and applying a sigma cut.
To determine the final sample, we repeat the matching process with the selected match criteria over an additional 1000 trials and select all stars that have successful matches in more than 95% of the trials. The result is a sample of 143 stars. The matched sources are shown in Figure 3. In the Appendix, we include similar figures for all LOS (Figure 6).
Figure 3. Results of matching observed radial velocities of young, massive stars (red dashed) with the radial velocities of distinct H i peaks along the same LOS (black). For each of the 1202 stars with radial and transverse velocity measurements, over 1000 trials, we sample random velocities within (gray) and find the 143 LOS that "match" (see Section 3.2). The ±1
are illustrated (gray shading), and the first moment (Figure 1(b)) is included (blue dotted line).
Download figure:
Standard image High-resolution imageAlthough the final sample represents only ∼12% of the initial sample of stars with both radial and transverse velocity constraints (1202 sources), this does not mean that the remaining 88% of the sample are not associated with H i gas. In fact, we observe that the majority of the initial sample of 1202 stars (∼80%) is associated with significant H i emission, which is defined as there being significant H i emission at the radial velocity of the star (). This indicates that most stars are embedded somewhere within the neutral gas distribution, if not associated with the strongest peaks. The final sample of 143 matched stars represents our highest confidence for tracing the H i kinematics with minimal contamination. We further observe that although relaxing the matching criteria (e.g., including wider uncertainty ranges, selecting different brightness temperature cuts to define H i components, etc.) changes the number and selection of stars that are included in the matched sample, it does not significantly affect the statistics of the sample (e.g., estimated contamination fraction) or the ensuing conclusions.
4. Results
In Figure 4(a), we display the residual PMs (i.e., with the COM PM from Kallivayalil et al. 2013 subtracted) for the 143 stars whose radial velocities match with significant H i peaks. The stars exhibit a radial velocity gradient in agreement with the H i radial velocity field (Figure 1(b)). However, the observed residual PMs do not exhibit a clear signature of rotation. For comparison, in Figure 4(b), we display the radial velocities and residual PMs predicted by the DT19 rotation model. As expected, the modeled residual PMs exhibit a clear signature of rotation.
Figure 4. (a) Radial velocities (color) and residual PMs (vectors) for the sample of 143 stars with significantly detected radial and transverse motions from spectroscopy (Evans & Howarth 2008; Lamb et al. 2016) and Gaia (Gaia Collaboration et al. 2018), whose radial velocities match significant H i components along the LOS (see Figure 3). All vectors are constructed to illustrate PMs over the next 7.2 Myr (arbitrarily, as in van der Marel & Sahlmann 2016). The systemic PM of the SMC (Kallivayalil et al. 2013) is included at top left (dotted), as is a scale bar to illustrate 100 km s−1 (dashed). (b) Same, for the predicted radial velocities and residual PMs (black vectors) from the SMC rotation model by DT19.
Download figure:
Standard image High-resolution imageComparing the panels of Figure 4, although the observed and modeled kinematics are not consistent throughout the SMC, the level of agreement may depend on location within the system. For example, a subset of stars have large observed tangential velocities (). These large velocities are likely dominated by dynamical effects (Oey et al. 2018). Stars in the the Wing region (southeast, see Figure 1(a)) are moving toward the LMC via the Magellanic Bridge as a result of their recent interaction history (Zivick et al. 2019). Furthermore, this net PM of the Wing relative to the Bar indicates that these two regions may be kinematically distinct (Oey et al. 2018). Although stars and gas throughout the SMC are affected at some level by the same dynamical interactions, in subsequent analysis, we will distinguish stars in the Wing and Bar to account for the possibility that stars in different regions are tracing distinct populations.
In Figure 5, we compare the observed and predicted radial (LOS) velocity (a), and the residual PMs to the west and north ((b) and (c), respectively) for the 143 matched stars. Stars located in Wing and Bar are distinguished from each other for clarity (see Figure 1(a) for definitions).
Figure 5. Comparison between radial velocity (vr) and PMs to the west and north ( are 1
and
are 1
Download figure:
Standard image High-resolution imageFigures 4 and 5 indicate that the observed kinematics are largely inconsistent with the predictions from the rotating disk model. To quantify the agreement, we use the K-sample A–D statistic. First, we bootstrap the sample with replacement over 1000 trials, and within each trial, we use a simple Monte Carlo exercise to draw random values of each velocity component within ±3
5. Discussion
Overall, the observed 3D kinematics of H i in the SMC traced by young massive stars are inconsistent with predictions from the rotating disk model. In the following, we discuss how to reconcile this result within the context of the SMC.
First, it is important to consider alternatives to the hypothesis that the SMC is well represented by a single rotating disk. As discussed by DT19, there are several caveats to the disk model, including the fact that the disk is not razor thin (as must be assumed to simplify the velocity component computation), the rotation is not necessarily perfectly circular, and the disk is not necessarily axisymmetric. In addition to these caveats, there is clear evidence that significant components of the system (e.g., the Wing) are dominated by radial motions due to the influence of the LMC (e.g., Zivick et al. 2019). The fact that the rotation model performs as well as it does at reproducing the radial velocity field of the gas suggests that the large scale motion is dominated by rotation and the discrepant stellar kinematics observed here may reflect only local, noncircular motions. Although a subset of our sample likely does trace local motions (e.g., effects of stellar feedback on small scales), the majority were selected to trace the dominant H i components. As seen in Figure 5(a), the observed radial velocities of stars matched to H i peaks trace distinct, coherent components in radial velocity—one at and one at
—which are found in both the Wing and the Bar. The intensity-weighted mean velocity along each LOS used to predict 3D gas motions can be insensitive to this multiplicity, especially if there are multiple components with similar intensities. To illustrate this effect, in Figure 3, we include the first moment for comparison with the H i peak velocities along each LOS. In addition, we note that although Gaia PMs in the SMC display a marginal signal of rotation (Gaia Collaboration et al. 2018, their Figure 16), it is inconsistent with the rotation inferred from the H i.
An alternative explanation for the discrepancy between the observed and model kinematics is that the SMC is not a single system but is rather composed of multiple structures along the LOS. As illustrated in Figures 2 and 5, the structure of the SMC is complex and exhibits multiple, significant H i components separated in radial velocity. This velocity structure has previously been attributed to separate subsystems within the SMC (e.g., the Mini Magellanic Cloud and the SMC Remnant; Mathewson et al. 1984), which may be separated in distance along the LOS by (Mathewson et al. 1986). Based on measurements of ultraviolet absorption lines toward stars in the SMC, the lower-velocity H i component appears to be sitting in front of the higher-velocity H i component along the LOS (Songaila et al. 1986; Wayte 1990; Danforth et al. 2002; Welty et al. 2012). This implies that the SMC is being ripped apart by its interactions with the LMC, which is in agreement with tidal models of the system (Murai & Fujimoto 1980), resulting in disparate remnant structures.
If the SMC is composed of separate subsystems at different distances, it is possible that the apparent rotation seen in the radial velocity field of the gas is dominated by rotation within one component over the other. However, exactly how we can reconcile potential rotation in this component alone with the nonrotating kinematics of the rest of the SMC is uncertain and is not incorporated into current theoretical models. In addition, it has been shown that models with purely radial motions can reproduce observed velocity fields just as well as disk models with purely rotational motion (e.g., Sylos Labini et al. 2019). Furthermore, the derived global properties of the system, including its dynamical center and mass—crucial input parameters to numerical models of the broader Magellanic System (Besla et al. 2012)—are based on fits to the global radial velocity field (e.g., including all velocities along each LOS; Figure 1). If only subsets of the SMC are actually rotating, these models (and the resulting parameters) would still need to be revisited.
6. Summary and Conclusions
In this paper, we present new constraints for the 3D kinematics of gas in the SMC. We trace the transverse motions of H i in the system for the first time using the observed proper motions of massive, young stars whose radial velocities match with significant H i components identified from new observations by the ASKAP telescope (McClure-Griffiths et al. 2018). We test the hypothesis that gas in the SMC is rotating in a disk by comparing the observed radial and transverse velocities of young stars with predictions from the latest model of the radial velocity field from DT19. Our main conclusions are the following:
- 1.
- 2.The 3D kinematics of the SMC are more complex than can be accurately inferred from the integrated radial velocity field alone.
- 3.Due to violent tidal and ram-pressure interactions with the LMC and the MW halo, it is likely that the SMC is composed of overlapping substructures, whose properties should be considered in future models of the system for extracting basic parameters, such as the COM, dynamical mass, and depth along the LOS.
Overall, this study represents one step in the process of unraveling the structure and kinematics of the SMC. In the near future, with the third data release from Gaia, we will have access to radial velocity measurements toward a new sample of stars throughout the SMC with which we will further dissect the kinematics of individual gaseous components. In addition, with resolved stellar photometry, we will constrain the distances to individual dusty components and complete the full 6D phase space of the system (3D positions and 3D velocities). As the Magellanic Clouds are the only systems for which we can achieve this level of resolution, this will not only inform our understanding of these systems but also our understanding of dwarf galaxy interactions throughout cosmic time.
We would like to thank the anonymous referee for thoughtful comments that have improved this manuscript. C.E.M. is supported by an NSF Astronomy and Astrophysics Postdoctoral Fellowship under award AST-1801471. N.Mc.-G. acknowledges funding from the Australian Research Council via grant FT150100024. This work took part under the program Milky Way-Gaia of the PSI2 project funded by the IDEX Paris-Saclay, ANR-11-IDEX-0003-02. This research has made use of NASA's Astrophysics Data System. 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.
Software: Astropy (Astropy Collaboration et al. 2013), NumPy (Van Der Walt et al. 2011), matplotlib (Hunter 2007), glue (Beaumont et al. 2015).
Appendix: Matching Radial Velocities
In this appendix, we provide supplemental figures and tables for illustrating the match process between the observed radial velocities of young, massive stars (Evans & Howarth 2008; Lamb et al. 2016) and the observed radial velocities of significant H i peaks along the same LOS. In Figure 6, we display plots for each of 1202 stars with significant radial velocities with corresponding proper motion measurements from Gaia DR2 (Gaia Collaboration et al. 2018). We include the H i spectrum () and the identified significant components along each LOS. The 143 matched LOS are identified by purple plot frames (and are also shown in Figure 3).
Figure 6.
Each panel compares the 21 cm spectrum and identified H i components (
; dotted green lines) with the observed stellar radial velocity (
; red dashed line). The uncertainties in vr,star are illustrated by 50 random draws from a normal distribution with
and
(gray). The 143 matches (see the text) are highlighted by bold, purple plot frames (also see Figure 3). If
is not visible within the panel, it falls outside the plotted velocity range. (An extended version of this figure is available.)
Download figure:
Standard image High-resolution imageIn Table 1, we include the coordinates, velocities, and PMs for the 143 H i-matched stars.
Table 1. Parameters for the 143 Stars Used in the Kinematic Comparison
R.A. | Decl. |
![]() |
![]() |
Type | References | ||||
---|---|---|---|---|---|---|---|---|---|
(deg) | (km s−1) | (mas yr−1) | |||||||
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
7.798 | −73.993 | 116.0 | 16.0 | B1-3 (V) | 2 | −0.656 | −1.550 | 0.127 | 0.107 |
8.410 | −73.973 | 160.0 | 9.0 | B1-3 (IV) | 2 | −0.772 | −1.580 | 0.125 | 0.092 |
8.412 | −74.064 | 160.0 | 9.0 | B9 (Ib) | 2 | −0.501 | −1.274 | 0.067 | 0.050 |
8.565 | −73.694 | 124.0 | 9.0 | B3-5 (III) | 2 | −0.834 | −1.402 | 0.136 | 0.107 |
8.793 | −74.122 | 160.0 | 9.0 | B9 (II) | 2 | −0.450 | −1.190 | 0.082 | 0.069 |
9.178 | −73.137 | 116.0 | 7.0 | B8 (II) | 2 | −0.628 | −1.088 | 0.056 | 0.049 |
9.178 | −73.794 | 105.0 | 6.0 | B8 (II) | 2 | −0.782 | −1.291 | 0.097 | 0.069 |
9.504 | −73.652 | 108.0 | 10.0 | B9 (Ib) | 2 | −0.593 | −1.219 | 0.047 | 0.043 |
9.536 | −73.516 | 113.0 | 9.0 | B9 (II) | 2 | −0.712 | −1.269 | 0.055 | 0.051 |
9.641 | −73.579 | 173.0 | 7.0 | B1-3 (II) | 2 | −0.418 | −1.128 | 0.068 | 0.056 |
Note. (1), (2): R.A., decl. coordinates; (3), (4): observed radial velocity and uncertainty; (5) stellar type; (6) radial velocity reference: 1 = Lamb et al. (2016) and 2 = Evans & Howarth (2008); (7)–(10): PMs and uncertainties in the west; (7), (9) north; and (8), (10) directions from Gaia Collaboration et al. (2018).
Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.
Download table as: DataTypeset image