Search

The internal velocity structure within dense gaseous cores plays a crucial role in providing the initial conditions for star formation in molecular clouds. However, the kinematic properties of dense gas at core scales (~0.01 - 0.1 pc) has not been extensively characterized because of instrument limitations until the unique capabilities of GBT-Argus became available. The ongoing GBT-Argus Large Program, Dynamics in Star-forming Cores (DiSCo) thus aims to investigate the origin and distribution of angular momenta of star-forming cores. DiSCo will survey all starless cores and Class 0 protostellar cores in the Perseus molecular complex down to ~0.01 pc scales with < 0.05 km/s velocity resolution using the dense gas tracer N$_2$H$^+$. Here, we present the first datasets from DiSCo toward the B1 and NGC 1333 regions in Perseus. Our results suggest that a dense core's internal velocity structure has little correlation with other core-scale properties, indicating these gas motions may be originated externally from cloud-scale turbulence. These first datasets also reaffirm the ability of GBT-Argus for studying dense core velocity structure and provided an empirical basis for future studies that address the angular momentum problem with a statistically broad sample., Comment: 17 pages, 12 figures, accepted by MNRAS

We report observations of the ground state transitions of $^{12}$CO, $^{13}$CO C$^{18}$O, HCN, and HCO$^+$ at 88-115 GHz in the inner region of the nearby galaxy IC 342. These data were obtained with the 16-pixel spectroscopic focal plane array Argus on the 100-m Robert C. Byrd Green Bank Telescope (GBT) at 6-9$^{\prime\prime}$ resolution. In the nuclear bar region, the intensity distributions of $^{12}$CO(1-0) and $^{13}$CO(1-0) emission trace moderate densities, and differ from the dense gas distributions sampled in C$^{18}$O(1-0), HCN(1-0), and HCO$^+$(1-0). We observe a constant HCN(1-0)-to-HCO$^+$(1-0) ratio of 1.2$\pm$0.1 across the whole $\sim$1 kpc bar. This indicates that HCN(1-0) and HCO$^+$(1-0) lines have intermediate optical depth, and that the corresponding $n_{\textrm{H}_2}$ of the gas producing the emission is of 10$^{4.5-6}$ cm$^{-3}$. We show that HCO$^+$(1-0) is thermalized and HCN(1-0) is close to thermalization. The very tight correlation between HCN(1-0) and HCO$^+$(1-0) intensities across the 1~kpc bar suggests that this ratio is more sensitive to the relative abundance of the two species than to the gas density. We confirm the angular offset ($\sim$10$^{\prime\prime}$) between the spatial distribution of molecular gas and the star formation sites. Finally, we find a breakdown of the $L_\textrm{IR}$-$L_\textrm{HCN}$ correlation at high spatial resolution due to the effect of incomplete sampling of star-forming regions by HCN emission in IC 342. The scatter of the $L_\textrm{IR}$-$L_\textrm{HCN}$ relation decreases as the spatial scale increases from 10$^{\prime\prime}$~to 30$^{\prime\prime}$ (170-510~pc), and is comparable to the scatter of the global relation at the scale of 340 pc., Comment: Accepted for publication in ApJ. 18 pages, 10 figures, and 4 tables. Comments are more than welcome!

