Source author record

Megan Reiter

Megan Reiter appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2022arXiv

Dynamics of young stellar clusters as planet forming environments

Most stars and thus most planetary systems do not form in isolation. The larger star-forming environment affects protoplanetary disks in multiple ways: gravitational interactions with other stars truncate disks and alter the architectures of exoplanet systems; external irradiation from nearby high-mass stars truncates disks and shortens their lifetimes; and remaining gas and dust in the environment affects dynamical evolution (if removed by feedback processes) and provides some shielding for disks from external irradiation. The dynamical evolution of the region regulates when and how long various feedback mechanisms impact protoplanetary disks. Density is a key parameter that regulates the intensity and duration of UV irradiation and the frequency of dynamical encounters. The evolution of larger star-forming complexes may also play an important role by mixing populations. Observations suggest that clusters are not a single-age population but multiple populations with small age differences which may be key to resolving several timescale issues (i.e., proplyd lifetimes, enrichment). In this review, we consider stellar clusters as the ecosystems in which most stars and therefore most planets form. We review recent observational and theoretical results and highlight upcoming contributions from facilities expected to begin observations in the next five years. Looking further ahead, we argue that the next frontier is large-scale surveys of low-mass stars in more distant high-mass star-forming regions. The future of ecosystem studies is bright as faint low-mass stars in more distant high-mass star-forming regions will be routinely observable in the era of Extremely Large Telescopes (ELTs).

preprint2021arXiv

Binary Formation in the Orion Nebula Cluster: Exploring the Sub-stellar Limit

We present results constraining the multiplicity of the very low mass stars and sub-stellar objects in the Orion Nebula Cluster (ONC). Our sample covers primary masses 0.012-0.1M$_{\odot}$ using archival Hubble Space Telescope data obtained with the Advanced Camera for Surveys using multiple filters. Studying the binary populations of clusters provides valuable constraints of how the birth environment affects binary formation and evolution. Prior surveys have shown that the binary populations of high-mass, high-density star clusters like the ONC may substantially differ from those in low-mass associations. Very low mass stellar and sub-stellar binaries at wide separations, $>$20AU, are statistically rare in the Galactic field and have been identified in stellar associations like Taurus-Auriga and Ophiuchus. They also may be susceptible to dynamical interactions, and their formation may be suppressed by feedback from on-going star formation. We implement a double point-spread function (PSF) fitting algorithm using empirical, position dependent PSF models to search for binary companions at projected separations $>$ 10 AU (25 mas). With this technique, we identify 7 very low mass binaries, 5 of which are new detections, resulting in a binary frequency of 12$^{+6}_{-3.2}$% over mass ratios of 0.5 - 1.0 and projected separations of 20 - 200 AU. We find an excess of very low mass binaries in the ONC compared to the Galactic field, with a probability of 10$^{-6}$ that the populations are statistically consistent. The sub-stellar population of the ONC may require further dynamical processing of the lowest binding energy binaries to resemble the field population.

preprint2020arXiv

Evolution of Stellar Feedback in HII Regions

Stellar feedback is needed to produce realistic giant molecular clouds (GMCs) and galaxies in simulations, but due to limited numerical resolution, feedback must be implemented using subgrid models. Observational work is an important means to test and anchor these models, but limited studies have assessed the relative dynamical role of multiple feedback modes, particularly at the earliest stages of expansion when HII regions are still deeply embedded. In this paper, we use multiwavelength (radio, infrared, and X-ray) data to measure the pressures associated with direct radiation ($P_{\rm dir}$), dust-processed radiation ($P_{\rm IR}$), photoionization heating ($P_{\rm HII}$), and shock-heating from stellar winds ($P_{\rm X}$) in a sample of 106 young, resolved HII regions with radii $\lesssim$0.5 pc to determine how stellar feedback drives their expansion. We find that the $P_{\rm IR}$ dominates in 84% of the regions and that the median $P_{\rm dir}$ and $P_{\rm HII}$ are smaller than the median $P_{\rm IR}$ by factors of $\approx 6$ and $\approx 9$, respectively. Based on the radial dependences of the pressure terms, we show that HII regions transition from $P_{\rm IR}$-dominated to $P_{\rm HII}$-dominated at radii of $\sim$3 pc. We find a median trapping factor of $f_{\rm trap} \sim$ 8 without any radial dependence for the sample, suggesting this value can be adopted in sub-grid feedback models. Moreover, we show that the total pressure is greater than the gravitational pressure in the majority of our sample, indicating that the feedback is sufficient to expel gas from the regions.

