Source author record

Paul F. Goldsmith

Paul F. Goldsmith 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

47works
5topics
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

47 published item(s)

preprint2023arXiv

Atomic oxygen abundance toward Sagittarius B2

A substantial fraction of oxygen in diffuse clouds is unaccounted for by observations and is postulated to be in an unknown refractory form, referred to as unidentified depleted oxygen (UDO), which, depending on the local gas density, may contribute up to 50% of the total oxygen content. Previous Infrared Space Observatory (ISO) observations suggest that a significant fraction of oxygen in even denser, translucent clouds may be in atomic form. We have analyzed velocity-resolved archival SOFIA observations of the 63 $μ$m fine-structure [O I] transition toward the high-mass star-forming region Sgr B2(M) in the Central Molecular Zone. The foreground spiral-arm clouds as well as the extended Sgr B2 envelope between the Sun and the background dust continuum source produce multiple [O i] absorption components, spectrally separated in velocity space. The gas-phase atomic oxygen column density in foreground clouds toward Sgr B2 is well correlated with the total hydrogen column density, with an average atomic oxygen abundance of $(2.51 \pm 0.69) \times 10^{-4}$ with respect to hydrogen nuclei. This value is in good agreement with the earlier ISO measurements on the same line of sight, and is about 35% lower than the total interstellar medium oxygen abundance in the low-density warm gas, as measured in the UV. We find no evidence that a significant fraction of the oxygen on the line of sight toward Sagittarius B2 is in the form of UDO.

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

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- VII. A catalogue of SiO clumps from ACA observations

To understand the nature of SiO emission, we conducted ACA observations of the SiO (2-1) lines toward 146 massive star-forming regions, as part of the ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS) survey. We detected SiO emission in 128 (87.7$\%$) sources and identified 171 SiO clumps, 105 of which are spatially separated from 3 mm continuum emission. A large amount of the SiO line profiles (60$\%$) are non-Gaussian. The velocity dispersion of the SiO lines ranges from 0.3 to 5.43 km s$^{-1}$. In 63 sources the SiO clumps are associated with H$_\rm{II}$ regions characterized by H40$α$ emission. We find that 68$\%$ (116) of the SiO clumps are associated with strong outflows. The median velocity dispersion of the SiO line for outflow sources and non-outflow sources is 1.91 km s$^{-1}$ and 0.99 km s$^{-1}$, respectively. These results indicate that outflow activities could be connected to strongly shocked gas. The velocity dispersion and [SiO]/[H$^{13}$CO$^+$] intensity ratio do not show any correlation with the dust temperature and particle number density of clumps. We find a positive correlation between the SiO line luminosity and the bolometric luminosity, implying stronger shock activities are associated with more luminous proto-clusters. The SiO clumps in associations with H$_\rm{II}$ regions were found to show a steeper feature in $L_\rm{sio}$/$L_\rm{bol}$. The SiO line luminosity and the fraction of shocked gas have no apparent evidence of correlation with the evolutionary stages traced by luminosity to mass ratio ($L_\rm{bol}/M$).

preprint2022arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- VIII. A search for hot cores by using C$_2$H$_5$CN, CH$_3$OCHO and CH$_3$OH lines

Hot cores characterized by rich lines of complex organic molecules are considered as ideal sites for investigating the physical and chemical environments of massive star formation. We present a search for hot cores by using typical nitrogen- and oxygen-bearing complex organic molecules (C$_2$H$_5$CN, CH$_3$OCHO and CH$_3$OH), based on ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS). The angular resolutions and line sensitivities of the ALMA observations are better than 2 arcsec and 10 mJy/beam, respectively. A total of 60 hot cores are identified with 45 being newly detected, in which the complex organic molecules have high gas temperatures ($>$ 100 K) and small source sizes ($<$ 0.1 pc). So far this is the largest sample of hot cores observed with similar angular resolution and spectral coverage. The observations have also shown nitrogen and oxygen differentiation in both line emission and gas distribution in 29 hot cores. Column densities of CH$_3$OH and CH$_3$OCHO increase as rotation temperatures rise. The column density of CH$_3$OCHO correlates tightly with that of CH$_3$OH. The pathways for production of different species are discussed. Based on the spatial position difference between hot cores and UC~H{\sc ii} regions, we conclude that 24 hot cores are externally heated while the other hot cores are internally heated. The observations presented here will potentially help establish a hot core template for studying massive star formation and astrochemistry.

preprint2022arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- XI. From inflow to infall in hub-filament systems

We investigate the presence of hub-filament systems in a large sample of 146 active proto-clusters, using H$^{13}$CO$^{+}$ J=1-0 molecular line data obtained from the ATOMS survey. We find that filaments are ubiquitous in proto-clusters, and hub-filament systems are very common from dense core scales ($\sim$0.1 pc) to clump/cloud scales ($\sim$1-10 pc). The proportion of proto-clusters containing hub-filament systems decreases with increasing dust temperature ($T_d$) and luminosity-to-mass ratios ($L/M$) of clumps, indicating that stellar feedback from H{\sc ii} regions gradually destroys the hub-filament systems as proto-clusters evolve. Clear velocity gradients are seen along the longest filaments with a mean velocity gradient of 8.71 km s$^{-1}$pc$^{-1}$ and a median velocity gradient of 5.54 km s$^{-1}$pc$^{-1}$. We find that velocity gradients are small for filament lengths larger than $\sim$1~pc, probably hinting at the existence of inertial inflows, although we cannot determine whether the latter are driven by large-scale turbulence or large-scale gravitational contraction. In contrast, velocity gradients below $\sim$1~pc dramatically increase as filament lengths decrease, indicating that the gravity of the hubs or cores starts to dominate gas infall at small scales. We suggest that self-similar hub-filament systems and filamentary accretion at all scales may play a key role in high-mass star formation.

preprint2022arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- XII: Fragmentation and multi-scale gas kinematics in protoclusters G12.42+0.50 and G19.88-0.53

We present new continuum and molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS) survey for the two protoclusters, G12.42+0.50 and G19.88-0.53. The 3 mm continuum maps reveal seven cores in each of the two globally contracting protoclusters. These cores satisfy the radius-mass relation and the surface mass density criteria for high-mass star formation. Similar to their natal clumps, the virial analysis of the cores suggests that they are undergoing gravitational collapse ($\rm α_{vir} << 2$). The clump to core scale fragmentation is investigated and the derived core masses and separations are found to be consistent with thermal Jeans fragmentation. We detect large-scale filamentary structures with velocity gradients and multiple outflows in both regions. Dendrogram analysis of the H$^{13}$CO$^{+}$ map identifies several branch and leaf structures with sizes $\sim$ 0.1 and 0.03 pc, respectively. The supersonic gas motion displayed by the branch structures is in agreement with the Larson power-law indicating that the gas kinematics at this spatial scale is driven by turbulence. The transition to transonic/subsonic gas motion is seen to occur at spatial scales of $\sim$0.1 pc indicating the dissipation of turbulence. In agreement with this, the leaf structures reveal gas motions that deviate from the slope of Larson's law. From the large-scale converging filaments to the collapsing cores, the gas dynamics in G12.42+0.50 and G19.88-0.53 show scale-dependent dominance of turbulence and gravity and the combination of these two driving mechanisms needs to be invoked to explain massive star formation in the protoclusters.

preprint2022arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions-IX. A pilot study towards IRDC G034.43+00.24 on multi-scale structures and gas kinematics

