Source author record

Doug Johnstone

Doug Johnstone 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)

preprint2025arXiv

FAUST XXX: Dust enhancement in the young binary L1551 IRS 5

Young binary stars with discs provide unique laboratories to study the earliest stages of planet formation in star-forming environments. The detection of substructure in discs around Class I protostars challenges current models of disc evolution, suggesting that planets may form earlier than previously expected ($<1$ Myr). In the context of the ALMA Large Program FAUST, we present observations of the circumbinary disc (CBD) around the young binary system L1551 IRS 5. The CBD exhibits two prominent over-densities in the continuum emission at the edge of the cavity, with the Northern over-density being about 20% brighter than the Southern one. By analysing the disc morphology and kinematics of L1551 IRS 5, we delineate dynamical constraints on the binary's orbital parameters. Additionally, we present 3D hydrodynamical models of the CBD to predict both the dust and the gas surface densities. Then, we compare the resulting synthetic observations with ALMA observations of the continuum emission at 1.3 mm and the C$^{18}$O line emission. Our analysis suggests that the density enhancements observed with ALMA in L1551 IRS 5 can be caused by interactions between the binary stars and the CBD, leading to dust concentration within the disc. We conclude that the observed over-density corresponds to a location where could potentially grow under favourable conditions.

preprint2024arXiv

JWST/NIRCam Imaging of Young Stellar Objects. I. Constraints on Planets Exterior to The Spiral Disk Around MWC 758

MWC 758 is a young star hosting a spiral protoplanetary disk. The spirals are likely companion-driven, and two previously-identified candidate companions have been identified -- one at the end the Southern spiral arm at ~0.6 arcsec, and one interior to the gap at ~0.1 arcsec. With JWST/NIRCam, we provide new images of the disk and constraints on planets exterior to ~1". We detect the two-armed spiral disk, a known background star, and a spatially resolved background galaxy, but no clear companions. The candidates that have been reported are at separations that are not probed by our data with sensitivity sufficient to detect them -- nevertheless, these observations place new limits on companions down to ~2 Jupiter-masses at ~150 au and ~0.5 Jupiter masses at ~600 au. Owing to the unprecedented sensitivity of JWST and youth of the target, these are among the deepest mass-detection limits yet obtained through direct imaging observations, and provide new insights into the system's dynamical nature.

preprint2024arXiv

JWST/NIRCam Imaging of Young Stellar Objects. II. Deep Constraints on Giant Planets and a Planet Candidate Outside of the Spiral Disk Around SAO 206462

We present JWST/NIRCam F187N, F200W, F405N and F410M direct imaging data of the disk surrounding SAO 206462. Previous images show a very structured disk, with a pair of spiral arms thought to be launched by one or more external perturbers. The spiral features are visible in three of the four filters, with the non-detection in F410M due to the large detector saturation radius. We detect with a signal-to-noise ratio of 4.4 a companion candidate (CC1) that, if on a coplanar circular orbit, would orbit SAO 206462 at a separation of $\sim300$ au, $2.25σ$ away from the predicted separation for the driver of the eastern spiral. According to the BEX models, CC1 has a mass of $M_\mathrm{CC1}=0.8\pm0.3~M_\mathrm{J}$. No other companion candidates were detected. At the location predicted by simulations of both spirals generated by a single massive companion, the NIRCam data exclude objects more massive than $\sim2.2~M_\mathrm{J}$ assuming the BEX evolutionary models. In terms of temperatures, the data are sensitive to objects with $T_{\text{eff}}\sim650-850$ K, when assuming planets emit like blackbodies ($R_\mathrm{p}$ between 1 and $3 R_\mathrm{J}$). From these results, we conclude that if the spirals are driven by gas giants, these must be either cold or embedded in circumplanetary material. In addition, the NIRCam data provide tight constraints on ongoing accretion processes. In the low extinction scenario we are sensitive to mass accretion rates of the order $\dot{M}\sim10^{-9} M_\mathrm{J}$ yr$^{-1}$. Thanks to the longer wavelengths used to search for emission lines, we reach unprecedented sensitivities to processes with $\dot{M}\sim10^{-7} M_\mathrm{J}$ yr$^{-1}$ even towards highly extincted environments ($A_\mathrm{V}\approx50$~mag).

preprint2024arXiv

The JCMT Transient Survey: Six-Year Summary of 450/850\,$μ$m Protostellar Variability and Calibration Pipeline Version 2.0

The JCMT Transient Survey has been monitoring eight Gould Belt low-mass star-forming regions since December 2015 and six somewhat more distant intermediate-mass star-forming regions since February 2020 with SCUBA-2 on the JCMT at \ShortS and \LongS and with an approximately monthly cadence. We introduce our Pipeline v2 relative calibration procedures for image alignment and flux calibration across epochs, improving on our previous Pipeline v1 by decreasing measurement uncertainties and providing additional robustness. These new techniques work at both \LongS and \ShortNS, where v1 only allowed investigation of the \LongS data. Pipeline v2 achieves better than $0.5^{\prime\prime}$ relative image alignment, less than a tenth of the submillimeter beam widths. The v2 relative flux calibration is found to be 1\% at \LongS and $<5$\% at \ShortNS. The improvement in the calibration is demonstrated by comparing the two pipelines over the first four years of the survey and recovering additional robust variables with v2. Using the full six years of the Gould Belt survey the number of robust variables increases by 50\,\%, and at \ShortS we identify four robust variables, all of which are also robust at \LongNS. The multi-wavelength light curves for these sources are investigated and found to be consistent with the variability being due to dust heating within the envelope in response to accretion luminosity changes from the central source.

