Source author record

James Di Francesco

James Di Francesco 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

58works
6topics
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

58 published item(s)

preprint2026arXiv

How To Use Thermal Dust Continuum Emission To Measure The Physical Properties Of Dusty Astrophysical Objects

Dust grains in the interstellar medium interact with photons across the electromagnetic spectrum. They are generally photon energy converters, absorbing short wavelength radiation and emitting long wavelength radiation. Sixty years ago in 1965, thermal emission from dust grains in the interstellar medium was discovered. This tutorial is a summary of the physics of thermal dust continuum emission and how to use observations of the intensity and flux density of dusty objects to calculate physical properties such as mass, column density, luminosity, dust temperature, and dust opacity spectral index. Equations are derived, when feasible, from first principles with all limits and assumptions explicitly stated. Properties of dust opacities appropriate for different astrophysical environments (e.g. diffuse ISM, dense cores, protoplanetary disks) are discussed and tabulated for the wavelengths of past, current, and future bolometer cameras. Corrections for observations at high redshift as well as the effects of telescope measurement limitations are derived. We also update the calculation of the mean molecular weight in different ISM environments and find that it is 1.404 per H atom, 2.809 per H2 molecule, and 2.351 per gas particle assuming protosolar metallicity and the latest values of the ISM gas phase abundances of metals.

preprint2022arXiv

A Survey of Deuterated Ammonia in the Cepheus Star-Forming Region L1251

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 $σ_{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$σ$ 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.

preprint2022arXiv

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): How do dense core properties affect the multiplicity of protostars?

During the transition phase from a prestellar to a protostellar cloud core, one or several protostars can form within a single gas core. The detailed physical processes of this transition, however, still remain unclear. We present 1.3 mm dust continuum and molecular line observations with the Atacama Large Millimeter/submillimeter Array (ALMA) toward 43 protostellar cores in the Orion Molecular Cloud Complex ($λ$ Orionis, Orion B, and Orion A) with an angular resolution of $\sim$ 0.35" ($\sim$ 140 au). In total, we detect 13 binary/multiple systems. We derive an overall multiplicity frequency (MF) of 28$\%$ $\pm$ 4$\%$ and a companion star fraction (CSF) of 51$\%$ $\pm$ 6$\%$, over a separation range of 300-8900 au. The median separation of companions is about 2100 au. The occurrence of stellar multiplicity may depend on the physical characteristics of the dense cores. Notably, those containing binary/multiple systems tend to show higher gas density and Mach number than cores forming single stars. The integral-shaped filament (ISF) of Orion A giant molecular cloud (GMC), which has the highest gas density and hosts high-mass star formation in its central region (the Orion Nebula cluster), shows the highest MF and CSF among the Orion GMCs. In contrast, the $λ$ Orionis Giant Molecular Cloud (GMC) has a lower MF and CSF than the Orion B and Orion A GMCs, indicating that feedback from HII regions may suppress the formation of multiple systems. We also find that the protostars comprising a binary/multiple system are usually at different evolutionary stages.

preprint2022arXiv

B-fields in Star-Forming Region Observations (BISTRO): Magnetic Fields in the Filamentary Structures of Serpens Main

We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~μ$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities.

preprint2022arXiv

Turbulence and Accretion: a High-resolution Study of the B5 Filaments

High-resolution observations of the Perseus B5 "core" have previously revealed that this subsonic region actually consists of several filaments that are likely in the process of forming a quadruple stellar system. Since subsonic filaments are thought to be produced at the $\sim 0.1$ pc sonic scale by turbulent compression, a detailed kinematic study is crucial to test such a scenario in the context of core and star formation. Here we present a detailed kinematic follow-up study of the B5 filaments at a 0.009 pc resolution using the VLA and GBT combined observations fitted with multi-component spectral models. Using precisely identified filament spines, we find a remarkable resemblance between the averaged width profiles of each filament and Plummer-like functions, with filaments possessing FWHM widths of $\sim 0.03$ pc. The velocity dispersion profiles of the filaments also show decreasing trends towards the filament spines. Moreover, the velocity gradient field in B5 appears to be locally well ordered ($\sim 0.04$ pc) but globally complex, with kinematic behaviors suggestive of inhomogeneous turbulent accretion onto filaments and longitudinal flows towards a local overdensity along one of the filaments.

preprint2021arXiv

Detection of a dense SiO jet in the evolved protostellar phase

Jets and outflows trace the accretion history of protostars. High-velocity molecular jets have been observed from several protostars in the early Class\,0 phase of star formation, detected with the high-density tracer SiO. Until now, no clear jet has been detected with SiO emission from isolated evolved Class\,I protostellar systems. We report a prominent dense SiO jet from a Class\,I source G205S3 (HOPS\,315: T$_{bol}$ $\sim$ 180 K, spectral index $\sim$ 0.417), with a moderately high mass-loss rate ($\sim$ 0.59 $\times$ 10$^{-6}$ M$_\odot$ yr$^{-1}$) estimated from CO emission. Together, these features suggest that G205S3 is still in a high accretion phase, similar to that expected of Class\,0 objects. We compare G205S3 to a representative Class\,0 system G206W2 (HOPS\,399) and literature Class\,0/I sources to explore the possible explanations behind the SiO emission seen at the later phase. We estimate a high inclination angle ($\sim$ 40$^\circ$) for G205S3 from CO emission, which may expose the infrared emission from the central core and mislead the spectral classification. However, the compact 1.3\,mm continuum, C$^{18}$O emission, location in the bolometric luminosity to sub-millimeter fluxes diagram, outflow force ($\sim$ 3.26 $\times$ 10$^{-5}$ M$_\odot$km s$^{-1}$/yr) are also analogous to that of Class\,I systems. We thus consider G205S3 to be at the very early phase of Class\,I, and in the late phase of ``high-accretion". The episodic ejection could be due to the presence of an unknown binary, a planetary companion, or dense clumps, where the required mass for such high accretion could be supplied by a massive circumbinary disk.

preprint2021arXiv

Molecular Cloud Cores with High Deuterium Fractions: Nobeyama Mapping Survey

We present the results of on-the-fly mapping observations of 44 fields containing 107 SCUBA-2 cores in the emission lines of molecules, N$_2$H$^+$, HC$_3$N, and CCS at 82$-$94 GHz using the Nobeyama 45-m telescope. This study aimed at investigating the physical properties of cores that show high deuterium fractions and might be close to the onset of star formation. We found that the distributions of the N$_2$H$^+$ and HC$_3$N line emissions are approximately similar to that of 850-$μ$m dust continuum emission, whereas the CCS line emission is often undetected or is distributed in a clumpy structure surrounding the peak position of the 850-$μ$m dust continuum emission. Occasionally (12%), we observe the CCS emission which is an early-type gas tracer toward the young stellar object, probably due to local high excitation. Evolution toward star formation does not immediately affect nonthermal velocity dispersion.

preprint2020arXiv

ALMA ACA and Nobeyama observations of two Orion cores in deuterated molecular lines

We mapped two molecular cloud cores in the Orion A cloud with the ALMA ACA 7-m Array and with the Nobeyama 45-m radio telescope. These cores have bright N$_2$D$^+$ emission in single-pointing observations with the Nobeyama 45-m radio telescope, have relatively high deuterium fraction, and are thought to be close to the onset of star formation. One is a star-forming core, and the other is starless. These cores are located along filaments observed in N$_2$H$^+$, and show narrow linewidths of 0.41 km s$^{-1}$ and 0.45 km s$^{-1}$ in N$_2$D$^+$, respectively, with the Nobeyama 45-m telescope. Both cores were detected with the ALMA ACA 7m Array in the continuum and molecular lines at Band 6. The starless core G211 shows clumpy structure with several sub-cores, which in turn show chemical differences. Also, the sub-cores in G211 have internal motions that are almost purely thermal. The starless sub-core G211D, in particular, shows a hint of the inverse P Cygni profile, suggesting infall motion. The star-forming core G210 shows an interesting spatial feature of two N$_2$D$^+$ peaks of similar intensity and radial velocity located symmetrically with respect to the single dust continuum peak. One interpretation is that the two N$_2$D$^+$ peaks represent an edge-on pseudo-disk. The CO outflow lobes, however, are not directed perpendicular to the line connecting both N$_2$D$^+$ peaks.

preprint2020arXiv

JCMT POL-2 and BISTRO Survey observations of magnetic fields in the L1689 molecular cloud

