Search

Understanding the chemical processes during starless core and prestellar core evolution is an important step in understanding the initial stages of star and disk formation. This project is a study of deuterated ammonia, o-NH$_2$D, in the L1251 star-forming region toward Cepheus. Twenty-two dense cores (twenty of which are starless or prestellar, and two of which have a protostar), previously identified by p-NH$_3$ (1,1) observations, were targeted with the 12m Arizona Radio Observatory telescope on Kitt Peak. o-NH$_2$D J$_{\rm{K_a} \rm{K_c}}^{\pm} =$ $1_{11}^{+} \rightarrow 1_{01}^{-}$ was detected in 13 (59\%) of the NH$_3$-detected cores with a median sensitivity of $\sigma_{T_{mb}} = 17$ mK. All cores detected in o-NH$_2$D at this sensitivity have p-NH$_3$ column densities $> 10^{14}$ cm$^{-2}$. The o-NH$_2$D column densities were calculated using the constant excitation temperature (CTEX) approximation while correcting for the filling fraction of the NH$_3$ source size. The median deuterium fraction was found to be 0.11 (including 3$\sigma$ upper limits). However, there are no strong, discernible trends in plots of deuterium fraction with any physical or evolutionary variables. If the cores in L1251 have similar initial chemical conditions, then this result is evidence of the cores physically evolving at different rates., Comment: 18 pages, 14 figures, MNRAS, in press

The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular clouds. We find that the results on relative alignment between cores and background magnetic field in synthetic observations slightly disagree with those measured in fully 3D simulation data, which is partly because cores identified in projected 2D maps tend to coexist within filamentary structures, while 3D cores are generally more rounded. In addition, we examine the progression of magnetic field from pc- to core-scale in the simulation, which is consistent with the anisotropic core formation model that gas preferably flow along the magnetic field toward dense cores. When comparing the observed cores identified from the GBT Ammonia Survey (GAS) and Planck polarization-inferred magnetic field orientations, we find that the relative core-field alignment has a regional dependence among different clouds. More specifically, we find that dense cores in the Taurus molecular cloud tend to align perpendicular to the background magnetic field, while those in Perseus and Ophiuchus tend to have random (Perseus) or slightly parallel (Ophiuchus) orientations with respect to the field. We argue that this feature of relative core-field orientation could be used to probe the relative significance of the magnetic field within the cloud., 18 pages, 11 figures, accepted for publication in MNRAS

Dynamical studies of dense structures within molecular clouds often conclude that the most massive clumps contain too little kinetic energy for virial equilibrium, unless they are magnetized to an unexpected degree. This raises questions about how such a state might arise, and how it might persist long enough to represent the population of massive clumps. In an effort to re-examine the origins of this conclusion, we use ammonia line data from the Green Bank Ammonia Survey and Planck-calibrated dust emission data from Herschel to estimate the masses and kinetic and gravitational energies for dense clumps in the Gould Belt clouds. We show that several types of systematic error can enhance the appearance of low kinetic-to-gravitational energy ratios: insufficient removal of foreground and background material; ignoring the kinetic energy associated with velocity differences across a resolved cloud; and over-correcting for stratification when evaluating the gravitational energy. Using an analysis designed to avoid these errors, we find that the most massive Gould Belt clumps harbor virial motions, rather than sub-virial ones. As a byproduct, we present a catalog of masses, energies, and virial energy ratios for 85 Gould Belt clumps., Comment: Submitted to ApJ