We present a comprehensive study of the gas kinematics associated with density structures at different spatial scales in the filamentary infrared dark cloud, G034.43+00.24 (G34). This study makes use of the H13CO+ (1-0) molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS) survey, which has spatial and velocity resolution of 0.04 pc and 0.2 km/s, respectively. Several tens of dendrogram structures have been extracted in the position-position-velocity space of H13CO+, which include 21 small-scale leaves and 20 larger-scale branches. Overall, their gas motions are supersonic but they exhibit the interesting behavior where leaves tend to be less dynamically supersonic than the branches. For the larger-scale, branch structures, the observed velocity-size relation (i.e., velocity variation/dispersion versus size) are seen to follow the Larson scaling exponent while the smaller-scale, leaf structures show a systematic deviation and display a steeper slope. We argue that the origin of the observed kinematics of the branch structures is likely to be a combination of turbulence and gravity-driven ordered gas flows. In comparison, gravity-driven chaotic gas motion is likely at the level of small-scale leaf structures. The results presented in our previous paper and this current follow-up study suggest that the main driving mechanism for mass accretion/inflow observed in G34 varies at different spatial scales. We therefore conclude that a scale-dependent combined effect of turbulence and gravity is essential to explain the star-formation processes in G34.

preprint2022arXiv

The Distribution of UV Radiation Field in the Molecular Clouds of Gould Belt

The distribution of ultraviolet (UV) radiation field provides critical constraints on the physical environments of molecular clouds. Within 1 kpc of our solar system and fostering protostars of different masses, the giant molecular clouds in the Gould Belt present an excellent opportunity to resolve the UV field structure in star forming regions. We performed spectral energy distribution (SED) fitting of the archival data from the Herschel Gould Belt Survey (HGBS). Dust radiative transfer analysis with the DUSTY code were applied to 23 regions in 14 molecular complexes of the Gould Belt, resulting in the spatial distribution of radiation field in these regions. For 10 of 15 regions with independent measurements of star formation rate, their star formation rate and UV radiation intensity largely conform to a linear correlation found in previous studies.

preprint2020arXiv

A SOFIA Survey of [CII] in the galaxy M51 II. [CII] and CO kinematics across spiral arms

We present the first complete, velocity-resolved [CII] 158um image of the M51 grand-design spiral galaxy, observed with the upGREAT instrument on SOFIA. [CII] is an important tracer of various phases of the interstellar medium (ISM), including ionized gas, neutral atomic, and diffuse molecular regions. We combine the [CII] data with HI, CO, 24um dust continuum, FUV, and near-infrared K-band observations to study the evolution of the ISM across M51's spiral arms in both position-position, and position-velocity space. Our data show strong velocity gradients in HI, 12CO, and [CII] at the locations of stellar arms (traced by K--band data) with a clear offset in position-velocity space between upstream molecular gas (traced by 12CO) and downstream star formation (traced by [CII]). We compare the observed position--velocity maps across spiral arms with synthetic observations from numerical simulations of galaxies with both dynamical and quasi-stationary steady spiral arms that predict both tangential and radial velocities at the location of spiral arms. We find that our observations, based on the observed velocity gradients and associated offset between CO and [CII], are consistent with the presence of shocks in spiral arms in the inner parts of M51 and in the arm connecting the companion galaxy, M51b, in the outer parts of M51.

preprint2020arXiv

Applications of Machine Learning Algorithms In Processing Terahertz Spectroscopic Data

We present the data reduction software and the distribution of Level 1 and Level 2 products of the Stratospheric Terahertz Observatory 2 (STO2). STO2, a balloon-borne Terahertz telescope, surveyed star-forming regions and the Galactic plane and produced approximately 300,000 spectra. The data are largely similar to spectra typically produced by single-dish radio telescopes. However, a fraction of the data contained rapidly varying fringe/baseline features and drift noise, which could not be adequately corrected using conventional data reduction software. To process the entire science data of the STO2 mission, we have adopted a new method to find proper off-source spectra to reduce large-amplitude fringes and new algorithms including Asymmetric Least Square (ALS), Independent Component Analysis (ICA), and Density-based spatial clustering of applications with noise (DBSCAN). The STO2 data reduction software efficiently reduced the amplitude of fringes from a few hundred to 10 K and resulted in baselines of amplitude down to a few K. The Level 1 products typically have the noise of a few K in [CII] spectra and ~1 K in [NII] spectra. Using a regridding algorithm, we made spectral maps of star-forming regions and the Galactic plane survey using an algorithm employing a Bessel-Gaussian kernel. Level 1 and 2 products are available to the astronomical community through the STO2 data server and the DataVerse. The software is also accessible to the public through Github. The detailed addresses are given in Section 4 of the paper on data distribution.

preprint2020arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- I. Survey description and a first look at G9.62+0.19

The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the "ATOMS" survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS $J = 2-1$, HCO$^+$ $J = 1-0$ and HCN $J = 1-0$, are found to show extended gas emission in low density regions within the clump; less than 25\% of their emission is from dense cores. SO, CH$_3$OH, H$^{13}$CN and HC$_3$N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over $\sim$1 pc), which may be caused by slow shocks from large--scale colliding flows or H{\sc ii} regions. Stellar feedback from an expanding H{\sc ii} region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analyzed with other survey data, e.g., "MALT90", "Orion B", "EMPIRE", "ALMA\_IMF", and "ALMAGAL", to deepen our understandings of "dense gas" star formation scaling relations and massive proto-cluster formation.

preprint2020arXiv

Distribution of Water Vapor in Molecular Clouds. II

The depth-dependent abundance of both gas-phase and solid-state water within dense, quiescent, molecular clouds is important to both the cloud chemistry and gas cooling. Where water is in the gas phase, it's free to participate in the network of ion-neutral reactions that lead to a host of oxygen-bearing molecules, and its many ortho and para energy levels make it an effective coolant for gas temperatures greater than 20K. Where water is abundant as ice on grain surfaces, and unavailable to cool the gas, significant amounts of oxygen are removed from the gas phase, suppressing the gas-phase chemical reactions that lead to a number of oxygen-bearing species, including O2. Models of FUV-illuminated clouds predict that the gas-phase water abundance peaks in the range Av ~3 and 8mag of the cloud surface, depending on the gas density and FUV field strength. Deeper within such clouds, water is predicted to exist mainly as ice on grain surfaces. More broadly, these models are used to analyze a variety of other regions, including outflow cavities associated with young stellar objects and the surface layers of protoplanetary disks. In this paper, we report the results of observational tests of FUV-illuminated cloud models toward the Orion Molecular Ridge and Cepheus B using data obtained from the Herschel Space Observatory and the Five College Radio Astronomy Observatory. Toward Orion, 2220 spatial positions were observed along the face-on Orion Ridge in the H2O 110-101 557GHz and NH3 J,K=1,0-0,0 572GHz lines. Toward Cepheus B, two strip scans were made in the same lines across the edge-on ionization front. These new observations demonstrate that gas-phase water exists primarily within a few magnitudes of dense cloud surfaces, strengthening the conclusions of an earlier study based on a much smaller data set, and indirectly supports the prediction that water ice is quite abundant in dense clouds.

preprint2020arXiv

Magnetic Field Morphology in Interstellar Clouds with the Velocity Gradient Technique

Magnetic fields, while ubiquitous in many astrophysical environments, are challenging to measure observationally. Based on the properties of anisotropy of eddies in magnetized turbulence, the Velocity Gradient Technique is a method synergistic to dust polarimetry that is capable of tracing plane-of-the-sky magnetic field, measuring the magnetization of interstellar media and estimating the fraction of gravitational collapsing gas in molecular clouds using spectral line observations. In this paper, we apply this technique to five low-mass star-forming molecular clouds in the Gould Belt and compare the results to the magnetic-field orientation obtained from polarized dust emission. We find the estimates of magnetic field orientations and magnetization for both methods are statistically similar. We estimate the fraction of collapsing gas in the selected clouds. By means of the Velocity Gradient Technique, we also present the plane-of-the-sky magnetic field orientation and magnetization of the Smith cloud, for which dust polarimetry data are unavailable.