preprint2020arXiv

Illuminating a tadpole's metamorphosis II: observing the on-going transformation with ALMA

We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations of the tadpole, a small globule in the Carina Nebula that hosts the HH 900 jet+outflow system. Our data include $^{12}$CO, $^{13}$CO, C$^{18}$O J=2-1, $^{13}$CO, C$^{18}$O J=3-2, and serendipitous detections of DCN J=3-2 and CS J=7-6. With angular resolution comparable to the Hubble Space Telescope (HST), our data reveal for the first time the bipolar molecular outflow in CO, seen only inside the globule, that is launched from the previously unseen jet-driving protostar (the HH 900 YSO). The biconical morphology joins smoothly with the externally irradiated outflow seen in ionized gas tracers outside the globule, tracing the overall morphology of a jet-driven molecular outflow. Continuum emission at the location of the HH 900 YSO appears to be slightly flattened perpendicular to outflow axis. Model fits to the continuum have a best-fit spectral index of $\sim 2$, suggesting cold dust and the onset of grain growth. In position-velocity space, $^{13}$CO and C$^{18}$O gas kinematics trace a C-shaped morphology, similar to infall profiles seen in other sources, although the global dynamical behaviour of the gas remains unclear. Line profiles of the CO isotopologues display features consistent with externally heated gas. We estimate a globule mass of $\sim 1.9$ M$_{\odot}$, indicating a remaining lifetime of $\sim 4$ Myr, assuming a constant photoevaporation rate. This long globule lifetime will shield the disk from external irradiation perhaps prolonging its life and enabling planet formation in regions where disks are typically rapidly destroyed.

preprint2020arXiv

Illuminating a tadpole's metamorphosis III: quantifying past and present environmental impact

We combine MUSE and ALMA observations with theoretical models to evaluate how a tadpole-shaped globule located in the Carina Nebula has been influenced by its environment. This globule is now relatively small (radius ~2500 au), hosts a protostellar jet+outflow (HH 900) and, with a blue-shifted velocity of ~10 km/s, is travelling faster than it should be if its kinematics were set by the turbulent velocity dispersion of the precursor cloud. Its outer layers are currently still subject to heating, but comparing the internal and external pressures implies that the globule is in a post-collapse phase. Intriguingly the outflow is bent, implying that the YSO responsible for launching it is comoving with the globule, which requires that the star formed after the globule was up to speed since otherwise it would have been left behind. We conclude that the most likely scenario is one in which the cloud was much larger before being subject to radiatively-driven implosion, which accelerated the globule to the high observed speeds under the photoevaporative rocket effect and triggered the formation of the star responsible for the outflow. The globule may now be in a quasi-steady state following collapse. Finally, the HH 900 YSO is likely $\gtrsim$1 M$_{\odot}$ and may be the only star forming in the globule. It may be that this process of triggered star formation has prevented the globule from fragmenting to form multiple stars (e.g., due to heating) and has produced a single higher mass star.

preprint2020arXiv

Observational constraints on the likelihood of $^{26}$Al in planet-forming environments

Recent work suggests that $^{26}$Al may determine the water budget in terrestrial exoplanets as its radioactive decay dehydrates planetesimals leading to rockier compositions. Here I consider the observed distribution of $^{26}$Al in the Galaxy and typical star-forming environments to estimate the likelihood of $^{26}$Al enrichment during planet formation. I do not assume Solar-System-specific constraints as I am interested in enrichment for exoplanets generally. Observations indicate that high-mass stars dominate the production of $^{26}$Al with nearly equal contributions from their winds and supernovae. $^{26}$Al abundances are comparable to those in the early Solar System in the high-mass star-forming regions where most stars (and thereby most planets) form. These high abundances appear to be maintained for a few Myr, much longer than the 0.7 Myr half-life. Observed bulk $^{26}$Al velocities are an order of magnitude slower than expected from winds and supernovae. These observations are at odds with typical model assumptions that $^{26}$Al is provided instantaneously by high velocity mass loss from supernovae and winds. Regular replenishment of $^{26}$Al especially when coupled with the small age differences that are common in high-mass star-forming complexes, may significantly increase the number of star/planet-forming systems exposed to $^{26}$Al. Exposure does not imply enrichment, but the order of magnitude slower velocity of $^{26}$Al may alter the fraction that is incorporated into planet-forming material. Together, this suggests that the conditions for rocky planet formation are not rare, nor are they ubiquitous, as small regions like Taurus that lack high-mass stars to produce $^{26}$Al may be less likely to form rocky planets. I conclude with suggested directions for future studies.