We present Herschel SPIRE and PACS maps of the Cepheus Flare clouds L1157, L1172, L1228, L1241, and L1251, observed by the Herschel Gould Belt Survey (HGBS) of nearby star-forming molecular clouds. Through modified blackbody fits to the SPIRE and PACS data, we determine typical cloud column densities of 0.5-1.0 $\times$ 10$^{21}$ cm$^{-2}$ and typical cloud temperatures of 14-15 K. Using the getsources identification algorithm, we extract 832 dense cores from the SPIRE and PACS data at 160-500 $\mu$m. From placement in a mass vs. size diagram, we consider 303 to be candidate prestellar cores, and 178 of these to be "robust" prestellar cores. From an independent extraction of sources at 70 $\mu$m, we consider 25 of the 832 dense cores to be protostellar. The distribution of background column densities coincident with candidate prestellar cores peaks at 2-4 $\times$ 10$^{21}$ cm$^{-2}$. About half of the candidate prestellar cores in Cepheus may have formed due to the widespread fragmentation expected to occur within filaments of "transcritical" line mass. The lognormal robust prestellar core mass function (CMF) drawn from all five Cepheus clouds peaks at 0.56 M$_{\odot}$ and has a width of $\sim$0.5 dex, similar to that of Aquila's CMF. Indeed, the width of Cepheus' aggregate CMF is similar to the stellar system Initial Mass Function (IMF). The similarity of CMF widths in different clouds and the system IMF suggests a common, possibly turbulent origin for seeding the fluctuations that evolve into prestellar cores and stars., Comment: 36 pages, 21 figures, 2 online text tables (not yet included here), accepted for publication in The Astrophysical Journal

When multiple star-forming gas structures overlap along the line-of-sight and emit optically thin emission at significantly different radial velocities, the emission can become non-Gaussian and often exhibits two distinct peaks. Traditional line-fitting techniques can fail to account adequately for these double-peaked profiles, providing inaccurate cloud kinematics measurements. We present a new method called Convnet Line-fitting Of Velocities in Emission-line Regions (CLOVER) for distinguishing between one-component, two-component, and noise-only emission lines using 1D convolutional neural networks trained with synthetic spectral cubes. CLOVER utilizes spatial information in spectral cubes by predicting on $3\times3$ pixel sub-cubes, using both the central pixel's spectrum and the average spectrum over the $3\times3$ grid as input. On an unseen set of 10,000 synthetic spectral cubes in each predicted class, CLOVER has classification accuracies of $\sim99\%$ for the one-component class and $\sim97\%$ for the two-component class. For the noise-only class, which is analogous to a signal-to-noise cutoff of four for traditional line-fitting methods, CLOVER has classification accuracy of $100\%$. CLOVER also has exceptional performance on real observations, correctly distinguishing between the three classes across a variety of star-forming regions. In addition, CLOVER quickly and accurately extracts kinematics directly from spectra identified as two-component class members. Moreover, we show that CLOVER is easily scalable to emission lines with hyperfine splitting, making it an attractive tool in the new era of large-scale NH$_3$ and N$_2$H$^+$ mapping surveys., Accepted for publication in ApJ

We perform a virial analysis of starless dense cores in three nearby star-forming regions : L1688 in Ophiuchus, NGC 1333 in Perseus, and B18 in Taurus. Our analysis takes advantage of comprehensive kinematic information for the dense gas in all of these regions made publicly available through the Green Bank Ammonia Survey Data Release 1, which used to estimate internal support against collapse. We combine this information with ancillary data used to estimate other important properties of the cores, including continuum data from the James Clerk Maxwell Telescope Gould Belt Survey for core identification, core masses, and core sizes. Additionally, we used \textit{Planck} and \textit{Herschel}-based column density maps for external cloud weight pressure, and Five College Radio Astronomy Observatory $^{13}$CO observations for external turbulent pressure. Our self-consistent analysis suggests that many dense cores in all three star-forming regions are not bound by gravity alone, but rather require additional pressure confinement to remain bound. Unlike a recent, similar study in Orion~A, we find that turbulent pressure represents a significant portion of the external pressure budget. Our broad conclusion emphasizing the importance of pressure confinement in dense core evolution, however, agrees with earlier work., 35 pages, 8 tables, and 14 figures consisting of 16 .pdf files. Accepted for publication in the Astrophysical Journal