preprint2020arXiv

Molecular Oxygen in the nearest QSO Mrk 231

We report the detection of an emission feature at the 12 sigma level with FWHM line width of about 450 km/s toward the nearest quasi-stellar object, QSO Mrk 231. Based on observations with the IRAM 30 m telescope and the NOEMA Interferometer, the 11-10 transition of molecular oxygen is the likely origin of line with rest frequency close to 118.75 GHz. The velocity of the O2 emission in Mrk 231 coincides with the red wing seen in CO emission, suggesting that it is associated with the outflowing molecular gas, located mainly at about ten kpc away from the central AGN. This first detection of extragalactic molecular oxygen provides an ideal tool to study AGN-driven molecular outflows on dynamic time scales of tens of Myr. O2 may be a significant coolant for molecular gas in such regions affected by AGN-driven outflows. New astrochemical models are needed to explain the implied high molecular oxygen abundance in such regions several kpc away from the center of galaxies.

preprint2020arXiv

Rotation of Two Micron All Sky Survey Clumps in Molecular Clouds

We have analyzed the rotational properties of 12 clumps using $^{13}$CO (1--0) and C$^{18}$O (1--0) maps of the Five College Radio Astronomy Observatory 13.7 m radio telescope. The clumps, located within molecular clouds, have radii ($R$) in the range of 0.06 -- 0.27\,pc. The direction of clump elongation is not correlated with the direction of the velocity gradient. We measured the specific angular momentum (J/M) to be between 0.0022 and 0.025 pc\,km\,s$^{-1}$ based on $^{13}$CO images, and between 0.0025 and 0.021 pc\,km\,s$^{-1}$ based on C$^{18}$O images. The consistency of $J/M$ based on different tracers indicates the $^{13}$CO and C$^{18}$O in dense clumps trace essentially the same material despite significantly different opacities. We also found that $J/M$ increases monotonically as a function of $R$ in power--law form, $J/M~\propto~R^{1.58~\pm~0.11}$. The ratio between rotation energy and gravitational energy, $β$, ranges from 0.0012 to 0.018. The small values of $β$ imply that rotation alone is not sufficient to support the clump against gravitational collapse.

preprint2019arXiv

ALMA observations reveal no preferred outflow--filament and outflow--magnetic field orientations

We present a statistical study on the orientation of outflows with respect to large-scale filaments and the magnetic fields. Although filaments are widely observed toward Galactic star-forming regions, the exact role of filaments in star formation is unclear. Studies toward low-mass star-forming regions revealed both preferred and random orientation of outflows respective to the filament long-axes, while outflows in massive star-forming regions mostly oriented perpendicular to the host filaments, and parallel to the magnetic fields at similar physical scales. Here, we explore outflows in a sample of 11 protoclusters in HII regions, a more evolved stage compared to IRDCs, using ALMA CO (3-2) line observations. We identify a total of 105 outflow lobes in these protoclusters. Among the 11 targets, 7 are embedded within parsec-scale filamentary structures detected in $^{13}$CO line and 870 $μm$ continuum emissions. The angles between outflow axes and corresponding filaments ($γ_\mathrm{Fil}$) do not show any hint of preferred orientations (i.e., orthogonal or parallel as inferred in numerical models) with respect to the position angle of the filaments. Identified outflow lobes are also not correlated with the magnetic fields and Galactic plane position angles. Outflows associated with filaments aligned along the large-scale magnetic fields are also randomly orientated. Our study presents the first statistical results of outflow orientation respective to large-scale filaments and magnetic fields in evolved massive star-forming regions. The random distribution suggests a lack of alignment of outflows with filaments, which may be a result of the evolutionary stage of the clusters.

preprint2018arXiv

An Ammonia Spectral Map of the L1495-B218 Filaments in the Taurus Molecular Cloud: II CCS & HC$_7$N Chemistry and Three Modes of Star Formation in the Filaments

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

preprint2016arXiv

Discovery of an extremely wide-angle bipolar outflow in AFGL 5142

Most bipolar outflows are associated with individual young stellar objects and have small opening angles. Here we report the discovery of an extremely wide-angle ($\sim$180$\arcdeg$) bipolar outflow ("EWBO") in a cluster forming region AFGL 5142 from low-velocity emission of the HCN (3-2) and HCO$^{+}$ (3-2) lines. This bipolar outflow is along a north-west to south-east direction with a line-of-sight flow velocity of about 3 km~s$^{-1}$ and is spatially connected to the high-velocity jet-like outflows. It seems to be a collection of low-velocity material entrained by the high-velocity outflows due to momentum feedback. The total ejected mass and mass loss rate due to both high velocity jet-like outflows and the "EWBO" are $\sim$24.5 M$_{\sun}$ and $\sim1.7\times10^{-3}$ M$_{\sun}$~yr$^{-1}$, respectively. Global collapse of the clump is revealed by the "blue profile" in the HCO$^{+}$ (1-0) line. A hierarchical network of filaments was identified in NH$_{3}$ (1,1) emission. Clear velocity gradients of the order of 10 km~s$^{-1}$~pc$^{-1}$ are found along filaments, indicating gas inflow along the filaments. The sum of the accretion rate along filaments and mass infall rate along the line of sight is $\sim$3.1$\times10^{-3}$ M$_{\sun}$~yr$^{-1}$, which exceeds the total mass loss rate, indicating that the central cluster is probably still gaining mass. The central cluster is highly fragmented and 22 condensations are identified in 1.1 mm continuum emission. The fragmentation process seems to be determined by thermal pressure and turbulence. The magnetic field may not play an important role in fragmentation.

preprint2016arXiv

Evolution of OH and CO-dark Molecular Gas Fraction Across a Molecular Cloud Boundary In Taurus

We present observations of 12CO J=1-0, 13CO J=1-0, HI, and all four ground-state transitions of the hydroxyl (OH) radical toward a sharp boundary region of the Taurus molecular cloud. Based on a PDR model that reproduces CO and [CI] emission from the same region, we modeled the three OH transitions, 1612, 1665, 1667 MHz successfully through escape probability non-LTE radiative transfer model calculations. We could not reproduce the 1720 MHz observations, due to un-modeled pumping mechanisms, of which the most likely candidate is a C-shock. The abundance of OH and CO-dark molecular gas (DMG) are well constrained. The OH abundance [OH]/[H2] decreases from 8*10-7 to 1*10-7 as Av increases from 0.4 to 2.7 mag, following an empirical law [OH]/[H2]= 1.5 * 10^{-7} + 9.0 * 10^{-7} * exp(-Av/0.81), which is higher than PDR model predictions for low extinction regions by a factor of 80. The overabundance of OH at extinctions at or below 1 mag is likely the result of a C-shock. The dark gas fraction (DGF, defined as fraction of molecular gas without detectable CO emission) decreases from 80% to 20%, following a gaussian profile DGF= 0.90 * exp(-( Av -0.79 )/0.71)^2) This trend of the DGF is consistent with our understanding that the DGF drops at low visual extinction due to photodissociation of H2 and drops at high visual extinction due to CO formation. The DGF peaks in the extinction range where H2 has already formed and achieved self-shielding but 12CO has not. Two narrow velocity components with a peak-to-peak spacing of ~ 1 km s-1 were clearly identified. Their relative intensity and variation in space and frequency suggest colliding streams or gas flows at the boundary region.

preprint2016arXiv

Geometry-Independent Determination of Radial Density Distributions in Molecular Cloud Cores and Other Astronomical Objects