preprint2024arXiv

The outflow of the protostar in B335: I

The isolated globule B335 contains a single, low luminosity Class 0 protostar associated with a bipolar nebula and outflow system seen nearly perpendicular to its axis. We observed the innermost regions of this outflow as part of JWST/NIRCam GTO program 1187, primarily intended for wide-field slitless spectroscopy of background stars behind the globule. We find a system of expanding shock fronts with kinematic ages of only a few decades emerging symmetrically from the position of the embedded protostar, which is not directly detected at NIRCam wavelengths. The innermost and youngest of the shock fronts studied here shows strong emission from CO. The next older shock front shows less CO and the third shock front shows only H_2 emission in our data. This third and most distant of these inner shock fronts shows substantial evolution of its shape since it was last observed with high spatial resolution in 1996 with Keck/NIRC. This may be evidence of a faster internal shock catching up with a slower one and of the two shocks merging.

preprint2022arXiv

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): A Hot Corino Survey toward Protostellar Cores in the Orion Cloud

The presence of complex organic molecules (COMs) in the interstellar medium (ISM) is of great interest since it may link to the origin and prevalence of life in the universe. Aiming to investigate the occurrence of COMs and their possible origins, we conducted a chemical census toward a sample of protostellar cores as part of the ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP) project. We report the detection of 11 hot corino sources, which exhibit compact emissions from warm and abundant COMs, among 56 Class 0/I protostellar cores. All the hot corino sources discovered are likely Class 0 and their sizes of the warm region ($>$ 100 K) are comparable to 100 au. The luminosity of the hot corino sources exhibits positive correlations with the total number of methanol and the extent of its emissions. Such correlations are consistent with the thermal desorption picture for the presence of hot corino and suggest that the lower luminosity (Class 0) sources likely have a smaller region with COMs emissions. With the same sample selection method and detection criteria being applied, the detection rates of the warm methanol in the Orion cloud (15/37) and the Perseus cloud (28/50) are statistically similar when the cloud distances and the limited sample size are considered. Observing the same set of COM transitions will bring a more informative comparison between the cloud properties.

preprint2022arXiv

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): Deriving Inclination Angle and Velocity of the Protostellar Jets from their SiO Knots

We have selected six sources (G209.55-19.68S2, G205.46-14.56S1$_{-}$A, G203.21-11.20W2, G191.90-11.21S, G205.46-14.56S3, and G206.93-16.61W2) from the Atacama Large Millimeter/submillimeter Array Survey of Orion Planck Galactic Cold Clumps (ALMASOP), in which these sources have been mapped in the CO (J=2-1), SiO (J=5-4), and C$^{18}$O (J=2-1) lines. These sources have high-velocity SiO jets surrounded by low-velocity CO outflows. The SiO jets consist of a chain of knots. These knots have been thought to be produced by semi-periodical variations in jet velocity. Therefore, we adopt a shock-forming model, which uses such variations to estimate the inclination angle and velocity of the jets. We also derive the inclination angle of the CO outflows using the wide-angle wind-driven shell model, and find it to be broadly consistent with that of the associated SiO jets. In addition, we apply this shock-forming model to another three protostellar sources with SiO jets in the literature -- HH 211, HH 212, and L1448C(N) -- and find that their inclination angle and jet velocity are consistent with those previously estimated from proper motion and radial velocity studies.

preprint2022arXiv

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): Evidence for a Molecular Jet Launched at an Unprecedented Early Phase of Protostellar evolution

Protostellar outflows and jets play a vital role in star formation as they carry away excess angular momentum from the inner disk surface, allowing the material to be transferred toward the central protostar. Theoretically, low velocity and poorly collimated outflows appear from the beginning of the collapse, at the first hydrostatic core (FHSC) stage. With growing protostellar core mass, high-density jets are launched which entrain an outflow from the infalling envelope. Until now, molecular jets have been observed at high velocity ($\gtrsim$ 100 km/s) in early Class\,0 protostars. We, for the first time, detect a dense molecular jet in SiO emission with small-velocity ($\sim$ 4.2 km\,s$^{-1}$, deprojected $\sim$ 24 km\,s$^{-1}$) from source G208.89-20.04Walma (hereafter, G208Walma) using ALMA Band\,6 observations. This object has some characteristics of FHSCs, such as a small outflow/jet velocity, extended 1.3\,mm continuum emission, and N$_2$D$^+$ line emission. Additional characteristics, however, are typical of early protostars: collimated outflow and SiO jet. The full extent of the outflow corresponds to a dynamical time scale of $\sim$ 930$^{+200}_{-100}$ years. The spectral energy distribution also suggests a very young source having an upper limit of T$_{bol}$ $\sim$ 31 K and L$_{bol}$ $\sim$ 0.8 L$_\sun$. We conclude that G208Walma is likely in the transition phase from FHSC to protostar, and the molecular jet has been launched within a few hundred years of initial collapse. Therefore, G208Walma may be the earliest object discovered in the protostellar phase with a molecular jet.

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

CCAT-prime: Design of the Mod-Cam receiver and 280 GHz MKID instrument module

