Massive Star-forming Regions Research Articles

Anatomy of the massive star-forming region S106

The central area (40″ × 40″) of the bipolar nebula S106 was mapped in the [O I] line at 63.2 μm (4.74 THz) with high angular (6″) and spectral (0.24 MHz) resolution, using the GREAT heterodyne receiver on board SOFIA. The spatial and spectral emission distribution of [O I] is compared to emission in the CO 16 →15, [C II] 158 μm, and CO 11 →10 lines, mm-molecular lines, and continuum. The [O I] emission is composed of several velocity components in the range from –30 to 25 km s−1. The high-velocity blue- and red-shifted emission (v = −30 to –9 km s−1 and 8 to 25 km s−1) can be explained as arising from accelerated photodissociated gas associated with a dark lane close to the massive binary system S106 IR, and from shocks caused by the stellar wind and/or a disk–envelope interaction. At velocities from –9 to –4 km s−1 and from 0.5 to 8 km s−1 line wings are observed in most of the lines that we attribute to cooling in photodissociation regions (PDRs) created by the ionizing radiation impinging on the cavity walls. The velocity range from –4 to 0.5 km s−1 is dominated by emission from the clumpy molecular cloud, and the [O I], [C II], and high-J CO lines are excited in PDRs on clump surfaces that are illuminated by the central stars. Modelling the line emission in the different velocity ranges with the KOSMA-τ code constrains a radiation field χ of a few times 104 and densities n of a few times 104 cm−3. Considering self-absorption of the [O I] line results in higher densities (up to 106 cm−3) only for the gas component seen at high blue- and red velocities. We thus confirm the scenario found in other studies that the emission of these lines can be explained by a two-phase PDR, but attribute the high-density gas to the high-velocity component only. The dark lane has a mass of ~275 M⊙ and shows a velocity difference of ~1.4 km s−1 along its projected length of ~1 pc, determined from H13CO+ 1 →0 mapping. Its nature depends on the geometry and can be interpreted as a massive accretion flow (infall rate of ~2.5 × 10−4 M⊙ yr−1), or the remains of it, linked to S106 IR/FIR. The most likely explanation is that the binary system is at a stage of its evolution where gas accretion is counteracted by the stellar winds and radiation, leading to the very complex observed spatial and kinematic emission distribution of the various tracers.

From Giant H ii regions and H ii galaxies to globular clusters and compact dwarf ellipticals

Massive starforming regions like Giant HII Regions (GHIIR) and HII Galaxies (HIIG) are emission line systems ionized by compact young massive star clusters (YMC) with masses ranging from $10^4$M$_\odot$ to $10^8$M$_\odot$. We model the photometric and dynamical evolution over a Hubble time of the massive gravitationally bound systems that populate the tight relation between absolute blue magnitude and velocity dispersion ($M_{B}-\sigma$) of GHIIR and HIIG and compare the resulting relation with that one of old stellar systems: globular clusters, elliptical galaxies, bulges of spirals. After 12~Gyr of evolution their position on the $\sigma$ vs. M$_B$ plane coincides -- depending on the initial mass -- either with the globular clusters for systems with initial mass $M < 10^6$M$_\odot$ or with a continuation of the ellipticals, bulges of spirals and ultracompact dwarfs for YMC with $M >10^6$M$_\odot$. The slope change in the $M_{B}-\sigma$ and $M_B$-size relations at cluster masses around $10^6$M$_\odot$ is due to the larger impact of the dynamical evolution on the lower mass clusters. We interpret our result as an indication that the YMC that ionize GHIIR and HIIG can evolve to form globular clusters and ultra compact dwarf ellipticals in about 12 Gyr so that present day globular clusters and ultra compact dwarf ellipticals may have formed in conditions similar to those observed in today GHIIR and HIIG.

A Submillimeter Burst of S255IR SMA1: The Rise and Fall of Its Luminosity

Abstract Temporal photometric variations at near-infrared to submillimeter wavelengths have been found in low-mass young stellar objects. These phenomena are generally interpreted as accretion events of star-disk systems with varying accretion rates. There is growing evidence suggesting that similar luminosity flaring also occurs in high-mass star/cluster-forming regions. We report in this Letter the rise and fall of the 900 μm continuum emission and the newly found 349.1 GHz methanol maser emission in the massive star-forming region S255IR SMA1 observed with the Submillimeter Array and the Atacama Large Millimeter/submillimeter Array. The level of flux variation at a factor of ∼2 at the submillimeter band and the relatively short 2 years duration of this burst suggest that the event is probably similar to those milder and more frequent minor bursts seen in 3D numerical simulations.