We present the observation and analysis of newly discovered coherent structures in the L1688 region of Ophiuchus and the B18 region of Taurus. Using data from the Green Bank Ammonia Survey, we identify regions of high density and near-constant, almost-thermal velocity dispersion. We reveal 18 coherent structures are revealed, 12 in L1688 and 6 in B18, each of which shows a sharp “transition to coherence” in velocity dispersion around its periphery. The identification of these structures provides a chance to statistically study the coherent structures in molecular clouds. The identified coherent structures have a typical radius of 0.04 pc and a typical mass of 0.4 M ☉, generally smaller than previously known coherent cores identified by Goodman et al., Caselli et al., and Pineda et al. We call these structures “droplets.” We find that, unlike previously known coherent cores, these structures are not virially bound by self-gravity and are instead predominantly confined by ambient pressure. The droplets have density profiles shallower than a critical Bonnor–Ebert sphere, and they have a velocity (V LSR) distribution consistent with the dense gas motions traced by NH3 emission. These results point to a potential formation mechanism through pressure compression and turbulent processes in the dense gas. We present a comparison with a magnetohydrodynamic simulation of a star-forming region, and we speculate on the relationship of droplets with larger, gravitationally bound coherent cores, as well as on the role that droplets and other coherent structures play in the star formation process.

We present initial results from the K-band focal plane array Examinations of Young STellar Object Natal Environments (KEYSTONE) survey, a large project on the 100-m Green Bank Telescope mapping ammonia emission across eleven giant molecular clouds at distances of $0.9-3.0$ kpc (Cygnus X North, Cygnus X South, M16, M17, MonR1, MonR2, NGC2264, NGC7538, Rosette, W3, and W48). This data release includes the NH$_3$ (1,1) and (2,2) maps for each cloud, which are modeled to produce maps of kinetic temperature, centroid velocity, velocity dispersion, and ammonia column density. Median cloud kinetic temperatures range from $11.4\pm2.2$ K in the coldest cloud (MonR1) to $23.0\pm6.5$ K in the warmest cloud (M17). Using dendrograms on the NH$_3$ (1,1) integrated intensity maps, we identify 856 dense gas clumps across the eleven clouds. Depending on the cloud observed, $40-100\%$ of the clumps are aligned spatially with filaments identified in H$_2$ column density maps derived from SED-fitting of dust continuum emission. A virial analysis reveals that 523 of the 835 clumps ($\sim63\%$) with mass estimates are bound by gravity alone. We find no significant difference between the virial parameter distributions for clumps aligned with the dust-continuum filaments and those unaligned with filaments. In some clouds, however, hubs or ridges of dense gas with unusually high mass and low virial parameters are located within a single filament or at the intersection of multiple filaments. These hubs and ridges tend to host water maser emission, multiple 70$\mu$m-detected protostars, and have masses and radii above an empirical threshold for forming massive stars., Comment: Accepted for publication in ApJ

Recent observations of global velocity gradients across and along molecular filaments have been interpreted as signs of gas accreting onto and along these filaments, potentially feeding star-forming cores and proto-clusters. The behavior of velocity gradients in filaments, however, has not been studied in detail, particularly on small scales (< 0.1 pc). In this paper, we present MUFASA, an efficient, robust, and automatic method to fit ammonia lines with multiple velocity components, generalizable to other molecular species. We also present CRISPy, a Python package to identify filament spines in 3D images (e.g., position-position-velocity cubes), along with a complementary technique to sort fitted velocity components into velocity-coherent filaments. In NGC 1333, we find a wealth of velocity gradient structures on a beam-resolved scale of ~0.05 pc. Interestingly, these local velocity gradients are not randomly oriented with respect to filament spines and their perpendicular, i.e., radial, component decreases in magnitude towards the spine for many filaments. Together with remarkably constant velocity gradients on larger scales along many filaments, these results suggest a scenario in which gas falling onto filaments is progressively damped and redirected to flow along these filaments., Accepted to ApJ