Mod-Cam is a first light and commissioning instrument for the CCAT-prime project's six-meter aperture Fred Young Submillimeter Telescope (FYST), currently under construction at 5600 m on Cerro Chajnantor in Chile's Atacama Desert. Prime-Cam, a first-generation science instrument for FYST, will deliver over ten times greater mapping speed than current and near-term facilities for unprecedented 280-850 GHz broadband and spectroscopic measurements with microwave kinetic inductance detectors (MKIDs). CCAT-prime will address a suite of science goals, from Big Bang cosmology to star formation and galaxy evolution over cosmic time. Mod-Cam deployment on FYST with a 280 GHz instrument module containing MKID arrays is planned for early science observations in 2024. Mod-Cam will be used to test instrument modules for Prime-Cam, which can house up to seven instrument modules. We discuss the design and status of the 0.9 m diameter, 1.8 m long Mod-Cam receiver and 40 cm diameter 280 GHz instrument module, with cold stages at 40 K, 4 K, 1 K, and 100 mK. We also describe the instrument module's cryogenic readout designs to enable the readout of more than 10,000 MKIDs across 18 networks.

preprint2022arXiv

CCAT-prime: The 850 GHz camera for Prime-Cam on FYST

The Fred Young Submillimeter Telescope (FYST) at the Cerro-Chajnantor Atacama Telescope prime (CCAT- prime) Facility will host Prime-Cam as a powerful, first generation camera with imaging polarimeters working at several wavelengths and spectroscopic instruments aimed at intensity mapping during the Epoch of Reionization. Here we introduce the 850 GHz (350 micron) instrument module. This will be the highest frequency module in Prime-Cam and the most novel for astronomical and cosmological surveys, taking full advantage of the atmospheric transparency at the high 5600 meter CCAT-prime siting on Cerro Chajnantor. With a 1.1 deg diameter field, the 850 GHz module will deploy ~40,000 Kinetic Inductance Detectors (KIDs) with Silicon platelet feedhorn coupling (both fabricated at NIST), and will provide unprecedented broadband intensity and polarization measurement capabilities. The 850 GHz module will be key to addressing pressing astrophysical questions regarding galaxy formation, Big Bang cosmology, and star formation within our own Galaxy. We present the motivation and overall design for the module, and initial laboratory characterization.

preprint2022arXiv

Dissecting the different components of the modest accretion bursts of the very young protostar HOPS 373

Observed changes in protostellar brightness can be complicated to interpret. In our JCMT~Transient monitoring survey, we discovered that a young binary protostar, HOPS 373, is undergoing a modest $30\%$ brightness increase at 850 $μ$m, caused by a factor of 1.8$-$3.3 enhancement in the accretion rate. The initial burst occurred over a few months, with a sharp rise and then shallower decay. A second rise occurred soon after the decay, and the source is still bright one year later. The mid-IR emission, the small-scale CO outflow mapped with ALMA, and the location of variable maser emission indicate that the variability is associated with the SW component. The near-infrared and NEOWISE $W1$ and $W2$ emission is located along the blueshifted CO outflow, spatially offset by $\sim3$ to $4^{\prime\prime}$ from the SW component. The $K$-band emission imaged by UKIRT shows a compact H$_2$ emission source at the edge of the outflow, with a tail tracing the outflow back to the source. The $W1$ emission, likely dominated by scattered light, brightens by 0.7 mag, consistent with expectations based on the sub-mm lightcurve. The signal of continuum variability in $K$-band and $W2$ is masked by stable H$_2$ emission, as seen in our Gemini/GNIRS spectrum, and perhaps by CO emission. These differences in emission sources complicate infrared searches for variability of the youngest protostars.

preprint2022arXiv

FAUST III. Misaligned rotations of the envelope, outflow, and disks in the multiple protostellar system of VLA 1623$-$2417

We report a study of the low-mass Class-0 multiple system VLA 1623AB in the Ophiuchus star-forming region, using H$^{13}$CO$^+$ ($J=3-2$), CS ($J=5-4$), and CCH ($N=3-2$) lines as part of the ALMA Large Program FAUST. The analysis of the velocity fields revealed the rotation motion in the envelope and the velocity gradients in the outflows (about 2000 au down to 50 au). We further investigated the rotation of the circum-binary VLA 1623A disk as well as the VLA 1623B disk. We found that the minor axis of the circum-binary disk of VLA 1623A is misaligned by about 12 degrees with respect to the large-scale outflow and the rotation axis of the envelope. In contrast, the minor axis of the circum-binary disk is parallel to the large-scale magnetic field according to previous dust polarization observations, suggesting that the misalignment may be caused by the different directions of the envelope rotation and the magnetic field. If the velocity gradient of the outflow is caused by rotation, the outflow has a constant angular momentum and the launching radius is estimated to be $5-16$ au, although it cannot be ruled out that the velocity gradient is driven by entrainments of the two high-velocity outflows. Furthermore, we detected for the first time a velocity gradient associated with rotation toward the VLA 16293B disk. The velocity gradient is opposite to the one from the large-scale envelope, outflow, and circum-binary disk. The origin of its opposite gradient is also discussed.

preprint2022arXiv

Gap Opening and Inner Disk Structure in the Strongly Accreting Transition Disk of DM Tau

Large inner dust gaps in transition disks are frequently posited as evidence of giant planets sculpting gas and dust in the disk, or the opening of a gap by photoevaporative winds. Although the former hypothesis is strongly supported by the observations of planets and deep depletions in gas within the gap some disks, many T Tauri stars hosting transition disks accrete at rates typical for an undepleted disk, raising the question of how gap opening occurs in these objects. We thus present an analysis of the structure of the transition disk around the T Tauri star DM Tau, which is strongly accreting ($\sim 10^{-8.3}~\mathrm{M}_\odot~ \mathrm{yr}^{-1}$) and turbulent ($α=0.078 \pm 0.02$). Using the DALI thermochemical code, we fit disk models to simultaneously reproduce the accretion rate, high level of turbulence, the gas traced by ALMA band 6 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O J=2--1 lines, and the observed dust emission from the mm continuum and spectral energy distribution. We find a shallow depletion in gas surface density of $\sim 10$ relative to the outer disk and a gas rich inner disk is consistent with the observations. The planet mass of $<1$ M$_\mathrm{Jup}$ implied by the gap depth is in tension with predictions for dust trapping in a highly viscous disk, which requires a more massive planet of of $\sim10$M$_\mathrm{Jup}$. Photoevaporative models including a dead zone can qualitatively reproduce some features of the DM Tau disk, but still struggle to explain the high accretion rates and the observed mm continuum flux.