Flux Density Variations at 3.6 cm in the Massive Star-forming Region W49A

Abstract A number of ultracompact H ii regions in Galactic star-forming environments have been observed to vary significantly in radio flux density on timescales of 10–20 years. Theory predicted that such variations should occur when the accretion flow that feeds a young massive star becomes unstable and clumpy. We have targeted the massive star-forming region W49A with the Karl G. Jansky Very Large Array for observations at 3.6 cm with the B-configuration at ∼0.″8 resolution, to compare to nearly identical observations taken almost 21 years earlier (2015 February and 1994 August). Most of the sources in the crowded field of ultracompact and hypercompact H ii regions exhibit no significant changes over this time period. However, one source, W49A/G2, decreased by 20% in peak intensity (from 71 ± 4 mJy/beam to 57 ± 3 mJy/beam), and 40% in integrated flux (from 0.109 ± 0.011 Jy to 0.067 ± 0.007 Jy), where we cite 5σ errors in peak intensity, and 10% errors in integrated flux. We present the radio images of the W49A region at the two epochs, the difference image that indicates the location of the flux density decrease, and discuss explanations for the flux density decrease near the position of W49A/G2.

The dense galactic environments of the Milky Way

AbstractStar formation takes place in the dense gas phase, and therefore a simple dense gas and star formation rate relation has been proposed. With the advent of multi-beam receivers, new observations show that the deviation from linear relations is possible. In addition, different dense gas tracers might also change significantly the measurement of dense gas mass and subsequently the relation between star formation rate and dense gas mass. We report the preliminary results the DEnse GAs in MAssive star-forming regions in the Milky Way (DEGAMA) survey that observed the dense gas toward a suite of well-characterized massive star-forming regions in the Milky Way. Using the resulting maps of HCO+ 1–0, HCN 1–0, CS 2–1, we discuss the current understanding of the dense gas phase where star formation takes place.

Revealing magnetic fields towards massive protostars: a multi-scale approach using masers and dust

AbstractMagnetic fields play a significant role during star formation processes, hindering the fragmentation and the collapse of the parental cloud, and affecting the accretion mechanisms and feedback phenomena. However, several questions still need to be addressed to clarify the importance of magnetic fields at the onset of high-mass star formation, such as how strong they are and at what evolutionary stage and spatial scales their action becomes relevant. Furthermore, the magnetic field parameters are still poorly constrained especially at small scales, i.e. few astronomical units from the central object, where the accretion disc and the base of the outflow are located. Thus we need to probe magnetic fields at different scales, at different evolutionary steps and possibly with different tracers. We show that the magnetic field morphology around high-mass protostars can be successfully traced at different scales by observing maser and dust polarised emission. A confirmation that they are effective tools is indeed provided by our recent results from 6.7 GHz MERLIN observations of the massive protostar IRAS 18089-1732, where we find that the small-scale magnetic field probed by methanol masers is consistent with the large-scale magnetic field probed by dust (Dall’Olio et al. 2017 A&amp;A 607, A111). Moreover we present results obtained from our ALMA Band 7 polarisation observations of G9.62+0.20, which is a massive star-forming region with a sequence of cores at different evolutionary stages (Dall’Olio et al. submitted to A&amp;A). In this region we resolve several protostellar cores embedded in a bright and dusty filamentary structure. The magnetic field morphology and strength in different cores is related to the evolutionary sequence of the star formation process which is occurring across the filament.

Magnetic fields and massive star formation