We present 850$μ$m polarization observations of the L1689 molecular cloud, part of the nearby Ophiuchus molecular cloud complex, taken with the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT). We observe three regions of L1689: the clump L1689N which houses the IRAS 16293-2422 protostellar system, the starless clump SMM-16, and the starless core L1689B. We use the Davis-Chandrasekhar-Fermi method to estimate plane-of-sky field strengths of $366\pm 55$ $μ$G in L1689N, $284\pm 34$ $μ$G in SMM-16, and $72\pm 33$ $μ$G in L1689B, for our fiducial value of dust opacity. These values indicate that all three regions are likely to be magnetically trans-critical with sub-Alfvénic turbulence. In all three regions, the inferred mean magnetic field direction is approximately perpendicular to the local filament direction identified in $Herschel$ Space Telescope observations. The core-scale field morphologies for L1689N and L1689B are consistent with the cloud-scale field morphology measured by the $Planck$ Space Observatory, suggesting that material can flow freely from large to small scales for these sources. Based on these magnetic field measurements, we posit that accretion from the cloud onto L1689N and L1689B may be magnetically regulated. However, in SMM-16, the clump-scale field is nearly perpendicular to the field seen on cloud scales by $Planck$, suggesting that it may be unable to efficiently accrete further material from its surroundings.

preprint2020arXiv

Protoplanetary disk masses in NGC 2024: Evidence for two populations

Protoplanetary disks in dense, massive star-forming regions are strongly affected by their environment. How this environmental impact changes over time is an important constraint on disk evolution and external photoevaporation models. We characterize the dust emission from 179 disks in the core of the young (0.5 Myr) NGC 2024 cluster. By studying how the disk mass varies within the cluster, and comparing these disks to those in other regions, we aim to determine how external photoevaporation influences disk properties over time. Using the Atacama Large Millimeter/submillimeter Array, a 2.9' x 2.9' mosaic centered on NGC 2024 FIR 3 was observed at 225 GHz with a resolution of 0.25'', or ~100 AU. The imaged region contains 179 disks identified at IR wavelengths, seven new disk candidates, and several protostars. The overall detection rate of disks is $32 \pm 4\%$. Few of the disks are resolved, with the exception of a giant (R = 300 AU) transition disk. Serendipitously, we observe a millimeter flare from an X-ray bright young stellar object (YSO), and resolve continuum emission from a Class 0 YSO in the FIR 3 core. Two distinct disk populations are present: a more massive one in the east, along the dense molecular ridge hosting the FIR 1-5 YSOs, with a detection rate of $45 \pm 7\%$. In the western population, towards IRS 1, only $15 \pm 4\%$ of disks are detected. NGC 2024 hosts two distinct disk populations. Disks along the dense molecular ridge are young (0.2 - 0.5 Myr) and partly shielded from the far ultraviolet radiation of IRS 2b; their masses are similar to isolated 1 - 3 Myr old SFRs. The western population is older and at lower extinctions, and may be affected by external photoevaporation from both IRS 1 and IRS 2b. However, it is possible these disks had lower masses to begin with.

preprint2020arXiv

Relative Alignment between Dense Molecular Cores and Ambient Magnetic Field: The Synergy of Numerical Models and Observations

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.

preprint2020arXiv

The Herschel Orion Protostar Survey: Far-Infrared Photometry and Colors of Protostars and Their Variations across Orion A and B

The degree to which the properties of protostars are affected by environment remains an open question. To investigate this, we look at the Orion A and B molecular clouds, home to most of the protostars within 500 pc. At ~400 pc, Orion is close enough to distinguish individual protostars across a range of environments in terms of both the stellar and gas projected densities. As part of the Herschel Orion Protostar Survey (HOPS), we used the Photodetector Array Camera and Spectrometer (PACS) to map 108 partially overlapping square fields with edge lengths of 5 arcmin or 8 arcmin and measure the 70 micron and 160 micron flux densities of 338 protostars within them. In this paper we examine how these flux densities and their ratio depend on evolutionary state and environment within the Orion complex. We show that Class 0 protostars occupy a region of the 70 micron flux density versus 160 micron to 70 micron flux density ratio diagram that is distinct from their more evolved counterparts. We then present evidence that the Integral-Shaped Filament (ISF) and Orion B contain protostars with more massive envelopes than those in the more sparsely populated LDN 1641 region. This can be interpreted as evidence for increasing star formation rates in the ISF and Orion B or as a tendency for more massive envelopes to be inherited from denser birth environments. We also provide technical details about the map-making and photometric procedures used in the HOPS program.

preprint2020arXiv

The Herschel view of the dense core population in the Ophiuchus molecular cloud

Herschel observations of nearby clouds in the Gould Belt support a paradigm for low-mass star formation, starting with the generation of molecular filaments, followed by filament fragmentation, and the concentration of mass into self-gravitating prestellar cores. With the unique far-infrared and submillimeter continuum imaging capabilities of the Herschel Space observatory, the closeby (d = 139 pc) Ophiuchus cloud was mapped at five wavelengths from 70 $μm$ to 500 $μm$ with the aim of providing a complete census of dense cores in this region, including unbound starless cores, bound prestellar cores, and protostellar cores. Taking advantage of the high dynamic range and multi-wavelength nature of the Herschel data, we used the multi-scale decomposition algorithms getsources and getfilaments to identify a complete sample of dense cores and filaments in the cloud and study their properties. The densest clouds of the Ophiuchus complex, L1688 and L1689, which thus far are only indirectly described as filamentary regions owing to the spatial distribution of their young stellar objects (YSOs), are confirmed to be dominated by filamentary structures. The tight correlation observed between prestellar cores and filamentary structures in L1688 and L1689 supports the view that solar-type star formation occurs primarily in dense filaments. While the sub clouds of the complex show disparities, L1689 being less efficient than L1688 at forming stars when considering their total mass budgets, both sub clouds share almost the same prestellar core formation efficiency in dense molecular gas. We also find evidence in the Herschel data for a remarkable concentric geometrical configuration in L1688 which is dominated by up to three arc-like compression fronts and presumably created by shockwave events emanating from the Sco OB2 association, including the neighboring massive (O9V) star $σ$ Sco.

preprint2020arXiv

Ubiquitous $\rm NH_3$ supersonic component in L1688 coherent cores

Context : Star formation takes place in cold dense cores in molecular clouds. Earlier observations have found that dense cores exhibit subsonic non-thermal velocity dispersions. In contrast, CO observations show that the ambient large-scale cloud is warmer and has supersonic velocity dispersions. Aims : We aim to study the ammonia ($\rm NH_3$) molecular line profiles with exquisite sensitivity towards the coherent cores in L1688 in order to study their kinematical properties in unprecedented detail. Methods : We used $\rm NH_3$ (1,1) and (2,2) data from the first data release (DR1) in the Green Bank Ammonia Survey (GAS). We first smoothed the data to a larger beam of 1' to obtain substantially more extended maps of velocity dispersion and kinetic temperature, compared to the DR1 maps. We then identified the coherent cores in the cloud and analysed the averaged line profiles towards the cores. Results : For the first time, we detected a faint (mean $\rm NH_3$(1,1) peak brightness $<$0.25 K in $T_{MB}$), supersonic component towards all the coherent cores in L1688. We fitted two components, one broad and one narrow, and derived the kinetic temperature and velocity dispersion of each component. The broad components towards all cores have supersonic linewidths ($\mathcal{M}_S \ge 1$). This component biases the estimate of the narrow dense core component's velocity dispersion by $\approx$28% and the kinetic temperature by $\approx$10%, on average, as compared to the results from single-component fits. Conclusions : Neglecting this ubiquitous presence of a broad component towards all coherent cores causes the typical single-component fit to overestimate the temperature and velocity dispersion. This affects the derived detailed physical structure and stability of the cores estimated from $\rm NH_3$ observations.

preprint2020arXiv

Velocity-coherent Filaments in NGC 1333: Evidence for Accretion Flow?

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.

preprint2019arXiv

Canada and the SKA from 2020-2030

This white paper submitted for the 2020 Canadian Long-Range Planning process (LRP2020) presents the prospects for Canada and the Square Kilometre Array (SKA) from 2020-2030, focussing on the first phase of the project (SKA1) scheduled to begin construction early in the next decade. SKA1 will make transformational advances in our understanding of the Universe across a wide range of fields, and Canadians are poised to play leadership roles in several. Canadian key SKA technologies will ensure a good return on capital investment in addition to strong scientific returns, positioning Canadian astronomy for future opportunities well beyond 2030. We therefore advocate for Canada's continued scientific and technological engagement in the SKA from 2020-2030 through participation in the construction and operations phases of SKA1.