We present a geometry-independent method for determining the shapes of radial volume density profiles of astronomical objects whose geometries are unknown, based on a single column density map. Such profiles are often critical to understand the physics and chemistry of molecular cloud cores, in which star formation takes place. The method presented here does not assume any geometry for the object being studied, thus removing a significant source of bias. Instead it exploits contour self-similarity in column density maps which appears to be common in data for astronomical objects. Our method may be applied to many types of astronomical objects and observable quantities so long as they satisfy a limited set of conditions which we describe in detail. We derive the method analytically, test it numerically, and illustrate its utility using 2MASS-derived dust extinction in molecular cloud cores. While not having made an extensive comparison of different density profiles, we find that the overall radial density distribution within molecular cloud cores is adequately described by an attenuated power law.

preprint2016arXiv

L1599B: Cloud Envelope and C+ Emission in a Region of Moderately Enhanced Radiation Field

We study the effects of an asymmetric radiation field on the properties of a molecular cloud envelope. We employ observations of carbon monoxide (12CO and 13CO), atomic carbon, ionized carbon, and atomic hydrogen to analyze the chemical and physical properties of the core and envelope of L1599B, a molecular cloud forming a portion of the ring at approximately 27 pc from the star Lambda Ori. The O III star provides an asymmetric radiation field that produces a moderate enhancement of the external radiation field. Observations of the [CII] fine structure line with the GREAT instrument on SOFIA indicate a significant enhanced emission on the side of the cloud facing the star, while the [Ci], 12CO and 13CO J = 1-0 and 2-1, and 12CO J = 3-2 data from the PMO and APEX telescopes suggest a relatively typical cloud interior. The atomic, ionic, and molecular line centroid velocities track each other very closely, and indicate that the cloud may be undergoing differential radial motion. The HI data from the Arecibo GALFA survey and the SOFIA/GREAT [CII] data do not suggest any systematic motion of the halo gas, relative to the dense central portion of the cloud traced by 12CO and 13CO.

preprint2016arXiv

The Magnetic Field of L1544: I. Near-Infrared Polarimetry and the Non-Uniform Envelope

The magnetic field (B-field) of the starless dark cloud L1544 has been studied using near-infrared (NIR) background starlight polarimetry (BSP) and archival data in order to characterize the properties of the plane-of-sky B-field. NIR linear polarization measurements of over 1,700 stars were obtained in the H-band and 201 of these were also measured in the K-band. The NIR BSP properties are correlated with reddening, as traced using the RJCE (H-M) method, and with thermal dust emission from the L1544 cloud and envelope seen in Herschel maps. The NIR polarization position angles change at the location of the cloud and exhibit their lowest dispersion of position angles there, offering strong evidence that NIR polarization traces the plane-of-sky B-field of L1544. In this paper, the uniformity of the plane-of-sky B-field in the envelope region of L1544 is quantitatively assessed. This allowed evaluating the approach of assuming uniform field geometry when measuring relative mass-to-flux ratios in the cloud envelope and core based on averaging of the envelope radio Zeeman observations, as in Crutcher et al. (2009). In L1544, the NIR BSP shows the envelope B-field to be significantly non-uniform and likely not suitable for averaging Zeeman properties without treating intrinsic variations. Deeper analyses of the NIR BSP and related data sets, including estimates of the B-field strength and testing how it varies with position and gas density, are the subjects of later papers in this series.

preprint2015arXiv

Herschel Galactic plane survey of [NII] fine structure emission

We present the first large scale high angular resolution survey of ionized nitrogen in the Galactic Plane through emission of its two fine structure transitions ([NII]) at 122 $μ$m and 205 $μ$m. The observations were largely obtained with the PACS instrument onboard the Herschel Space Observatory. The lines-of-sight were in the Galactic plane, following those of the Herschel OTKP project GOT C+. Both lines are reliably detected at the 10$^{-8}$ - 10$^{-7}$ $W$m$^{-2}$sr$^{-1}$ level over the range -60$^{o}$ $\leq$ $l$ $\leq$ 60$^{o}$. The $rms$ of the intensity among the 25 PACS spaxels of a given pointing is typically less than one third of the mean intensity, showing that the emission is extended. [NII] is produced in gas in which hydrogen is ionized, and collisional excitation is by electrons. The ratio of the two fine structure transitions provides a direct measurement of the electron density, yielding $n(e)$ largely in the range 10 to 50 cm$^{-3}$ with an average value of 29 cm$^{-3}$ and N$^+$ column densities 10$^{16}$ to 10$^{17}$ cm$^{-2}$. [NII] emission is highly correlated with that of [CII], and we calculate that between 1/3 and 1/2 of the [CII] emission is associated with the ionized gas. The relatively high electron densities indicate that the source of the [NII] emission is not the Warm Ionized Medium (WIM), which has electron densities more than 100 times smaller. Possible origins of the observed [NII] include the ionized surfaces of dense atomic and molecular clouds, the extended low density envelopes of HII regions, and low-filling factor high-density fluctuations of the WIM.

preprint2015arXiv

Herschel observations of interstellar chloronium. II - Detections toward G29.96-0.02, W49N, W51, and W3(OH), and determinations of the ortho-to-para and $^{35}$Cl/$^{37}$Cl isotopic ratios

We report additional detections of the chloronium molecular ion, H$_2$Cl$^+$, toward four bright submillimeter continuum sources: G29.96, W49N, W51, and W3(OH). With the use of the HIFI instrument on the Herschel Space Observatory, we observed the $2_{12}-1_{01}$ transition of ortho-H$_2^{35}$Cl$^+$ at 781.627 GHz in absorption toward all four sources. Much of the detected absorption arises in diffuse foreground clouds that are unassociated with the background continuum sources and in which our best estimates of the $N({\rm H_2Cl^+})/N({\rm H})$ ratio lie in the range $(0.9 - 4.8) \times 10^{-9}$. These chloronium abundances relative to atomic hydrogen can exceed the predictions of current astrochemical models by up to a factor of 5. Toward W49N, we have also detected the $2_{12}-1_{01}$ transition of ortho-H$_2^{37}$Cl$^+$ at 780.053 GHz and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ at 485.418 GHz. These observations imply $\rm H_2^{35}Cl^+/H_2^{37}Cl^+$ column density ratios that are consistent with the solar system $^{35}$Cl/$^{37}$Cl isotopic ratio of 3.1, and chloronium ortho-to-para ratios consistent with 3, the ratio of spin statistical weights.

preprint2015arXiv

Outflows and Bubbles in Taurus: Star-formation Feedback Sufficient to Maintain Turbulence

We have identified outflows and bubbles in the Taurus molecular cloud based on the $\sim 100$ deg$^2$ Five College Radio Astronomy Observatory $^{12}$CO(1-0) and $^{13}$CO(1-0) maps and the Spitzer young stellar object catalogs. In the main 44 deg$^2$ area of Taurus we found 55 outflows, of which 31 were previously unknown. We also found 37 bubbles in the entire 100 deg$^2$ area of Taurus, all of which had not been found before. The total kinetic energy of the identified outflows is estimated to be $\bf \sim 3.9 \times 10^{45}$ erg, which is \textbf{1\%} of the cloud turbulent energy. The total kinetic energy of the detected bubbles is estimated to be $\sim 9.2 \times 10^{46}$ erg, which is 29\% of the turbulent energy of Taurus. The energy injection rate from outflows is $\bf \sim 1.3 \times 10^{33}~\rm erg\ s^{-1}$, \textbf{0.4 - 2 times} the dissipation rate of the cloud turbulence. The energy injection rate from bubbles is $\sim 6.4 \times 10^{33}$ erg s$^{-1}$, \textbf{2 - 10 times} the turbulent dissipation rate of the cloud. The gravitational binding energy of the cloud is $\bf \sim 1.5 \times 10^{48}$ {\bf erg}, \textbf{385} and 16 times the energy of outflows and bubbles, respectively. We conclude that neither outflows nor bubbles can \textbf{provide enough energy to balance the overall gravitational binding energy and the turbulent energy of Taurus. However,} in the current epoch, stellar feedback is sufficient to maintain the observed turbulence in Taurus.