AbstractMassive stars ( ${\rm{M}} &gt; \,8{M_ \odot }$ ) often form in parsec-scale molecular clumps that collapse and fragment, leading to the birth of a cluster of stellar objects. The role of magnetic fields during the formation of massive dense cores is still not clear. The steady improvement in sensitivity of (sub)millimeter interferometers over the past decade enabled observations of dust polarization of large samples of massive star formation regions. We carried out a polarimetric survey with the Submillimeter Array of 14 massive star forming clumps in continuum emission at a wavelength of 0.89 mm. This unprecedentedly large sample of massive star forming regions observed by a submillimeter interferometer before the advent of ALMA revealed compelling evidence of strong magnetic influence on the gas dynamics from 1 pc to 0.1 pc scales. We found that the magnetic fields in dense cores tend to be either parallel or perpendicular to the mean magnetic fields in their parental molecular clumps. Furthermore, the main axis of protostellar outflows does not appear to be aligned with the mean magnetic fields in the dense core where outflows are launched. These findings suggest that from 1 pc to 0.1 pc scales, magnetic fields are dynamically important in the collapse of clumps and the formation of dense cores. From the dense core scale to the accretion disk scale of ∼102 au, however, gravity and angular momentum appear to be more dominant relative to the magnetic field.

The birth environment of the solar system constrained by the relative abundances of the solar radionuclides

AbstractThe relative abundances of the radionuclides in the solar system at the time of its birth are crucial arbiters for competing hypotheses regarding the birth environment of the Sun. The presence of short-lived radionuclides, as evidenced by their decay products in meteorites, has been used to suggest that particular, sometimes exotic, stellar sources were proximal to the Sun’s birth environment. The recent confirmation of neutron star - neutron star (NS-NS) mergers and associated kilonovae as potentially dominant sources of r-process nuclides can be tested in the case of the solar birth environment using the relative abundances of the longer-lived nuclides. Critical analysis of the 15 radionuclides and their stable partners for which abundances and production ratios are well known suggests that the Sun formed in a typical massive star-forming region (SFR). The apparent overabundances of short-lived radionuclides (e.g. 26Al, 41Ca, 36Cl) in the early solar system appears to be an artifact of a heretofore under-appreciation for the important influences of enrichment by Wolf-Rayet winds in SFRs. The long-lived nuclides (e.g. 238U, 244Pu, 247Cr, 129I) are consistent with an average time interval between production events of 108 years, seemingly too short to be the products of NS-NS mergers alone. The relative abundances of all of these nuclides can be explained by their mean decay lifetimes and an average residence time in the ISM of ∼200 Myr. This residence time evidenced by the radionuclides is consistent with the average lifetime of dust in the ISM and the timescale for converting molecular cloud mass to stars.

3D shape of Orion A with Gaia DR2 An informed view on Star Formation Rates and Efficiencies

AbstractThe giant molecular cloud Orion A is the closest massive star-forming region to earth (d ∼ 400 pc). It contains the rich Orion Nebula Cluster (ONC) in the North, and low-mass star-forming regions (L1641, L1647) to the South. To get a better understanding of the differences in star formation activity, we perform an analysis of the gas mass distribution and star formation rate across the cloud. We find that the gas is roughly uniformly distributed, while, oddly, the ONC region produced about a factor of ten more stars compared to the rest of the cloud. For a better interpretation of this phenomenon, we use Gaia DR2 parallaxes, to analyse distances of young stellar objects, using them as proxy for cloud distances. We find that the ONC region indeed lies at about 400 pc while the low-mass star-forming parts are inclined about 70∘ from the plane of the sky reaching until ∼470 pc. With this we estimate that Orion A is an about 90 pc long filamentary cloud (about twice as long as previously assumed), with its “Head” (the ONC region) being “bent” and oriented towards the galactic mid-plane. This striking new view allows us to perform a more robust analysis of this important star-forming region in the future.

An Isothermal Outflow in High-mass Star-forming Region G240.31+0.07

Abstract We present Atacama Pathfinder EXperiment (APEX) observations toward the massive star-forming region G240.31+0.07 in the CO J = 3–2, 6–5, and 7–6 lines. We detect a parsec-sized, bipolar, and high-velocity outflow in all the lines, which allow us, in combination with the existing CO J = 2–1 data, to perform a multi-line analysis of physical conditions of the outflowing gas. The CO 7–6/6–5, 6–5/3–2, and 6–5/2–1 ratios are found to be nearly constant over a velocity range of ∼5–25 km s−1 for both blueshifted and redshifted lobes. We carry out rotation diagram and large velocity gradient calculations of the four lines and find that the outflow is approximately isothermal with a gas temperature of ∼50 K and that the CO column density clearly decreases with the outflow velocity. If the CO abundance and the velocity gradient do not vary much, the decreasing CO column density indicates a decline in the outflow gas density with velocity. By comparing with theoretical models of outflow driving mechanisms, our observations and calculations suggest that the massive outflow in G240.31+0.07 is being driven by a wide-angle wind and further support a disk-mediated accretion at play for the formation of the central high-mass star.