preprint2019arXiv

CLOVER: Convnet Line-fitting Of Velocities in Emission-line Regions

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.

preprint2019arXiv

Development Plans for the Atacama Large Millimeter/submillimeter Array (ALMA)

(abridged) The Atacama Large Millimeter/submillimeter Array (ALMA) was the top-ranked priority for a new ground-based facility in the 2000 Canadian Long Range Plan. Ten years later, at the time of LRP2010, ALMA construction was well underway, with first science observations anticipated for 2011. In the past 8 years, ALMA has proved itself to be a high-impact, high-demand observatory, with record numbers of proposals submitted to the annual calls and large numbers of highly cited scientific papers across fields from protoplanetary disks to high-redshift galaxies and quasars. The LRP2010 ALMA white paper laid out 8 specific metrics that could be used to judge the success of Canada's participation in ALMA. Among these metrics were publications (number; impact), collaborations (international; multi-wavelength), and student training. To call out one particular metric, Canadians are making excellent use of ALMA in training graduate students and postdocs: as of June 2018, 12 of 23 Canadian first-author papers were led by a graduate student, and a further 4 papers were led by postdocs. All 8 metrics argue for Canada's involvement in ALMA over the past decade to be judged a success. The successful achievement of these wide-ranging goals argues strongly for Canada's continuing participation in ALMA over the next decade and beyond. Looking forward, our community needs to: (1) maintain Canadian access to ALMA and our competitiveness in using ALMA; (2) preserve full Canadian funding for our share of ALMA operations; (3) identify components of ALMA development in which Canada can play a significant role, including stimulating expertise in submillimetre instrumentation to capitalize on future opportunities; and (4) keep Canadians fully trained and engaged in ALMA, as new capabilities become available, reaching the widest possible community of potential users.

preprint2019arXiv

The Formation of Stars -- From Filaments to Cores to Protostars and Protoplanetry Disks

Star formation involves the flow of gas and dust within molecular clouds into protostars and young stellar objects (YSOs) due to gravity. Along the way, these flows are shaped significantly by many other mechanisms, including pressure, turbulent motions, magnetic fields, stellar feedback, jets, and angular momentum. How all these mechanisms interact nonlinearly with each other on various length scales leads to the formation and evolution of substructures within clouds, including filaments, clumps, cores, disks, outflows, the protostars/YSOs themselves, and planets. In this white paper, prepared for the 2020 Long Range Plan panel which will recommend Canada's future directions for astronomy, we describe the observational and theoretical leadership in the star formation field that Canada's vibrant community has demonstrated over the past decade. Drawing from this extensive background, we identify five key questions that must be addressed for further progress to be made in understanding star formation in the next decade. Addressing these questions will improve our understanding of the dynamics of the dense gas and the role of the magnetic field in star formation, the optical properties of the dust used to trace mass and magnetic fields, the sources of variability in star-forming objects on short timescales, and the physical processes that specifically promote the clustering of stars. We further highlight key facilities in which Canada should become involved to continue making progress in this field. Single-dish facilities we recommend include LSST, trans-atmospheric far-infrared telescopes like BLAST-TNG and SPICA, and ground-based telescopes like JCMT, GBT, and CCAT-p. Interferometric facilities we recommend include ALMA, ngVLA, and SKA1.

preprint2019arXiv

The Next Generation Very Large Array

The next generation Very Large Array (ngVLA) is a transformational radio observatory being designed by the U.S. National Radio Astronomy Observatory (NRAO). It will provide order of magnitude improvements in sensitivity, resolution, and uv coverage over the current Jansky Very Large Array (VLA) at ~1.2-50 GHz and extend the frequency range up to 70-115 GHz. This document is a white paper written by members of the Canadian community for the 2020 Long Range Plan panel, which will be making recommendations on Canada's future directions in astronomy. Since Canadians have been historically major users of the VLA and have been valued partners with NRAO for ALMA, Canada's participation in ngVLA is welcome. Canadians have been actually involved in ngVLA discussions for the past five years, and have played leadership roles in the ngVLA Science and Technical Advisory Councils. Canadian technologies are also very attractive for the ngVLA, in particular our designs for radio antennas, receivers, correlates, and data archives, and our industrial capacities to realize them. Indeed, the Canadian designs for the ngVLA antennas and correlator/beamformer are presently the baseline models for the project. Given the size of Canada's radio community and earlier use of the VLA (and ALMA), we recommend Canadian participation in the ngVLA at the 7% level. Such participation would be significant enough to allow Canadian leadership in gVLA's construction and usage. Canada's participation in ngVLA should not preclude its participation in SKA; access to both facilities is necessary to meet Canada's radio astronomy needs. Indeed, ngVLA will fill the gap between those radio frequencies observable with the SKA and ALMA at high sensitivities and resolutions. Canada's partnership in ngVLA will give it access to cutting-edge facilities together covering approximately three orders of magnitude in frequency.

preprint2019arXiv

The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars I. Identifying and Characterizing the Protostellar Content of the OMC2-FIR4 and OMC2-FIR3 Regions

We present ALMA (0.87~mm) and VLA (9~mm) observations toward OMC2-FIR4 and OMC2-FIR3 within the Orion integral-shaped filament that are thought to be the nearest regions of intermediate mass star formation. We characterize the continuum sources within these regions on $\sim$40~AU (0\farcs1) scales and associated molecular line emission at a factor of $\sim$30 better resolution than previous observations at similar wavelengths. We identify six compact continuum sources within OMC2-FIR4, four in OMC2-FIR3, and one additional source just outside OMC2-FIR4. This continuum emission is tracing the inner envelope and/or disk emission on less than 100~AU scales. HOPS-108 is the only protostar in OMC2-FIR4 that exhibits emission from high-excitation transitions of complex organic molecules (e.g., methanol and other lines) coincident with the continuum emission. HOPS-370 in OMC2-FIR3 with L~$\sim$~360~\lsun, also exhibits emission from high-excitation methanol and other lines. The methanol emission toward these two protostars is indicative of temperatures high enough to thermally evaporate methanol from icy dust grains; overall these protostars have characteristics similar to hot corinos. We do not identify a clear outflow from HOPS-108 in \twco, but find evidence of interaction between the outflow/jet from HOPS-370 and the OMC2-FIR4 region. The multitude of observational constraints indicate that HOPS-108 is likely a low to intermediate-mass protostar in its main mass accretion phase and it is the most luminous protostar in OMC2-FIR4. The high resolution data presented here are essential for disentangling the embedded protostars from their surrounding dusty environments and characterizing them.

preprint2016arXiv

ALMA Observations of the Debris Disk of Solar Analogue Tau Ceti

We present 1.3 mm observations of the Sun-like star $τ$ Ceti with the Atacama Large Millimeter/submillimeter Array (ALMA) that probe angular scales of $\sim1$'' (4 AU). This first interferometric image of the $τ$ Ceti system, which hosts both a debris disk and possible multiplanet system, shows emission from a nearly face-on belt of cold dust with a position angle of $90^\circ$ surrounding an unresolved central source at the stellar position. To characterize this emission structure, we fit parametric models to the millimeter visibilities. The resulting best-fit model yields an inner belt edge of $6.2^{+9.8}_{-4.6}$ AU, consistent with inferences from lower resolution, far-infrared Herschel observations. While the limited data at sufficiently short baselines preclude us from placing stronger constraints on the belt properties and its relation to the proposed five planet system, the observations do provide a strong lower limit on the fractional width of the belt, $ΔR/R > 0.75$ with $99\%$ confidence. This fractional width is more similar to broad disks such as HD 107146 than narrow belts such as the Kuiper Belt and Fomalhaut. The unresolved central source has a higher flux density than the predicted flux of the stellar photosphere at 1.3 mm. Given previous measurements of an excess by a factor of $\sim2$ at 8.7 mm, this emission is likely due to a hot stellar chromosphere.

preprint2016arXiv

An ALMA Search for Substructure, Fragmentation, and Hidden Protostars in Starless Cores in Chamaeleon I