preprint2019arXiv

A Search for Intermediate Separation Low Mass Binaries in the Orion Nebula Cluster

We present the results of a binary population study in the Orion Nebula Cluster (ONC) using archival Hubble Space Telescope (HST) data obtained with the Advanced Camera for Surveys (ACS) in Johnson V filter (HST Proposal 10246, PI M. Robberto). Young clusters and associations hold clues to the origin and properties of multiple star systems. Binaries with separations $< 100 $ AU are useful as tracers of the initial binary population since they are not as likely to be destroyed through dynamical interactions. Low mass, low stellar density star-forming regions such as Taurus-Auriga, reveal an excess of multiples compared to the Galactic Field. Studying the binary population of higher mass, higher stellar density star-forming regions like the ONC provides useful information concerning the origin of the Galactic Field star population. In this survey, we characterize the previously unexplored (and incomplete) separation parameter space of binaries in the ONC (15 - 160 AU) by fitting a double-PSF model built from empirical PSFs. We identified 14 candidate binaries (11 new detections) and find that 8$_{-2\%}^{+4\%}$ of our observed sample are in binary systems, complete over mass ratios and separations of 0.6 $< $ q $< $ 1.0 and 30 $< $ a $< $ 160 AU. This is consistent with the Galactic Field M-dwarf population over the same parameter ranges, 6.5\% $\pm$ 3\%. Therefore, high mass star forming regions like the ONC would not require further dynamical evolution for their binary population to resemble the Galactic Field, as some models have hypothesized for young clusters.

preprint2016arXiv

[Fe II] jets from intermediate-mass protostars in Carina

We present new HST/WFC3-IR narrowband [Fe II] images of protostellar jets in the Carina Nebula. Combined with 5 previously published sources, we have a sample of 18 jets and 2 HH objects. All of the jets we targeted with WFC3 show bright infrared [Fe II] emission, and a few H$α$ candidate jets are confirmed as collimated outflows based on the morphology of their [Fe II] emission. Continuum-subtracted images clearly separate jet emission from the adjacent ionization front, providing a better tracer of the collimated jet than H$α$ and allowing us to connect these jets with their embedded driving sources. The [Fe II] 1.64 $μ$m/H$α$ flux ratio measured in the jets is $\gtrsim 5$ times larger than in the adjacent ionization fronts. The low-ionization jet core requires high densities to shield Fe$^+$ against further ionization by the FUV radiation from O-type stars in the H II region. High jet densities imply high mass-loss rates, consistent with the intermediate-mass driving sources we identify for 13 jets. The remaining jets emerge from opaque globules that obscure emission from the protostar. In many respects, the HH jets in Carina look like a scaled-up version of the jets driven by low-mass protostars. Altogether, these observations suggest that [Fe II] emission is a reliable tracer of dense, irradiated jets driven by intermediate-mass protostars. We argue that highly collimated outflows are common to more massive protostars, and that they suggest the outflow physics inferred for low-mass stars formation scales up to at least $\sim8$ M$_{\odot}$.

preprint2016arXiv

Ancient eruptions of Eta Carinae: A tale written in proper motions

We analyze eight epochs of Hubble Space Telescope H$α$+[N II] imaging of Eta Carinae's outer ejecta. Proper motions of nearly 800 knots reveal that the detected ejecta are divided into three apparent age groups, dating to around 1250 A.D., to around 1550 A.D., and to during or shortly before the Great Eruption of the 1840s. Ejecta from these groups reside in different locations and provide a firm constraint that Eta Car experienced multiple major eruptions prior to the 19th century. The 1250 and 1550 events did not share the same axisymmetry as the Homunculus; the 1250 event was particularly asymmetric, even one-sided. In addition, the ejecta in the S ridge, which have been associated with the Great Eruption, appear to predate the ejection of the Homunculus by several decades. We detect essentially ballistic expansion across multiple epochs. We find no evidence for large-scale deceleration of the observed knots that could power the soft X-ray shell by plowing into surrounding material, suggesting that the observed X-rays arise instead from fast, rarefied ejecta from the 1840s overtaking the older dense knots. Early deceleration and subsequent coasting cannot explain the origin of the older outer ejecta---significant episodic mass loss prior to the 19th century is required. The timescale and geometry of the past eruptions provide important constraints for any theoretical physical mechanisms driving Eta Car's behavior. Non-repeating mechanisms such as the merger of a close binary in a triple system would require additional complexities to explain the observations.