preprint2015arXiv

Planck cold clumps in the $λ$ Orionis complex: I. Discovery of an extremely young Class 0 protostellar object and a proto-brown dwarf candidate in a bright rimmed clump PGCC G192.32-11.88

We are performing a series of observations with ground-based telescopes toward Planck Galactic cold clumps (PGCCs) in the $λ$ Orionis complex in order to systematically investigate the effects of stellar feedback. In the particular case of PGCC G192.32-11.88, we discovered an extremely young Class 0 protostellar object (G192N) and a proto-brown dwarf candidate (G192S). G192N and G192S are located in a gravitationally bound bright-rimmed clump. The velocity and temperature gradients seen in line emission of CO isotopologues indicate that PGCC G192.32-11.88 is externally heated and compressed. G192N probably has the lowest bolometric luminosity ($\sim0.8$ L$_{\sun}$) and accretion rate (6.3$\times10^{-7}$ M$_{\sun}$~yr$^{-1}$) when compared with other young Class 0 sources (e.g. PACS Bright Red sources (PBRs)) in the Orion complex. It has slightly larger internal luminosity ($0.21\pm0.01$ L$_{\sun}$) and outflow velocity ($\sim$14 km~s$^{-1}$) than the predictions of first hydrostatic cores (FHSCs). G192N might be among the youngest Class 0 sources, which are slightly more evolved than a FHSC. Considering its low internal luminosity ($0.08\pm0.01$ L$_{\odot}$) and accretion rate (2.8$\times10^{-8}$ M$_{\sun}$~yr$^{-1}$), G192S is an ideal proto-brown dwarf candidate. The star formation efficiency ($\sim$0.3\%-0.4\%) and core formation efficiency ($\sim$1\%) in PGCC G192.32-11.88 are significantly smaller than in other giant molecular clouds or filaments, indicating that the star formation therein is greatly suppressed due to stellar feedback.

preprint2014arXiv

A Herschel [CII] Galactic plane survey III: [CII] as a tracer of star formation

We study the relationship between the [CII] emission and the star formation rate (SFR) in the Galactic plane and separate the relationship of different ISM phases to the SFR. We compare these relationships to those in external galaxies and local clouds, allowing examinations of these relationships over a wide range of physical scales. We compare the distribution of the [CII] emission, with its different contributing ISM phases, as a function of Galactocentric distance with the SFR derived from radio continuum observations. We also compare the SFR with the surface density distribution of atomic and molecular gas, including the CO-dark H2 component. The [CII] and SFR are well correlated at Galactic scales with a relationship that is in general agreement with that found for external galaxies. By combining [CII] and SFR data points in the Galactic plane with those in external galaxies and nearby star forming regions, we find that a single scaling relationship between the [CII] luminosity and SFR applies over six orders of magnitude. The [CII] emission from different ISM phases are each correlated with the SFR, but only the combined emission shows a slope that is consistent with extragalactic observations. These ISM components have roughly comparable contributions to the Galactic [CII] luminosity: dense PDRs (30%), cold HI (25%), CO-dark H2 (25%), and ionized gas (20%). The SFR-gas surface density relationship shows a steeper slope compared to that observed in galaxies, but one that it is consistent with those seen in nearby clouds. The different slope is a result of the use of a constant CO-to-H2 conversion factor in the extragalactic studies, which in turn is related to the assumption of constant metallicity in galaxies. We find a linear correlation between the SFR surface density and that of the dense molecular gas.

preprint2014arXiv

Airships: A New Horizon for Science

The "Airships: A New Horizon for Science" study at the Keck Institute for Space Studies investigated the potential of a variety of airships currently operable or under development to serve as observatories and science instrumentation platforms for a range of space, atmospheric, and Earth science. The participants represent a diverse cross-section of the aerospace sector, NASA, and academia. Over the last two decades, there has been wide interest in developing a high altitude, stratospheric lighter-than-air (LTA) airship that could maneuver and remain in a desired geographic position (i.e., "station-keeping") for weeks, months or even years. Our study found considerable scientific value in both low altitude (< 40 kft) and high altitude (> 60 kft) airships across a wide spectrum of space, atmospheric, and Earth science programs. Over the course of the study period, we identified stratospheric tethered aerostats as a viable alternative to airships where station-keeping was valued over maneuverability. By opening up the sky and Earth's stratospheric horizon in affordable ways with long-term flexibility, airships allow us to push technology and science forward in a project-rich environment that complements existing space observatories as well as aircraft and high-altitude balloon missions.

preprint2014arXiv

Herschel HIFI observations of O$_2$ toward Orion: special conditions for shock enhanced emission

We report observations of molecular oxygen (O$_2$) rotational transitions at 487 GHz, 774 GHz, and 1121 GHz toward Orion Peak A. The O2 lines at 487 GHz and 774 GHz are detected at velocities of 10-12 km/s with line widths 3 km/s; however, the transition at 1121 GHz is not detected. The observed line characteristics, combined with the results of earlier observations, suggest that the region responsible for the O$_2$ emission is 9" (6e16 cm) in size, and is located close to the H2 Peak 1position (where vibrationally-excited H$_2$ emission peaks), and not at Peak A, 23" away. The peak O2 column density is 1.1e18/cm2. The line velocity is close to that of 621 GHz water maser emission found in this portion of the Orion Molecular Cloud, and having a shock with velocity vector lying nearly in the plane of the sky is consistent with producing maximum maser gain along the line-of-sight. The enhanced O$_2$ abundance compared to that generally found in dense interstellar clouds can be explained by passage of a low-velocity C-shock through a clump with preshock density 2e4/cm3, if a reasonable flux of UV radiation is present. The postshock O$_2$ can explain the emission from the source if its line of sight dimension is ~10 times larger than its size on the plane of the sky. The special geometry and conditions required may explain why O$_2$ emission has not been detected in the cores of other massive star-forming molecular clouds.

preprint2014arXiv

Magnetic Fields in High-Mass Infrared Dark Clouds

High-mass Stars are cosmic engines known to dominate the energetics in the Milky Way and other galaxies. However, their formation is still not well understood. Massive, cold, dense clouds, often appearing as Infrared Dark Clouds (IRDCs), are the nurseries of massive stars. No measurements of magnetic fields in IRDCs in a state prior to the onset of high-mass star formation (HMSF) have previously been available, and prevailing HMSF theories do not consider strong magnetic fields. Here, we report observations of magnetic fields in two of the most massive IRDCs in the Milky Way. We show that IRDCs G11.11-0.12 and G0.253+0.016 are strongly magnetized and that the strong magnetic field is as important as turbulence and gravity for HMSF. The main dense filament in G11.11-0.12 is perpendicular to the magnetic field, while the lower density filament merging onto the main filament is parallel to the magnetic field. The implied magnetic field is strong enough to suppress fragmentation sufficiently to allow HMSF. Other mechanisms reducing fragmentation, such as the entrapment of heating from young stars via high mass surface densities, are not required to facilitate HMSF.

preprint2014arXiv

Water deuterium fractionation in the high-mass star-forming region G34.26+0.15 based on Herschel/HIFI data