preprint2022arXiv

The Atacama Large Aperture Submillimeter Telescope: Key science drivers

The Atacama Large Aperture Submillimeter Telescope (AtLAST) is a concept for a 50m class single-dish telescope that will provide high sensitivity, fast mapping of the (sub-)millimeter sky. Expected to be powered by renewable energy sources, and to be constructed in the Atacama desert in the 2030s, AtLAST's suite of up to six state-of-the-art instruments will take advantage of its large field of view and high throughput to deliver efficient continuum and spectroscopic observations of the faint, large-scale emission that eludes current facilities. Here we present the key science drivers for the telescope characteristics, and discuss constraints that the transformational science goals place on future instrumentation.

preprint2022arXiv

Two Rings and a Marginally Resolved, 5 AU, Disk Around LkCa 15 Identified Via Near Infrared Sparse Aperture Masking Interferometry

Sparse aperture masking interferometry (SAM) is a high resolution observing technique that allows for imaging at and beyond a telescope's diffraction limit. The technique is ideal for searching for stellar companions at small separations from their host star; however, previous analysis of SAM observations of young stars surrounded by dusty disks have had difficulties disentangling planet and extended disk emission. We analyse VLT/SPHERE-IRDIS SAM observations of the transition disk LkCa\,15, model the extended disk emission, probe for planets at small separations, and improve contrast limits for planets. We fit geometrical models directly to the interferometric observables and recover previously observed extended disk emission. We use dynamic nested sampling to estimate uncertainties on our model parameters and to calculate evidences to perform model comparison. We compare our extended disk emission models against point source models to robustly conclude that the system is dominated by extended emission within 50 au. We report detections of two previously observed asymmetric rings at $\sim$17 au and $\sim$45 au. The peak brightness location of each model ring is consistent with the previous observations. We also, for the first time with imaging, robustly recover an elliptical Gaussian inner disk, previously inferred via SED fitting. This inner disk has a FWHM of ~5 au and a similar inclination and orientation as the outer rings. Finally, we recover no clear evidence for candidate planets. By modelling the extended disk emission, we are able to place a lower limit on the near infrared companion contrast of at least 1000.

preprint2021arXiv

CCAT-prime: Designs and status of the first light 280 GHz MKID array and Mod-Cam receiver

The CCAT-prime project's first light array will be deployed in Mod-Cam, a single-module testbed and first light cryostat, on the Fred Young Submillimeter Telescope (FYST) in Chile's high Atacama desert in late 2022. FYST is a six-meter aperture telescope being built on Cerro Chajnantor at an elevation of 5600 meters to observe at millimeter and submillimeter wavelengths.1 Mod-Cam will pave the way for Prime-Cam, the primary first generation instrument, which will house up to seven instrument modules to simultaneously observe the sky and study a diverse set of science goals from monitoring protostars to probing distant galaxy clusters and characterizing the cosmic microwave background (CMB). At least one feedhorn-coupled array of microwave kinetic inductance detectors (MKIDs) centered on 280 GHz will be included in Mod-Cam at first light, with additional instrument modules to be deployed along with Prime-Cam in stages. The first 280 GHz detector array was fabricated by the Quantum Sensors Group at NIST in Boulder, CO and includes 3,456 polarization-sensitive MKIDs. Current mechanical designs allow for up to three hexagonal arrays to be placed in each single instrument module. We present details on this first light detector array, including mechanical designs and cold readout plans, as well as introducing Mod-Cam as both a testbed and predecessor to Prime-Cam.

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

Quantifying Variability of YSOs in the Mid-IR Over Six Years with NEOWISE

Variability in young stellar objects (YSOs) can be caused by various time-dependent phenomena associated with star formation, including accretion rates, geometric changes in the circumstellar disks, stochastic hydromagnetic interactions between stellar surfaces and inner disk edges, reconnections within the stellar magnetosphere, and hot/cold spots on stellar surfaces. We uncover and characterize $\sim$1700 variables from a sample of $\sim$5400 YSOs in nearby low-mass star-forming regions using mid-IR light curves obtained from the 6.5-years NEOWISE All Sky Survey. The mid-IR variability traces a wide range of dynamical, physical, and geometrical phenomenon. We classify six types of YSO mid-IR variability based on their light curves: secular variability ($Linear, Curved, Periodic$) and stochastic variability ($Burst, Drop, Irregular$). YSOs in earlier evolutionary stages have higher fractions of variables and higher amplitudes for the variability, with the recurrence timescale of FUor-type outbursts (defined here as $Δ$W1 or $Δ$W2 $>1$ mag followed by inspection of candidates) of $\sim$1000 years in the early embedded protostellar phase. Known eruptive young stars and subluminous objects show fractions of variables similar to the fraction ($\sim55\%$) found in typical protostars, suggesting that these two distinct types are not distinct in variability over the 6.5-year timescale. Along with brightness variability, we also find a diverse range of secular color variations, which can be attributed to a competitive interplay between the variable accretion luminosity of the central source and the variable extinction by material associated with the accretion process.