The TOP-SCOPE Survey of PGCCs: PMO and SCUBA-2 Observations of 64 PGCCs in the Second Galactic Quadrant

Abstract In order to understand the initial conditions and early evolution of star formation in a wide range of Galactic environments, we carried out an investigation of 64 Planck Galactic cold clumps (PGCCs) in the second quadrant of the Milky Way. Using the 13CO and C18O J = 1–0 lines and 850 μm continuum observations, we investigated cloud fragmentation and evolution associated with star formation. We extracted 468 clumps and 117 cores from the 13CO line and 850 μm continuum maps, respectively. We made use of the Bayesian distance calculator and derived the distances of all 64 PGCCs. We found that in general, the mass–size plane follows a relation of m ∼ r 1.67. At a given scale, the masses of our objects are around 1/10 of that of typical Galactic massive star-forming regions. Analysis of the clump and core masses, virial parameters, densities, and mass–size relation suggests that the PGCCs in our sample have a low core formation efficiency (∼3.0%), and most PGCCs are likely low-mass star-forming candidates. Statistical study indicates that the 850 μm cores are more turbulent, more optically thick, and denser than the 13CO clumps for star formation candidates, suggesting that the 850 μm cores are likely more appropriate future star formation candidates than the 13CO clumps.

Simultaneous Survey of Water and Class I Methanol Masers toward Red MSX Sources

Abstract We report simultaneous single-dish surveys of 22 GHz H2O and 44 and 95 GHz class I CH3OH masers toward 299 Red Midcourse Space Experiment Sources in the protostellar stage. The detection rates are 45% at 22 GHz, 28% at 44 GHz, and 23% at 95 GHz. There are 15, 53, and 51 new discoveries at 22, 44, and 95 GHz, respectively. We detect high-velocity (&gt;30 km s−1) features in 27 H2O maser sources. The 95 GHz maser emission is detected only in 44 GHz maser sources. The two transitions show strong correlations in the peak velocity, peak flux density, and isotropic maser luminosity, indicating that they are likely generated in the same sites by the same mechanisms. The 44 GHz masers have much narrower distributions than 22 GHz masers in the relative peak velocity and velocity range, while 6.7 GHz class II CH3OH masers have distributions intermediate between the two. The maser luminosity significantly correlates with the parental clump mass, while it correlates well with the bolometric luminosity of the central protostar only when data of the low-mass regime from the literature are added. Comparison with the results of previous maser surveys toward massive star-forming regions suggests that the detection rates of 22 and 44 GHz masers tend to increase as the central objects evolve. This is contrary to the trends found in low- and intermediate-mass star-forming regions. Thus, the occurrence of both masers might depend on the surrounding environments as well as on the evolution of the central object.

Radiative transfer modelling of W33A MM1: 3D structure and dynamics of a complex massive star-forming region

We present a composite model and radiative transfer simulations of the massive star forming core W33A MM1. The model was tailored to reproduce the complex features observed with ALMA at $\approx 0.2$ arcsec resolution in CH$_3$CN and dust emission. The MM1 core is fragmented into six compact sources coexisting within $\sim 1000$ au. In our models, three of these compact sources are better represented as disc-envelope systems around a central (proto)star, two as envelopes with a central object, and one as a pure envelope. The model of the most prominent object (Main) contains the most massive (proto)star ($M_\star\approx7~M_\odot$) and disc+envelope ($M_\mathrm{gas}\approx0.4~M_\odot$), and is the most luminous ($L_\mathrm{Main} \sim 10^4~L_\odot$). The model discs are small (a few hundred au) for all sources. The composite model shows that the elongated spiral-like feature converging to the MM1 core can be convincingly interpreted as a filamentary accretion flow that feeds the rising stellar system. The kinematics of this filament is reproduced by a parabolic trajectory with focus at the center of mass of the region. Radial collapse and fragmentation within this filament, as well as smaller filamentary flows between pairs of sources are proposed to exist. Our modelling supports an interpretation where what was once considered as a single massive star with a $\sim 10^3$ au disc and envelope, is instead a forming stellar association which appears to be virialized and to form several low-mass stars per high-mass object.