Understanding water deuterium fractionation is important for constraining the mechanisms of water formation in interstellar clouds. Observations of HDO and H$_2^{18}$O transitions were carried out towards the high-mass star-forming region G34.26+0.15 with the HIFI instrument onboard the Herschel Space Observatory, as well as with ground-based single-dish telescopes. Ten HDO lines and three H$_2^{18}$O lines covering a broad range of upper energy levels (22-204 K) were detected. We used a non-LTE 1D analysis to determine the HDO/H$_2$O ratio as a function of radius in the envelope. Models with different water abundance distributions were considered in order to reproduce the observed line profiles. The HDO/H$_2$O ratio is found to be lower in the hot core ($\sim$3.5 $\times$ 10$^{-4}$ - 7.5 $\times$ 10$^{-4}$) than in the colder envelope ($\sim$1.0 $\times$ 10$^{-3}$ - 2.2 $\times$ 10$^{-3}$). This is the first time that a radial variation of the HDO/H$_2$O ratio has been found to occur in a high-mass source. The chemical evolution of this source was modeled as a function of its radius and the observations are relatively well reproduced. The comparison between the chemical model and the observations leads to an age of $\sim$10$^5$ years after the infrared dark cloud stage.

preprint2013arXiv

A Herschel [CII] Galactic plane survey I: the global distribution of ISM gas components

[Abridged] The [CII] 158um line is an important tool for understanding the life cycle of interstellar matter. Ionized carbon is present in a variety of phases of the interstellar medium, including the diffuse ionized medium, warm and cold atomic clouds, clouds in transition from atomic to molecular, and dense and warm photon dominated regions (PDRs). The Galactic Observations of Terahertz C+ (GOTC+) project surveys the [CII] line over the entire Galactic disk with velocity-resolved observations using the Herschel/HIFI instrument. We present the first longitude-velocity maps of the [CII] emission for Galactic latitudes b=0deg, +-0.5deg, and +-1.0deg. [CII] emission is mostly associated with spiral arms, mainly emerging from Galactocentric distances between 4 and 10 kpc. We estimate that most of the observed [CII] emission is produced by dense PDRs (47%), with smaller contributions from CO-dark H2 gas (28%), cold atomic gas (21%), and ionized gas (4%). Atomic gas inside the Solar radius is mostly in the form of cold neutral medium (CNM), while the warm neutral medium (WNM) gas dominates the outer galaxy. The average fraction of CNM relative to total atomic gas is 43%. We find that the warm and diffuse CO-dark H2 is distributed over a larger range of Galactocentric distances (4-11 kpc) than the cold and dense H2 gas traced by 12CO and 13CO (4-8kpc). The fraction of CO-dark H2 to total H2 increases with Galactocentric distance, ranging from 20% at 4 kpc to 80% at 10 kpc. On average, CO-dark H2 accounts for 30% of the molecular mass of the Milky Way. When the CO-dark H2 component is included, the radial distribution of the CO-to-H2 conversion factor is steeper than that when only molecular gas traced by CO is considered. Most of the observed [CII] emission emerging from dense PDRs is associated with modest far-ultraviolet fields in the range chi0~1-30.

preprint2013arXiv

Alignment Between Flattened Protostellar Infall Envelopes and Ambient Magnetic Fields

We present 350 micron polarization observations of four low-mass cores containing Class 0 protostars: L483, L1157, L1448-IRS2, and Serp-FIR1. This is the second paper in a larger survey aimed at testing magnetically regulated models for core-collapse. One key prediction of these models is that the mean magnetic field in a core should be aligned with the symmetry axis (minor axis) of the flattened YSO inner envelope (aka pseudodisk). Furthermore, the field should exhibit a pinched or hour-glass shaped morphology as gravity drags the field inward towards the central protostar. We combine our results for the four cores with results for three similar cores that were published in the first paper from our survey. An analysis of the 350 micron polarization data for the seven cores yields evidence of a positive correlation between mean field direction and pseudodisk symmetry axis. Our rough estimate for the probability of obtaining by pure chance a correlation as strong as the one we found is about 5%. In addition, we combine together data for multiple cores to create a source-averaged magnetic field map having improved signal-to-noise ratio, and this map shows good agreement between mean field direction and pseudodisk axis (they are within 15 degrees). We also see hints of a magnetic pinch in the source-averaged map. We conclude that core-scale magnetic fields appear to be strong enough to guide gas infall, as predicted by the magnetically regulated models. Finally, we find evidence of a positive correlation between core magnetic field direction and bipolar outflow axis.

preprint2013arXiv

Deep observations of O2 toward a low-mass protostar with Herschel-HIFI

According to traditional gas-phase chemical models, O2 should be abundant in molecular clouds, but until recently, attempts to detect interstellar O2 line emission with ground- and space-based observatories have failed. Following the multi-line detections of O2 with low abundances in the Orion and rho Oph A molecular clouds with Herschel, it is important to investigate other environments, and we here quantify the O2 abundance near a solar-mass protostar. Observations of O2, at 487 GHz toward a deeply embedded low-mass Class 0 protostar, NGC 1333-IRAS 4A, are presented, using the HIFI instrument on the Herschel Space Observatory. Complementary data of the chemically related NO and CO molecules are obtained as well. The high spectral resolution data are analysed using radiative transfer models to infer column densities and abundances, and are tested directly against full gas-grain chemical models. The deep HIFI spectrum fails to show O2 at the velocity of the dense protostellar envelope, implying one of the lowest abundance upper limits of O2/H2 at <6x10^-9 (3 sigma). However, a tentative (4.5 sigma) detection of O2 is seen at the velocity of the surrounding NGC 1333 molecular cloud, shifted by 1 km/s relative to the protostar. For the protostellar envelope, pure gas-phase models and gas-grain chemical models require a long pre-collapse phase (~0.7-1x10^6 years), during which atomic and molecular oxygen are frozen out onto dust grains and fully converted to H2O, to avoid overproduction of O2 in the dense envelope. The same model also reproduces the limits on the chemically related NO molecule. The tentative detection of O2 in the surrounding cloud is consistent with a low-density PDR model with small changes in reaction rates. The low O2 abundance in the collapsing envelope around a low-mass protostar suggests that the gas and ice entering protoplanetary disks is very poor in O2.

preprint2013arXiv

Diffuse Molecular Cloud Densities from UV Measurements of CO Absorption

We use UV measurements of interstellar CO towards nearby stars to calculate the density in the diffuse molecular clouds containing the molecules responsible for the observed absorption. Chemical models and recent calculations of the excitation rate coefficients indicate that the regions in which CO is found have hydrogen predominantly in molecular form. We carry out statistical equilibrium calculations using CO-H2 collision rates to solve for the H2 density in the observed sources without including effects of radiative trapping. We have assumed kinetic temperatures of 50 K and 100 K, finding this choice to make relatively little difference to the lowest transition. For the sources having T_ex(1-0) only, for which we could determine upper and lower density limits, we find <n(H2)> = 49 cm-3. While we can find a consistent density range for a good fraction of the sources having either two or three values of the excitation temperature, there is a suggestion that the higher-J transitions are sampling clouds or regions within diffuse molecular cloud material that have higher densities than the material sampled by the J = 1-0 transition. The assumed kinetic temperature and derived H2 density are anticorrelated when the J = 2-1 transition data, the J = 3-2 transition data, or both are included. For sources with either two or three values of the excitation temperature, we find average values of the midpoint of the density range that is consistent with all of the observations equal to 68 cm-3 for T_k = 100 K and 92 cm-3 for T_k = 50 K. The data for this set of sources imply that diffuse molecular clouds are characterized by an average thermal pressure between 4600 and 6800 Kcm-3.

preprint2013arXiv

Low Virial Parameters in Molecular Clouds: Implications for High Mass Star Formation and Magnetic Fields