We present a new upper limit on the cosmic molecular gas density at $z=2.4-3.4$ obtained using the first year of observations from the CO Mapping Array Project (COMAP). COMAP data cubes are stacked on the 3D positions of 243 quasars selected from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) catalog, yielding a 95% upper limit for flux from CO(1-0) line emission of 0.129 Jy km/s. Depending on the balance of the emission between the quasar host and its environment, this value can be interpreted as an average CO line luminosity $L'_\mathrm{CO}$ of eBOSS quasars of $\leq 1.26\times10^{11}$ K km pc$^2$ s$^{-1}$, or an average molecular gas density $\rho_\mathrm{H_2}$ in regions of the universe containing a quasar of $\leq 1.52\times10^8$ M$_\odot$ cMpc$^{-3}$. The $L'_\mathrm{CO}$ upper limit falls among CO line luminosities obtained from individually-targeted quasars in the COMAP redshift range, and the $\rho_\mathrm{H_2}$ value is comparable to upper limits obtained from other Line Intensity Mapping (LIM) surveys and their joint analyses. Further, we forecast the values obtainable with the COMAP/eBOSS stack after the full 5-year COMAP Pathfinder survey. We predict that a detection is probable with this method, depending on the CO properties of the quasar sample. Based on the achieved sensitivity, we believe that this technique of stacking LIM data on the positions of traditional galaxy or quasar catalogs is extremely promising, both as a technique for investigating large galaxy catalogs efficiently at high redshift and as a technique for bolstering the sensitivity of LIM experiments, even with a fraction of their total expected survey data., Comment: 26 pages, 22 figures; Version 2 is as accepted by ApJ

We report on results from GBT/ARGUS $^{12}$CO(1-0) observations for the giant low surface brightness galaxy Malin 1, which allow us to determine an upper limit for its CO mass, and hence its molecular gas mass and molecular gas mass surface density $\Sigma_{H_2}$. Although we performed very deep observations through 17 hours on source integration time, reaching a noise level of $\sim 0.2$ mK (T$^{*}_{A}$) with a corresponding extended source CO limit (3$\sigma$) of 0.09 K km s$^{-1}$, 19 times more sensitive than previous works, we do not detect the $^{12}$CO(1-0) emission line. However, the observations allow us to estimate an upper limit (3$\sigma$) for the CO mass of about $7.4 \times 10^9$ M$_\odot$ for the extended emission, and $1.4 \times 10^8$ M$_\odot$ for the central part of the galaxy. With these figures we conclude that the molecular gas surface density is lower than 0.3 M$_\odot$ pc$^{-2}$, and the corresponding molecular to atomic gas mass ratio is lower than 0.13. The evidence suggests a quite different physical conditions for the interstellar medium in Malin 1 compared to that of normal, high surface brightness spirals. This, in one way to another, keeps an usual molecular gas tracer as CO hidden from our observations, in spite of the diverse stellar and structural properties of Malin 1 observed by several authors since more than 30 years., Comment: This is the final (proof-corrected) version, published in ApJ Letters

We introduce COMAP-EoR, the next generation of the Carbon Monoxide Mapping Array Project aimed at extending CO intensity mapping to the Epoch of Reionization. COMAP-EoR supplements the existing 30 GHz COMAP Pathfinder with two additional 30 GHz instruments and a new 16 GHz receiver. This combination of frequencies will be able to simultaneously map CO(1--0) and CO(2--1) at reionization redshifts ($z\sim5-8$) in addition to providing a significant boost to the $z\sim3$ sensitivity of the Pathfinder. We examine a set of existing models of the EoR CO signal, and find power spectra spanning several orders of magnitude, highlighting our extreme ignorance about this period of cosmic history and the value of the COMAP-EoR measurement. We carry out the most detailed forecast to date of an intensity mapping cross-correlation, and find that five out of the six models we consider yield signal to noise ratios (S/N) $\gtrsim20$ for COMAP-EoR, with the brightest reaching a S/N above 400. We show that, for these models, COMAP-EoR can make a detailed measurement of the cosmic molecular gas history from $z\sim2-8$, as well as probe the population of faint, star-forming galaxies predicted by these models to be undetectable by traditional surveys. We show that, for the single model that does not predict numerous faint emitters, a COMAP-EoR-type measurement is required to rule out their existence. We briefly explore prospects for a third-generation Expanded Reionization Array (COMAP-ERA) capable of detecting the faintest models and characterizing the brightest signals in extreme detail., Comment: Paper 7 of 7 in series. 19 pages, 10 figures, to be submitted to ApJ