preprint2020arXiv

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): I. Detection of New Hot Corinos with ACA

We report the detection of four new hot corino sources, G211.47-19.27S, G208.68-19.20N1, G210.49-19.79W and G192.12-11.10 from a survey study of Planck Galactic Cold Clumps in the Orion Molecular Cloud Complex with the Atacama Compact Array (ACA). Three sources had been identified as low mass Class 0 protostars in the Herschel Orion Protostar Survey (HOPS). One source in the lambda Orionis region is firstly reported as a protostellar core. We have observed abundant complex organic molecules (COMs), primarily methanol but also other oxygen-bearing COMs (in G211.47-19.27S and G208.68-19.20N1) and the molecule of prebiotic interest NH2CHO (in G211.47-19.27S), signifying the presence of hot corinos. While our spatial resolution is not sufficient for resolving most of the molecular emission structure, the large linewidth and high rotational temperature of COMs suggest that they likely reside in the hotter and innermost region immediately surrounding the protostar. In G211.47-19.27S, the D/H ratio of methanol ([CH2DOH]/[CH3OH]) and the 12C/13C ratio of methanol ([CH3OH]/[13CH3OH]) are comparable to those of other hot corinos. Hydrocarbons and long carbon-chain molecules such as c-C3H2 and HCCCN are also detected in the four sources, likely tracing the outer and cooler molecular envelopes.

preprint2020arXiv

Dual-Wavelength ALMA Observations of Dust Rings in Protoplanetary Disks

We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations for three protoplanetary disks in Taurus at 2.9\,mm and comparisons with previous 1.3\,mm data both at an angular resolution of $\sim0.''1$ (15\,au for the distance of Taurus). In the single-ring disk DS Tau, double-ring disk GO Tau, and multiple-ring disk DL Tau, the same rings are detected at both wavelengths, with radial locations spanning from 50 to 120\,au. To quantify the dust emission morphology, the observed visibilities are modeled with a parametric prescription for the radial intensity profile. The disk outer radii, taken as 95\% of the total flux encircled in the model intensity profiles, are consistent at both wavelengths for the three disks. Dust evolution models show that dust trapping in local pressure maxima in the outer disk could explain the observed patterns. Dust rings are mostly unresolved. The marginally resolved ring in DS Tau shows a tentatively narrower ring at the longer wavelength, an observational feature expected from efficient dust trapping. The spectral index ($α_{\rm mm}$) increases outward and exhibits local minima that correspond to the peaks of dust rings, indicative of the changes in grain properties across the disks. The low optical depths ($τ\sim$0.1--0.2 at 2.9\,mm and 0.2--0.4 at 1.3\,mm) in the dust rings suggest that grains in the rings may have grown to millimeter sizes. The ubiquitous dust rings in protoplanetary disks modify the overall dynamics and evolution of dust grains, likely paving the way towards the new generation of planet formation.

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

Radiative Transfer modeling of EC 53: An Episodically Accreting Class I Young Stellar Object

In the episodic accretion scenario, a large fraction of the protostellar mass accretes during repeated and large bursts of accretion. Since outbursts on protostars are typically identified at specific wavelengths, interpreting these outbursts requires converting this change in flux to a change in total luminosity. The Class I young stellar object EC 53 in the Serpens Main cloud has undergone repeated increases in brightness at 850 $μ$m that are likely caused by bursts of accretion. In this study, we perform two- and three-dimensional continuum radiative transfer modeling to quantify the internal luminosity rise in EC 53 that corresponds to the factor of $\sim$1.5 enhancement in flux at 850 $μ$m. We model the spectral energy distribution and radial intensity profile in both the quiescent and outburst phases. The internal luminosity in the outburst phase is $\sim 3.3$ times brighter than the luminosity in the quiescent phase. The radial intensity profile analysis demonstrates that the detected sub-mm flux variation of EC 53 comes from the heated envelope by the accretion burst. We also find that the role of external heating of the EC 53 envelope by the interstellar radiation field is insignificant.

preprint2020arXiv

Steady Wind-Blown Cavities within Infalling Rotating Envelopes: Application to the Broad Velocity Component in Young Protostars

Wind-driven outflows are observed around a broad range of accreting objects throughout the Universe, ranging from forming low-mass stars to super-massive black holes. We study the interaction between a central isotropic wind and an infalling, rotating, envelope, determining the steady-state cavity shape formed at their interface under the assumption of weak mixing. The shape of the resulting wind-blown cavity is elongated and self-similar, with a physical size determined by the ratio between wind ram pressure and envelope thermal pressure. We compute the growth of a warm turbulent mixing-layer between the shocked wind and the deflected envelope, and calculate the resultant broad line profile, under the assumption of a linear (Couette-type) velocity profile across the layer. We then test our model against the warm broad velocity component observed in CO $J$=16--15 by Herschel/HIFI in the protostar Serpens-Main SMM1. Given independent observational constraints on the temperature and density of the dust envelope around SMM1, we find an excellent match to all its observed properties (line profile, momentum, temperature) and to the SMM1 outflow cavity width for a physically reasonable set of parameters: a ratio of wind to infall mass-flux $\simeq 4\%$, a wind speed $v_{\rm w} \simeq 30$ km/s, an interstellar abundance of CO and H$_2$, and a turbulent entrainment efficiency consistent with laboratory experiments. The inferred ratio of ejection to disk accretion rate, $\simeq 6-20\%$, is in agreement with current disk wind theories. Thus, the model provides a new framework to reconcile the modest outflow cavity widths in protostars with the large observed flow velocities. Being self-similar, it is applicable over a broader range of astrophysical contexts as well.