preprint2015arXiv

A Spitzer/IRAC Characterization of Galactic AGB and RSG Stars

We present new Spitzer/IRAC observations of 55 dusty Long Period Variables (LPVs, 48 AGB and 6 RSG stars) in the Galaxy that have different chemistry, variability type, and mass-loss rate. O-rich AGB stars (including intrinsic S-type) tend to have redder [3.6]-[8.0] colors than carbon stars for a given [3.6]-[4.5] color due to silicate features increasing the flux in the 8.0 μm IRAC band. For colors including the 5.8 μm band, carbon stars separate into two distinct sequences, likely due to a variable photospheric C$_3$ feature that is only visible in relatively unobscured, low mass-loss rate sources. Semiregular variables tend to have smaller IR excess in [3.6]-[8.0] color than Miras, consistent with the hypothesis that semiregular variables lose mass discontinuously. Miras have redder colors for longer periods while semiregular variables do not. Galactic AGB stars follow the period-luminosity sequences found for the Magellanic Clouds. Mira variables fall along the fundamental pulsation sequence, while semiregular variables are mostly on overtone sequences. We also derive a relationship between mass-loss rate and [3.6]-[8.0] color. The fits are similar in shape to those found by other authors for AGBs in the LMC, but discrepant in overall normalization, likely due to different assumptions in the models used to derive mass-loss rates. We find that IR colors are not unique discriminators of chemical type, suggesting caution when using color selection techniques to infer the chemical composition of AGB dust returned to the ISM.

preprint2015arXiv

Disentangling the outflow and protostars in HH 900 in the Carina Nebula

HH 900 is a peculiar protostellar outflow emerging from a small, tadpole-shaped globule in the Carina nebula. Previous Hα imaging with HST/ACS showed an ionized outflow with a wide opening angle that is distinct from the highly collimated structures typically seen in protostellar jets. We present new narrowband near-IR [Fe II] images taken with the Wide Field Camera 3 on the Hubble Space Telescope that reveal a remarkably different structure than Hα. In contrast to the unusual broad Hα outflow, the [Fe II] emission traces a symmetric, collimated bipolar jet with the morphology and kinematics that are more typical of protostellar jets. In addition, new Gemini adaptive optics images reveal near-IR H$_2$ emission coincident with the Hα emission, but not the [Fe II]. Spectra of these three components trace three separate and distinct velocity components: (1) H$_2$ from the slow, entrained molecular gas, (2) Hα from the ionized skin of the accelerating outflow sheath, and (3) [Fe II] from the fast, dense, and collimated protostellar jet itself. Together, these data require a driving source inside the dark globule that remains undetected behind a large column density of material. In contrast, Hα and H$_2$ emission trace the broad outflow of material entrained by the jet, which is irradiated outside the globule. As it get dissociated and ionized, it remains visible for only a short time after it is dragged into the H II region.

preprint2015arXiv

HH 666: Different kinematics from Hα and [Fe II] emission provide a missing link between jets and outflows

HH 666 is an externally irradiated protostellar outflow in the Carina Nebula for which we present new near-IR [Fe II] spectra obtained with the FIRE spectrograph at Magellan Observatory. Earlier Hα and near-IR [Fe II] imaging revealed that the two emission lines trace substantially different morphologies in the inner ~40" of the outflow. Hα traces a broad cocoon that surrounds the collimated [Fe II] jet that extends throughout the parent dust pillar. New spectra show that this discrepancy extends to their kinematics. Near-IR [Fe II] emission traces steady, fast velocities of +/- 200 km/s from the eastern and western limbs of the jet. We compare this to a previously published Hα spectrum that reveals a Hubble-flow velocity structure near the jet-driving source. New, second-epoch HST/ACS Hα images reveal the lateral spreading of the Hα outflow lobe away from the jet axis. Hα proper motions also indicate a sudden increase in the mass-loss rate ~1000 yr ago, while steady [Fe II] emission throughout the inner jet suggests that the burst is ongoing. An accretion burst sustained for ~1000 yr is an order of magnitude longer than expected for FU Orionis outbursts, but represents only a small fraction of the total age of the HH 666 outflow. Altogether, available data suggests that [Fe II] traces the highly collimated protostellar jet while Hα traces the entrained and irradiated outflow. HH 666 appears to be a missing link between bare jets seen in H II regions and entrained molecular outflows seen from embedded protostars in more quiescent regions.