We present the current state of models for the $z\sim3$ carbon monoxide (CO) line-intensity signal targeted by the CO Mapping Array Project (COMAP) Pathfinder in the context of its early science results. Our fiducial model, relating dark matter halo properties to CO luminosities, informs parameter priors with empirical models of the galaxy-halo connection and previous CO(1-0) observations. The Pathfinder early science data spanning wavenumbers $k=0.051$-$0.62\,$Mpc$^{-1}$ represent the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum. Our 95% upper limit on the redshift-space clustering amplitude $A_{\rm clust}\lesssim70\,\mu$K$^2$ greatly improves on the indirect upper limit of $420\,\mu$K$^2$ reported from the CO Power Spectrum Survey (COPSS) measurement at $k\sim1\,$Mpc$^{-1}$. The COMAP limit excludes a subset of models from previous literature, and constrains interpretation of the COPSS results, demonstrating the complementary nature of COMAP and interferometric CO surveys. Using line bias expectations from our priors, we also constrain the squared mean line intensity-bias product, $\langle{Tb}\rangle^2\lesssim50\,\mu$K$^2$, and the cosmic molecular gas density, $\rho_\text{H2}<2.5\times10^8\,M_\odot\,$Mpc$^{-3}$ (95% upper limits). Based on early instrument performance and our current CO signal estimates, we forecast that the five-year Pathfinder campaign will detect the CO power spectrum with overall signal-to-noise of 9-17. Between then and now, we also expect to detect the CO-galaxy cross-spectrum using overlapping galaxy survey data, enabling enhanced inferences of cosmic star-formation and galaxy-evolution history., Comment: Paper 5 of 7 in series. 17 pages + appendix and bibliography (30 pages total); 15 figures, 6 tables; accepted for publication in ApJ; v3 reflects the accepted version with minor changes and additions to text

We present early results from the COMAP Galactic Plane Survey conducted between June 2019 and April 2021, spanning $20^\circ<\ell<40^\circ$ in Galactic longitude and $|b|<1.\!\!^{\circ}5$ in Galactic latitude with an angular resolution of $4.5^{\prime}$. The full survey will span $\ell \sim 20^{\circ}$- $220^{\circ}$ and will be the first large-scale radio continuum survey at $30$ GHz with sub-degree resolution. We present initial results from the first part of the survey, including diffuse emission and spectral energy distributions (SEDs) of HII regions and supernova remnants. Using low and high frequency surveys to constrain free-free and thermal dust emission contributions, we find evidence of excess flux density at $30\,$GHz in six regions that we interpret as anomalous microwave emission. Furthermore we model UCHII contributions using data from the $5\,$GHz CORNISH catalogue and reject this as the cause of the $30\,$GHz excess. Six known supernova remnants (SNR) are detected at $30\,$GHz, and we measure spectral indices consistent with the literature or show evidence of steepening. The flux density of the SNR W44 at $30\,$GHz is consistent with a power-law extrapolation from lower frequencies with no indication of spectral steepening in contrast with recent results from the Sardinia Radio Telescope. We also extract five hydrogen radio recombination lines to map the warm ionized gas, which can be used to estimate electron temperatures or to constrain continuum free-free emission. The full COMAP Galactic plane survey, to be released in 2023/2024, will be an invaluable resource for Galactic astrophysics., Comment: Paper 6 of 7 in series. 28 pages, 10 figures, submitted to ApJ

We present the power spectrum methodology used for the first-season COMAP analysis, and assess the quality of the current data set. The main results are derived through the Feed-feed Pseudo-Cross-Spectrum (FPXS) method, which is a robust estimator with respect to both noise modeling errors and experimental systematics. We use effective transfer functions to take into account the effects of instrumental beam smoothing and various filter operations applied during the low-level data processing. The power spectra estimated in this way have allowed us to identify a systematic error associated with one of our two scanning strategies, believed to be due to residual ground or atmospheric contamination. We omit these data from our analysis and no longer use this scanning technique for observations. We present the power spectra from our first season of observing and demonstrate that the uncertainties are integrating as expected for uncorrelated noise, with any residual systematics suppressed to a level below the noise. Using the FPXS method, and combining data on scales $k=0.051-0.62 \,\mathrm{Mpc}^{-1}$ we estimate $P_\mathrm{CO}(k) = -2.7 \pm 1.7 \times 10^4\mu\textrm{K}^2\mathrm{Mpc}^3$, the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum in the literature., Comment: Paper 4 of 7 in series. 18 pages, 11 figures, as accepted in ApJ