We present an Atacama Large Millimeter/submillimeter Array (ALMA) 106 GHz (Band 3) continuum survey of the complete population of dense cores in the Chamaeleon I molecular cloud. We detect a total of 24 continuum sources in 19 different target fields. All previously known Class 0 and Class I protostars in Chamaeleon I are detected, whereas all of the 56 starless cores in our sample are undetected. We show that the Spitzer+Herschel census of protostars in Chamaeleon I is complete, with the rate at which protostellar cores have been misclassified as starless cores calculated as <1/56, or < 2%. We use synthetic observations to show that starless cores collapsing following the turbulent fragmentation scenario are detectable by our ALMA observations when their central densities exceed ~10^8 cm^-3, with the exact density dependent on the viewing geometry. Bonnor-Ebert spheres, on the other hand, remain undetected to central densities at least as high as 10^10 cm^-3. Our starless core non-detections are used to infer that either the star formation rate is declining in Chamaeleon I and most of the starless cores are not collapsing, matching the findings of previous studies, or that the evolution of starless cores are more accurately described by models that develop less substructure than predicted by the turbulent fragmentation scenario, such as Bonnor-Ebert spheres. We outline future work necessary to distinguish between these two possibilities.

preprint2016arXiv

Astrochemical Properties of Planck Cold Clumps

We observed thirteen Planck cold clumps with the James Clerk Maxwell Telescope/SCUBA-2 and with the Nobeyama 45 m radio telescope. The N$_2$H$^+$ distribution obtained with the Nobeyama telescope is quite similar to SCUBA-2 dust distribution. The 82 GHz HC$_3$N, 82 GHz CCS, and 94 GHz CCS emission are often distributed differently with respect to the N$_2$H$^+$ emission. The CCS emission, which is known to be abundant in starless molecular cloud cores, is often very clumpy in the observed targets. We made deep single-pointing observations in DNC, HN$^{13}$C, N$_2$D$^+$, cyclic-C$_3$H$_2$ toward nine clumps. The detection rate of N$_2$D$^+$ is 50\%. Furthermore, we observed the NH$_3$ emission toward 15 Planck cold clumps to estimate the kinetic temperature, and confirmed that most of targets are cold ($\lesssim$ 20 K). In two of the starless clumps observe, the CCS emission is distributed as it surrounds the N$_2$H$^+$ core (chemically evolved gas), which resembles the case of L1544, a prestellar core showing collapse. In addition, we detected both DNC and N$_2$D$^+$. These two clumps are most likely on the verge of star formation. We introduce the Chemical Evolution Factor (CEF) for starless cores to describe the chemical evolutionary stage, and analyze the observed Planck cold clumps.

preprint2016arXiv

Infall/Expansion Velocities in the Low-Mass Dense Cores L492, L694-2, and L1521F: Dependence on Position and Molecular Tracer

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.

preprint2016arXiv

The JCMT and Herschel Gould Belt Surveys: A comparison of SCUBA-2 and Herschel data of dense cores in the Taurus dark cloud L1495

We present a comparison of SCUBA-2 850-$μ$m and Herschel 70--500-$μ$m observations of the L1495 filament in the Taurus Molecular Cloud with the goal of characterising the SCUBA-2 Gould Belt Survey (GBS) data set. We identify and characterise starless cores in three data sets: SCUBA-2 850-$μ$m, Herschel 250-$μ$m, and Herschel 250-$μ$m spatially filtered to mimic the SCUBA-2 data. SCUBA-2 detects only the highest-surface-brightness sources, principally detecting protostellar sources and starless cores embedded in filaments, while Herschel is sensitive to most of the cloud structure, including extended low-surface-brightness emission. Herschel detects considerably more sources than SCUBA-2 even after spatial filtering. We investigate which properties of a starless core detected by Herschel determine its detectability by SCUBA-2, and find that they are the core's temperature and column density (for given dust properties). For similar-temperature cores, such as those seen in L1495, the surface brightnesses of the cores are determined by their column densities, with the highest-column-density cores being detected by SCUBA-2. For roughly spherical geometries, column density corresponds to volume density, and so SCUBA-2 selects the densest cores from a population at a given temperature. This selection effect, which we quantify as a function of distance, makes SCUBA-2 ideal for identifying those cores in Herschel catalogues that are closest to forming stars. Our results can now be used by anyone wishing to use the SCUBA-2 GBS data set.

preprint2015arXiv

Akari, SCUBA2 and Herschel data of pre-stellar cores

We show Akari data, Herschel data and data from the SCUBA2 camera on JCMT, of molecular clouds. We focus on pre-stellar cores within the clouds. We present Akari data of the L1147-1157 ring in Cepheus and show how the data indicate that the cores are being externally heated. We present SCUBA2 and Herschel data of the Ophiuchus region and show how the environment is also affecting core evolution in this region. We discuss the effects of the magnetic field in the Lupus I region, and how this lends support to a model for the formation and evolution of cores in filamentary molecular clouds.

preprint2015arXiv

ALMA Observations of the Largest Proto-Planetary Disk in the Orion Nebula, 114-426: A CO Silhouette

We present ALMA observations of the largest protoplanetary disk in the Orion Nebula, 114-426. Detectable 345 GHz (856 micron) dust continuum is produced only in the 350 AU central region of the ~1000 AU diameter silhouette seen against the bright H-alpha background in HST images. Assuming optically thin dust emission at 345 GHz, a gas-to-dust ratio of 100, and a grain temperature of 20 K, the disk gas-mass is estimated to be 3.1 +/- 0.6 Jupiter masses. If most solids and ices have have been incorporated into large grains, however, this value is a lower limit. The disk is not detected in dense-gas tracers such as HCO+ J=4-3, HCN J=4-3, or CS =7-6. These results may indicate that the 114-426 disk is evolved and depleted in some light organic compounds found in molecular clouds. The CO J=3-2 line is seen in absorption against the bright 50 to 80 K background of the Orion A molecular cloud over the full spatial extent and a little beyond the dust continuum emission. The CO absorption reaches a depth of 27 K below the background CO emission at VLSR ~6.7 km/s about 0.52 arcseconds (210 AU) northeast and 12 K below the background CO emission at VLSR ~ 9.7 km/s about 0.34 arcseconds (140 AU) southwest of the suspected location of the central star, implying that the embedded star has a mass less than 1 Solar mass .

preprint2015arXiv

Detection of Infall in the Protostar B335 with ALMA

Observations of the isolated globule B335 with ALMA have yielded absorption features against the continuum that are redshifted from the systemic velocity in both HCN and HCO$^+$ lines. These features provide unambiguous evidence for infall toward a central luminosity source. Previously developed models of inside-out collapse can match the observed line profiles of HCN and HCO$^+$ averaged over the central 50 AU. At the new distance of 100 pc, the inferred infall radius is 0.012 pc, the mass infall rate is $3 \times 10^{-6}$ solar masses per year, the age is 50,000 years, and the accumulated mass in the central zone is 0.15 solar masses, most of which must be in the star or in parts of a disk that are opaque at 0.8 mm. The continuum detection indicates an optically thin mass (gas and dust) of only $7.5\times 10^{-4}$ solar masses in the central region, consistent with only a very small disk mass.

preprint2015arXiv

Infrared and Radio observations of a small group of protostellar objects in the molecular core, L1251-C

We present a multi-wavelength observational study of a low-mass star-forming region, L1251-C, with observational results at wavelengths from the near-infrared to the millimeter. Spitzer Space Telescope observations confirmed that IRAS 22343+7501 is a small group of protostellar objects. The extended emission to east-west direction with its intensity peak at the center of L1251A has been detected at 350 and 850 micron with the CSO and JCMT telescopes, tracing dense envelope materials around L1251A. The single-dish data from the KVN and TRAO telescopes show inconsistencies between the intensity peaks of several molecular line emission and that of the continuum emission, suggesting complex distributions of molecular abundances around L1251A. The SMA interferometer data, however, show intensity peaks of CO 2-1 and 13CO 2-1 located at the position of IRS 1, which is both the brightest source in IRAC image and the weakest source in the 1.3 mm dust continuum map. IRS 1 is the strongest candidate for the driving source of the newly detected compact CO 2-1 outflow. Over the whole region (14' by 14') of L125l-C, 3 Class I and 16 Class II sources have been detected, including three YSOs in L1251A. A comparison with the average projected distance among 19 YSOs in L1251-C and that among 3 YSOs in L1251A suggests L1251-C is an example of low-mass cluster formation, where protostellar objects are forming in a small group.

preprint2015arXiv

Next Generation Very Large Array Memo No. 6, Science Working Group 1: The Cradle of Life