preprint2020arXiv

The Relationship between Mid-Infrared and Sub-Millimetre Variability of Deeply Embedded Protostars

We study the relationship between the mid-infrared and sub-mm variability of deeply embedded protostars using the multi-epoch data from the Wide Infrared Survey Explorer ($WISE$/NEOWISE) and the ongoing James Clerk Maxwell Telescope (JCMT) transient survey. Our search for signs of stochastic (random) and/or secular (roughly monotonic in time) variability in a sample of 59 young stellar objects (YSOs) revealed that 35 are variable in at least one of the two surveys. This variability is dominated by secular changes. Of those objects with secular variability, 14 objects ($22\%$ of the sample) show correlated secular variability over mid-IR and sub-mm wavelengths. Variable accretion is the likely mechanism responsible for this type of variability. Fluxes of YSOs that vary in both wavelengths follow a relation of $\log_{10} F_{4.6}(t)=η\log_{10} F_{850}(t)$ between the mid-IR and sub-mm, with $η=5.53\pm0.29$. This relationship arises from the fact that sub-mm fluxes respond to the dust temperature in the larger envelope whereas the mid-IR emissivity is more directly proportional to the accretion luminosity. The exact scaling relation, however, depends on the structure of the envelope, the importance of viscous heating in the disc, and dust opacity laws.

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

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 circumstellar environment around the embedded protostar EC 53

EC53 is an embedded protostar with quasi-periodic emission in the near-IR and sub-mm. We use ALMA high-resolution observations of continuum and molecular line emission to describe the circumstellar environment of EC 53. The continuum image reveals a disk with a flux that suggests a mass of 0.075 Msun, much less than the estimated mass in the envelope, and an in-band spectral index that indicates grain growth to centimeter sizes. Molecular lines trace the outflow cavity walls, infalling and rotating envelope, and/or the Keplerian disk. The rotation profile of the C17O 3-2 line emission cannot isolate the Keplerian motion clearly although the lower limit of the protostellar mass can be calculated as 0.3 +- 0.1 Msun if the Keplerian motion is adopted. The weak CH3OH emission, which is anti-correlated with the HCO+ 4-3 line emission, indicates that the water snow line is more extended than what expected from the current luminosity, attesting to bygone outburst events. The extended snow line may persist for longer at the disk surface because the lower density increases the freeze-out timescale of methanol and water.

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.

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

The Eruption of the Candidate Young Star ASASSN-15qi

Outbursts on young stars are usually interpreted as accretion bursts caused by instabilities in the disk or the star-disk connection. However, some protostellar outbursts may not fit into this framework. In this paper, we analyze optical and near-infrared spectra and photometry to characterize the 2015 outburst of the probable young star ASASSN-15qi. The $\sim 3.5$ mag brightening in the $V$ band was sudden, with an unresolved rise time of less than one day. The outburst decayed exponentially by 1 mag for 6 days and then gradually back to the pre-outburst level after 200 days. The outburst is dominated by emission from $\sim10,000$ K gas. An explosive release of energy accelerated matter from the star in all directions, seen in a spectacular cool, spherical wind with a maximum velocity of 1000 km/s. The wind and hot gas both disappeared as the outburst faded and the source returned to its quiescent F-star spectrum. Nebulosity near the star brightened with a delay of 10-20 days. Fluorescent excitation of H$_2$ is detected in emission from vibrational levels as high as $v=11$, also with a possible time delay in flux increase. The mid-infrared spectral energy distribution does not indicate the presence of warm dust emission, although the optical photospheric absorption and CO overtone emission could be related to a gaseous disk. Archival photometry reveals a prior outburst in 1976. Although we speculate about possible causes for this outburst, none of the explanations are compelling.

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

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

Mid-J CO observations of Perseus B1-East 5: evidence for turbulent dissipation via low-velocity shocks

Giant molecular clouds contain supersonic turbulence and magnetohydrodynamic simulations predict that this turbulence should decay rapidly. Such turbulent dissipation has the potential to create a warm (T ~100 K) gas component within a molecular cloud. We present observations of the CO J = 5-4 and 6-5 transitions, taken with the Herschel Space Observatory, towards the Perseus B1-East 5 region. We combine these new observations with archival measurements of lower rotational transitions and fit photodissociation region models to the data. We show that Perseus B1-E5 has an anomalously large CO J = 6-5 integrated intensity, consistent with a warm gas component existing within the region. This excess emission is consistent with predictions for shock heating due to the dissipation of turbulence in low velocity shocks with the shocks having a volume filling factor of 0.15 per cent. We find that B1-E has a turbulent energy dissipation rate of 3.5 x 10$^{32}$ erg / s and a dissipation time-scale that is only a factor of 3 larger than the flow crossing time-scale.

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

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

Kompaneets Model Fitting of the Orion-Eridanus Superbubble

Winds and supernovae from OB associations create large cavities in the interstellar medium referred to as superbubbles. The Orion molecular clouds are the nearest high mass star-forming region and have created a highly elongated, 20 degree x 45 degree, superbubble. We fit Kompaneets models to the Orion-Eridanus superbubble and find that a model where the Eridanus side of the superbubble is oriented away from the Sun provides a marginal fit. Because this model requires an unusually small scaleheight of 40 pc and has the superbubble inclined 35 degrees from the normal to the Galactic plane, we propose that this model should be treated as a general framework for modelling the Orion-Eridanus superbubble, with a secondary physical mechanism not included in the Kompaneets model required to fully account for the orientation and elongation of the superbubble.

preprint2014arXiv