preprint2014arXiv

Kinematics of powerful jets from intermediate-mass protostars in the Carina nebula

We present measurements of proper motions and radial velocities of four powerful Herbig-Haro (HH) jets in the Carina nebula: HH 666, HH 901, HH 902, and HH 1066. Two epochs of Hubble Space Telescope (HST) imaging separated by a time baseline of ~4.4 years provide proper motions that allow us to measure the transverse velocities of the jets, while ground-based spectra sample their Doppler velocities. Together these yield full three-dimensional space velocities. Aside from HH 666, their identification as outflows was previously inferred only from morphology in images. Proper motions now show decisively that these objects are indeed jets, and confirm that the intermediate-mass protostars identified as the candidate driving sources for HH 666 and HH 1066 are indeed the origin of these outflows. The appearance of two new knots in the HH 1066 jet suggest recent (~35 yr) changes in the accretion rate, underscoring the variable nature of accretion and outflow in the formation of intermediate-mass stars. In fact, kinematics and mass-ejection histories for all the jets suggest highly episodic mass loss, and point toward pronounced accretion fluctuations. Overall, we measure velocities similar to those found for low-mass protostars. However, the HH jets in Carina have higher densities and are more massive than their low-mass counterparts. Coarse estimates suggest that the heavy jets of intermediate-mass protostars can compete with or even exceed inject ~10 or more times the cumulative momentum injection of lower-mass protostars.

preprint2013arXiv

HD 106906 b: A planetary-mass companion outside a massive debris disk

We report the discovery of a planetary-mass companion, HD 106906 b, with the new Magellan Adaptive Optics (MagAO) + Clio2 system. The companion is detected with Clio2 in three bands: $J$, $K_S$, and $L^\prime$, and lies at a projected separation of 7.1'' (650 AU). It is confirmed to be comoving with its $13\pm2$ Myr-old F5 host using Hubble Space Telescope/Advanced Camera for Surveys astrometry over a time baseline of 8.3 yr. DUSTY and COND evolutionary models predict the companion's luminosity corresponds to a mass of $11\pm2 M_{Jup}$, making it one of the most widely separated planetary-mass companions known. We classify its Magellan/Folded-Port InfraRed Echellette $J/H/K$ spectrum as L$2.5\pm1$; the triangular $H$-band morphology suggests an intermediate surface gravity. HD 106906 A, a pre-main-sequence Lower Centaurus Crux member, was initially targeted because it hosts a massive debris disk detected via infrared excess emission in unresolved Spitzer imaging and spectroscopy. The disk emission is best fit by a single component at 95 K, corresponding to an inner edge of 15-20 AU and an outer edge of up to 120 AU. If the companion is on an eccentric ($e>0.65$) orbit, it could be interacting with the outer edge of the disk. Close-in, planet-like formation followed by scattering to the current location would likely disrupt the disk and is disfavored. Furthermore, we find no additional companions, though we could detect similar-mass objects at projected separations $>35$ AU. In situ formation in a binary-star-like process is more probable, although the companion-to-primary mass ratio, at $<1\%$, is unusually small.

preprint2013arXiv

HST/WFC3 Imaging of Protostellar Jets in Carina: [Fe II] Emission Tracing Massive Jets from Intermediate Mass Protostars

We present narrowband WFC3-UVIS and -IR images of four externally irradiated protostellar jets in the Carina nebula: HH666, HH901, HH902, and HH1066. These massive jets are unusual because they are bathed in UV radiation from dozens of nearby O stars, but despite the strong incident ionizing radiation, portions of the jet remain neutral. Near-IR [Fe II] images reveal dense, neutral gas that was not seen in previous studies of H-alpha emission. We show that near-IR [Fe II] emitting gas must be self-shielded from Lyman continuum photons, regardless of its excitation mechanism (shocks, FUV radiation, or both). High densities are required for the survival of Fe+ amid the strong Lyman continuum from Tr14, raising estimates of the mass-loss rates by an order of magnitude. Higher jet mass-loss rates require higher accretion rates onto their driving protostars, implying that these jets are driven by intermediate-mass (IM; ~2-8 Msun) stars. Indeed, the IR driving sources of two of these outflows have luminosities that require IM protostars. All four of these HH jets are highly collimated, with opening angles of only a few degrees, similar to those observed in low-mass protostars. We propose that these jets reflect essentially the same outflow phenomenon seen in wide-angle molecular outflows associated with IM and high-mass protostars, but that the collimated atomic jet core is irradiated and rendered observable in Carina's harsh radiative environment. In more quiescent environments, this atomic core remains invisible, and outflows traced by shock-excited molecules in the outflow cavity give the impression that these outflows have a wider opening angle. Thus, the externally irradiated jets in Carina constitute a new view of collimated jets from IM protostars, and offer strong additional evidence that stars up to at least ~8 Msun form by the same accretion mechanisms as low-mass stars.