This paper discusses compelling science cases for a future long-baseline interferometer operating at millimeter and centimeter wavelengths, like the proposed Next Generation Vary Large Array (ngVLA). We report on the activities of the Cradle of Life science working group, which focused on the formation of low- and high-mass stars, the formation of planets and evolution of protoplanetary disks, the physical and compositional study of Solar System bodies, and the possible detection of radio signals from extraterrestrial civilizations. We propose 19 scientific projects based on the current specification of the ngVLA. Five of them are highlighted as possible Key Science Projects: (1) Resolving the density structure and dynamics of the youngest HII regions and high-mass protostellar jets, (2) Unveiling binary/multiple protostars at higher resolution, (3) Mapping planet formation regions in nearby disks on scales down to 1 AU, (4) Studying the formation of complex molecules, and (5) Deep atmospheric mapping of giant planets in the Solar System. For each of these projects, we discuss the scientific importance and feasibility. The results presented here should be considered as the beginning of a more in-depth analysis of the science enabled by such a facility, and are by no means complete or exhaustive.

preprint2015arXiv

Protoplanetary Disk Masses in the Young NGC 2024 Cluster

We present the results from a Submillimeter Array survey of the 887 micron continuum emission from the protoplanetary disks around 95 young stars in the young cluster NGC 2024. Emission was detected from 22 infrared sources, with flux densities from ~5 to 330 mJy; upper limits (at 3sigma) for the other 73 sources range from 3 to 24 mJy. For standard assumptions, the corresponding disk masses range from ~0.003 to 0.2Msolar, with upper limits at 0.002--0.01Msolar. The NGC 2024 sample has a slightly more populated tail at the high end of its disk mass distribution compared to other clusters, but without more information on the nature of the sample hosts it remains unclear if this difference is statistically significant or a superficial selection effect. Unlike in the Orion Trapezium, there is no evidence for a disk mass dependence on the (projected) separation from the massive star IRS2b in the NGC 2024 cluster. We suggest that this is due to either the cluster youth or a comparatively weaker photoionizing radiation field.

preprint2015arXiv

The AU Mic Debris Disk: far-infrared and submillimeter resolved imaging

We present far-infrared and submillimeter maps from the Herschel Space Observatory and the James Clerk Maxwell Telescope of the debris disk host star AU Microscopii. Disk emission is detected at 70, 160, 250, 350, 450, 500 and 850 micron. The disk is resolved at 70, 160 and 450 micron. In addition to the planetesimal belt, we detect thermal emission from AU Mic's halo for the first time. In contrast to the scattered light images, no asymmetries are evident in the disk. The fractional luminosity of the disk is $3.9 \times 10^{-4}$ and its mm-grain dust mass is 0.01 MEarth (+/- 20%). We create a simple spatial model that reconciles the disk SED as a blackbody of 53 +/- 2 K (a composite of 39 and 50 K components) and the presence of small (non-blackbody) grains which populate the extended halo. The best fit model is consistent with the "birth ring" model explored in earlier works, i.e., an edge-on dust belt extending from 8.8-40 AU, but with an additional halo component with an $r^{-1.5}$ surface density profile extending to the limits of sensitivity (140 AU). We confirm that AU Mic does not exert enough radiation force to blow out grains. For stellar mass loss rates of 10-100x solar, compact (zero porosity) grains can only be removed if they are very small, consistently with previous work, if the porosity is 0.9, then grains approaching 0.1 micron can be removed via corpuscular forces (i.e., the stellar wind).

preprint2015arXiv

The Kinematic and Chemical Properties of a Potential Core-Forming Clump: Perseus B1-E

We present 13CO and C18O (1-0), (2-1), and (3-2) maps towards the core-forming Perseus B1-E clump using observations from the James Clerk Maxwell Telescope (JCMT), Submillimeter Telescope (SMT) of the Arizona Radio Observatory, and IRAM 30 m telescope. We find that the 13CO and C18O line emission both have very complex velocity structures, indicative of multiple velocity components within the ambient gas. The (1-0) transitions reveal a radial velocity gradient across B1-E of 1 km/s/pc that increases from north-west to south-east, whereas the majority of the Perseus cloud has a radial velocity gradient increasing from south-west to north-east. In contrast, we see no evidence of a velocity gradient associated with the denser Herschel-identified substructures in B1-E. Additionally, the denser substructures have much lower systemic motions than the ambient clump material, which indicates that they are likely decoupled from the large-scale gas. Nevertheless, these substructures themselves have broad line widths (0.4 km/s) similar to that of the C18O gas in the clump, which suggests they inherited their kinematic properties from the larger-scale, moderately dense gas. Finally, we find evidence of C18O depletion only toward one substructure, B1-E2, which is also the only object with narrow (transonic) line widths. We suggest that as prestellar cores form, their chemical and kinematic properties are linked in evolution, such that these objects must first dissipate their turbulence before they deplete in CO.

preprint2015arXiv

Young Stellar Objects in the Gould Belt

We present the full catalog of Young Stellar Objects (YSOs) identified in the 18 molecular clouds surveyed by the Spitzer Space Telescope "cores to disks" (c2d) and "Gould Belt" (GB) Legacy surveys. Using standard techniques developed by the c2d project, we identify 3239 candidate YSOs in the 18 clouds, 2966 of which survive visual inspection and form our final catalog of YSOs in the Gould Belt. We compile extinction corrected SEDs for all 2966 YSOs and calculate and tabulate the infrared spectral index, bolometric luminosity, and bolometric temperature for each object. We find that 326 (11%), 210 (7%), 1248 (42%), and 1182 (40%) are classified as Class 0+I, Flat-spectrum, Class II, and Class III, respectively, and show that the Class III sample suffers from an overall contamination rate by background AGB stars between 25% and 90%. Adopting standard assumptions, we derive durations of 0.40-0.78 Myr for Class 0+I YSOs and 0.26-0.50 Myr for Flat-spectrum YSOs, where the ranges encompass uncertainties in the adopted assumptions. Including information from (sub)millimeter wavelengths, one-third of the Class 0+I sample is classified as Class 0, leading to durations of 0.13-0.26 Myr (Class 0) and 0.27-0.52 Myr (Class I). We revisit infrared color-color diagrams used in the literature to classify YSOs and propose minor revisions to classification boundaries in these diagrams. Finally, we show that the bolometric temperature is a poor discriminator between Class II and Class III YSOs.

preprint2014arXiv

ALMA observations of a misaligned binary protoplanetary disk system in Orion

We present ALMA observations of a wide binary system in Orion, with projected separation 440 AU, in which we detect submillimeter emission from the protoplanetary disks around each star. Both disks appear moderately massive and have strong line emission in CO 3-2, HCO+ 4-3, and HCN 3-2. In addition, CS 7-6 is detected in one disk. The line-to-continuum ratios are similar for the two disks in each of the lines. From the resolved velocity gradients across each disk, we constrain the masses of the central stars, and show consistency with optical-infrared spectroscopy, both indicative of a high mass ratio ~9. The small difference between the systemic velocities indicates that the binary orbital plane is close to face-on. The angle between the projected disk rotation axes is very high, ~72 degrees, showing that the system did not form from a single massive disk or a rigidly rotating cloud core. This finding, which adds to related evidence from disk geometries in other systems, protostellar outflows, stellar rotation, and similar recent ALMA results, demonstrates that turbulence or dynamical interactions act on small scales well below that of molecular cores during the early stages of star formation.

preprint2014arXiv

ALMA Observations of the Orion Proplyds

We present ALMA observations of protoplanetary disks ("proplyds") in the Orion Nebula Cluster. We imaged 5 individual fields at 856um containing 22 HST-identified proplyds and detected 21 of them. Eight of those disks were detected for the first time at submillimeter wavelengths, including the most prominent, well-known proplyd in the entire Orion Nebula, 114-426. Thermal dust emission in excess of any free-free component was measured in all but one of the detected disks, and ranged between 1-163 mJy, with resulting disk masses of 0.3-79 Mjup. An additional 26 stars with no prior evidence of associated disks in HST observations were also imaged within the 5 fields, but only 2 were detected. The disk mass upper limits for the undetected targets, which include OB stars, theta1Ori C and theta1Ori F, range from 0.1-0.6 Mjup. Combining these ALMA data with previous SMA observations, we find a lack of massive (>3 Mjup) disks in the extreme-UV dominated region of Orion, within 0.03 pc of O-star theta1Ori C. At larger separations from theta1Ori C, in the far-UV dominated region, there is a wide range of disk masses, similar to what is found in low-mass star forming regions. Taken together, these results suggest that a rapid dissipation of disk masses likely inhibits potential planet formation in the extreme-UV dominated regions of OB associations, but leaves disks in the far-UV dominated regions relatively unaffected.