Observatory/data centre partnerships and the VO-centric archive: The JCMT Science Archive experience

We present, as a case study, a description of the partnership between an observatory (JCMT) and a data centre (CADC) that led to the development of the JCMT Science Archive (JSA). The JSA is a successful example of a service designed to use Virtual Observatory (VO) technologies from the start. We describe the motivation, process and lessons learned from this approach.

preprint2014arXiv

Synthetic Observations of the Evolution of Starless Cores in a Molecular Cloud Simulation: Comparisons with JCMT Data and Predictions for ALMA

Interpreting the nature of starless cores has been a prominent goal in star formation for many years. In order to characterise the evolutionary stages of these objects, we perform synthetic observations of a numerical simulation of a turbulent molecular cloud. We find that nearly all cores that we detect are associated with filaments and eventually form protostars. We conclude that observed starless cores which appear Jeans unstable are only marginally larger than their respective Jeans masses (within a factor of 3). We note single dish observations such as those performed with the JCMT appear to miss significant core structure on small scales due to beam averaging. Finally, we predict that interferometric observations with ALMA Cycle 1 will resolve the important small scale structure, which has so far been missed by mm-wavelength observations.

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.

preprint2014arXiv

The Origin of Ionized Filaments Within the Orion-Eridanus Superbubble

The Orion-Eridanus superbubble, formed by the nearby Orion high mass star-forming region, contains multiple bright H$α$ filaments on the Eridanus side of the superbubble. We examine the implications of the H$α$ brightnesses and sizes of these filaments, the Eridanus filaments. We find that either the filaments must be highly elongated along the line of sight or they cannot be equilibrium structures illuminated solely by the Orion star-forming region. The Eridanus filaments may, instead, have formed when the Orion-Eridanus superbubble encountered and compressed a pre-existing, ionized gas cloud, such that the filaments are now out of equilibrium and slowly recombining.

preprint2013arXiv

Continuum Variability of Deeply Embedded Protostars as a Probe of Envelope Structure

Stars may be assembled in large growth spurts, however the evidence for this hypothesis is circumstantial. Directly studying the accretion at the earliest phases of stellar growth is challenging because young stars are deeply embedded in optically thick envelopes, which have spectral energy distributions that peak in the far-IR, where observations are difficult. In this paper, we consider the feasibility of detecting accretion outbursts from these younger stars by investigating the timescales for how the protostellar envelope responds to changes in the emission properties of the central source. The envelope heats up in response to an outburst, brightening at all wavelengths and with the emission peak moving to shorter wavelengths. The timescale for this change depends on the time for dust grains to heat and re-emit photons and the time required for the energy to escape the inner, optically-thick portion of the envelope. We find that the dust response time is much shorter than the photon propagation time and thus the timescale over which the emission varies is set by time delays imposed by geometry. These times are hours to days near the peak of the spectral energy distribution and weeks to months in the sub-mm. The ideal location to quickly detect continuum variability is therefore in the mid- to far-IR, near the peak of the spectral energy distribution, where the change in emission amplitude is largest. Searching for variability in sub-mm continuum emission is also feasible, though with a longer time separation and a weaker relationship between the amount of detected emission amplitude and change in central source luminosity. Such observations would constrain accretion histories of protostars and would help to trace the disk/envelope instabilities that lead to stellar growth.

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

The Star Formation Rate of Molecular Clouds

We review recent advances in the analytical and numerical modeling of the star formation rate in molecular clouds and discuss the available observational constraints. We focus on molecular clouds as the fundamental star formation sites, rather than on the larger-scale processes that form the clouds and set their properties. Molecular clouds are shaped into a complex filamentary structure by supersonic turbulence, with only a small fraction of the cloud mass channeled into collapsing protostars over a free-fall time of the system. In recent years, the physics of supersonic turbulence has been widely explored with computer simulations, leading to statistical models of this fragmentation process, and to the prediction of the star formation rate as a function of fundamental physical parameters of molecular clouds, such as the virial parameter, the rms Mach number, the compressive fraction of the turbulence driver, and the ratio of gas to magnetic pressure. Infrared space telescopes, as well as ground-based observatories have provided unprecedented probes of the filamentary structure of molecular clouds and the location of forming stars within them.

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

Aspect Ratio Dependence of the Free-Fall Time for Non-Spherical Symmetries

We investigate the collapse of non-spherical substructures, such as sheets and filaments, which are ubiquitous in molecular clouds. Such non-spherical substructures collapse homologously in their interiors but are influenced by an edge effect that causes their edges to be preferentially accelerated. We analytically compute the homologous collapse timescales of the interiors of uniform-density, self-gravitating filaments and find that the homologous collapse timescale scales linearly with the aspect ratio. The characteristic timescale for an edge driven collapse mode in a filament, however, is shown to have a square root dependence on the aspect ratio. For both filaments and circular sheets, we find that selective edge acceleration becomes more important with increasing aspect ratio. In general, we find that lower dimensional objects and objects with larger aspect ratios have longer collapse timescales. We show that estimates for star formation rates, based upon gas densities, can be overestimated by an order of magnitude if the geometry of a cloud is not taken into account.

preprint2012arXiv

Filamentary Star Formation: Observing the Evolution toward Flattened Envelopes