We describe the first season COMAP analysis pipeline that converts raw detector readouts to calibrated sky maps. This pipeline implements four main steps: gain calibration, filtering, data selection, and map-making. Absolute gain calibration relies on a combination of instrumental and astrophysical sources, while relative gain calibration exploits real-time total-power variations. High efficiency filtering is achieved through spectroscopic common-mode rejection within and across receivers, resulting in nearly uncorrelated white noise within single-frequency channels. Consequently, near-optimal but biased maps are produced by binning the filtered time stream into pixelized maps; the corresponding signal bias transfer function is estimated through simulations. Data selection is performed automatically through a series of goodness-of-fit statistics, including $\chi^2$ and multi-scale correlation tests. Applying this pipeline to the first-season COMAP data, we produce a dataset with very low levels of correlated noise. We find that one of our two scanning strategies (the Lissajous type) is sensitive to residual instrumental systematics. As a result, we no longer use this type of scan and exclude data taken this way from our Season 1 power spectrum estimates. We perform a careful analysis of our data processing and observing efficiencies and take account of planned improvements to estimate our future performance. Power spectrum results derived from the first-season COMAP maps are presented and discussed in companion papers., Comment: Paper 3 of 7 in series. 26 pages, 23 figures, submitted to ApJ

The CO Mapping Array Project (COMAP) aims to use line intensity mapping of carbon monoxide (CO) to trace the distribution and global properties of galaxies over cosmic time, back to the Epoch of Reionization (EoR). To validate the technologies and techniques needed for this goal, a Pathfinder instrument has been constructed and fielded. Sensitive to CO(1-0) emission from $z=2.4$-$3.4$ and a fainter contribution from CO(2-1) at $z=6$-8, the Pathfinder is surveying $12$ deg$^2$ in a 5-year observing campaign to detect the CO signal from $z\sim3$. Using data from the first 13 months of observing, we estimate $P_\mathrm{CO}(k) = -2.7 \pm 1.7 \times 10^4\mu\mathrm{K}^2 \mathrm{Mpc}^3$ on scales $k=0.051-0.62 \mathrm{Mpc}^{-1}$ - the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum. Based on these observations alone, we obtain a constraint on the amplitude of the clustering component (the squared mean CO line temperature-bias product) of $\langle Tb\rangle^2<49$ $\mu$K$^2$ - nearly an order-of-magnitude improvement on the previous best measurement. These constraints allow us to rule out two models from the literature. We forecast a detection of the power spectrum after 5 years with signal-to-noise ratio (S/N) 9-17. Cross-correlation with an overlapping galaxy survey will yield a detection of the CO-galaxy power spectrum with S/N of 19. We are also conducting a 30 GHz survey of the Galactic plane and present a preliminary map. Looking to the future of COMAP, we examine the prospects for future phases of the experiment to detect and characterize the CO signal from the EoR., Comment: Paper 1 of 7 in series. 18 pages, 16 figures, submitted to ApJ

Line intensity mapping (LIM) is a new technique for tracing the global properties of galaxies over cosmic time. Detection of the very faint signals from redshifted carbon monoxide (CO), a tracer of star formation, pushes the limits of what is feasible with a total-power instrument. The CO Mapping Project (COMAP) Pathfinder is a first-generation instrument aiming to prove the concept and develop the technology for future experiments, as well as delivering early science products. With 19 receiver channels in a hexagonal focal plane arrangement on a 10.4 m antenna, and an instantaneous 26-34 GHz frequency range with 2 MHz resolution, it is ideally suited to measuring CO($J$=1-0) from $z\sim3$. In this paper we discuss strategies for designing and building the Pathfinder and the challenges that were encountered. The design of the instrument prioritized LIM requirements over those of ancillary science. After a couple of years of operation, the instrument is well understood, and the first year of data is already yielding useful science results. Experience with this Pathfinder will drive the design of the next generations of experiments., Comment: Paper 2 of 7 in series. 27 pages, 28 figures, submitted to ApJ