preprint2014arXiv

Evidence for Large Grains in the Star-forming Filament OMC-2/3

We present a new 3.3 mm continuum map of the OMC-2/3 region. When paired with previously published maps of 1.2mm continuum and NH3-derived temperature, we derive the emissivity spectral index of dust emission in this region, tracking its changes across the filament and cores. We find that the median value of the emissivity spectral index is 0.9, much shallower than previous estimates in other nearby molecular clouds. We find no significant difference between the emissivity spectral index of dust in the OMC-2/3 filament and the starless or protostellar cores. Furthermore, the temperature and emissivity spectral index, beta, are anti-correlated at the 4 sigma level. The low values of the emissivity spectral index found in OMC-2/3 can be explained by the presence of millimeter-sized dust grains in the dense regions of the filaments to which these maps are most sensitive. Alternatively, a shallow dust emissivity spectral index may indicate non-powerlaw spectral energy distributions, significant free-free emission, or anomalous microwave emission. We discuss the possible implications of millimeter-sized dust grains compared to the alternatives.

preprint2014arXiv

The disk around the brown dwarf KPNO Tau 3

We present submillimeter observations of the young brown dwarfs KPNO Tau 1, KPNO Tau 3, and KPNO Tau 6 at 450 micron and 850 micron taken with the Submillimeter Common-User Bolometer Array on the James Clerke Maxwell Telescope. KPNO Tau 3 and KPNO Tau 6 have been previously identified as Class II objects hosting accretion disks, whereas KPNO Tau 1 has been identified as a Class III object and shows no evidence of circumsubstellar material. Our 3 sigma detection of cold dust around KPNO Tau 3 implies a total disk mass of (4.0 +/- 1.1) x 10^{-4} Msolar (assuming a gas to dust ratio of 100:1). We place tight constraints on any disks around KPNO Tau 1 or KPNO Tau 6 of <2.1 x 10^{-4} Msolar and <2.7 x 10^{-4} Msolar, respectively. Modeling the spectral energy distribution of KPNO Tau 3 and its disk suggests the disk properties (geometry, dust mass, and grain size distribution) are consistent with observations of other brown dwarf disks and low-mass T-Tauri stars. In particular, the disk-to-host mass ratio for KPNO Tau 3 is congruent with the scenario that at least some brown dwarfs form via the same mechanism as low-mass stars.

preprint2014arXiv

The JCMT Legacy Survey of the Gould Belt: a molecular line study of the Ophiuchus molecular cloud

CO, $^{13}$CO and C$^{18}$O ${\it J}$ = 3--2 observations are presented of the Ophiuchus molecular cloud. The $^{13}$CO and C$^{18}$O emission is dominated by the Oph A clump, and the Oph B1, B2, C, E, F and J regions. The optically thin(ner) C$^{18}$O line is used as a column density tracer, from which the gravitational binding energy is estimated to be $4.5 \times 10^{39}$ J (2282 $M_\odot$ km$^2$ s$^{-2}$). The turbulent kinetic energy is $6.3 \times 10^{38}$ J (320 $M_\odot$ km$^2$ s$^{-2}$), or 7 times less than this, and therefore the Oph cloud as a whole is gravitationally bound. Thirty protostars were searched for high velocity gas, with eight showing outflows, and twenty more having evidence of high velocity gas along their lines-of-sight. The total outflow kinetic energy is $1.3 \times 10^{38}$ J (67 $M_\odot$ km$^2$ s$^{-2}$), corresponding to 21$\%$ of the cloud's turbulent kinetic energy. Although turbulent injection by outflows is significant, but does ${\it not}$ appear to be the dominant source of turbulence in the cloud. 105 dense molecular clumplets were identified, which had radii $\sim$ 0.01--0.05 pc, virial masses $\sim$ 0.1--12 $M_\odot$, luminosities $\sim$ 0.001--0.1 K~km s$^{-1}$ pc$^{-2}$, and excitation temperatures $\sim$ 10--50K. These are consistent with the standard GMC based size-line width relationships, showing that the scaling laws extend down to size scales of hundredths of a parsec, and to sub solar-mass condensations. There is however no compelling evidence that the majority of clumplets are undergoing free-fall collapse, nor that they are pressure confined.

preprint2013arXiv

A Herschel and APEX Census of the Reddest Sources in Orion: Searching for the Youngest Protostars

We perform a census of the reddest, and potentially youngest, protostars in the Orion molecular clouds using data obtained with the PACS instrument onboard the Herschel Space Observatory and the LABOCA and SABOCA instruments on APEX as part of the Herschel Orion Protostar Survey (HOPS). A total of 55 new protostar candidates are detected at 70 um and 160 um that are either too faint (m24 > 7 mag) to be reliably classified as protostars or undetected in the Spitzer/MIPS 24 um band. We find that the 11 reddest protostar candidates with log (lambda F_lambda 70) / (lambda F_lambda 24) > 1.65 are free of contamination and can thus be reliably explained as protostars. The remaining 44 sources have less extreme 70/24 colors, fainter 70 um fluxes, and higher levels of contamination. Taking the previously known sample of Spitzer protostars and the new sample together, we find 18 sources that have log (lambda F_lambda 70) / (lambda F_lambda 24) > 1.65; we name these sources "PACS Bright Red sources", or PBRs. Our analysis reveals that the PBRs sample is composed of Class 0 like sources characterized by very red SEDs (T_bol < 45 K) and large values of sub-millimeter fluxes (L_smm/L_bol > 0.6%). Modified black-body fits to the SEDs provide lower limits to the envelope masses of 0.2 M_sun to 2 M_sun and luminosities of 0.7 L_sun to 10 L_sun. Based on these properties, and a comparison of the SEDs with radiative transfer models of protostars, we conclude that the PBRs are most likely extreme Class 0 objects distinguished by higher than typical envelope densities and hence, high mass infall rates.

preprint2013arXiv

A Roadmap for Canadian Submillimetre Astronomy

We survey the present landscape in submillimetre astronomy for Canada and describe a plan for continued engagement in observational facilities to ~2020. Building on Canada's decadal Long Range Plan process, we emphasize that continued involvement in a large, single-dish facility is crucial given Canada's substantial investment in ALMA and numerous PI-led submillimetre experiments. In particular, we recommend: i) an extension of Canadian participation in the JCMT until at least the unique JCMT Legacy Survey program is able to realize the full scientific potential provided by the world-leading SCUBA-2 instrument; and ii) involvement of the entire Canadian community in CCAT, with a large enough share in the partnership for Canadian astronomers to participate at all levels of the facility. We further recommend continued participation in ALMA development, involvement in many focused PI-led submillimetre experiments, and partnership in SPICA.

preprint2013arXiv

Correlating Infall with Deuterium Fractionation in Dense Cores

We present a survey of HCO+ (3-2) observations pointed towards dense cores with previous measurements of N(N2D+)/N(N2H+). Of the 26 cores in this survey, five show the spectroscopic signature of outward motion, nine exhibit neither inward nor outward motion, eleven appear to be infalling, and one is not detected. We compare the degree of deuterium fractionation with infall velocities calculated from the HCO+ spectra and find that those cores with [D]/[H] > 0.1 are more likely to have the signature of inward motions than cores with smaller [D]/[H] ratios. Infall motions are also much more common in cores with masses exceeding their thermal Jeans masses. The fastest infall velocity measured belongs to one of the two protostellar cores in our survey, L1521F, and the observed motions are typically on the order of the sound speed.

preprint2013arXiv

Does a prestellar core always become protostellar? Tracing the evolution of cores from the prestellar to protostellar phase

Recently, a subset of starless cores whose thermal Jeans mass is apparently overwhelmed by the mass of the core has been identified, e.g., the core {\small L183}. In literature, massive cores such as this one are often referred to as "super-Jeans cores". As starless cores are perhaps on the cusp of forming stars, a study of their dynamics will improve our understanding of the transition from the prestellar to the protostellar phase. In the present work we use non-magnetic polytropes belonging originally to the family of the Isothermal sphere. For the purpose, perturbations were applied to individual polytropes, first by replacing the isothermal gas with a gas that was cold near the centre of the polytrope and relatively warm in the outer regions, and second, through a slight compression of the polytrope by raising the external confining pressure. Using this latter configuration we identify thermodynamic conditions under which a core is likely to remain starless. In fact, we also argue that the attribute "super-Jeans" is subjective and that these cores do not formally violate the Jeans stability criterion. On the basis of our test results we suggest that gas temperature in a star-forming cloud is crucial towards the formation and evolution of a core. Simulations in this work were performed using the particle-based Smoothed Particle Hydrodynamics algorithm. However, to establish numerical convergence of the results we suggest similar tests with a grid-scheme, such as the Adaptive mesh refinement.