OH as an Alternate Tracer for Molecular Gas: Excitation Temperatures of the OH 18 cm Main Lines in W5

Abstract We present excitation temperatures T ex for the OH 18 cm main lines at 1665 and 1667 MHz measured directly in front of the W5 star-forming region, using observations from the Green Bank Telescope and the Very Large Array. We find unequivocally that T ex at 1665 MHz is greater than T ex at 1667 MHz. Our method exploits variations in the continuum emission from W5, and the fact that the continuum brightness temperatures T C in this nebula are close to the excitation temperatures of the OH lines in the foreground gas. The result is that an OH line can appear in emission in one location and in absorption in a neighboring location, and the value of T C where the profiles switch from emission to absorption indicates T ex. Absolute measurements of T ex for the main lines were subject to greater uncertainty because of unknown effects of geometry of the OH features. We also employed the traditional “expected profile” method for comparison with our “continuum background” method and found that the continuum background method provided more precise results and was the one to definitively show the T ex difference. Our best estimate values are K, K, and K. The T ex values we have measured for the ISM in front of W5 are similar to those found in the quiescent ISM, indicating that proximity to massive star-forming regions does not generally result in widespread anomalous excitation of OH emission.

A masing event in NGC 6334I: contemporaneous flaring of hydroxyl, methanol, and water masers

As a product of the maser monitoring program with the 26m telescope of the Hartebeesthoek Radio Astronomy Observatory (HartRAO), we present an unprecedented, contemporaneous flaring event of 10 maser transitions in hydroxyl, methanol, and water that began in 2015 January in the massive star-forming region NGC6334I in the velocity range -10 to -2 km/s. The 6.7 GHz methanol and 22.2 GHz water masers began flaring within 22 days of each other, while the 12.2 GHz methanol and 1665 MHz hydroxyl masers flared 80 and 113 days later respectively. The 1665 MHz, 6.7 GHz, and 22.2 GHz masers have all remained in their flared state for nearly 3 years. The brightest flaring components increased by factors of 66, 21, 26, and 20 in the 12.2 and 6.7 GHz methanol, 1665 MHz hydroxyl and 22.2 GHz water maser transitions respectively; some weaker components increased by up to a factor of 145. We also report new maser emission in the 1720, 6031, and 6035 MHz OH lines and the 23.1 GHz methanol line, along with the detection of only the fifth 4660 MHz OH maser. We note the correlation of this event with the extraordinary (sub)millimeter continuum outburst from the massive protostellar system NGC6334I-MM1 and discuss the implications of the observed time lags between different maser velocity components on the nature of the outburst. Finally, we identify two earlier epoch maser flaring events likely associated with this object, which suggest a recurring accretive phenomenon that generates powerful radiative outbursts.

Nitrogen fractionation in high-mass star-forming cores across the Galaxy

The fractionation of nitrogen (N) in star-forming regions is a poorly understood process. To put more stringent observational constraints on the N-fractionation, we have observed with the IRAM-30m telescope a large sample of 66 cores in massive star-forming regions. We targeted the (1-0) rotational transition of HN$^{13}$C, HC$^{15}$N, H$^{13}$CN and HC$^{15}$N, and derived the $^{14}$N/$^{15}$N ratio for both HCN and HNC. We have completed this sample with that already observed by Colzi et al. (2018), and thus analysed a total sample of 87 sources. The $^{14}$N/$^{15}$N ratios are distributed around the Proto-Solar Nebula value with a lower limit near the terrestrial atmosphere value ($\sim$272). We have also derived the $^{14}$N/$^{15}$N ratio as a function of the Galactocentric distance and deduced a linear trend based on unprecedented statistics. The Galactocentric dependences that we have found are consistent, in the slope, with past works but we have found a new local $^{14}$N/$^{15}$N value of $\sim$400, i.e. closer to the Prosolar Nebula value. A second analysis was done, and a parabolic Galactocentric trend was found. Comparison with Galactic chemical evolution models shows that the slope until 8 kpc is consistent with the linear analysis, while the flattening trend above 8 kpc is well reproduced by the parabolic analysis.

Astrochemical evolution along star formation: Overview of the IRAM Large Program ASAI.