Whether or not molecular clouds and embedded cloud fragments are stable against collapse is of utmost importance for the study of the star formation process. Only "supercritical" cloud fragments are able to collapse and form stars. The virial parameter, alpha=M_vir/M, which compares the virial to the actual mass, provides one way to gauge stability against collapse. Supercritical cloud fragments are characterized by alpha<2, as indicated by a comprehensive stability analysis considering perturbations in pressure and density gradients. Past research has suggested that virial parameters alpha>2 prevail in clouds. This would suggest that collapse towards star formation is a gradual and relatively slow process, and that magnetic fields are not needed to explain the observed cloud structure. Here, we review a range of very recent observational studies that derive virial parameters <<2 and compile a catalogue of 1325 virial parameter estimates. Low values of alpha are in particular observed for regions of high mass star formation (HMSF). These observations may argue for a more rapid and violent evolution during collapse. This would enable "competitive accretion" in HMSF, constrain some models of "monolithic collapse", and might explain the absence of high--mass starless cores. Alternatively, the data could point at the presence of significant magnetic fields ~1 mG at high gas densities. We examine to what extent the derived observational properties might be biased by observational or theoretical uncertainties. For a wide range of reasonable parameters, our conclusions appear to be robust with respect to such biases.

preprint2012arXiv

Collisional Excitation of the [CII] Fine Structure Transition in Interstellar Clouds

We analyze the collisional excitation of the 158 micron (1900.5 GHz) fine structure transition of ionized carbon (C+) in terms of line intensities produced by simple cloud models. The single C+ fine structure transition is a very important coolant of the atomic interstellar medium and of photon dominated regions in which carbon is partially or completely in ionized form. The [CII] line is widely used as a tracer of star formation in the Milky Way and other galaxies. Excitation of the [CII] fine structure transition can be via collisions with hydrogen molecules, atoms, and electrons. Velocity-resolved observations of [CII] have become possible with the HIFI instrument on Herschel and the GREAT instrument on SOFIA. Analysis of these observations is complicated by the fact that it is difficult to determine the optical depth of the [CII] line due to the relative weakness and blending of the components of the analogous transition of 13C$+. We discuss the excitation and radiative transition of the [CII] line, deriving analytic results for several limiting cases and carry out numerical solutions using a large velocity gradient model for a more inclusive analysis. We show that for antenna temperatures up to 1/3 of the brightness temperature of the gas kinetic temperature, the antenna temperature is linearly proportional to the column density of C+ irrespective of the optical depth of the transition, which can be referred to as the effectively optically thin (EOT) approximation. We review the critical densities for excitation of the [CII] line by various collision partners. We briefly analyze C+ absorption and conclude with a discussion of C+ cooling and how the considerations for line intensities affect the behavior of this important coolant of the ISM.

preprint2012arXiv

Herschel observations of interstellar chloronium

Using the Herschel Space Observatory's Heterodyne Instrument for the Far-Infrared (HIFI), we have observed para-chloronium (H2Cl+) toward six sources in the Galaxy. We detected interstellar chloronium absorption in foreground molecular clouds along the sight-lines to the bright submillimeter continuum sources Sgr A (+50 km/s cloud) and W31C. Both the para-H2-35Cl+ and para-H2-37Cl+ isotopologues were detected, through observations of their 1(11)-0(00) transitions at rest frequencies of 485.42 and 484.23 GHz, respectively. For an assumed ortho-to-para ratio of 3, the observed optical depths imply that chloronium accounts for ~ 4 - 12% of chlorine nuclei in the gas phase. We detected interstellar chloronium emission from two sources in the Orion Molecular Cloud 1: the Orion Bar photodissociation region and the Orion South condensation. For an assumed ortho-to-para ratio of 3 for chloronium, the observed emission line fluxes imply total beam-averaged column densities of ~ 2.0E+13 cm-2 and ~ 1.2E+13 cm-2, respectively, for chloronium in these two sources. We obtained upper limits on the para-H2-35Cl+ line strengths toward H2 Peak 1 in the Orion Molecular cloud and toward the massive young star AFGL 2591. The chloronium abundances inferred in this study are typically at least a factor ~10 larger than the predictions of steady-state theoretical models for the chemistry of interstellar molecules containing chlorine. Several explanations for this discrepancy were investigated, but none has proven satisfactory, and thus the large observed abundances of chloronium remain puzzling.

preprint2012arXiv

Herschel Search for O2 Toward the Orion Bar

We report the results of a search for molecular oxygen (O2) toward the Orion Bar, a prominent photodissociation region at the southern edge of the HII region created by the luminous Trapezium stars. We observed the spectral region around the frequency of the O2 N_J = 3_3 - 1_2 transition at 487 GHz and the 5_4 - 3_4 transition at 774 GHz using the Heterodyne Instrument for the Far Infrared on the Herschel Space Observatory. Neither line was detected, but the 3sigma upper limits established here translate to a total line-of-sight O2 column density < 1.5 10^16 cm^-2 for an emitting region whose temperature is between 30K and 250 K, or < 1 10^16 cm^-2 if the O2 emitting region is primarily at a temperature of ~< 100 K. Because the Orion Bar is oriented nearly edge-on relative to our line of sight, the observed column density is enhanced by a factor estimated to be between 4 and 20 relative to the face-on value. Our upper limits imply that the face-on O2 column density is less than 4 10^15 cm^-2, a value that is below, and possibly well below, model predictions for gas with a density of 10^4 - 10^5 cm^-3 exposed to a far ultraviolet flux 10^4 times the local value, conditions inferred from previous observations of the Orion Bar. The discrepancy might be resolved if: (1) the adsorption energy of O atoms to ice is greater than 800 K; (2) the total face-on Av of the Bar is less than required for O2 to reach peak abundance; (3) the O2 emission arises within dense clumps with a small beam filling factor; or, (4) the face-on depth into the Bar where O2 reaches its peak abundance, which is density dependent, corresponds to a sky position different from that sampled by our Herschel beams.

preprint2012arXiv

Is the Taurus B213 Region a True Filament?: Observations of Multiple Cyanoacetylene Transitions

We have obtained spectra of the J=2-1 and J=10-9 transitions of cyanoacetylene (\hc3n) toward a collection of positions in the most prominent filament, B213, in the Taurus molecular cloud. The analysis of the excitation conditions of these transitions reveals an average gas H$_2$ volume density of $(1.8\pm 0.7) \times10^{4} $ \cc. Based on column density derived from 2MASS and this volume density, the line of sight dimension of the high density portion of B213 is found to be $\simeq$ 0.12 pc, which is comparable to the smaller projected dimension and much smaller than the elongated dimension of B213 ($\sim$2.4 pc). B213 is thus likely a true cylinder--like filament rather than a sheet seen edge-on. The line width and velocity gradient seen in \hc3n are also consistent with Taurus B213 being a self-gravitating filament in the early stage of either fragmentation and/or collape.

preprint2011arXiv

Herschel Measurements of Molecular Oxygen in Orion