preprint2013arXiv

From Filamentary Networks to Dense Cores in Molecular Clouds: Toward a New Paradigm for Star Formation

Recent studies of the nearest star-forming clouds of the Galaxy at submillimeter wavelengths with the Herschel Space Observatory have provided us with unprecedented images of the initial and boundary conditions of the star formation process. The Herschel results emphasize the role of interstellar filaments in the star formation process and connect remarkably well with nearly a decade's worth of numerical simulations and theory that have consistently shown that the ISM should be highly filamentary on all scales and star formation is intimately related to self-gravitating filaments. In this review, we trace how the apparent complexity of cloud structure and star formation is governed by relatively simple universal processes - from filamentary clumps to galactic scales. We emphasize two crucial and complementary aspects: (i) the key observational results obtained with Herschel over the past three years, along with relevant new results obtained from the ground on the kinematics of interstellar structures, and (ii) the key existing theoretical models and the many numerical simulations of interstellar cloud structure and star formation. We then synthesize a comprehensive physical picture that arises from the confrontation of these observations and simulations.

preprint2013arXiv

Water and Methanol Maser Survey of Protostars in the Orion Molecular Cloud Complex

The results of a maser survey toward ninety-nine protostars in the Orion molecular cloud complex are presented. The target sources are low-mass protostars identified from infrared observations. Single-dish observations were carried out in the water maser line at 22 GHz and the methanol class I maser lines at 44, 95, and 133 GHz. Most of the detected sources were mapped to determine the source positions. Five water maser sources were detected, and they are excited by HH 1-2 VLA 3, HH 1-2 VLA 1, L1641N MM1/3, NGC 2071 IRS 1/3, and an object in the OMC 3 region. The water masers showed significant variability in intensity and velocity with time scales of a month or shorter. Four methanol emission sources were detected, and those in the OMC 2 FIR 3/4 and L1641N MM1/3 regions are probably masers. The methanol emission from the other two sources in the NGC 2071 IRS 1-3 and V380 Ori NE regions are probably thermal. For the water masers, the number of detections per protostar in the survey region is about 2%, which suggests that the water masers of low-mass protostars are rarely detectable. The methanol class I maser of low-mass protostars is an even rarer phenomenon, with a detection rate much smaller than 1%.

preprint2012arXiv

A First Look at the Auriga-California Giant Molecular Cloud With Herschel and the CSO: Census of the Young Stellar Objects and the Dense Gas

We have mapped the Auriga/California molecular cloud with the Herschel PACS and SPIRE cameras and the Bolocam 1.1 mm camera on the Caltech Submillimeter Observatory (CSO) with the eventual goal of quantifying the star formation and cloud structure in this Giant Molecular Cloud (GMC) that is comparable in size and mass to the Orion GMC, but which appears to be forming far fewer stars. We have tabulated 60 compact 70/160um sources that are likely pre-main-sequence objects and correlated those with Spitzer and WISE mid-IR sources. At 1.1 mm we find 18 cold, compact sources and discuss their properties. The most important result from this part of our study is that we find a modest number of additional compact young objects beyond those identified at shorter wavelengths with Spitzer. We also describe the dust column density and temperature structure derived from our photometric maps. The column density peaks at a few x 10^22 cm^-2 (N_H2) and is distributed in a clear filamentary structure along which nearly all the pre-main-sequence objects are found. We compare the YSO surface density to the gas column density and find a strong non-linear correlation between them. The dust temperature in the densest parts of the filaments drops to ~10K from values ~ 14--15K in the low density parts of the cloud. We also derive the cumulative mass fraction and probability density function of material in the cloud which we compare with similar data on other star-forming clouds.

preprint2012arXiv

A High Resolution Study of the HI-H2 Transition across the Perseus Molecular Cloud

To investigate the fundamental principles of H2 formation in a giant molecular cloud (GMC), we derive the HI and H2 surface density (Sigma_HI and Sigma_H2) images of the Perseus molecular cloud on sub-pc scales (~0.4 pc). We use the far-infrared data from the Improved Reprocessing of the IRAS Survey and the V-band extinction image provided by the COMPLETE Survey to estimate the dust column density image of Perseus. In combination with the HI data from the Galactic Arecibo L-band Feed Array HI Survey and an estimate of the local dust-to-gas ratio, we then derive the Sigma_H2 distribution across Perseus. We find a relatively uniform Sigma_HI ~ 6-8 Msun pc^-2 for both dark and star-forming regions, suggesting a minimum HI surface density required to shield H2 against photodissociation. As a result, a remarkably tight and consistent relation is found between Sigma_H2/Sigma_HI and Sigma_HI+Sigma_H2. The transition between the HI- and H2-dominated regions occurs at N(HI)+2N(H2) ~ (8-14) x 10^20 cm^-2. Our findings are consistent with predictions for H2 formation in equilibrium, suggesting that turbulence may not be of primary importance for H2 formation. However, the importance of a warm neutral medium for H2 shielding, an internal radiation field, and the timescale of H2 formation still remain as open questions. We also compare H2 and CO distributions and estimate the fraction of "CO-dark" gas, f_DG ~ 0.3. While significant spatial variations of f_DG are found, we do not find a clear correlation with the mean V-band extinction.

preprint2012arXiv

An Observed Lack of Substructure in Starless Cores II: Super-Jeans Cores

We present SMA and CARMA continuum and spectral line observations of five dense cores located in the Perseus and Ophiuchus molecular clouds whose masses exceed their thermal Jeans masses. Three of these cores have previously been identified as being starless and two have been classified as being possibly protostellar. We find that one core is certainly protostellar. The other four cores, however, are starless and undetected in both C18O (2-1) and 1.3mm continuum emission. These four starless cores have flat density profiles out to at least 0.006 pc, which is typical for starless cores in general. Density profiles predicted by some collapse models, especially in the early stages of infall, are consistent with our observations. Archival data reveal that these starless cores have significant non-thermal support against collapse, although they may still be unstable.

preprint2012arXiv

Changes Of Dust Opacity With Density in the Orion A Molecular Cloud

We have studied the opacity of dust grains at submillimeter wavelengths by estimating the optical depth from imaging at 160, 250, 350, and 500 um from the Herschel Gould Belt Survey and comparing this to a column density obtained from the 2MASS-derived color excess E(J-Ks). Our main goal was to investigate the spatial variations of the opacity due to "big" grains over a variety of environmental conditions and thereby quantify how emission properties of the dust change with column (and volume) density. The central and southern areas of the Orion A molecular cloud examined here, with NH ranging from 1.5X10^21 cm^-2 to 50X10^21 cm^-2, are well suited to this approach. We fit the multi-frequency Herschel spectral energy distributions (SEDs) of each pixel with a modified blackbody to obtain the temperature, T, and optical depth, τ(1200), at a fiducial frequency of 1200 GHz (250 um). Using a calibration of NH/E(J-Ks)for the interstellar medium (ISM) we obtained the opacity (dust emission cross-section per H nucleon), σ_e(1200), for every pixel. From a value of ~ 1X10^-25 cm^2 H^-1 at the lowest column densities that is typical of the high latitude diffuse ISM, σ_e(1200) increases as NH^0.28 over the range studied. This is suggestive of grain evolution. Integrating the SEDs over frequency, we also calculated the specific power P (emission power per H) for the big grains. In low column density regions where dust clouds are optically thin to the interstellar radiation field (ISRF), P is typically 3.7 X 10^-31 W H^-1, again close to that in the high latitude diffuse ISM. However, we find evidence for a decrease of P in high column density regions, which would be a natural outcome of attenuation of the ISRF that heats the grains, and for localized increases for dust illuminated by nearby stars or embedded protostars.

preprint2011arXiv

Filaments and ridges in Vela C revealed by Herschel: from low-mass to high-mass star-forming sites