Filamentary structures are ubiquitous from large-scale molecular clouds (few parsecs) to small-scale circumstellar envelopes around Class 0 sources (~1000 AU to ~0.1 pc). In particular, recent observations with the Herschel Space Observatory emphasize the importance of large-scale filaments (few parsecs) and star formation. The small-scale flattened envelopes around Class 0 sources are reminiscent of the large-scale filaments. We propose an observationally derived scenario for filamentary star formation that describes the evolution of filaments as part of the process for formation of cores and circumstellar envelopes. If such a scenario is correct, small-scale filamentary structures (0.1 pc in length) with higher densities embedded in starless cores should exist, although to date almost all the interferometers have failed to observe such structures. We perform synthetic observations of filaments at the prestellar stage by modeling the known Class 0 flattened envelope in L1157 using both the Combined Array for Research in Millimeter-wave Astronomy (CARMA) and the Atacama Large Millimeter/Submillimeter Array (ALMA). We show that with reasonable estimates for the column density through the flattened envelope, the CARMA D-array at 3mm wavelengths is not able to detect such filamentary structure, so previous studies would not have detected them. However, the substructures may be detected with CARMA D+E array at 3 mm and CARMA E array at 1 mm as a result of more appropriate resolution and sensitivity. ALMA is also capable of detecting the substructures and showing the structures in detail compared to the CARMA results with its unprecedented sensitivity. Such detection will confirm the new proposed paradigm of non-spherical star formation.

preprint2012arXiv

The Small-Scale Physical Structure and Fragmentation Difference of Two Embedded Intermediate Mass Protostars in Orion

Intermediate mass protostars, the bridge between the very common solar-like protostars and the more massive, but rarer, O and B stars, can only be studied at high physical spatial resolutions in a handful of clouds. In this paper we present and analyze the continuum results from an observing campaign at the Submillimeter Array targeting two well-studied intermediate mass protostars in Orion, NGC 2071 and L1641 S3 MMS 1. The extended SMA (eSMA) probes structure at angular resolutions up to 0.2", revealing protostellar disks on scales of 200 AU. Continuum flux measurements on these scales indicate that a significant amount of mass, a few tens of M{\odot}, are present. Envelope, stellar, and disk masses are derived using both compact, extended and eSMA configurations and compared against SED-fitting models. We hypothesize that fragmentation into three components occurred within NGC 2071 at an early time, when the envelopes were less than 10% of their current masses, e.g. < 0.5 M{\odot}. No fragmentation occurred for L1641 S3 MMS 1. For NGC 2071 evidence is given that the bulk of the envelope material currently around each source was accreted after the initial fragmentation. In addition, about 30% of the total core mass is not yet associated to one of the three sources. A global accretion model is favored and a potential accretion history of NGC 2071 is presented. It is shown that the relatively low level of fragmentation in NGC 2071 was stifled compared to the expected fragmentation from a Jean's argument.

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.

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

An Observed Lack of Substructure in Starless Cores

In this paper we present the results of a high resolution (5") CARMA and SZA survey of the 3mm continuum emission from 11 of the brightest (at 1.1mm) starless cores in the Perseus molecular cloud. We detect 2 of the 11 cores, both of which are composed of single structures, and the median 3 sigma upper limit for the non-detections is 0.2 M_sun in a 5" beam. These results are consisent with, and as stringent as, the low detection rate of compact 3mm continuum emission in dense cores in Perseus reported by Olmi et al. (2005). From the non-detection of multiple components in any of the eleven cores we conclude that starless core mass functions derived from bolometer maps at resolutions from 10"-30" (e.g. with MAMBO, SCUBA or Bolocam) are unlikely to be significantly biased by the blending of lower mass cores with small separations. These observations provide additional evidence that the majority of starless cores in Perseus have inner density profiles shallower than r^-2.

preprint2010arXiv

Dense Gas Tracers in Perseus: Relating the N2H+, NH3, and Dust Continuum Properties of Pre- and Proto-Stellar Cores

We investigate 35 pre-stellar cores and 36 proto-stellar cores in the Perseus molecular cloud. We find a very tight correlation between the physical parameters describing the N2H+ and NH3 gas. Both the velocity centroids and the line widths of N2H+ and NH3 correlate much better than either species correlates with CO, as expected if the nitrogen-bearing species are probing primarily the dense core gas where the CO has been depleted. We also find a tight correlation in the inferred abundance ratio between N2H+ and para-NH3 across all cores, with N(p-NH3)/N(N2H+)= 22 +/- 10. We find a mild correlation between NH3 (and N2H+) column density and the (sub)millimeter dust continuum derived H2 column density for pre-stellar cores, N(p-NH3)/N(H2) ~ 10e-8, but do not find a fixed ratio for proto-stellar cores. The observations suggest that in the Perseus molecular cloud the formation and destruction mechanisms for the two nitrogen-bearing species are similar, regardless of the physical conditions in the dense core gas. While the equivalence of N2H+ and NH3 as powerful tracers of dense gas is validated, the lack of correspondence between these species and the (sub)millimeter dust continuum observations for proto-stellar cores is disconcerting and presently unexplained.

preprint2010arXiv

The Dynamics of Dense Cores in the Perseus Molecular Cloud II: The Relationship Between Dense Cores and the Cloud

We utilize the extensive datasets available for the Perseus molecular cloud to analyze the relationship between the kinematics of small-scale dense cores and the larger structures in which they are embedded. The kinematic measures presented here can be used in conjunction with those discussed in our previous work as strong observational constraints that numerical simulations (or analytic models) of star formation should match. We find that dense cores have small motions with respect to the 13CO gas, about one third of the 13CO velocity dispersion along the same line of sight. Within each extinction region, the core-to-core velocity dispersion is about half of the total (13CO) velocity dispersion seen in the region. Large-scale velocity gradients account for roughly half of the total velocity dispersion in each region, similar to what is predicted from large-scale turbulent modes following a power spectrum of P(k) ~ k^{-4}.