We report observations of three rotational transitions of molecular oxygen (O2) in emission from the H2 Peak 1 position of vibrationally excited molecular hydrogen in Orion. We observed the 487 GHz, 774 GHz, and 1121 GHz lines using HIFI on the Herschel Space Observatory, having velocities of 11 km s-1 to 12 km s-1 and widths of 3 km s-1. The beam-averaged column density is N(O2) = 6.5\times1016 cm-2, and assuming that the source has an equal beam filling factor for all transitions (beam widths 44, 28, and 19"), the relative line intensities imply a kinetic temperature between 65 K and 120 K. The fractional abundance of O2 relative to H2 is 0.3 - 7.3\times10-6. The unusual velocity suggests an association with a ~ 5" diameter source, denoted Peak A, the Western Clump, or MF4. The mass of this source is ~ 10 M\odot and the dust temperature is \geq 150 K. Our preferred explanation of the enhanced O2 abundance is that dust grains in this region are sufficiently warm (T \geq 100 K) to desorb water ice and thus keep a significant fraction of elemental oxygen in the gas phase, with a significant fraction as O2. For this small source, the line ratios require a temperature \geq 180 K. The inferred O2 column density \simeq 5\times1018 cm-2 can be produced in Peak A, having N(H2) \simeq 4\times1024 cm-2. An alternative mechanism is a low-velocity (10 to 15 km s-1) C-shock, which can produce N(O2) up to 1017 cm-2.

preprint2011arXiv

The Magnetic Field in Taurus Probed by Infrared Polarization

We present maps of the plane-of-sky magnetic field within two regions of the Taurus molecular cloud: one in the dense core L1495/B213 filament, the other in a diffuse region to the west. The field is measured from the polarization of background starlight seen through the cloud. In total, we measured 287 high-quality near-infrared polarization vectors in these regions. In L1495/B213, the percent polarization increases with column density up to Av ~ 9 mag, the limits of our data. The Radiative Torques model for grain alignment can explain this behavior, but models that invoke turbulence are inconsistent with the data. We also combine our data with published optical and near-infrared polarization measurements in Taurus. Using this large sample, we estimate the strength of the plane-of-sky component of the magnetic field in nine subregions. This estimation is done with two different techniques that use the observed dispersion in polarization angles. Our values range from 5-82 microgauss and tend to be higher in denser regions. In all subregions, the critical index of the mass-to-magnetic flux ratio is sub-unity, implying that Taurus is magnetically supported on large scales (~2 pc). Within the region observed, the B213 filament makes a sharp turn to the north and the direction of the magnetic field also takes a sharp turn, switching from being perpendicular to the filament to becoming parallel. This behavior can be understood if we are observing the rim of a bubble. We argue that it has resulted from a supernova remnant associated with a recently discovered nearby gamma-ray pulsar.

preprint2010arXiv

A Sample of [CII] Clouds Tracing Dense Clouds in Weak FUV Fields observed by Herschel

The [CII] fine--structure line at 158um is an excellent tracer of the warm diffuse gas in the ISM and the interfaces between molecular clouds and their surrounding atomic and ionized envelopes. Here we present the initial results from Galactic Observations of Terahertz C+ (GOTC+), a Herschel Key Project devoted to study the [CII] fine structure emission in the galactic plane using the HIFI instrument. We use the [CII] emission together with observations of CO as a probe to understand the effects of newly--formed stars on their interstellar environment and characterize the physical and chemical state of the star-forming gas. We collected data along 16 lines--of--sight passing near star forming regions in the inner Galaxy near longitudes 330 degrees and 20 degrees. We identify fifty-eight [CII] components that are associated with high--column density molecular clouds as traced by 13CO emission. We combine [CII], 12CO, and 13CO observations to derive the physical conditions of the [CII]--emitting regions in our sample of high--column density clouds based on comparison with results from a grid of Photon Dominated Region (PDR) models. From this unbiased sample, our results suggest that most of [CII] emission originates from clouds with H2 volume densities between 10e3.5 and 10e5.5 cm^-3 and weak FUV strength (CHI_0=1-10). We find two regions where our analysis suggests high densities >10e5 cm^-3 and strong FUV fields (CHI=10e4-10e6), likely associated with massive star formation. We suggest that [CII] emission in conjunction with CO isotopes is a good tool to differentiate between regions of massive star formation (high densities/strong FUV fields) and regions that are distant from massive stars (lower densities/weaker FUV fields) along the line--of--sight

preprint2010arXiv

Interstellar OH+, H2O+ and H3O+ along the sight-line to G10.6-0.4

We report the detection of absorption lines by the reactive ions OH+, H2O+ and H3O+ along the line of sight to the submillimeter continuum source G10.6$-$0.4 (W31C). We used the Herschel HIFI instrument in dual beam switch mode to observe the ground state rotational transitions of OH+ at 971 GHz, H2O+ at 1115 and 607 GHz, and H3O+ at 984 GHz. The resultant spectra show deep absorption over a broad velocity range that originates in the interstellar matter along the line of sight to G10.6$-$0.4 as well as in the molecular gas directly associated with that source. The OH+ spectrum reaches saturation over most velocities corresponding to the foreground gas, while the opacity of the H2O+ lines remains lower than 1 in the same velocity range, and the H3O+ line shows only weak absorption. For LSR velocities between 7 and 50 kms$^{-1}$ we estimate total column densities of $N$(OH+) $> 2.5 \times 10^{14}$ cm$^{-2}$, $N$(H2O+) $\sim 6 \times 10^{13}$ cm$^{-2}$ and $N$(H3O+) $\sim 4.0 \times 10^{13}$ cm$^{-2}$. These detections confirm the role of O$^+$ and OH$^+$ in initiating the oxygen chemistry in diffuse molecular gas and strengthen our understanding of the gas phase production of water. The high ratio of the OH+ by the H2O+ column density implies that these species predominantly trace low-density gas with a small fraction of hydrogen in molecular form.

preprint2010arXiv

The relation between gas and dust in the Taurus Molecular Cloud

(abridged) We report a study of the relation between dust and gas over a 100deg^2 area in the Taurus molecular cloud. We compare the H2 column density derived from dust extinction with the CO column density derived from the 12CO and 13CO J= 1-0 lines. We derive the visual extinction from reddening determined from 2MASS data. The comparison is done at an angular size of 200", corresponding to 0.14pc at a distance of 140pc. We find that the relation between visual extinction Av and N(CO) is linear between Av~3 and 10 mag in the region associated with the B213--L1495 filament. In other regions the linear relation is flattened for Av > 4 mag. We find that the presence of temperature gradients in the molecular gas affects the determination of N(CO) by ~30--70% with the largest difference occurring at large column densities. Adding a correction for this effect and accounting for the observed relation between the column density of CO and CO2 ices and Av, we find a linear relationship between the column of carbon monoxide and dust for observed visual extinctions up to the maximum value in our data 23mag. We have used these data to study a sample of dense cores in Taurus. Fitting an analytical column density profile to these cores we derive an average volume density of about 1.4e4 cm^-3 and a CO depletion age of about 4.2e5 years. We estimate the H2 mass of Taurus to be about 1.5e4 M_sun, independently derived from the Av and N(CO) maps. We derive a CO integrated intensity to H2 conversion factor of about 2.1e20 cm^-2 (K km/s)^-1, which applies even in the region where the [CO]/[H_2] ratio is reduced by up to two orders of magnitude. The distribution of column densities in our Taurus maps resembles a log--normal function but shows tails at large and low column densities.

preprint2009arXiv

A Flexible Quasioptical Input System for a Submillimeter Muliobject Spectrometer

We present a conceptual design for the input optical system for a multi-object spectrometer operating at submillimeter wavelengths. The Mirror MOS is based on a sequence of mirrors that enables low-loss propagation of beams from selected positions distributed throughout the focal plane to the spectroscopic receiver inputs. This approach should be useful for observations of sources which have a relatively low density on the sky, for which it is inefficient to use a traditional array receiver with uniformly spaced, relatively closely packed beams. Our concept is based on assigning a patrol region to each of the receivers, which have inputs distributed over the focal plane of the telescope. The input to each receiver can be positioned at any point within this patrol region. This approach, with only 4 reflections, offers very low loss. The Gaussian beam optical system can be designed to produce frequency-independent illumination of the telescope, which is an important advantage for broadband systems such those required for determination of redshifts of submillimeter galaxies.