We present the first Herschel PACS and SPIRE results of the Vela C molecular complex in the far-infrared and submillimetre regimes at 70, 160, 250, 350, and 500 um, spanning the peak of emission of cold prestellar or protostellar cores. Column density and multi-resolution analysis (MRA) differentiates the Vela C complex into five distinct sub-regions. Each sub-region displays differences in their column density and temperature probability distribution functions (PDFs), in particular, the PDFs of the `Centre-Ridge' and `South-Nest' sub-regions appear in stark contrast to each other. The Centre-Ridge displays a bimodal temperature PDF representative of hot gas surrounding the HII region RCW 36 and the cold neighbouring filaments, whilst the South-Nest is dominated by cold filamentary structure. The column density PDF of the Centre-Ridge is flatter than the South-Nest, with a high column density tail, consistent with formation through large-scale flows, and regulation by self-gravity. At small to intermediate scales MRA indicates the Centre-Ridge to be twice as concentrated as the South-Nest, whilst on larger scales, a greater portion of the gas in the South-Nest is dominated by turbulence than in the Centre-Ridge. In Vela C, high-mass stars appear to be preferentially forming in ridges, i.e., dominant high column density filaments.

preprint2011arXiv

How Starless Are Starless Cores?

In this paper, we present the results of CARMA continuum and spectral line observations of the dense core Per-Bolo 45. Although this core has previously been classified as starless, we find evidence for an outflow and conclude that Per-Bolo 45 is actually an embedded, low-luminosity protostar. We discuss the impact of newly discovered, low-luminosity, embedded objects in the Perseus molecular cloud on starless core and protostar lifetimes. We estimate that the starless core lifetime has been overestimated by 4-18% and the Class 0/I protostellar lifetime has been underestimated by 5-20%. Given the relatively large systematic uncertainties involved in these calculations, variations on the order of 10% do not significantly change either core lifetimes or the expected protostellar luminosity function. Finally, we suggest that high resolution (sub)millimeter surveys of known cores lacking near-infrared and mid-infrared emission are necessary to make an accurate census of starless cores.

preprint2011arXiv

Initial conditions for star formation in clusters: physical and kinematical structure of the starless core Oph A-N6

We present high spatial (<300 AU) and spectral (0.07 km/s) resolution Submillimeter Array observations of the dense starless cluster core Oph A-N6, in the 1 mm dust continuum and the 3-2 line of N2H+ and N2D+. The dust continuum observations reveal a compact source not seen in single-dish observations, of size ~1000 AU and mass 0.005-0.01 M\odot. The combined line and single-dish observations reveal a core of size 3000 \times 1400 AU elongated in a NW-SE direction, with almost no variation in either line width or line center velocity across the map, and very small non-thermal motions. The deuterium fraction has a peak value of ~0.15 and is >0.05 over much of the core. The N2H+ column density profile across the major axis of Oph A-N6 is well represented by an isothermal cylinder, with temperature 20 K, peak density 7.1 \times 10^6 cm^{-3}, and N2H+ abundance 2.7 \times 10^{-10}. The mass of Oph A-N6 is estimated to be 0.29 M\odot, compared to a value of 0.18 M\odot from the isothermal cylinder analysis, and 0.63 M\odot for the critical mass for fragmentation of an isothermal cylinder. Compared to isolated low-mass cores, Oph A-N6 shows similar narrow line widths and small velocity variation, with a deuterium fraction similar to "evolved" dense cores. It is significantly smaller than isolated cores, with larger peak column and volume density. The available evidence suggests Oph A-N6 has formed through the fragmentation of the Oph A filament and is the precursor to a low-mass star. The dust continuum emission suggests it may already have begun to form a star.

preprint2010arXiv

"Starless" Super-Jeans Cores in Four Gould Belt Clouds

From a survey of 729 cores based on JCMT/SCUBA data, we present an analysis of 17 candidate starless cores with masses that exceed their stable Jeans masses. We re-examine the classification of these super-Jeans cores using Spitzer maps and find that 3 are re-classified as protostellar, 11 have ambiguous emission near the core positions, and 3 appear to be genuinely starless. We suggest the 3 starless and 11 undetermined super-Jeans cores represent excellent targets for future observational and computational study to understand the evolution of dense cores and the process of star formation.

preprint2010arXiv

IRAS 22198+6336: Discovery of an Intermediate-Mass Hot Core

We present new SMA and PdBI observations of the intermediate-mass object IRAS 22198+6336 in the millimeter continuum and in several molecular line transitions. The millimeter continuum emission reveals a strong and compact source with a mass of around 5 Msun and with properties of Class 0 objects. CO emission shows an outflow with a quadrupolar morphology centered on the position of the dust condensation. The CO outflow emission seems to come from two distinct outflows, one of them associated with SiO outflow emission. A large set of molecular lines has been detected toward a compact dense core clearly coincident with the compact millimeter source, and showing a velocity gradient perpendicular to the outflow traced by CO and SiO. The chemically rich spectrum and the rotational temperatures derived from CH$_3$CN and CH$_3$OH (100-150 K) indicate that IRAS 22198+6336 is harbouring one the few intermediate-mass hot cores known at present.

preprint2010arXiv

The Submillimetre Universe

Submillimetre continuum radiation allows us to probe cold objects, particularly the earliest, dusty phases of star formation, high-redshift galaxies and circumstellar disks. The submillimetre window gives a unique view of the physical and dynamical conditions in the neutral and molecular interstellar medium. In the next decade a combination of wide-field surveys with single-dish telescopes and targeted follow-up with ALMA and other facilities should enable rapid progress in answering questions about the origins of planetary systems, stars and galaxies.

preprint2007arXiv

PROSAC: A Submillimeter Array Survey of Low-Mass Protostars. I. Overview of Program: Envelopes, Disks, Outflows and Hot Cores

This paper presents a large spectral line and continuum survey of 8 deeply embedded, low-mass protostellar cores using the Submillimeter Array. Each source was observed in high excitation lines of some of the most common molecular species, CO, HCO+, CS, SO, H2CO, CH3OH and SiO. Line emission from 11 species originating from warm and dense gas have been imaged at high angular resolution (1-3"; typically 200-600 AU) together with continuum emission at 230 GHz (1.3 mm) and 345 GHz (0.8 mm). Compact continuum emission is observed for all sources which likely originates in marginally optically thick circumstellar disks, with typical lower limits to their masses of 0.1 M_sun (1-10% of the masses of their envelopes) and having a dust opacity law with beta approximately 1. Prominent outflows are present in CO 2-1 observations in all sources: the most diffuse outflows are found in the sources with the lowest ratios of disk-to-envelope mass, and it is suggested that these sources are in a phase where accretion of matter from the envelope has almost finished and the remainder of the envelope material is being dispersed by the outflows. Other characteristic dynamical signatures are found with inverse P Cygni profiles indicative of infalling motions seen in the 13CO 2-1 lines toward NGC1333-IRAS4A and -IRAS4B. Outflow-induced shocks are present on all scales in the protostellar environments and are most clearly traced by the emission of CH3OH in NGC1333-IRAS4A and -IRAS4B. These observations suggest that the emission of CH3OH and H2CO from these proposed "hot corinos" are related to the shocks caused by the protostellar outflows. Only one source, NGC1333-IRAS2A, has evidence for hot, compact CH3OH emission coincident with the embedded protostar.

preprint2006arXiv

The Spitzer c2d Survey of Nearby Dense Cores: III: Low Mass Star Formation in a Small Group, L1251B

We present a comprehensive study of a low-mass star-forming region,L1251B, at wavelengths from the near-infrared to the millimeter. L1251B, where only one protostar, IRAS 22376+7455, was known previously, is confirmed to be a small group of protostars based on observations with the Spitzer Space Telescope. The most luminous source of L1251B is located 5" north of the IRAS position. A near-infrared bipolar nebula, which is not associated with the brightest object and is located at the southeast corner of L1251B, has been detected in the IRAC bands. OVRO and SMA interferometric observations indicate that the brightest source and the bipolar nebula source in the IRAC bands are deeply embedded disk sources.Submillimeter continuum observations with single-dish telescopes and the SMA interferometric observations suggest two possible prestellar objects with very high column densities. Outside of the small group, many young stellar object candidates have been detected over a larger region of 12' x 12'. Extended emission to the east of L1251B has been detected at 850 micron; this "east core" may be a site for future star formation since no point source has been detected with IRAC or MIPS. This region is therefore a possible example of low-mass cluster formation, where a small group of pre- and protostellar objects (L1251B) is currently forming, alongside a large starless core (the east core).