In the absence of the two emission lines H$\alpha$ and [NII] (6584\AA) in a BPT diagram, we show that other spectral information is sufficiently informative to distinguish AGN galaxies from star-forming galaxies. We use pattern recognition methods and a sample of galaxy spectra from the Sloan Digital Sky Survey (SDSS) to show that, in this survey, the flux and equivalent width of [OIII] (5007\AA) and H$\beta$, along with the 4000\AA-break, can be used to classify galaxies in a BPT diagram. This method provides a higher accuracy of predictions than those which use stellar mass and [OIII]/H$\beta$. First, we use BPT diagrams and various physical parameters to re-classify the galaxies. Next, using confusion matrices, we determine the `correctly' predicted classes as well as confused cases. In this way, we investigate the effect of each parameter in the confusion matrices and rank the physical parameters used in the discrimination of the different classes. We show that in this survey, for example, $\rm{g - r}$ colour can provide the same accuracy as galaxy stellar mass to predict whether or not a galaxy hosts an AGN. Finally, with the same information, we also rank the parameters involved in the discrimination of Seyfert and LINER galaxies., Comment: Accepted for publication in MNRAS. 12 pages, 14 figures

The standard model of prestellar core collapse suggests that this process works from the inside and moves outwards, with the fastest motions at the center. The relative abundances of many molecules also vary within cores, with certain molecules found only in specific regions characterized by narrow ranges of temperature and density. These characteristics lead to the hypothesis that the observed infall speeds in starless cores depend on both the position of the observations and the molecular tracer chosen. By measuring line emission at multiple positions across a core using an array of tracer molecules, one can determine whether these theoretical dependencies match observational evidence. Although surveys of infall motions in dense cores have been carried out for years, very few surveys have been awarded enough time to map infall across cores using multiple spectral line observations. To fill this gap, we present IRAM 30m maps of N2H (1-0), DCO(2-1), DCO(3-2) and HCO(3-2) emission towards two prestellar cores (L1544 and L694) and one protostellar core (L1521F). We find that the measured infall velocity varies as a function of position across each core and varies with the choice of molecular line, likely as a result of radial variations in core chemistry and dynamics. Subject headings: astronomy, astrophysics, star formation, dense cores, molecular clouds

We use Green Bank Ammonia Survey observations of NH$_3$ (1,1) and (2,2) emission with 32'' FWHM resolution from a ~ 10 pc$^{2}$ portion of the Cepheus-L1251 molecular cloud to identify hierarchical dense gas structures. Our dendrogram analysis of the NH$_3$ data results in 22 top-level structures, which reside within 13 lower-level, parent structures. The structures are compact (0.01 pc $\lesssim R_{eff} \lesssim$ 0.1 pc) and are spatially correlated with the highest H$_2$ column density portions of the cloud. We also compare the ammonia data to a catalog of dense cores identified by higher-resolution (18.2'' FWHM) Herschel Space Observatory observations of dust continuum emission from Cepheus-L1251. Maps of kinetic gas temperature, velocity dispersion, and NH$_3$ column density, derived from detailed modeling of the NH$_3$ data, are used to investigate the stability and chemistry of the ammonia-identified and Herschel-identified structures. We show that the dust and dense gas in the structures have similar temperatures, with median $T_{dust}$ and $T_K$ measurements of 11.7 $\pm$ 1.1 K and 10.3 $\pm$ 2.0 K, respectively. Based on a virial analysis, we find that the ammonia-identified structures are gravitationally dominated, yet may be in or near a state of virial equilibrium. Meanwhile, the majority of the Herschel-identified dense cores appear to be not bound by their own gravity and instead confined by external pressure. CCS $(2_0-1_0)$ and HC$_5$N $(9-8)$ emission from the region reveal broader line widths and centroid velocity offsets when compared to the NH$_3$ (1,1) emission in some cases, likely due to these carbon-based molecules tracing the turbulent outer layers of the dense cores., Accepted for publication in ApJ

We present observations of the Cepheus Flare obtained as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Legacy Survey (GBLS) with the SCUBA-2 instrument. We produce a catalogue of sources found by SCUBA-2, and separate these into starless cores and protostars. We determine masses and densities for each of our sources, using source temperatures determined by the Herschel Gould Belt Survey. We compare the properties of starless cores in four different molecular clouds: L1147/58, L1172/74, L1251 and L1228. We find that the core mass functions for each region typically show shallower-than-Salpeter behaviour. We find that L1147/58 and L1228 have a high ratio of starless cores to Class II protostars, while L1251 and L1174 have a low ratio, consistent with the latter regions being more active sites of current star formation, while the former are forming stars less actively. We determine that if modelled as thermally supported Bonnor–Ebert spheres, most of our cores have stable configurations accessible to them. We estimate the external pressures on our cores using archival 13CO velocity dispersion measurements and find that our cores are typically pressure confined, rather than gravitationally bound. We perform a virial analysis on our cores, and find that they typically cannot be supported against collapse by internal thermal energy alone, due primarily to the measured external pressures. This suggests that the dominant mode of internal support in starless cores in the Cepheus Flare is either non-thermal motions or internal magnetic fields., Monthly Notices of the Royal Astronomical Society, 464 (4), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

Understanding the early stages of star formation is a research field of ongoing development, both theoretically and observationally. In this context, molecular data have been continuously providing observational constraints on the gas dynamics at different excitation conditions and depths in the sources. We have investigated the Barnard 59 core, the only active site of star formation in the Pipe Nebula, to achieve a comprehensive view of the kinematic properties of the source. These information were derived by simultaneously fitting ammonia inversion transition lines (1,1) and (2,2). Our analysis unveils the imprint of protostellar feedback, such as increasing line widths, temperature and turbulent motions in our molecular data. Combined with complementary observations of dust thermal emission, we estimate that the core is gravitationally bound following a virial analysis. If the core is not contracting, another source of internal pressure, most likely the magnetic field, is supporting it against gravitational collapse and limits its star formation efficiency., 18 pages, 18 figures

Although surveys of infall motions in dense cores have been carried out for years, few surveys have focused on mapping infall across cores using multiple spectral line observations. To fill this gap, we present IRAM 30-m Telescope maps of N2H+(1-0), DCO+(2-1), DCO+(3-2), and HCO+(3-2) emission towards two prestellar cores (L492 and L694-2) and one protostellar core (L1521F). We find that the measured infall velocity varies with position across each core and choice of molecular line, likely as a result of radial variations in core chemistry and dynamics. Line-of-sight infall speeds estimated from DCO+(2-1) line profiles can decrease by 40-50 m/s when observing at a radial offset >= 0.04 pc from the core's dust continuum emission peak. Median infall speeds calculated from all observed positions across a core can also vary by as much as 65 m/s depending on the transition. These results show that while single-pointing, single-transition surveys of core infall velocities may be good indicators of whether a core is either contracting or expanding, the magnitude of the velocities they measure are significantly impacted by the choice of molecular line, proximity to the core center, and core evolutionary state., Accepted for publication in ApJ

We present deep radio continuum observations of the star-forming core of the Serpens South Infrared Dark Cloud with the Karl G. Jansky Very Large Array (VLA). Observations were conducted in two bands centered at 7.25 GHz (4.14 cm) and 4.75 GHz (6.31 cm) with an rms of 8.5 and 11.1 microJy/beam, respectively. We also use 2MASS, Spitzer and Herschel data to put our radio observations in the context of young stellar populations characterized by near and far infrared observations. Within a 5 arcmin x 5 arcmin region of interest around the central cluster, we detect roughly eighteen radio sources, seven of which we determine are protostellar in nature due to their radio spectral indices and their association with infrared sources. We find evidence for a previously undetected embedded Class 0 protostar and reaffirm Class 0 protostellar classifications determined by previous millimeter wavelength continuum studies. We use our infrared data to derive mid-infrared luminosities for three of our protostellar sources and find relative agreement between the known YSO radio luminosity vs bolometric luminosity correlation. Lastly, we marginally detect an additional six radio sources at the 2-3 sigma level that lie within two arcseconds of infrared YSO candidates, providing motivation for higher sensitivity studies to clarify the nature of these sources and further probe embedded and/or low luminosity YSOs in Serpens South.