Dense cores are the final place where turbulence is dissipated. It has been proposed from theoretical arguments that the non-thermal velocity dispersion should be narrower both for molecular ions (compared to neutrals) and for transitions with higher critical densities. To test these hypotheses, we compare the velocity dispersion of N$_2$H$^+$ (1--0) (n$_{\rm crit}$ = $6\times10^4$ cm$^{-3}) and NH$_3$ (n$_{\rm crit}=2\times10^3$ cm$^{-3}), in the dense core Barnard 5. We analyse well resolved and high signal-to-noise observations of NH$_3$ (1,1) and (2,2) obtained with combining GBT and VLA data, and N$_2$H$^+$ (1--0) obtained with GBT Argus, which present a similar morphology. % Surprisingly, the non-thermal velocity dispersion of the ion is systematically higher than that of the neutral by 20\%. The derived sonic Mach number, $\mathcal{M}_s = \sigma_{\rm NT}/c_s$, has peak values $\mathcal{M}_{s, {\rm N_2H^+}} = 0.59$ and $\mathcal{M}_{s, {\rm NH}_3} = 0.48$ for N$_2$H$^+$ and NH$_3$, respectively. % This observed difference may indicate that the magnetic field even deep within the dense core is still oscillating, as it is in the turbulent region outside the core. The ions should be more strongly dynamically coupled to this oscillating field than the neutrals, thus accounting for their broader linewidth. If corroborated by further observations, this finding would shed additional light on the transition to quiescence in dense cores., Comment: Accepted for publication in ApJ. 14 pages, 11 Figures

We present the first results of high-spectral resolution (0.023 km/s) N$_2$H$^+$ observations of dense gas dynamics at core scales (~0.01 pc) using the recently commissioned Argus instrument on the Green Bank Telescope (GBT). While the fitted linear velocity gradients across the cores measured in our targets nicely agree with the well-known power-law correlation between the specific angular momentum and core size, it is unclear if the observed gradients represent core-scale rotation. In addition, our Argus data reveal detailed and intriguing gas structures in position-velocity (PV) space for all 5 targets studied in this project, which could suggest that the velocity gradients previously observed in many dense cores actually originate from large-scale turbulence or convergent flow compression instead of rigid-body rotation. We also note that there are targets in this study with their star-forming disks nearly perpendicular to the local velocity gradients, which, assuming the velocity gradient represents the direction of rotation, is opposite to what is described by the classical theory of star formation. This provides important insight on the transport of angular momentum within star-forming cores, which is a critical topic on studying protostellar disk formation., Comment: 14 pages, 10 figures, accepted by MNRAS

We combine recent simulation work on the SFR--[C II] correlation at high redshift with empirical modeling of the galaxy--halo connection (via UniverseMachine) to forecast [C II] auto power spectra from $z\sim4$ to $z\sim8$. We compare these to sensitivities realistically expected from various instruments expected to come on-line in the next decade. If the predictions of our model are correct, [C II] should be detectable up to $z\sim6$ in this generation of surveys, but detecting [C II] past the end of reionization will require a generational leap in line-intensity survey capabilities., Comment: 11 pages + appendix (16 total), 3 figures, 5 tables; accepted for publication in ApJ