preprint2011arXiv

Evidence for Inflow in High-Mass Star-Forming Clumps

We analyze the HCO+ 3-2 and H13CO+ 3-2 line profiles of 27 high-mass star-forming regions to identify asymmetries that are suggestive of mass inflow. Three quantitative measures of line asymmetry are used to indicate whether a line profile is blue, red or neither - the ratio of the temperature of the blue and red peaks, the line skew and the dimensionless parameter delta_v. We find nine HCO+ 3-2 line profiles with a significant blue asymmetry and four with significant red asymmetric profiles. Comparing our HCO+ 3-2 results to HCN 3-2 observations from Wu et al. (2003, 2010), we find that eight of the blue and three of red have profiles with the same asymmetry in HCN. The eight sources with blue asymmetries in both tracers are considered strong candidates for inflow. Quantitative measures of the asymmetry (e.g. delta_v) tend to be larger for HCN. This, combined with possible HCO+ abundance enhancements in outflows, suggests that HCN may be a better tracer of inflow. Understanding the behavior of common molecular tracers like HCO+ in clumps of different masses is important for properly analyzing the line profiles seen in a sample of sources representing a broad range of clump masses. Such studies will soon be possible with the large number of sources with possible self-absorption seen in spectroscopic follow-up observations of clumps identified in the Bolocam Galactic Plane Survey.

preprint2011arXiv

Testing the global star formation relation: An HCO+ (3-2) mapping study of Red MSX sources in the Bolocam Galactic Plane Survey

We present an analysis of the relation between the star formation rate (SFR) and mass of dense gas in Galactic clumps and nearby galaxies. Using the bolometric luminosity as a measure of SFR and the molecular line luminosity of HCO+ (3-2) as a measure of dense gas mass, we find that the relation between SFR and M_{dense} is approximately linear. This is similar to published results derived using HCN (1-0) as a dense gas tracer. HCO+ (3-2) and HCN (1-0) have similar conditions for excitation. Our work includes 16 Galactic clumps that are in both the Bolocam Galactic Plane Survey and the Red MSX Survey, 27 water maser sources from the literature, and the aforementioned HCN (1-0) data. Our results agree qualitatively with predictions of recent theoretical models which state that the nature of the relation should depend on how the critical density of the tracer compares with the mean density of the gas.

preprint2011arXiv

The Physical Properties of High-Mass Star-Forming Clumps: A Systematic Comparison of Molecular Tracers

We present observations of HCO+ and H^13CO+, N2H+, HCS+, HCN and HN^13C, SO and ^34SO, CCH, SO_2, and CH_3OH-E towards a sample of 27 high-mass clumps coincident with water maser emission. All transitions are observed with or convolved to nearly identical resolution (30"), allowing for inter-comparison of the clump properties derived from the mapped transitions. We find N2H+ emission is spatially differentiated compared to the dust and the other molecules towards a few very luminous cores (10 of 27) and the N2H+ integrated intensity does not correlate well with dust continuum flux. We calculate the effective excitation density, n_eff, the density required to excite a 1 K line in T_kin=20 K gas for each molecular tracer. The intensity of molecular tracers with larger effective excitation densities (n_eff > 10^5 cm^-3) appear to correlate more strongly with the submillimeter dust continuum intensity. The median sizes of the clumps are anti-correlated with the n_eff of the tracers (which span more than three orders of magnitude). Virial mass is not correlated with n_eff, especially where the lines are optically thick as the linewidths may be broadened significantly by non-virial motions. The median mass surface density and median volume density of the clumps is correlated with n_eff indicating the importance of understanding the excitation conditions of the molecular tracer when deriving the average properties of an ensemble of cores.