Evidence is mounting that the small bodies of our Solar System, such as comets and asteroids, have at least partially inherited their chemical composition from the first phases of the Solar System formation. It then appears that the molecular complexity of these small bodies is most likely related to the earliest stages of star formation. It is therefore important to characterize and to understand how the chemical evolution changes with solar-type protostellar evolution. We present here the Large Program "Astrochemical Surveys At IRAM" (ASAI). Its goal is to carry out unbiased millimeter line surveys between 80 and 272 GHz of a sample of ten template sources, which fully cover the first stages of the formation process of solar-type stars, from prestellar cores to the late protostellar phase. In this article, we present an overview of the surveys and results obtained from the analysis of the 3 mm band observations. The number of detected main isotopic species barely varies with the evolutionary stage and is found to be very similar to that of massive star-forming regions. The molecular content in O- and C- bearing species allows us to define two chemical classes of envelopes, whose composition is dominated by either a) a rich content in O-rich complex organic molecules, associated with hot corino sources, or b) a rich content in hydrocarbons, typical of Warm Carbon Chain Chemistry sources. Overall, a high chemical richness is found to be present already in the initial phases of solar-type star formation.

The Massive Star-forming Regions Omnibus X-ray Catalog, Second Installment

Abstract We present the second installment of the Massive Star-forming Regions (MSFRs) Omnibus X-ray Catalog (MOXC2), a compilation of X-ray point sources detected in Chandra/ACIS observations of 16 Galactic MSFRs and surrounding fields. MOXC2 includes 13 ACIS mosaics, three containing a pair of unrelated MSFRs at different distances, with a total catalog of 18,396 point sources. The MSFRs sampled range over distances of 1.3 kpc to 6 kpc and populations varying from single massive protostars to the most massive Young Massive Cluster known in the Galaxy. By carefully detecting and removing X-ray point sources down to the faintest statistically significant limit, we facilitate the study of the remaining unresolved X-ray emission. Through comparison with mid-infrared images that trace photon-dominated regions and ionization fronts, we see that the unresolved X-ray emission is due primarily to hot plasmas threading these MSFRs, the result of feedback from the winds and supernovae of massive stars. The 16 MSFRs studied in MOXC2 more than double the MOXC1 sample, broadening the parameter space of ACIS MSFR explorations and expanding Chandra's substantial contribution to contemporary star formation science.

On the formation of runaway stars BN and x in the Orion Nebula Cluster

We explore scenarios for the dynamical ejection of stars BN and x from source I in the Kleinmann-Low nebula of the Orion Nebula Cluster (ONC), which is important because it is the closest region of massive star formation. This ejection would cause source I to become a close binary or a merger product of two stars. We thus consider binary-binary encounters as the mechanism to produce this event. By running a large suite of N-body simulations, we find that it is nearly impossible to match the observations when using the commonly adopted masses for the participants, especially a source I mass of 7 M⊙. The only way to recreate the event is if source I is more massive, that is, ~20 M⊙. However, even in this case, the likelihood of reproducing the observed system is low. We discuss the implications of these results for understanding this important star-forming region.

Infall and outflow motions towards a sample of massive star-forming regions from the RMS survey

We present the results of an outflow and infall survey towards a distance limited sample of 31 massive star forming regions drawn from the RMS survey. The presence of young, active outflows is identified from SiO (8-7) emission and the infall dynamics are explored using HCO$^+$/H$^{13}$CO$^+$ (4-3) emission. We investigate if the infall and outflow parameters vary with source properties, exploring whether regions hosting potentially young active outflows show similarities or differences with regions harbouring more evolved, possibly momentum driven, "fossil" outflows. SiO emission is detected towards approximately 46% of the sources. When considering sources with and without an SiO detection (i.e. potentially active and fossil outflows respectively), only the $^{12}$CO outflow velocity shows a significant difference between samples, indicating SiO is more prevalent towards sources with higher outflow velocities. Furthermore, we find the SiO luminosity increases as a function of the Herschel 70$\mu$m to WISE 22$\mu$m flux ratio, suggesting the production of SiO is prevalent in younger, more embedded regions. Similarly, we find tentative evidence that sources with an SiO detection have a smaller bolometric luminosity-to-mass ratio, indicating SiO (8-7) emission is associated with potentially younger regions. We do not find a prevalence towards sources displaying signatures of infall in our sample. However, the higher energy HCO$^+$ transitions may not be the best suited tracer of infall at this spatial resolution in these regions.