We present deep CCS and HC$_7$N observations of the L1495-B218 filaments in the Taurus molecular cloud obtained using the K-band focal plane array on the 100m Green Bank Telescope. We observed the L1495-B218 filaments in CCS $J_N$ = 2$_1$$-$1$_0$ and HC$_7$N $J$ = 21$-$20 with a spectral resolution of 0.038 km s$^{-1}$ and an angular resolution of 31$''$. We observed strong CCS emission in both evolved and young regions and weak emission in two evolved regions. HC$_7$N emission is observed only in L1495A-N and L1521D. We find that CCS and HC$_7$N intensity peaks do not coincide with NH$_3$ or dust continuum intensity peaks. We also find that the fractional abundance of CCS does not show a clear correlation with the dynamical evolutionary stage of dense cores. Our findings and chemical modeling indicate that the fractional abundances of CCS and HC$_7$N are sensitive to the initial gas-phase C/O ratio, and they are good tracers of young condensed gas only when the initial C/O is close to solar value. Kinematic analysis using multiple lines including NH$_3$, HC$_7$N, CCS, CO, HCN, \& HCO$^+$ suggests that there may be three different star formation modes in the L1495-B218 filaments. At the hub of the filaments, L1495A/B7N has formed a stellar cluster with large-scale inward flows (fast mode), while L1521D, a core embedded in a filament, is slowly contracting due to its self-gravity (slow mode). There is also one isolated core that appears to be marginally stable and may undergo quasi-static evolution (isolated mode)., Comment: ApJ accepted

Line-intensity mapping (LIM or IM) is an emerging field of observational work, with strong potential to fit into a larger effort to probe large-scale structure and small-scale astrophysical phenomena using multiple complementary tracers. Taking full advantage of such complementarity means, in part, undertaking line-intensity surveys with galaxy surveys in mind. We consider the potential for detection of a cross-correlation signal between COMAP and blind surveys based on photometric redshifts (as in COSMOS) or based on spectroscopic data (as with the HETDEX survey of Lyman-$\alpha$ emitters). We find that obtaining $\sigma_z/(1+z)\lesssim0.003$ accuracy in redshifts and $\gtrsim10^{-4}$ sources per Mpc$^3$ with spectroscopic redshift determination should enable a CO-galaxy cross spectrum detection significance at least twice that of the CO auto spectrum. Either a future targeted spectroscopic survey or a blind survey like HETDEX may be able to meet both of these requirements., Comment: 19 pages + appendix (31 pages total), 16 figures, 6 tables; accepted for publication in ApJ

Line-intensity mapping surveys probe large-scale structure through spatial variations in molecular line emission from a population of unresolved cosmological sources. Future such surveys of carbon monoxide line emission, specifically the CO(1-0) line, face potential contamination from a disjoint population of sources emitting in a hydrogen cyanide emission line, HCN(1-0). This paper explores the potential range of the strength of HCN emission and its effect on the CO auto power spectrum, using simulations with an empirical model of the CO/HCN--halo connection. We find that effects on the observed CO power spectrum depend on modeling assumptions but are very small for our fiducial model based on our understanding of the galaxy--halo connection, with the bias in overall CO detection significance due to HCN expected to be less than 1%., Comment: 13 pages + appendix (16 pages total), 8 figures, 5 tables; revised in response to referee, switched to AASTeX 6.0; accepted for publication in ApJ

Intensity mapping, which images a single spectral line from unresolved galaxies across cosmological volumes, is a promising technique for probing the early universe. Here we present predictions for the intensity map and power spectrum of the CO(1-0) line from galaxies at z ~ 2.4 -2.8, based on a parameterized model for the galaxy-halo connection, and demonstrate the extent to which properties of high-redshift galaxies can be directly inferred from such observations. We find that our fiducial prediction should be detectable by a realistic experiment. Motivated by significant modeling uncertainties, we demonstrate the effect on the power spectrum of varying each parameter in our model. Using simulated observations, we infer constraints on our model parameter space with an MCMC procedure, and show corresponding constraints on the LIR-LCO relation and the CO luminosity function. These constraints would be complementary to current high-redshift galaxy observations, which can detect the brightest galaxies but not complete samples from the faint end of the luminosity function. By probing these populations in aggregate, CO intensity mapping could be a valuable tool for probing molecular gas and its relation to star formation in high-redshift galaxies., Comment: Updated to match ApJ accepted version, with corrections to figure/table formatting (supersedes the published version, wherever they differ). 14 pages + appendix (22 pages total), 12 figures, 3 tables