We study the formation, evolution and collapse of dense cores by tracking density structures in a magnetohydrodynamic (MHD) simulation. We identify cores using the dendrogram algorithm and utilize machine learning techniques, including principal component analysis (PCA) and the k-means clustering algorithm to analyze the full density and velocity dispersion profiles of these cores. We find that there exists an evolutionary sequence consisting of three distinct phases: i) the formation of turbulent density structures (Phase I), ii) the dissipation of turbulence and the formation of coherent cores (Phase II), and iii) the transition to protostellar cores through gravitational collapse (Phase III). In dynamically evolving molecular clouds, the existence of these three phases corresponds to the coexistence of three populations of cores with distinct physical properties. The prestellar and protostellar cores frequently analyzed in previous studies of observations and simulations belong to the last phase in this evolutionary picture. We derive typical lifetimes of 1.4$\pm$1.0$\times$10$^5$ yr, 3.3$\pm$1.4$\times$10$^5$ yr and 3.3$\pm$1.4$\times$10$^5$ yr, respectively for Phase I, II and III. We find that cores can form from both converging flows and filament fragmentation and that cores may form both inside and outside the filaments. We then compare our results to previous observations of coherent cores and provide suggestions for future observations to study cores belonging to the three phases.

We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particularly efficient. Methanol and sulphur monoxide are most abundant in a narrow zone that follows the eastern side of the core. This side is sheltered from the stronger external radiation field coming from the west. We show that photodissociation on the illuminated side can give rise to an asymmetric methanol distribution, but that the stark contrast observed in H-MM1 is hard to explain without assuming enhanced desorption on the shaded side. The region of the brightest emission has a wavy structure that rolls up at one end. This is the signature of Kelvin-Helmholtz instability occurring in sheared flows. We suggest that in this zone, methanol and sulphur are released as a result of grain-grain collisions induced by shear vorticity.

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

For the past 150 years, the prevailing view of the local Interstellar Medium (ISM) was based on a peculiarity known as the Gould's Belt, an expanding ring of young stars, gas, and dust, tilted about 20$^\circ$ to the Galactic plane. Still, the physical relation between local gas clouds has remained practically unknown because the distance accuracy to clouds is of the same order or larger than their sizes. With the advent of large photometric surveys and the Gaia satellite astrometric survey this situation has changed. Here we report the 3-D structure of all local cloud complexes. We find a narrow and coherent 2.7 kpc arrangement of dense gas in the Solar neighborhood that contains many of the clouds thought to be associated with the Gould Belt. This finding is inconsistent with the notion that these clouds are part of a ring, disputing the Gould Belt model. The new structure comprises the majority of nearby star-forming regions, has an aspect ratio of about 1:20, and contains about 3 million solar masses of gas. Remarkably, the new structure appears to be undulating and its 3-D distribution is well described by a damped sinusoidal wave on the plane of the Milky Way, with an average period of about 2 kpc and a maximum amplitude of about 160 pc. Our results represent a first step in the revision of the local gas distribution and Galactic structure and offer a new, broader context to studies on the transformation of molecular gas into stars.

If dark matter is efficiently captured by a planet, energy released in its annihilation can exceed that planet's total heat output. Building on prior work, we treat Earth's composition and dark matter capture in detail and present improved limits on dark matter-nucleon scattering cross sections for dark matter masses ranging from 0.1 to $10^{10}$ GeV. We also extend Earth limits by applying the same treatment to Mars. The scope of dark matter models considered is expanded to include spin-dependent nuclear interactions including isospin-independent, proton only, and neutron only interactions. We find that Earth and Mars heating bounds are alleviated for dark matter s-wave self-annihilation cross sections $\lesssim 10^{-44}~{\rm cm^2}$.

Accurate distances to local molecular clouds are critical for understanding the star and planet formation process, yet distance measurements are often obtained inhomogeneously on a cloud-by-cloud basis. We have recently developed a method which combines stellar photometric data with Gaia DR2 parallax measurements in a Bayesian framework to infer the distances of nearby dust clouds to a typical accuracy of $\sim5\%$. After refining the technique to target lower latitudes and incorporating deep optical data from DECam in the southern Galactic plane, we have derived a catalog of distances to molecular clouds in Reipurth (2008, Star Formation Handbook, vols I and II) which contains a large fraction of the molecular material in the solar neighborhood. Comparison with distances derived from maser parallax measurements towards the same clouds shows our method produces consistent distances with $\lesssim10\%$ scatter for clouds across our entire distance spectrum (150 pc $-$ 2.5 kpc). We hope this catalog of homogeneous distances will serve as a baseline for future work.

We present and release the full dataset for the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. This survey used the Submillimeter Array (SMA) to image the 74 known protostars within the Perseus molecular cloud. The SMA was used in two array configurations to capture outflows for scales $>$30$^{\prime\prime}$ ($>$9000 au) and to probe scales down to $\sim$1$^{\prime\prime}$ ($\sim$300 au). The protostars were observed with the 1.3 mm and 850 $\mu$m receivers simultaneously to detect continuum at both wavelengths and molecular line emission from CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1), N$_2$D$^+$(3-2), CO(3-2), HCO$^+$(4-3), and H$^{13}$CO$^+$(4-3). Some of the observations also used the SMA's recently upgraded correlator, SWARM, whose broader bandwidth allowed for several more spectral lines to be observed (e.g., SO, H$_2$CO, DCO$^+$, DCN, CS, CN). Of the main continuum and spectral tracers observed, 84% of the images and cubes had emission detected. The median C$^{18}$O(2-1) linewidth is $\sim$1.0 km s$^{-1}$, which is slightly higher than those measured with single-dish telescopes at scales of 3000-20000 au. Of the 74 targets, six are suggested to be first hydrostatic core candidates, and we suggest that L1451-mm is the best candidate. We question a previous continuum detection toward L1448 IRS2E. In the SVS13 system, SVS13A certainly appears to be the most evolved source, while SVS13C appears to be hotter and more evolved than SVS13B. The MASSES survey is the largest publicly available interferometric continuum and spectral line protostellar survey to date, and is largely unbiased as it only targets protostars in Perseus. All visibility ($uv$) data and imaged data are publicly available at https://dataverse.harvard.edu/dataverse/full_MASSES/.

Using a population of large-scale filaments extracted from an AREPO simulation of a Milky Way-like galaxy, we seek to understand the extent to which observed large-scale filament properties (with lengths $\gtrsim 100$ pc) can be explained by galactic dynamics alone. From an observer's perspective in the disk of the galaxy, we identify filaments forming purely due to galactic dynamics, without the effects of feedback or local self-gravity. We find that large-scale Galactic filaments are intrinsically rare, and we estimate that at maximum approximately one filament per $\rm kpc^{2}$ should be identified in projection, when viewed from the direction of our Sun in the Milky Way. In this idealized scenario, we find filaments in both the arm and interarm regions, and hypothesize that the former may be due to gas compression in the spiral-potential wells, with the latter due to differential rotation. Using the same analysis pipeline applied previously to observations, we analyze the physical properties of large-scale Galactic filaments, and quantify their sensitivity to projection effects and galactic environment (i.e. whether they lie in the arm or interarm regions). We find that observed "Giant Molecular Filaments" are consistent with being non-self-gravitating structures dominated by galactic dynamics. Straighter, narrower, and denser "Bone-like" filaments, like the paradigmatic Nessie filament, have similar column densities, velocity gradients, and Galactic plane heights ($z\approx$ 0 pc) to those in our simple model, but additional physical effects (such as feedback and self-gravity) must be invoked to explain their lengths and widths.

We present an analysis of the internal velocity structures of the newly identified sub-0.1 pc coherent structures, droplets, in L1688 and B18. By fitting 2D linear velocity fields to the observed maps of velocity centroids, we determine the magnitudes of linear velocity gradients and examine the potential rotational motions that could lead to the observed velocity gradients. The results show that the droplets follow the same power-law relation between the velocity gradient and size found for larger-scale dense cores. Assuming that rotational motion giving rise to the observed velocity gradient in each core is a solid-body rotation of a rotating body with a uniform density, we derive the "net rotational motions" of the droplets. We find a ratio between rotational and gravitational energies, $\beta$, of $\sim 0.046$ for the droplets, and when including both droplets and larger-scale dense cores, we find $\beta \sim 0.039$. We then examine the alignment between the velocity gradient and the major axis of each droplet, using methods adapted from the histogram of relative orientations (HRO) introduced by Soler et al. (2013). We find no definitive correlation between the directions of velocity gradients and the elongations of the cores. Lastly, we discuss physical processes other than rotation that may give rise to the observed velocity field.

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

We present a uniform catalog of accurate distances to local molecular clouds informed by the Gaia DR2 data release. Our methodology builds on that of Schlafly et al. (2014). First, we infer the distance and extinction to stars along sightlines towards the clouds using optical and near-infrared photometry. When available, we incorporate knowledge of the stellar distances obtained from Gaia DR2 parallax measurements. We model these per-star distance-extinction estimates as being caused by a dust screen with a 2-D morphology derived from Planck at an unknown distance, which we then fit for using a nested sampling algorithm. We provide updated distances to the Schlafly et al. (2014) sightlines towards the Dame et al. (2001) and Magnani et al. (1985) clouds, finding good agreement with the earlier work. For a subset of 27 clouds, we construct interactive pixelated distance maps to further study detailed cloud structure, and find several clouds which display clear distance gradients and/or are comprised of multiple components. We use these maps to determine robust average distances to these clouds. The characteristic combined uncertainty on our distances is approximately 5-6%, though this can be higher for clouds at farther distances, due to the limitations of our single-cloud model.

We compare the magnetic field orientation for the young giant molecular cloud Vela C inferred from 500-$\mu$m polarization maps made with the BLASTPol balloon-borne polarimeter to the orientation of structures in the integrated line emission maps from Mopra observations. Averaging over the entire cloud we find that elongated structures in integrated line-intensity, or zeroth-moment maps, for low density tracers such as $^{12}$CO and $^{13}$CO $J$ $\rightarrow$ 1 - 0 are statistically more likely to align parallel to the magnetic field, while intermediate or high density tracers show (on average) a tendency for alignment perpendicular to the magnetic field. This observation agrees with previous studies of the change in relative orientation with column density in Vela C, and supports a model where the magnetic field is strong enough to have influenced the formation of dense gas structures within Vela C. The transition from parallel to no preferred/perpendicular orientation appears to happen between the densities traced by $^{13}$CO and by C$^{18}$O $J$ $\rightarrow$ 1 - 0. Using RADEX radiative transfer models to estimate the characteristic number density traced by each molecular line we find that the transition occurs at a molecular hydrogen number density of approximately $10^3$ cm$^{-3}$. We also see that the Centre-Ridge (the highest column density and most active star-forming region within Vela C) appears to have a transition at a lower number density, suggesting that this may depend on the evolutionary state of the cloud.

Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing.

In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors.

We introduce the histogram of oriented gradients (HOG), a tool developed for machine vision that we propose as a new metric for the systematic characterization of observations of atomic and molecular gas and the study of molecular cloud formation models. In essence, the HOG technique takes as input extended spectral-line observations from two tracers and provides an estimate of their spatial correlation across velocity channels. We characterize HOG using synthetic observations of HI and $^{13}$CO(J=1-0) emission from numerical simulations of MHD turbulence leading to the formation of molecular gas after the collision of two atomic clouds. We find a significant spatial correlation between the two tracers in velocity channels where $v_{HI}\approx v_{^{13}CO}$, independent of the orientation of the collision with respect to the line of sight. We use HOG to investigate the spatial correlation of the HI, from the THOR survey, and the $^{13}$CO(J=1-0) emission, from the GRS, toward the portion of the Galactic plane 33.75$\lt l\lt$35.25$^{o}$ and $|b|\lt$1.25$^{o}$. We find a significant spatial correlation between the tracers in extended portions of the studied region. Although some of the regions with high spatial correlation are associated with HI self-absorption features, suggesting that it is produced by the cold atomic gas, the correlation is not exclusive to this kind of region. The HOG results also indicate significant differences between individual regions: some show spatial correlation in channels around $v_{HI}\approx v_{^{13}CO}$ while others present this correlation in velocity channels separated by a few km/s. We associate these velocity offsets to the effect of feedback and to the presence of physical conditions that are not included in the atomic-cloud-collision simulations, such as more general magnetic field configurations, shear, and global gas infall.

We present combined observations of the NH3 (J,K) = (1,1) and (2,2) inversion transitions towards OMC1 in Orion A obtained by the Karl G. Jansky Very Large Array (VLA) and the 100 m Robert C. Byrd Green Bank Telescope (GBT). With an angular resolution of 6" (0.01 pc), these observations reveal with unprecedented detail the complex filamentary structure extending north of the active Orion BN/KL region in a field covering 6' x 7'. We find a 0.012 pc wide filament within OMC1, with an aspect ratio of ~37:1, that was missed in previous studies. Its orientation is directly compared to the relative orientation of the magnetic field from the James Clerk Maxwell Telescope BISTRO survey in Orion A. We find a small deviation of ~11 deg between the mean orientation of the filament and the magnetic field, suggesting that they are almost parallel to one another. The filament's column density is estimated to be 2-3 orders of magnitude larger than the filaments studied with Herschel and is possibly self-gravitating given the low values of turbulence found. We further produce maps of the gas kinematics by forward modeling the hyperfine structure of the NH3 (J,K) = (1,1) and (2,2) lines. The resulting distribution of velocity dispersions peaks at ~0.5 km/s, close to the subsonic regime of the gas. This value is about 0.2 km/s smaller than previously measured in single-dish observations of the same region, suggesting that higher angular and spectral resolution observations will identify even lower velocity dispersions that might reach the subsonic turbulence regime in dense gas filaments.

We present a new technique to determine distances to major star-forming regions across the Perseus Molecular Cloud, using a combination of stellar photometry, astrometric data, and $\rm ^{12} CO$ spectral-line maps. Incorporating the Gaia DR2 parallax measurements when available, we start by inferring the distance and reddening to stars from their Pan-STARRS1 and 2MASS photometry, based on a technique presented in Green et al. 2014; Green et al. 2015 and implemented in their 3D "Bayestar" dust map of three-quarters of the sky. We then refine the Green et al. technique by using the velocity slices of a CO spectral cube as dust templates and modeling the cumulative distribution of dust along the line of sight towards these stars as a linear combination of the emission in the slices. Using a nested sampling algorithm, we fit these per-star distance-reddening measurements to find the distances to the CO velocity slices towards each star-forming region. This results in distance estimates explicitly tied to the velocity structure of the molecular gas. We determine distances to the B5, IC348, B1, NGC1333, L1448, and L1451 star-forming regions and find that individual clouds are located between $\approx 275-300$ pc, with typical combined uncertainties of $\approx 5\%$. We find that the velocity gradient across Perseus corresponds to a distance gradient of about 25 pc, with the eastern portion of the cloud farther away than the western portion. We determine an average distance to the complex of $294\pm 17$ pc, about 60 pc higher than the distance derived to the western portion of the cloud using parallax measurements of water masers associated with young stellar objects. The method we present is not limited to the Perseus Complex, but may be applied anywhere on the sky with adequate CO data in the pursuit of more accurate 3D maps of molecular clouds in the solar neighborhood and beyond.

We use the $\mathit{Gaia}$ DR2 distances of about 700 mid-infrared selected young stellar objects in the benchmark giant molecular cloud Orion A to infer its 3D shape and orientation. We find that Orion A is not the fairly straight filamentary cloud that we see in (2D) projection, but instead a cometary-like cloud oriented toward the Galactic plane, with two distinct components: a denser and enhanced star-forming (bent) Head, and a lower density and star-formation quieter $\sim$75 pc long Tail. The true extent of Orion A is not the projected $\sim$40 pc but $\sim$90 pc, making it by far the largest molecular cloud in the local neighborhood. Its aspect ratio ($\sim$30:1) and high column-density fraction ($\sim45\%$) make it similar to large-scale Milky Way filaments ("bones"), despite its distance to the galactic mid-plane being an order of magnitude larger than typically found for these structures.

The past year has witnessed discovery of the first identified counterparts to a gravitational wave transient (GW 170817A) and a very high-energy neutrino (IceCube-170922A). These source identifications, and ensuing detailed studies, have realized longstanding dreams of astronomers and physicists to routinely carry out observations of cosmic sources by other than electromagnetic means, and inaugurated the era of "multi-messenger" astronomy. While this new era promises extraordinary physical insights into the universe, it brings with it new challenges, including: highly heterogeneous, high-volume, high-velocity datasets; globe-spanning cross-disciplinary teams of researchers, regularly brought together into transient collaborations; an extraordinary breadth and depth of domain-specific knowledge and computing resources required to anticipate, model, and interpret observations; and the routine need for adaptive, distributed, rapid-response observing campaigns to fully exploit the scientific potential of each source. We argue, therefore, that the time is ripe for the community to conceive and propose an Institute for Multi-Messenger Astrophysics that would coordinate its resources in a sustained and strategic fashion to efficiently address these challenges, while simultaneously serving as a center for education and key supporting activities. In this fashion, we can prepare now to realize the bright future that we see, beyond, through these newly opened windows onto the universe.

We present the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey, which uses the Submillimeter Array (SMA) interferometer to map the continuum and molecular lines for all 74 known Class 0/I protostellar systems in the Perseus molecular cloud. The primary goal of the survey is to observe an unbiased sample of young protostars in a single molecular cloud so that we can characterize the evolution of protostars. This paper releases the MASSES 1.3 mm data from the subcompact configuration ($\sim$4$^{\prime\prime}$ or $\sim$1000 au resolution), which is the SMA's most compact array configuration. We release both $uv$ visibility data and imaged data for the spectral lines CO(2-1),$^{13}$CO(2-1), C$^{18}$O(2-1), and N$_2$D$^+$(3-2), as well as for the 1.3 mm continuum. We identify the tracers that are detected toward each source. We also show example images of continuum and CO(2-1) outflows, analyze C$^{18}$O(2-1) spectra, and present data from the SVS 13 star-forming region. The calculated envelope masses from the continuum show a decreasing trend with bolometric temperature (a proxy for age). Typical C$^{18}$O(2-1) linewidths are 1.45 km s$^{-1}$, which is higher than the C$^{18}$O linewidths detected toward Perseus filaments and cores. We find that N$_2$D$^+$(3-2) is significantly more likely to be detected toward younger protostars. We show that the protostars in SVS 13 are contained within filamentary structures as traced by C$^{18}$O(2-1) and N$_2$D$^+$(3-2). We also present the locations of SVS 13A's high velocity (absolute line-of-sight velocities $>$150 km s$^{-1}$) red and blue outflow components. Data can be downloaded from https://dataverse.harvard.edu/dataverse/MASSES .

As the complexity and volume of datasets have increased along with the capabilities of modular, open-source, easy-to-implement, visualization tools, scientists' need for, and appreciation of, data visualization has risen too. Until recently, scientists thought of the "explanatory" graphics created at a research project's conclusion as "pretty pictures" needed only for journal publication or public outreach. The plots and displays produced during a research project - often intended only for experts - were thought of as a separate category, what we here call "exploratory" visualization. In this view, discovery comes from exploratory visualization, and explanatory visualization is just for communication. Our aim in this paper is to spark conversation amongst scientists, computer scientists, outreach professionals, educators, and graphics and perception experts about how to foster flexible data visualization practices that can facilitate discovery and communication at the same time. We present an example of a new finding made using the glue visualization environment to demonstrate how the border between explanatory and exploratory visualization is easily traversed. The linked-view principles as well as the actual code in glue are easily adapted to astronomy, medicine, and geographical information science - all fields where combining, visualizing, and analyzing several high-dimensional datasets yields insight. Whether or not scientists can use such a flexible "undisciplined" environment to its fullest potential without special training remains to be seen. We conclude with suggestions for improving the training of scientists in visualization practices, and of computer scientists in the iterative, non-workflow-like, ways in which modern science is carried out.

The characterization of our Galaxy's longest filamentary gas features has been the subject of several studies in recent years, producing not only a sizeable sample of large-scale filaments, but also confusion as to whether all these features (e.g. "Bones", "Giant Molecular Filaments") are the same. They are not. We undertake the first standardized analysis of the physical properties ($\rm{H_2}\;$column densities, dust temperatures, morphologies, radial column density profiles) and kinematics of large-scale filaments in the literature. We expand and improve upon prior analyses by using the same data sets, techniques, and spiral arm models to disentangle the filaments' inherent properties from selection criteria and methodology. Our results suggest that the myriad filament finding techniques are uncovering different physical structures, with length ($11-269\,\rm\,pc$), width ($1-40\,\rm\,pc$), mass ($\rm3\times10^3\; M_\odot-1.1\times10^{6}\;M_\odot$), aspect ratio (3:1-117:1), and high column density fraction ($0.2-100\%$) varying by over an order of magnitude across the sample of 45 filaments. We develop a radial profile fitting code, $\texttt{RadFil}$, which is publicly available. We also perform a $\textit{position-position-velocity}$ ($\textit{p-p-v}$) analysis on a subsample and find that while $60-70\%$ lie spatially in the plane of the Galaxy, only $30-45\%$ concurrently exhibit spatial $\textit{and}$ kinematic proximity to spiral arms. In a parameter space defined by aspect ratio, dust temperature, and column density, we broadly distinguish three filament categories, which could indicate different formation mechanisms or histories. Highly elongated "Bone-like" filaments show the most potential for tracing gross spiral structure (e.g. arms, spurs), while other categories could be large concentrations of molecular gas (GMCs, core complexes).

We present the first results from a new, high resolution, $^{12}$CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) molecular line survey of the Orion A cloud, hereafter referred to as the CARMA-NRO Orion Survey. CARMA observations have been combined with single-dish data from the Nobeyama 45m telescope to provide extended images at about 0.01 pc resolution, with a dynamic range of approximately 1200 in spatial scale. Here we describe the practical details of the data combination in uv space, including flux scale matching, the conversion of single dish data to visibilities, and joint deconvolution of single dish and interferometric data. A $\Delta$-variance analysis indicates that no artifacts are caused by combining data from the two instruments. Initial analysis of the data cubes, including moment maps, average spectra, channel maps, position-velocity diagrams, excitation temperature, column density, and line ratio maps provides evidence of complex and interesting structures such as filaments, bipolar outflows, shells, bubbles, and photo-eroded pillars. The implications for star formation processes are profound and follow-up scientific studies by the CARMA-NRO Orion team are now underway. We plan to make all the data products described here generally accessible; some are already available at https://dataverse.harvard.edu/dataverse/CARMA-NRO-Orion

The American Astronomical Society's WorldWide Telescope (WWT) project enables terabytes of astronomical images, data, and stories to be viewed and shared among researchers, exhibited in science museums, projected into full-dome immersive planetariums and virtual reality headsets, and taught in classrooms from middle school to college levels. We review the WWT ecosystem, how WWT has been used in the astronomical community, and comment on future directions.

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

We use gas temperature and velocity dispersion data from the Green Bank Ammonia Survey and core masses and sizes from the James Clerk Maxwell Telescope Gould Belt Survey to estimate the virial states of dense cores within the Orion A molecular cloud. Surprisingly, we find that almost none of the dense cores are sufficiently massive to be bound when considering only the balance between self-gravity and the thermal and non-thermal motions present in the dense gas. Including the additional pressure binding imposed by the weight of the ambient molecular cloud material and additional smaller pressure terms, however, suggests that most of the dense cores are pressure confined.

The column density probability distribution function (N-PDF) of GMCs has been used as a diagnostic of star formation. Simulations and analytic predictions have suggested the N-PDF is composed of a low density lognormal component and a high density power-law component, tracing turbulence and gravitational collapse, respectively. In this paper, we study how various properties of the true 2D column density distribution create the shape, or "anatomy" of the PDF. We test our ideas and analytic approaches using both a real, observed, PDF based on Herschel observations of dust emission as well as a simulation that uses the ENZO code. Using a dendrogram analysis, we examine the three main components of the N-PDF: the lognormal component, the power-law component, and the transition point between these two components. We find that the power-law component of an N-PDF is the summation of N-PDFs of power-law substructures identified by the dendrogram algorithm. We also find that the analytic solution to the transition point between lognormal and power-law components proposed by Burkhart, Stalpes & Collins (2017) is applicable when tested on observations and simulations, within the uncertainties. We reconfirm and extend the results of Lombardi, Alves & Lada (2015), which stated that the lognormal component of the N-PDF is difficult to constrain due to the artificial choice of the map area. Based on the resulting anatomy of the N-PDF, we suggest avoiding analyzing the column density structures of a star forming region based solely on fits to the lognormal component of an N-PDF. We also suggest applying the N-PDF analysis in combination with the dendrogram algorithm, to obtain a more complete picture of the global and local environments and their effects on the density structures.

We present new Submillimeter Array (SMA) observations of CO(2-1) outflows toward young, embedded protostars in the Perseus molecular cloud as part of the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. For 57 Perseus protostars, we characterize the orientation of the outflow angles and compare them with the orientation of the local filaments as derived from $Herschel$ observations. We find that the relative angles between outflows and filaments are inconsistent with purely parallel or purely perpendicular distributions. Instead, the observed distribution of outflow-filament angles are more consistent with either randomly aligned angles or a mix of projected parallel and perpendicular angles. A mix of parallel and perpendicular angles requires perpendicular alignment to be more common by a factor of $\sim$3. Our results show that the observed distributions probably hold regardless of the protostar's multiplicity, age, or the host core's opacity. These observations indicate that the angular momentum axis of a protostar may be independent of the large-scale structure. We discuss the significance of independent protostellar rotation axes in the general picture of filament-based star formation.

We report ALMA observations of polarized dust emission from the protostellar source Ser-emb 8 at a linear resolution of 140 au. Assuming models of dust-grain alignment hold, the observed polarization pattern gives a projected view of the magnetic field structure in this source. Contrary to expectations based on models of strongly magnetized star formation, the magnetic field in Ser-emb 8 does not exhibit an hourglass morphology. Combining the new ALMA data with previous observational studies, we can connect magnetic field structure from protostellar core (~80,000 au) to disk (~100 au) scales. We compare our observations with four magnetohydrodynamic gravo-turbulence simulations made with the AREPO code that have initial conditions ranging from super-Alfv\'enic (weakly magnetized) to sub-Alfv\'enic (strongly magnetized). These simulations achieve the spatial dynamic range necessary to resolve the collapse of protostars from the parsec scale of star-forming clouds down to the ~100 au scale probed by ALMA. Only in the very strongly magnetized simulation do we see both the preservation of the field direction from cloud to disk scales and an hourglass-shaped field at < 1000 au scales. We conduct an analysis of the relative orientation of the magnetic field and the density structure in both the Ser-emb 8 ALMA observations and the synthetic observations of the four AREPO simulations. We conclude that the Ser-emb 8 data are most similar to the weakly magnetized simulations, which exhibit random alignment, in contrast to the strongly magnetized simulation, where the magnetic field plays a role in shaping the density structure in the source. In the weak-field case, it is turbulence -- not the magnetic field -- that shapes the material that forms the protostar, highlighting the dominant role that turbulence can play across many orders of magnitude in spatial scale.

There is a considerable deficiency in the number of known supernova remnants (SNRs) in the Galaxy compared to that expected. Searches for extended low-surface brightness radio sources may find new Galactic SNRs, but confusion with the much larger population of HII regions makes identifying such features challenging. SNRs can, however, be separated from HII regions using their significantly lower mid-infrared (MIR) to radio continuum intensity ratios. We use the combination of high-resolution 1-2 GHz continuum data from The HI, OH, Recombination line survey of the Milky Way (THOR) and lower-resolution VLA 1.4 GHz Galactic Plane Survey (VGPS) continuum data, together with MIR data from the Spitzer GLIMPSE, Spitzer MIPSGAL, and WISE surveys to identify SNR candidates. To ensure that the candidates are not being confused with HII regions, we exclude radio continuum sources from the WISE Catalog of Galactic HII Regions, which contains all known and candidate H II regions in the Galaxy. We locate 76 new Galactic SNR candidates in the THOR and VGPS combined survey area of 67.4deg>l>17.5deg, |b|<1.25deg and measure the radio flux density for 52 previously-known SNRs. The candidate SNRs have a similar spatial distribution to the known SNRs, although we note a large number of new candidates near l=30deg, the tangent point of the Scutum spiral arm. The candidates are on average smaller in angle compared to the known regions, 6.4'+/-4.7' versus 11.0'+/-7.8', and have lower integrated flux densities. If the 76 candidates are confirmed as true SNRs, for example using radio polarization measurements or by deriving radio spectral indices, this would more than double the number of known Galactic SNRs in the survey area. This large increase would still, however, leave a discrepancy between the known and expected SNR populations of about a factor of two.

We present an overview of the first data release (DR1) and first-look science from the Green Bank Ammonia Survey (GAS). GAS is a Large Program at the Green Bank Telescope to map all Gould Belt star-forming regions with $A_V \gtrsim 7$ mag visible from the northern hemisphere in emission from NH$_3$ and other key molecular tracers. This first release includes the data for four regions in Gould Belt clouds: B18 in Taurus, NGC 1333 in Perseus, L1688 in Ophiuchus, and Orion A North in Orion. We compare the NH$_3$ emission to dust continuum emission from Herschel, and find that the two tracers correspond closely. NH$_3$ is present in over 60\% of lines-of-sight with $A_V \gtrsim 7$ mag in three of the four DR1 regions, in agreement with expectations from previous observations. The sole exception is B18, where NH$_3$ is detected toward ~ 40\% of lines-of-sight with $A_V \gtrsim 7$ mag. Moreover, we find that the NH$_3$ emission is generally extended beyond the typical 0.1 pc length scales of dense cores. We produce maps of the gas kinematics, temperature, and NH$_3$ column densities through forward modeling of the hyperfine structure of the NH$_3$ (1,1) and (2,2) lines. We show that the NH$_3$ velocity dispersion, ${\sigma}_v$, and gas kinetic temperature, $T_K$, vary systematically between the regions included in this release, with an increase in both the mean value and spread of ${\sigma}_v$ and $T_K$ with increasing star formation activity. The data presented in this paper are publicly available.

A well-instrumented Green Bank Telescope (GBT) operating at high frequency represents a unique scientific resource for the US community. As a filled-aperture, 100m-diameter telescope, the GBT is ideally suited to fast mapping of extended, low surface brightness emission with excellent instantaneous frequency coverage. This capability makes the GBT a key facility for a range of cutting edge science described in this document, only possible at these frequencies. We note that the ability to perform the necessary observations is unique and highly complementary to the capabilities offered by interferometers, and should be preserved. We argue that rather than divesting from this exceptional resource, it makes sense for the US community to invest moderately to maintain GBT operations and to instrument it in an optimal manner, enabling it to become an extraordinary complement to existing and future radio interferometers. Adequately instrumented, the GBT would be a pillar for 20-115 GHz science in the US and the world.

Recently, Goodman et al. (2014) argued that the very long, very thin infrared dark cloud "Nessie" lies directly in the Galactic mid-plane and runs along the Scutum-Centaurus arm in position-position-velocity ($p-p-v$) space as traced by lower density $\rm {CO}$ and higher density ${\rm NH}_3$ gas. Nessie was presented as the first "bone" of the Milky Way, an extraordinarily long, thin, high-contrast filament that can be used to map our Galaxy's "skeleton." Here, we present evidence for additional bones in the Milky Way Galaxy, arguing that Nessie is not a curiosity but one of several filaments that could potentially trace Galactic structure. Our ten bone candidates are all long, filamentary, mid-infrared extinction features which lie parallel to, and no more than 20 pc from, the physical Galactic mid-plane. We use $\rm {CO}$, ${\rm N}_2{\rm H}^+$, $\rm {HCO}^+$, and ${\rm NH}_3$ radial velocity data to establish the three-dimensional location of the candidates in ${\it p-p-v}$ space. Of the ten candidates, six also: have a projected aspect ratio of $\geqq$50:1; run along, or extremely close to, the Scutum-Centaurus arm in ${\it p-p-v}$ space; and exhibit no abrupt shifts in velocity. The evidence presented here suggests that these candidates are marking the locations of significant spiral features, with the bone called filament 5 ("BC_18.88-0.09") being a close analog to Nessie in the Northern Sky. As molecular spectral-line and extinction maps cover more of the sky at increasing resolution and sensitivity, it should be possible to find more bones in future studies.

Star formation in molecular clouds can be triggered by the dynamical action of winds from massive stars. Furthermore, X-ray and UV fluxes from massive stars can influence the life time of surrounding circumstellar disks. We present the results of a 53 ks XMM-Newton observation centered on the Rho Ophiuchi A+B binary system. Rho Ophiiuchi lies in the center of a ring of dust, likely formed by the action of its winds. This region is different from the dense core of the cloud (L1688 Core F) where star formation is at work. X-rays are detected from Rho Ophiuchi as well as a group of surrounding X-ray sources. We detected 89 X-ray sources, 47 of them have at least one counterpart in 2MASS + All-WISE catalogs. Based on IR and X-ray properties, we can distinguish between young stellar objects (YSOs) belonging to the cloud and background objects. Among the cloud members, we detect 3 debris disk objects and 22 disk-less / Class III young stars. We show that these stars have ages in $5-10$ Myr, and are significantly older than the YSOs in L1688. We speculate that they are the result of an early burst of star formation in the cloud. An X-ray energy of $\ge5\times10^{44}$ ergs has been injected into the surrounding medium during the past $5$ Myr, we discuss the effects of such energy budget in relation to the cloud properties and dynamics.

We investigate the alignment between outflow axes in nine of the youngest binary/multiple systems in the Perseus Molecular Cloud. These systems have typical member spacing larger than 1000 AU. For outflow identification, we use 12CO(2-1) and 12CO(3-2) data from a large survey with the Sub- millimeter Array: Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES). The distribution of outflow orientations in the binary pairs is consistent with random or preferentially anti-aligned distributions, demonstrating that these outflows are misaligned. This result suggests that these systems are possibly formed in environments where the distribution of angular momentum is complex and disordered, and these systems do not come from the same co-rotating structures or from an initial cloud with aligned vectors of angular momentum.

The all-Galaxy CO survey of Dame, Hartmann, & Thaddeus (2001) is by far the most uniform, large-scale Galactic CO survey. Using a dendrogram-based decomposition of this survey, we present a catalog of 1064 massive molecular clouds throughout the Galactic plane. This catalog contains $2.5 \times 10^8$ solar masses, or $25^{+10.7}_{-5.8} \%$ of the Milky Way's estimated H$_2$ mass. We track clouds in some spiral arms through multiple quadrants. The power index of Larson's first law, the size-linewidth relation, is consistent with 0.5 in all regions - possibly due to an observational bias - but clouds in the inner Galaxy systematically have significantly (~ 30%) higher linewidths at a given size, indicating that their linewidths are set in part by Galactic environment. The mass functions of clouds in the inner Galaxy versus the outer Galaxy are both qualitatively and quantitatively distinct. The inner Galaxy mass spectrum is best described by a truncated power-law with a power index of $\gamma=-1.6\pm0.1$ and an upper truncation mass $M_0 = (1.0 \pm 0.2) \times 10^7 M_\odot$, while the outer Galaxy mass spectrum is better described by a non-truncating power law with $\gamma=-2.2\pm0.1$ and an upper mass $M_0 = (1.5 \pm 0.5) \times 10^6 M_\odot$, indicating that the inner Galaxy is able to form and host substantially more massive GMCs than the outer Galaxy. Additionally, we have simulated how the Milky Way would appear in CO from extragalactic perspectives, for comparison with CO maps of other galaxies.

Orion A hosts the nearest massive star factory, thus offering a unique opportunity to resolve the processes connected with the formation of both low- and high-mass stars. Here we present the most detailed and sensitive near-infrared (NIR) observations of the entire molecular cloud to date. With the unique combination of high image quality, survey coverage, and sensitivity, our NIR survey of Orion A aims at establishing a solid empirical foundation for further studies of this important cloud. In this first paper we present the observations, data reduction, and source catalog generation. To demonstrate the data quality, we present a first application of our catalog to estimate the number of stars currently forming inside Orion A and to verify the existence of a more evolved young foreground population. We used the European Southern Observatory's (ESO) Visible and Infrared Survey Telescope for Astronomy (VISTA) to survey the entire Orion A molecular cloud in the NIR $J, H$, and $K_S$ bands, covering a total of $\sim$18.3 deg$^2$. We implemented all data reduction recipes independently of the ESO pipeline. Estimates of the young populations toward Orion A are derived via the $K_S$-band luminosity function. Our catalog (799995 sources) increases the source counts compared to the Two Micron All Sky Survey by about an order of magnitude. The 90% completeness limits are 20.4, 19.9, and 19.0 mag in $J, H$, and $K_S$, respectively. The reduced images have 20% better resolution on average compared to pipeline products. We find between 2300 and 3000 embedded objects in Orion A and confirm that there is an extended foreground population above the Galactic field, in agreement with previous work. The Orion A VISTA catalog represents the most detailed NIR view of the nearest massive star-forming region and provides a fundamental basis for future studies of star formation processes toward Orion.

We present continuum and molecular line observations at 230 GHz and 345 GHz from the Sub-millimeter Array (SMA) toward three protostars in the Perseus L1448N region. The data are from the large project "Mass Assembly of Stellar Systems and their Evolution with the SMA" (MASSES). Three dust continuum sources, Source B, Source NW, and Source A, are detected at both frequencies. These sources have corresponding emission peaks in C18O (J=2-1), 13CO (J=2-1), and HCO+ (J=4-3), and have offsets with N2D+ (J=3-2) peaks. High angular resolution data from a complimentary continuum survey with the Karl G. Jansky Very Large Array show that Source B is associated with three 8 mm continuum objects, Source NW with two, and Source A remains single. These results suggest that multiplicity in L1448N exists at different spatial scales from a few thousand AU to < 100 AU. Velocity gradients in each source obtained from two-dimensional fits to the SMA C18O emission are found to be perpendicular to within 20 degrees of the outflow directions as revealed by 12CO (J=2-1). We have observed that Sources B and NW with multiplicity have higher densities than Source A without multiplicity. This suggests that thermal Jeans fragmentation can be relevant in the fragmentation process. However, we have not observed a difference in the ratio between rotational and gravitational energy between sources with and without multiplicity. We also have not observed a trend between non-thermal velocity dispersions and the level of fragmentation. Our study has provided the first direct and comprehensive comparison between multiplicity and core properties in low-mass protostars, although based on small number statistics.

It is now well-established that the elemental abundance patterns of stars holds key clues not only to their formation but also to the assembly histories of galaxies. One of the most exciting possibilities is the use of stellar abundance patterns as "chemical tags" to identify stars that were born in the same molecular cloud. In this paper we assess the prospects of chemical tagging as a function of several key underlying parameters. We show that in the fiducial case of $10^4$ distinct cells in chemical space and $10^5-10^6$ stars in the survey, one can expect to detect $\sim 10^2-10^3$ groups that are $\ge 5\sigma$ overdensities in the chemical space. However, we find that even very large overdensities in chemical space do not guarantee that the overdensity is due to a single set of stars from a common birth cloud. In fact, for our fiducial model parameters, the typical $5\sigma$ overdensity is comprised of stars from a wide range of clusters with the most dominant cluster contributing only 25% of the stars. The most important factors limiting the identification of disrupted clusters via chemical tagging are the number of chemical cells in the chemical space and the survey sampling rate of the underlying stellar population. Both of these factors can be improved through strategic observational plans. While recovering individual clusters through chemical tagging may prove challenging, we show, in agreement with previous work, that different CMFs imprint different degrees of clumpiness in chemical space. These differences provide the opportunity to statistically reconstruct the slope and high mass cutoff of CMF and its evolution through cosmic time.

The very long and thin infrared dark cloud "Nessie" is even longer than had been previously claimed, and an analysis of its Galactic location suggests that it lies directly in the Milky Way's mid-plane, tracing out a highly elongated bone-like feature within the prominent Scutum-Centaurus spiral arm. Re-analysis of mid-infrared imagery from the Spitzer Space Telescope shows that this IRDC is at least 2, and possibly as many as 5 times longer than had originally been claimed by Nessie's discoverers (Jackson et al. 2010); its aspect ratio is therefore at least 300:1, and possibly as large as 800:1. A careful accounting for both the Sun's offset from the Galactic plane ($\sim 25$ pc) and the Galactic center's offset from the $(l^{II},b^{II})=(0,0)$ position shows that the latitude of the true Galactic mid-plane at the 3.1 kpc distance to the Scutum-Centaurus Arm is not $b=0$, but instead closer to $b=-0.4$, which is the latitude of Nessie to within a few pc. An analysis of the radial velocities of low-density (CO) and high-density (${\rm NH}_3$) gas associated with the Nessie dust feature suggests that Nessie runs along the Scutum-Centaurus Arm in position-position-velocity space, which means it likely forms a dense `spine' of the arm in real space as well. The Scutum-Centaurus arm is the closest major spiral arm to the Sun toward the inner Galaxy, and, at the longitude of Nessie, it is almost perpendicular to our line of sight, making Nessie the easiest feature to see as a shadow elongated along the Galactic Plane from our location. Future high-resolution dust mapping and molecular line observations of the harder-to-find Galactic "bones" should allow us to exploit the Sun's position above the plane to gain a (very foreshortened) view "from above" of the Milky Way's structure.

X-rays from massive stars are ubiquitous yet not clearly understood. In an XMM-Newton observation devoted to observe the first site of star formation in the $\rho$ Ophiuchi dark cloud, we detect smoothly variable X-ray emission from the B2IV+B2V system of $\rho$ Ophiuchi. Tentatively we assign the emission to the primary component. The light curve of the pn camera shows a first phase of low, almost steady rate, then a rise phase of duration of 10 ks, followed by a high rate phase. The variability is seen primarily in the band 1.0-8.0 keV while little variability is detected below 1 keV. The spectral analysis of the three phases reveals the presence of a hot component at 3.0 keV that adds up to two relatively cold components at 0.9 keV and 2.2 keV. We explain the smooth variability with the emergence of an extended active region on the surface of the primary star due to its fast rotation (v $sin~i \sim315$ km/s). We estimate that the region has diameter in the range $0.5-0.6$ R$_*$. The hard X-ray emission and its variability hint a magnetic origin, as suggested for few other late-O$-$early-B type stars. We also discuss an alternative explanation based on the emergence from occultation of a young (5-10 Myr) low mass companion bright and hot in X-rays.

We present Brut, an algorithm to identify bubbles in infrared images of the Galactic midplane. Brut is based on the Random Forest algorithm, and uses bubbles identified by >35,000 citizen scientists from the Milky Way Project to discover the identifying characteristics of bubbles in images from the Spitzer Space Telescope. We demonstrate that Brut's ability to identify bubbles is comparable to expert astronomers. We use Brut to re-assess the bubbles in the Milky Way Project catalog, and find that 10-30% of the objects in this catalog are non-bubble interlopers. Relative to these interlopers, high-reliability bubbles are more confined to the mid plane, and display a stronger excess of Young Stellar Objects along and within bubble rims. Furthermore, Brut is able to discover bubbles missed by previous searches -- particularly bubbles near bright sources which have low contrast relative to their surroundings. Brut demonstrates the synergies that exist between citizen scientists, professional scientists, and machine learning techniques. In cases where "untrained" citizens can identify patterns that machines cannot detect without training, machine learning algorithms like Brut can use the output of citizen science projects as input training sets, offering tremendous opportunities to speed the pace of scientific discovery. A hybrid model of machine learning combined with crowdsourced training data from citizen scientists can not only classify large quantities of data, but also address the weakness of each approach if deployed alone.

The question whether magnetic fields play an important role in the processes of molecular cloud and star formation has been debated for decades. Recent observations have revealed a simple picture that may help illuminate these questions: magnetic fields have a tendency to preserve their orientation at all scales that have been probed - from 100-pc scale inter-cloud media down to sub-pc scale cloud cores. This ordered morphology has implications for the way in which self-gravity and turbulence interact with magnetic fields: both gravitational contraction and turbulent velocities should be anisotropic, due to the influence of dynamically important magnetic fields. Such anisotropy is now observed. Here we review these recent observations and discuss how they can improve our understanding of cloud/star formation.

This article offers a short guide to the steps scientists can take to ensure that their data and associated analyses continue to be of value and to be recognized. In just the past few years, hundreds of scholarly papers and reports have been written on questions of data sharing, data provenance, research reproducibility, licensing, attribution, privacy, and more, but our goal here is not to review that literature. Instead, we present a short guide intended for researchers who want to know why it is important to "care for and feed" data, with some practical advice on how to do that.

Educational technology has attained significant importance as a mechanism for supporting experiential learning of science concepts. However, the growth of this mechanism is limited by the significant time and technical expertise needed to develop such products, particularly in specialized fields of science. We sought to test whether interactive, educational, online software modules can be developed effectively by students as a curriculum component of an advanced science course. We discuss a set of fifteen such modules developed by Harvard University graduate students to demonstrate various concepts related to astronomy and physics. Their successful development of these modules demonstrates that online software tools for education and outreach on specialized topics can be produced while simultaneously fulfilling project-based learning objectives. We describe a set of technologies suitable for module development and present in detail four examples of modules developed by the students. We offer recommendations for incorporating educational software development within a graduate curriculum and conclude by discussing the relevance of this novel approach to new online learning environments like edX

The physical properties of molecular clouds are often measured using spectral-line observations, which provide the only probes of the clouds' velocity structure. It is hard, though, to assess whether and to what extent intensity features in position-position-velocity (PPV) space correspond to "real" density structures in position-position-position (PPP) space. In this paper, we create synthetic molecular cloud spectral-line maps of simulated molecular clouds, and present a new technique for measuring the reality of individual PPV structures. Our procedure projects density structures identified in PPP space into corresponding intensity structures in PPV space and then measures the geometric overlap of the projected structures with structures identified from the synthetic observation. The fractional overlap between a PPP and PPV structure quantifies how well the synthetic observation recovers information about the 3D structure. Applying this machinery to a set of synthetic observations of CO isotopes, we measure how well spectral-line measurements recover mass, size, velocity dispersion, and virial parameter for a simulated star-forming region. By disabling various steps of our analysis, we investigate how much opacity, chemistry, and gravity affect measurements of physical properties extracted from PPV cubes. For the simulations used here, our results suggest that superposition induces a ~40% uncertainty in masses, sizes, and velocity dispersions derived from 13CO. The virial parameter is most affected by superposition, such that estimates of the virial parameter derived from PPV and PPP information typically disagree by a factor of ~2. This uncertainty makes it particularly difficult to judge whether gravitational or kinetic energy dominate a given region, since the majority of virial parameter measurements fall within a factor of 2 of the equipartition level alpha ~ 2.

Magnetohydrodynamic turbulence is able to create hierarchical structures in the interstellar medium that are correlated on a wide range of scales via the energy cascade. We use hierarchical tree diagrams known as dendrograms to characterize structures in synthetic Position-Position-Velocity (PPV) emission cubes of optically thin isothermal magnetohydrodynamic turbulence. We show that the structures and degree of hierarchy observed in PPV space are related to the physics of the gas, i.e. self-gravity and the global sonic and Alfvenic Mach number. Simulations with higher Alfvenic Mach number, self-gravity and supersonic flows display enhanced hierarchical structure. We observed a strong sonic and Alfvenic dependency when we apply the the statistical moments (i.e. mean, variance, skewness, kurtosis) to the dendrogram distribution. Larger magnetic field and sonic Mach number correspond to larger values of the moments. Application of the dendrogram to 3D density cubes, also known as Position-Position-Position cubes (PPP), reveals that the dominant emission contours in PPP and PPV are related for supersonic gas but not for subsonic. We also explore the effects of smoothing, thermal broadening and velocity resolution on the dendrograms in order to make our study more applicable to observational data. These results all point to hierarchical tree diagrams as being a promising additional tool for studying ISM turbulence and star forming regions in the direction of obtaining information on the degree of self-gravity, the Mach numbers and the complicated relationship between PPV and PPP.

We investigate the shape of the extinction law in two 1-degree square fields of the Perseus Molecular Cloud complex. We combine deep red-optical (r, i, and z-band) observations obtained using Megacam on the MMT with UKIDSS near-infrared (J, H, and K-band) data to measure the colours of background stars. We develop a new hierarchical Bayesian statistical model, including measurement error, intrinsic colour variation, spectral type, and dust reddening, to simultaneously infer parameters for individual stars and characteristics of the population. We implement an efficient MCMC algorithm utilising generalised Gibbs sampling to compute coherent probabilistic inferences. We find a strong correlation between the extinction (Av) and the slope of the extinction law (parameterized by Rv). Because the majority of the extinction toward our stars comes from the Perseus molecular cloud, we interpret this correlation as evidence of grain growth at moderate optical depths. The extinction law changes from the diffuse value of Rv = 3 to the dense cloud value of Rv = 5 as the column density rises from Av = 2 mags to Av = 10 mags. This relationship is similar for the two regions in our study, despite their different physical conditions, suggesting that dust grain growth is a fairly universal process.

Astronomical researchers often think of analysis and visualization as separate tasks. In the case of high-dimensional data sets, though, interactive exploratory data visualization can give far more insight than an approach where data processing and statistical analysis are followed, rather than accompanied, by visualization. This paper attempts to charts a course toward "linked view" systems, where multiple views of high-dimensional data sets update live as a researcher selects, highlights, or otherwise manipulates, one of several open views. For example, imagine a researcher looking at a 3D volume visualization of simulated or observed data, and simultaneously viewing statistical displays of the data set's properties (such as an x-y plot of temperature vs. velocity, or a histogram of vorticities). Then, imagine that when the researcher selects an interesting group of points in any one of these displays, that the same points become a highlighted subset in all other open displays. Selections can be graphical or algorithmic, and they can be combined, and saved. For tabular (ASCII) data, this kind of analysis has long been possible, even though it has been under-used in Astronomy. The bigger issue for Astronomy and several other "high-dimensional" fields is the need systems that allow full integration of images and data cubes within a linked-view environment. The paper concludes its history and analysis of the present situation with suggestions that look toward cooperatively-developed open-source modular software as a way to create an evolving, flexible, high-dimensional, linked-view visualization environment useful in astrophysical research.

We present a hierarchical Bayesian method for fitting infrared spectral energy distributions (SEDs) of dust emission to observed fluxes. Under the standard assumption of optically thin single temperature (T) sources the dust SED as represented by a power--law modified black body is subject to a strong degeneracy between T and the spectral index beta. The traditional non-hierarchical approaches, typically based on chi-square minimization, are severely limited by this degeneracy, as it produces an artificial anti-correlation between T and beta even with modest levels of observational noise. The hierarchical Bayesian method rigorously and self-consistently treats measurement uncertainties, including calibration and noise, resulting in more precise SED fits. As a result, the Bayesian fits do not produce any spurious anti-correlations between the SED parameters due to measurement uncertainty. We demonstrate that the Bayesian method is substantially more accurate than the chi-square fit in recovering the SED parameters, as well as the correlations between them. As an illustration, we apply our method to Herschel and sub millimeter ground-based observations of the star-forming Bok globule CB244. This source is a small, nearby molecular cloud containing a single low-mass protostar and a starless core. We find that T and beta are weakly positively correlated -- in contradiction with the chi-square fits, which indicate a T-beta anti-correlation from the same data-set. Additionally, in comparison to the chi-square fits the Bayesian SED parameter estimates exhibit a reduced range in values.

Despite over 30 years of study, the mass-area relationship within and among clouds is still poorly understood both observationally and theoretically. Modern extinction datasets should have sufficient resolution and dynamic range to characterize this relationship for nearby molecular clouds, although recent papers using extinction data seem to yield different interpretations regarding the nature and universality of this aspect of cloud structure. In this paper we try to unify these various results and interpretations by accounting for the different ways cloud properties are measured and analyzed. We interpret the mass-area relationship in terms of the column density distribution function and its possible variation within and among clouds. We quantitatively characterize regional variations in the column density PDF. We show that structures both within and among clouds possess the same degree of "universality", in that their PDF means do not systematically scale with structure size. Because of this, mass scales linearly with area.

The WorldWide Telescope computer program, released to researchers and the public as a free resource in 2008 by Microsoft Research, has changed the way the ever-growing Universe of online astronomical data is viewed and understood. The WWT program can be thought of as a scriptable, interactive, richly visual browser of the multi-wavelength Sky as we see it from Earth, and of the Universe as we would travel within it. In its web API format, WWT is being used as a service to display professional research data. In its desktop format, WWT works in concert (thanks to SAMP and other IVOA standards) with more traditional research applications such as ds9, Aladin and TOPCAT. The WWT Ambassadors Program (founded in 2009) recruits and trains astrophysically-literate volunteers (including retirees) who use WWT as a teaching tool in online, classroom, and informal educational settings. Early quantitative studies of WWTA indicate that student experiences with WWT enhance science learning dramatically. Thanks to the wealth of data it can access, and the growing number of services to which it connects, WWT is now a key linking technology in the Seamless Astronomy environment we seek to offer researchers, teachers, and students alike.

Determining the initial stellar multiplicity is a challenging problem since protostars are faint and deeply embedded at early times; once formed, multiple protostellar systems may significantly dynamically evolve before they are optically revealed. Interferometers such as CARMA and ALMA make it possible to probe the scales at which turbulent fragmentation occurs in dust continuum emission, potentially constraining early stellar multiplicity. In this Letter we present synthetic observations of starless and protostellar cores undergoing fragmentation on scales of a few thousand AU to produce wide binary systems. We show that interferometric observations of starless cores by CARMA should be predominantly featureless at early stages, although wide protostellar companions should be apparent. The enhanced capabilities of ALMA improve the detection of core morphology so that it may be possible to detect substructure at earlier times. In either case, spatial filtering from interferometry reduces the observed core substructure and often eradicates traces of existing filamentary morphology on scales down to 0.025 pc. However, some missing structure may be recaptured by combining data from the ALMA full science and Atacama compact arrays.

The ADS All-Sky Survey (ADSASS) is an ongoing effort aimed at turning the NASA Astrophysics Data System (ADS), widely known for its unrivaled value as a literature resource for astronomers, into a data resource. The ADS is not a data repository per se, but it implicitly contains valuable holdings of astronomical data, in the form of images, tables and object references contained within articles. The objective of the ADSASS effort is to extract these data and make them discoverable and available through existing data viewers. The resulting ADSASS data layer promises to greatly enhance workflows and enable new research by tying astronomical literature and data assets into one resource.

We present results from three-dimensional, self-gravitating, radiation-hydrodynamic simulations of low-mass protostellar outflows. We construct synthetic observations in 12CO in order to compare with observed outflows and evaluate the effects of beam resolution and outflow orientation on inferred outflow properties. To facilitate the comparison, we develop a quantitative prescription for measuring outflow opening angles. Using this prescription, we demonstrate that, in both simulations and synthetic observations, outflow opening angles broaden with time similarly to observed outflows. However, the interaction between the outflowing gas and the turbulent core envelope produces significant asymmetry between the red and blue shifted outflow lobes. We find that applying a velocity cutoff may result in outflow masses that are underestimated by a factor 5 or more, and masses derived from optically thick CO emission further underpredict the mass of the high-velocity gas by a factor of 5-10. Derived excitation temperatures indicate that outflowing gas is hotter than the ambient gas with temperature rising over time, which is in agreement with the simulation gas temperatures. However, excitation temperatures are otherwise not well correlated with the actual gas temperature.

We present a study on the shells (and bubbles) in the Perseus molecular cloud using the COMPLETE survey large-scale 12CO(1-0) and 13CO(1-0) maps. The twelve shells reported here are spread throughout most of the Perseus cloud and have circular or arc-like morphologies with a range in radius of about 0.1 to 3 pc. Most of them have not been detected before most likely as maps of the region lacked the coverage and resolution needed to distinguish them. The majority of the shells are coincident with infrared nebulosity of similar shape and have a candidate powering source near the center. We suggest they are formed by the interaction of spherical or very wide-angle winds powered by young stars inside or near the Perseus molecular cloud -a cloud that is commonly considered to be mostly forming low-mass stars. Two of the twelve shells are powered by high-mass stars close to the cloud, while the others appear to be powered by low or intermediate mass stars in the cloud. We argue that winds with a mass loss rate of about 10^-8 to 10^-6 M_sun/yr are required to produce the observed shells. Our estimates indicate that the energy input rate from these stellar winds is similar to the turbulence dissipation rate. We conclude that in Perseus the total energy input from both collimated protostellar outflows and powerful spherical winds from young stars is sufficient to maintain the turbulence in the molecular cloud. Large scale molecular line and IR continuum maps of a sample of clouds will help determine the frequency of this phenomenon in other star forming regions.

We present the detection of a dust continuum source at 3-mm (CARMA) and 1.3-mm (SMA), and 12CO(2-1) emission (SMA) towards the L1451-mm dense core. These detections suggest a compact object and an outflow where no point source at mid-infrared wavelengths is detected using Spitzer. An upper limit for the dense core bolometric luminosity of 0.05 Lsun is obtained. By modeling the broadband SED and the continuum interferometric visibilities simultaneously, we confirm that a central source of heating is needed to explain the observations. This modeling also shows that the data can be well fitted by a dense core with a YSO and disk, or by a dense core with a central First Hydrostatic Core (FHSC). Unfortunately, we are not able to decide between these two models, which produce similar fits. We also detect 12CO(2-1) emission with red- and blue-shifted emission suggesting the presence of a slow and poorly collimated outflow, in opposition to what is usually found towards young stellar objects but in agreement with prediction from simulations of a FHSC. This presents the best candidate, so far, for a FHSC, an object that has been identified in simulations of collapsing dense cores. Whatever the true nature of the central object in L1451-mm, this core presents an excellent laboratory to study the earliest phases of low-mass star formation.

The VAO (Virtual Astronomical Observatory) Science Council (VAO-SC) met on July 27-28, 2011 at the Harvard-Smithsonian Center for Astrophysics in Cambridge MA, to review the VAO performance during its first year of operations. In this meeting the VAO demonstrated the new tools for astronomers that are being released in September 2011 and presented plans for the second year of activities, resulting from studies conducted during the first year. This document contains the recommendations of the VAO-SC for the second year of activity of the VAO.

We apply Support Vector Machines -- a machine learning algorithm -- to the task of classifying structures in the Interstellar Medium. As a case study, we present a position-position velocity data cube of 12 CO J=3--2 emission towards G16.05-0.57, a supernova remnant that lies behind the M17 molecular cloud. Despite the fact that these two objects partially overlap in position-position-velocity space, the two structures can easily be distinguished by eye based on their distinct morphologies. The Support Vector Machine algorithm is able to infer these morphological distinctions, and associate individual pixels with each object at >90% accuracy. This case study suggests that similar techniques may be applicable to classifying other structures in the ISM -- a task that has thus far proven difficult to automate.

We review an approach to observation-theory comparisons we call "Taste-Testing." In this approach, synthetic observations are made of numerical simulations, and then both real and synthetic observations are "tasted" (compared) using a variety of statistical tests. We first lay out arguments for bringing theory to observational space rather than observations to theory space. Next, we explain that generating synthetic observations is only a step along the way to the quantitative, statistical, taste tests that offer the most insight. We offer a set of examples focused on polarimetry, scattering and emission by dust, and spectral-line mapping in starforming regions. We conclude with a discussion of the connection between statistical tests used to date and the physics we seek to understand. In particular, we suggest that the "lognormal" nature of molecular clouds can be created by the interaction of many random processes, as can the lognormal nature of the IMF, so that the fact that both the "Clump Mass Function" (CMF) and IMF appear lognormal does not necessarily imply a direct relationship between them.

We present a ~6.5'x8' Expanded Very Large Array (EVLA) mosaic observations of the NH3 (1,1) emission in the Barnard 5 region in Perseus, with an angular resolution of 6". This map covers the coherent region, where the dense gas presents subsonic non-thermal motions (as seen from single dish observations with the Green Bank Telescope, GBT). The combined EVLA and GBT observations reveal, for the first time, a striking filamentary structure (20" wide or 5,000 AU at the distance of Perseus) in this low-mass star forming region. The integrated intensity profile of this structure is consistent with models of an isothermal filament in hydrostatic equilibrium. The observed separation between the B5-IRS1 young stellar object (YSO), in the central region of the core, and the northern starless condensation matches the Jeans length of the dense gas. This suggests that the dense gas in the coherent region is fragmenting. The region observed displays a narrow velocity dispersion, where most of the gas shows evidence for subsonic turbulence, and where little spatial variations are present. It is only close to the YSO where an increase in the velocity dispersion is found, but still displaying subsonic non-thermal motions

We present a study of dense structures in the L1495 filament in the Taurus Molecular Cloud and examine its star-forming properties. In particular we construct a dust extinction map of the filament using deep near-infrared observations, exposing its small-scale structure in unprecedented detail. The filament shows highly fragmented substructures and a high mass-per-length value of M_line=17 Msun/pc, reflecting star-forming potential in all parts of it. However, a part of the filament, namely B211, is remarkably devoid of young stellar objects. We argue that in this region the initial filament collapse and fragmentation is still taking place and star formation is yet to occur. In the star-forming part of the filament, we identify 39 cores with masses from 0.4...10 Msun and preferred separations in agreement with the local Jeans length. Most of these cores exceed the Bonnor-Ebert critical mass, and are therefore likely to collapse and form stars. The Dense Core Mass Function follows a power law with exponent Gamma=1.2, a form commonly observed in star-forming regions.

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}.

We present two Very Large Array observations of the pre-main-sequence star PV Cephei, taken with a separation of 10.5 years. These data show that the proper motions of this star are $\mu_\alpha \cos \delta = +10.9 \pm 3.0$ mas yr$^{-1}$; $\mu_\delta = +0.2 \pm 1.8$ mas yr$^{-1}$, very similar to those -previously known- of HD 200775, the B2Ve star that dominates the illumination of the nearby reflection nebula NGC 7023. This result suggests that PV Cephei is not a rapidly moving run-away star as suggested by previous studies. The large velocity of PV Cephei had been inferred from the systematic eastward displacement of the bisectors of successive pairs of Herbig Haro knots along its flow. These systematic shifts might instead result from an intrinsic dissymmetry in the ejection mechanisms, or from an asymmetric distribution of the circumstellar material.

We study the rotational properties of magnetized and self-gravitating molecular cloud cores formed in 2 very high resolution 3D molecular cloud simulations.The simulations have been performed using the code RAMSES at an effective resolution of 4096^3.One simulation represents a mildly magnetically-supercritical cloud and the other a strongly magnetically-supercritical cloud.A noticeable difference between the 2 simulations is the core formation efficiency (CFE) of the high density cores.In the strongly supercritical simulations the CFE is ~33 % per free-fall time of the cloud tff,cl, whereas in the mildly supercritical simulations this value goes down to ~6%/tff,cl. A comparison of the intrinsic specific angular momentum j3D distributions of the cores with the distribitions of j2D derived using synthetic 2D velocity maps of the cores,shows that the synthetic observations tend to overestimate the true value of j by a factor of ~10.The origin of this discrepancy lies in the fact that contrary to the intrinsic determination which sums up the individual gas parcels contributions to j, the determination of j using the observational procedure which is based on a measurement on the global velocity gradient under the hypothesis of uniform rotation smoothes out the complex fluctuations present in the 3D velocity field. Our results provide a natural explanation for the discrepancy by a factor ~10 observed between the intrinsic 3D distributions of j and the corresponding distributions derived in real observations.We suggest that measurements of j which are based on the measurement of the observed global velocity gradients may need to be reduced by a factor of ~10 in order to derive a more accurate estimate of j in the cores.

Recommendations of the VAO-Science Council following the meeting of March 26-27, 2010. Refer to the webpage http://www.aui.edu/vao.php

We present a study on the impact of molecular outflows in the Perseus molecular cloud complex using the COMPLETE survey large-scale 12CO(1-0) and 13CO(1-0) maps. We used three-dimensional isosurface models generated in RA-DEC-Velocity space to visualize the maps. This rendering of the molecular line data allowed for a rapid and efficient way to search for molecular outflows over a large (~ 16 sq. deg.) area. Our outflow-searching technique detected previously known molecular outflows as well as new candidate outflows. Most of these new outflow-related high-velocity features lie in regions that have been poorly studied before. These new outflow candidates more than double the amount of outflow mass, momentum, and kinetic energy in the Perseus cloud complex. Our results indicate that outflows have significant impact on the environment immediately surrounding localized regions of active star formation, but lack the energy needed to feed the observed turbulence in the entire Perseus complex. This implies that other energy sources, in addition to protostellar outflows, are responsible for turbulence on a global cloud scale in Perseus. We studied the impact of outflows in six regions with active star formation within Perseus of sizes in the range of 1 to 4 pc. We find that outflows have enough power to maintain the turbulence in these regions and enough momentum to disperse and unbind some mass from them. We found no correlation between outflow strength and star formation efficiency for the six different regions we studied, contrary to results of recent numerical simulations. The low fraction of gas that potentially could be ejected due to outflows suggests that additional mechanisms other than cloud dispersal by outflows are needed to explain low star formation efficiencies in clusters.

We measure the mass and size of cloud fragments in several molecular clouds continuously over a wide range of spatial scales (0.05 < r / pc < 3). Based on the recently developed "dendrogram-technique", this characterizes dense cores as well as the enveloping clouds. "Larson's 3rd Law" of constant column density, m(r) = C*r^2, is not well suited to describe the derived mass-size data. Solar neighborhood clouds not forming massive stars (< 10 M_sun; Pipe Nebula, Taurus, Perseus, and Ophiuchus) obey m(r) < 870 M_sun (r / pc)^1.33 . In contrast to this, clouds forming massive stars (Orion A, G10.15$-$0.34, G11.11$-$0.12) do exceed the aforementioned relation. Thus, this limiting mass-size relation may approximate a threshold for the formation of massive stars. Across all clouds, cluster-forming cloud fragments are found to be---at given radius---more massive than fragments devoid of clusters. The cluster-bearing fragments are found to roughly obey a mass-size law m = C*r^1.27 (where the exponent is highly uncertain in any given cloud, but is certainly smaller than 1.5).

We present NH3 observations of the B5 region in Perseus obtained with the Green Bank Telescope (GBT). The map covers a region large enough (~11'x14') that it contains the entire dense core observed in previous dust continuum surveys. The dense gas traced by NH3(1,1) covers a much larger area than the dust continuum features found in bolometer observations. The velocity dispersion in the central region of the core is small, presenting subsonic non-thermal motions which are independent of scale. However, it is thanks to the coverage and high sensitivity of the observations that we present the detection, for the first time, of the transition between the coherent core and the dense but more turbulent gas surrounding it. This transition is sharp, increasing the velocity dispersion by a factor of 2 in less than 0.04 pc (the 31" beam size at the distance of Perseus, ~250 pc). The change in velocity dispersion at the transition is ~3 km/s/pc. The existence of the transition provides a natural definition of dense core: the region with nearly-constant subsonic non-thermal velocity dispersion. From the analysis presented here we can not confirm nor rule out a corresponding sharp density transition.

Power law mass-size and linewidth-size correlations, two of "Larson's laws," are often studied to assess the dynamical state of clumps within molecular clouds. Using the result of a hydrodynamic simulation of a molecular cloud, we investigate how geometric projection may affect the derived Larson relationships. We find that large scale structures in the column density map have similar masses and sizes to those in the 3D simulation (PPP). Smaller scale clumps in the column density map are measured to be more massive than the PPP clumps, due to the projection of all emitting gas along lines of sight. Further, due to projection effects, structures in a synthetic spectral observation (PPV) may not necessarily correlate with physical structures in the simulation. In considering the turbulent velocities only, the linewidth-size relationship in the PPV cube is appreciably different from that measured from the simulation. Including thermal pressure in the simulated linewidths imposes a minimum linewidth, which results in a better agreement in the slopes of the linewidth-size relationships, though there are still discrepancies in the offsets, as well as considerable scatter. Employing commonly used assumptions in a virial analysis, we find similarities in the computed virial parameters of the structures in the PPV and PPP cubes. However, due to the discrepancies in the linewidth- and mass- size relationships in the PPP and PPV cubes, we caution that applying a virial analysis to observed clouds may be misleading due to geometric projection effects. We speculate that consideration of physical processes beyond kinetic and gravitational pressure would be required for accurately assessing whether complex clouds, such as those with highly filamentary structure, are bound.

We use a new contour-based map analysis technique to measure the mass and size of molecular cloud fragments continuously over a wide range of spatial scales (0.05 < r / pc < 10), i.e., from the scale of dense cores to those of entire clouds. The present paper presents the method via a detailed exploration of the Perseus Molecular Cloud. Dust extinction and emission data are combined to yield reliable scale-dependent measurements of mass. This scale-independent analysis approach is useful for several reasons. First, it provides a more comprehensive characterization of a map (i.e., not biased towards a particular spatial scale). Such a lack of bias is extremely useful for the joint analysis of many data sets taken with different spatial resolution. This includes comparisons between different cloud complexes. Second, the multi-scale mass-size data constitutes a unique resource to derive slopes of mass-size laws (via power-law fits). Such slopes provide singular constraints on large-scale density gradients in clouds.

The ability to represent scientific data and concepts visually is becoming increasingly important due to the unprecedented exponential growth of computational power during the present digital age. The data sets and simulations scientists in all fields can now create are literally thousands of times as large as those created just 20 years ago. Historically successful methods for data visualization can, and should, be applied to today's huge data sets, but new approaches, also enabled by technology, are needed as well. Increasingly, "modular craftsmanship" will be applied, as relevant functionality from the graphically and technically best tools for a job are combined as-needed, without low-level programming.

In this paper we present a dust emission map of the starless core TMC-1C taken at 2100 microns. Along with maps at 160, 450, 850 and 1200 microns, we study the dust emissivity spectral index from the (sub)millimeter spectral energy distribution, and find that it is close to the typically assumed value of beta = 2. We also map the dust temperature and column density in TMC-1C, and find that at the position of the dust peak (A_V ~ 50), the line-of-sight-averaged temperature is ~7 K. Employing simple Monte Carlo modeling, we show that the data are consistent with a constant value for the emissivity spectral index over the whole map of TMC-1C.

One of the key problems in star formation research is to determine the role of magnetic fields. Starting from the atomic inter-cloud medium (ICM) which has density nH ~ 1 per cubic cm, gas must accumulate from a volume several hundred pc across in order to form a typical molecular cloud. Star formation usually occurs in cloud cores, which have linear sizes below 1 pc and densities nH2 > 10^5 per cubic cm. With current technologies, it is hard to probe magnetic fields at scales lying between the accumulation length and the size of cloud cores, a range corresponds to many levels of turbulent eddy cascade, and many orders of magnitude of density amplification. For field directions detected from the two extremes, however, we show here that a significant correlation is found. Comparing this result with molecular cloud simulations, only the sub-Alfvenic cases result in field orientations consistent with our observations.

We study the mass spectrum of sub-structures in the Perseus Molecular Cloud Complex traced by 13CO (1-0), finding that $dN/dM\propto M^{-2.4}$ for the standard Clumpfind parameters. This result does not agree with the classical $dN/dM\propto M^{-1.6}$. To understand this discrepancy we study the robustness of the mass spectrum derived using the Clumpfind algorithm. Both 2D and 3D Clumpfind versions are tested, using 850 $\mu$m dust emission and 13CO spectral-line observations of Perseus, respectively. The effect of varying threshold is not important, but varying stepsize produces a different effect for 2D and 3D cases. In the 2D case, where emission is relatively isolated (associated with only the densest peaks in the cloud), the mass spectrum variability is negligible compared to the mass function fit uncertainties. In the 3D case, however, where the 13CO emission traces the bulk of the molecular cloud, the number of clumps and the derived mass spectrum are highly correlated with the stepsize used. The distinction between "2D" and "3D" here is more importantly also a distinction between "sparse" and "crowded" emission. In any "crowded" case, Clumpfind should not be used blindly to derive mass functions. Clumpfind's output in the "crowded" case can still offer a statistical description of emission useful in inter-comparisons, but the clump-list should not be treated as a robust region decomposition suitable to generate a physically-meaningful mass function. We conclude that the 13CO mass spectrum depends on the observations resolution, due to the hierarchical structure of MC.

This White Paper, submitted to the National Academy of Sciences' Astro2010 Decadal Review Committee, focuses on 2 central themes in the study of young brown dwarfs -- their formation mechanism and disk characteristics -- which are of direct relevance to fundamental questions of stellar and planetary origins and properties.

We investigate the effect of line of sight temperature variations and noise on two commonly used methods to determine dust properties from dust continuum observations of dense cores. One method employs a direct fit to a modified blackbody SED; the other involves a comparison of flux ratios to an analytical prediction. Fitting fluxes near the SED peak produces inaccurate temperature and dust spectral index estimates due to the line of sight temperature (and density) variations. Longer wavelength fluxes in the Rayleigh-Jeans part of the spectrum (>~ 600 micron for typical cores) may more accurately recover the spectral index, but both methods are very sensitive to noise. The temperature estimate approaches the density weighted temperature, or "column temperature," of the source as short wavelength fluxes are excluded. An inverse temperature - spectral index correlation naturally results from SED fitting, due to the inaccurate isothermal assumption, as well as noise uncertainties. We show that above some "threshold" temperature, the temperatures estimated through the flux ratio method can be highly inaccurate. In general, observations with widely separated wavelengths, and including shorter wavelengths, result in higher threshold temperatures; such observations thus allow for more accurate temperature estimates of sources with temperatures less than the threshold temperature. When only three fluxes are available, a constrained fit, where the spectral index is fixed, produces less scatter in the temperature estimate when compared to the estimate from the flux ratio method.

We investigate how uncertainties in flux measurements affect the results from modified blackbody SED fits. We show that an inverse correlation between the dust temperature T and spectral index (beta) naturally arises from least squares fits due to the uncertainties, even for sources with a single T and beta. Fitting SEDs to noisy fluxes solely in the Rayleigh-Jeans regime produces unreliable T and beta estimates. Thus, for long wavelength observations (lambda >~ 200 micron), or for warm sources (T >~ 60 K), it becomes difficult to distinguish sources with different temperatures. We assess the role of noise in recent observational results that indicate an inverse and continuously varying T - beta relation. Though an inverse and continuous T - beta correlation may be a physical property of dust in the ISM, we find that the observed inverse correlation may be primarily due to noise.

We present the chemistry, temperature, and dynamical state of a sample of 193 dense cores or core candidates in the Perseus Molecular cloud and compare the properties of cores associated with young stars and clusters with those which are not. The combination of our NH3 and CCS observations with previous millimeter, sub-millimeter, and Spitzer data available for this cloud enable us both to determine core properties precisely and to accurately classify cores as starless or protostellar. The properties of cores in different cluster environments and before-and-after star formation provide important constraints on simulations of star-formation, particularly under the paradigm that the essence of star formation is set by the turbulent formation of prestellar cores. We separate the influence of stellar content from that of cluster environment and find that cores within clusters have (1) higher kinetic temperatures and (2) lower fractional abundances of CCS and NH3. Cores associated with protostars have (1) slightly higher kinetic temperatures (2) higher NH3 excitation temperatures), (3) are at higher column density, have (4) slightly more non-thermal/turbulent NH3 linewidths, have (5) higher masses and have (6) lower fractional abundance of CCS. We find that neither cluster environment nor protostellar content makes a significant difference to the dynamical state of cores as estimated by the virial parameter -- most cores in each category are gravitationally bound. Overall, cluster environment and protostellar content have a smaller influence on the properties of the cores than is typically assumed, and the variation within categories is larger than the differences between categories.

We use the COMPLETE Survey's observations of the Perseus star-forming region to assess and intercompare three methods for measuring column density in molecular clouds: extinction mapping (NIR); thermal emission mapping (FIR); and mapping the intensity of CO isotopologues. The structures shown by all three tracers are morphologically similar, but important differences exist. Dust-based measures give similar, log-normal, distributions for the full Perseus region, once careful calibration corrections are made. We also compare dust- and gas-based column density distributions for physically-meaningful sub-regions of Perseus, and we find significant variations in the distributions for those regions. Even though we have used 12CO data to estimate excitation temperatures, and we have corrected for opacity, the 13CO maps seem unable to give column distributions that consistently resemble those from dust measures. We have edited out the effects of the shell around the B-star HD 278942. In that shell's interior and in the parts where it overlaps the molecular cloud, there appears to be a dearth of 13CO, likely due either to 13CO not yet having had time to form in this young structure, and/or destruction of 13CO in the molecular cloud. We conclude that the use of either dust or gas measures of column density without extreme attention to calibration and artifacts is more perilous than even experts might normally admit. And, the use of 13CO to trace total column density in detail, even after proper calibration, is unavoidably limited in utility due to threshold, depletion, and opacity effects. If one's main aim is to map column density, then dust extinction seems the best probe. Linear fits amongst column density tracers are given, quantifying the inherent uncertainties in using one tracer (when compared with others). [abridged]

Using far-infrared emission maps taken by IRAS and Spitzer and a near-infrared extinction map derived from 2MASS data, we have made dust temperature and column density maps of the Perseus molecular cloud. We show that the emission from transiently heated very small grains and the big grain dust emissivity vary as a function of extinction and dust temperature, with higher dust emissivities for colder grains. This variable emissivity can not be explained by temperature gradients along the line of sight or by noise in the emission maps, but is consistent with grain growth in the higher density and lower temperature regions. By accounting for the variations in the dust emissivity and VSG emission, we are able to map the temperature and column density of a nearby molecular cloud with better accuracy than has previously been possible.

We demonstrate the utility of dendrograms at representing the essential features of the hierarchical structure of the isosurfaces for molecular line data cubes. The dendrogram of a data cube is an abstraction of the changing topology of the isosurfaces as a function of contour level. The ability to track hierarchical structure over a range of scales makes this analysis philosophically different from local segmentation algorithms like CLUMPFIND. Points in the dendrogram structure correspond to specific volumes in data cubes defined by their bounding isosurfaces. We further refine the technique by measuring the properties associated with each isosurface in the analysis allowing for a multiscale calculation of molecular gas properties. Using COMPLETE 13CO(1-0) data from the L1448 region in Perseus and mock observations of a simulated data cube, we identify regions that have a significant contribution by self-gravity to their energetics on a range of scales. We find evidence for self-gravitation on all spatial scales in L1448 though not in all regions. In the simulated observations, nearly all of the emission is found in objects that would be self-gravitating if gravity were included in the simulation. We reconstruct the size-line width relationship within the data cube using the dendrogram-derived properties and find it follows the standard relation: s_v ~ R^0.58. Finally, we show that constructing the dendrogram of CO J=1-0 emission from the Orion-Monoceros region allows for the identification of giant molecular clouds in a blended molecular line data set using only a physically motivated definition (self-gravitating clouds with masses 5x10^4 Msun.

We use the COMPLETE data to derive new calibrations of the X-factor and the 13CO abundance within Perseus. We divide Perseus into six sub-regions. The standard X factor, X=N(H2)/W(12CO), is derived both for the whole Perseus Complex and for each of the six sub-regions with values consistent with previous estimates. The X factor is heavily affected by the saturation of the emission above AV~4 mag, and variations are found between regions. We derive linear fits to relate W(12CO) and AV using only points below 4 mag of extinction, this yields a better estimation of the AV than the X-factor. We derive linear relations of W(13CO), N(13CO) and W(C18O) with AV . The extinction threshold above which 13CO(1-0) and C18O(1-0) are detected is about 1 mag larger than previous estimates. 12CO(1-0) and 13CO(1-0) lines saturate above 4 and 5 mag, respectively, whereas C18O(1-0) never saturates (up to 10 mag). Approximately 60% of the positions with 12CO emission have sub-thermally excited lines, and almost all positions have 12CO excitation temperatures below the dust temperature. Using the Meudon PDR code we find that 12CO and 13CO emission can be explained by uniform slab models with densities ranging between about 10^3 and 10^4 cm-3. Local variations in the volume density and non-thermal motions (linked to different star formation activity) can explain the observations. Higher densities are needed to reproduce CO data toward active star forming sites, where the larger internal motions driven by the young protostars allow more photons from the embedded high density cores to escape the cloud. In the most quiescent region, the 12CO and 13CO emission appears to arise from an almost uniform thin layer of molecular material at densities around 10^4 cm-3.

We present ammonia observations of 193 dense cores and core candidates in the Perseus molecular cloud made using the Robert F. Byrd Green Bank Telescope. We simultaneously observed the NH3(1,1), NH3(2,2), CCS (2_1 -> 1_0) and CC34S (2_1 -> 1_0) transitions near 23 GHz for each of the targets with a spectral resolution of dv ~ 0.024 km/s. We find ammonia emission associated with nearly all of the (sub)millimeter sources as well as at several positions with no associated continuum emission. For each detection, we have measured physical properties by fitting a simple model to every spectral line simultaneously. Where appropriate, we have refined the model by accounting for low optical depths, multiple components along the line of sight and imperfect coupling to the GBT beam. For the cores in Perseus, we find a typical kinetic temperature of T=11 K, a typical column density of N(NH3)~ 10^14.5 /cm^2 and velocity dispersions ranging from sigma_v = 0.07 km/s to 0.7 km/s. However, many cores with velocity dispersions > 0.2 km/s show evidence for multiple velocity components along the line of sight.

In studies of star-forming regions, near-infrared excess (NIRX) sources--objects with intrinsic colors redder than normal stars--constitute both signal (young stars) and noise (e.g. background galaxies). We hunt down (identify) galaxies using near-infrared observations in the Perseus star-forming region by combining structural information, colors, and number density estimates. Galaxies at moderate redshifts (z = 0.1 - 0.5) have colors similar to young stellar objects (YSOs) at both near- and mid-infrared (e.g. Spitzer) wavelengths, which limits our ability to identify YSOs from colors alone. Structural information from high-quality near-infrared observations allows us to better separate YSOs from galaxies, rejecting 2/5 of the YSO candidates identified from Spitzer observations of our regions and potentially extending the YSO luminosity function below K of 15 magnitudes where galaxy contamination dominates. Once they are identified we use galaxies as valuable extra signal for making extinction maps of molecular clouds. Our new iterative procedure: the Galaxies Near Infrared Color Excess method Revisited (GNICER), uses the mean colors of galaxies as a function of magnitude to include them in extinction maps in an unbiased way. GNICER increases the number of background sources used to probe the structure of a cloud, decreasing the noise and increasing the resolution of extinction maps made far from the galactic plane.

We have mapped the starless core TMC-1C in a variety of molecular lines with the IRAM 30m telescope. High density tracers show clear signs of self-absorption and sub-sonic infall asymmetries are present in N2H+ (1-0) and DCO+ (2-1) lines. The inward velocity profile in N2H+ (1-0) is extended over a region of about 7,000 AU in radius around the dust continuum peak, which is the most extended ``infalling'' region observed in a starless core with this tracer. The kinetic temperature (~12 K) measured from C17O and C18O suggests that their emission comes from a shell outside the colder interior traced by the mm continuum dust. The C18O (2-1) excitation temperature drops from 12 K to ~10 K away from the center. This is consistent with a volume density drop of the gas traced by the C18O lines, from ~4x10^4 cm^-3 towards the dust peak to ~6x10^3 cm^-3 at a projected distance from the dust peak of 80" (or 11,000 AU). The column density implied by the gas and dust show similar N2H+ and CO depletion factors (f_D < 6). This can be explained with a simple scenario in which: (i) the TMC-1C core is embedded in a relatively dense environment (H2 ~10^4 cm^-3), where CO is mostly in the gas phase and the N2H+ abundance had time to reach equilibrium values; (ii) the surrounding material (rich in CO and N2H+) is accreting onto the dense core nucleus; (iii) TMC-1C is older than 3x10^5 yr, to account for the observed abundance of N2H+ across the core (~10^-10 w.r.t. H2); and (iv) the core nucleus is either much younger (~10^4 yr) or ``undepleted'' material from the surrounding envelope has fallen towards it in the past 10,000 yr.

We have mapped the central 10'x10' of the dense core TMC-1C at 450, 850 and 1200 microns using SCUBA on the James Clerk Maxwell Telescope and MAMBO on the IRAM 30m telescope. We show that although one can, in principle, use images at these wavelengths to map the emissivity spectral index, temperature and column density independently, noise and calibration errors would have to be less than 2% to accurately derive these three quantities from a set of three emission maps. Because our data are not this free of errors, we use our emission maps to fit the dust temperature and column density assuming a constant value of the emissivity spectral index and explore the effects of noise on the derived physical parameters. We find that the derived extinction values for TMC-1C are large for a starless core (80 mag Av), and the derived temperatures are low (6 K) in the densest regions of the core, using our derived value of beta = 1.8.

The AstroMed project at Harvard University's Initiative in Innovative Computing (IIC) is working on improved visualization and data sharing solutions applicable to the fields of both astronomy and medicine. The current focus is on the application of medical imaging visualization and analysis techniques to three-dimensional astronomical data. The 3D Slicer and OsiriX medical imaging tools have been used to make isosurface and volumetric models in RA-DEC-velocity space of the Perseus star forming region from the COMPLETE Survey of Star Forming Region's spectral line maps. 3D Slicer, a brain imaging and visualization computer application developed at Brigham and Women's Hospital's Surgical Planning Lab, is capable of displaying volumes (i.e. data cubes), displaying slices in any direction through the volume, generating 3D isosurface models from the volume which can be viewed and rotated in 3D space, and making 3D models of label maps (for example CLUMPFIND output). OsiriX is able to generate volumetric models from data cubes and allows the user in real time to change the displayed intensity level, crop the models without losing the data, manipulate the model and viewing angle, and use a variety of projections. In applying 3D Slicer to 12CO and 13CO spectral line data cubes of Perseus, the visualization allowed for a rapid review of over 8 square degrees and 150,000 spectra, and the cataloging of 217 high velocity points. These points were further investigated in half of Perseus and all known outflows were detected, and 20 points were identified in these regions as possibly being associated with undocumented outflows. All IIC developed tools, as well as 3D Slicer and OsiriX, are freely available.

Almost all young stars are found in multiple systems. This suggests that protostellar cores almost always fragment into multiple objects. The observed properties of multiple systems such as their separation distribution and mass ratios provide strong constraints on star formation theories. We review the observed properties of young and old multiple systems and find that the multiplicity of stars changes. Such an evolution is probably due to (a) the dynamical decay of small-N systems and/or (b) the destruction of multiple systems within dense clusters. We review simulations of the fragmentation of rotating and turbulent molecular cores. Such models almost always produce multiple systems, however the properties of those systems do not match observations at all well. Magnetic fields appear to supress fragmentation, prehaps suggesting that they are not dynamically important in the formation of multiple systems. We finish by discussing possible reasons why theory fails to match observation, and the future prospects for this field.

We present an overview of data available for the Ophiuchus and Perseus molecular clouds from ``Phase I'' of the COMPLETE Survey of Star-Forming Regions. This survey provides a range of data complementary to the Spitzer Legacy Program ``From Molecular Cores to Planet Forming Disks.'' Phase I includes: Extinction maps derived from 2MASS near-infrared data using the NICER algorithm; extinction and temperature maps derived from IRAS 60 and 100um emission; HI maps of atomic gas; 12CO and 13CO maps of molecular gas; and submillimetre continuum images of emission from dust in dense cores. Not unexpectedly, the morphology of the regions appears quite different depending on the column-density tracer which is used, with IRAS tracing mainly warmer dust and CO being biased by chemical, excitation and optical depth effects. Histograms of column-density distribution are presented, showing that extinction as derived from 2MASS/NICER gives the closest match to a log-normal distribution as is predicted by numerical simulations. All the data presented in this paper, and links to more detailed publications on their implications are publically available at the COMPLETE website.

We have used a numerical simulation of a turbulent cloud to synthesize maps of the thermal emission from dust at a variety of far-IR and sub-mm wavelengths. The average column density and external radiation field in the simulation is well matched to clouds such as Perseus and Ophiuchus. We use pairs of single-wavelength emission maps to derive the dust color temperature and column density, and we compare the derived column densities with the true column density. We demonstrate that longer wavelength emission maps yield less biased estimates of column density than maps made towards the peak of the dust emission spectrum. We compare the scatter in the derived column density with the observed scatter in Perseus and Ophiuchus. We find that while in Perseus all of the observed scatter in the emission-derived versus the extinction-derived column density can be attributed to the flawed assumption of isothermal dust along each line of sight, in Ophiuchus there is additional scatter above what can be explained by the isothermal assumption. Our results imply that variations in dust emission properties within a molecular cloud are not necessarily a major source of uncertainty in column density measurements.

The Perseus molecular cloud complex is a ~30pc long chain of molecular clouds most well-known for the two star-forming clusters NGC1333 and IC348 and the well-studied outflow source in B5. However, when studied at mid- to far-infrared wavelengths the region is dominated by a ~10pc diameter shell of warm dust, likely generated by an HII region caused by the early B-star HD278942. Using a revised calibration technique the COMPLETE team has produced high-sensitivity temperature and column-density maps of the Perseus region from IRAS Sky Survey Atlas (ISSA) 60 and 100um data. In this paper, we combine the ISSA based dust-emission maps with other observations collected as part of the COMPLETE Survey, along with archival H-alpha and MSX observations. Molecular line observations from FCRAO and extinction maps constructed by applying the NICER method to the 2MASS catalog provide independent estimates of the ``true'' column-density of the shell. H-alpha emission in the region of the shell confirms that it is most likely an HII region located behind the cloud complex, and 8um data from MSX indicates that the shell may be interacting with the cloud. Finally, the two polarisation components previously seen towards background stars in the region can be explained by the association of the stronger component with the shell. If confirmed, this would be the first observation of a parsec-scale swept-up magnetic field.

We present new deep near-infrared images of dark clouds in the Perseus molecular complex. These images show beautiful extended emission which we model as scattered ambient starlight and name ``cloudshine''. The brightness and color variation of cloudshine complicates the production of extinction maps, the best tracer of column density in clouds. However, since the profile of reflected light is essentially a function of mass distribution, cloudshine provides a new way to study the structure of dark clouds. Previous work has used optical scattered light to study the density profile of tenuous clouds; extending this technique into the infrared provides a high-resolution view into the interiors of very dense clouds, bypassing the complexities of using thermal dust emission, which is biased by grain temperature, or molecular tracers, which have complicated depletion patterns. As new wide-field infrared cameras are used to study star-forming regions at greater depth, cloudshine will be widely observed and should be seen as a new high-resolution tool, rather than an inconvenience.

We have created new dust temperature and column density maps of Perseus, Ophiuchus and Serpens using 60 and 100 micron data from the IRIS recalibration of IRAS data. We describe an optimized method for finding the dust temperature, emissivity spectral index, and optical depth using optical and NIR extinction maps. The creation of these temperature and extinction maps (covering tens of square degrees of molecular clouds) is one of the first results from the ongoing COordinated Molecular Probe Line Extinction Thermal Emission (COMPLETE) Survey of Star Forming regions. However, while the extinctions derived from the IRIS emission maps are globally accurate, we warn that FIR emission is not a good proxy for extinction on the scale of one pixel (5').

3D Slicer, a brain imaging and visualization computer application developed at Brigham and Women's Hospital's Surgical Planning Lab, was used to visualize the IC 348 star-forming complex in Perseus in RA-DEC-Velocity space. This project is part of a collaboration between Harvard University's Astronomy Department and Medical School, and serves as a proof-of-concept technology demonstration for Harvard's Institute for Innovative Computing (IIC). 3D Slicer is capable of displaying volumes (data cubes), slices in any direction through the volume, 3D models generated from the volume, and 3D models of label maps. The 3D spectral line maps of IC 348 were created using 13CO and C18O data collected at the FCRAO (Five College Radio Astronomy Observatory). The 3D visualization makes the identification of the cloud's inner dense cores and velocity structure components easier than current conventional astronomical methods. It is planned for 3D Slicer to be enhanced with astrophysics-specific features resulting in an astronomical version to be available online through the IIC and in conjunction with the National Virtual Observatory (NVO).

We have mapped the central 10'x10' of the dense core TMC-1C at 450 and 850 microns using SCUBA on the James Clerk Maxwell Telescope. The unusually high quality of the 450 micron map allows us to make a detailed analysis of the temperature and column density profiles of the core. We find that the dust temperature at the center of TMC-1C is 7 K, rising to 11 K at the edges. We discuss the possibility and effects of a variable emissivity spectral index on the derived mass profile. The low dust temperature of TMC-1C results in a high derived mass for the core, significantly larger than the virial mass estimated from the linewidth of the N2H+ (1-0) transition. This result is valid within a wide range of dust properties and ellipticities of the core. The N2H+ (1-0) spectra, taken with the IRAM 30m telescope, show signs of self-absorption, which provide evidence of sub-sonic infall motions. The derived density profile and infall velocity is compared to the predictions of several popular star formation models, and the Bonnor-Ebert model isthe best fit analytic model.

At age two, humans are just learning to speak, and they see the world in an unbiased way. At age two, the COMPLETE (COordinated Molecular Probe Line Extinction Thermal Emission) Survey of Star-Forming Regions is just beginning to hint at what it will do for our world view, and it is unbiased as well. The coordinated observations that comprise COMPLETE have already shown us: 1) that large (~10pc) scale bubbles created by winds from evolved stars dominate maps of extended dust emission even in relatively low-mass star-forming regions (e.g. Perseus and Ophiuchus); 2) that "giant" (>1 pc-scale) outflows from young stars can be easily detected by applying statistical tools to an unbiased CO survey (e.g. we have doubled the number of known outflows in Perseus using the Spectral Correlation Function); and 3) that many apparent high column-density "clouds" seen in projection are actually the superposition of many velocity features (likely at different distances) along the line of sight. The upcoming detailed analysis of these new observations will quantify the degree to which the apparently ubiquitous spherical and collimated outflows can sculpt and stir molecular clouds. Taken together with recent numerical simulations' results and measurements of high velocities for some young stars, the COMPLETE Survey also raises, and should soon answer, questions about how, when, and how far stars move away from their birthplaces as functions of time.

Three independent lines of evidence imply that the young star PV Ceph is moving at roughly 20 km/s through the interstellar medium. The first, and strongest, suggestion of motion comes from the geometry of the HH knots in the "giant" Herbig-Haro (HH) flow associated with PV Ceph. Bisectors of lines drawn between pairs of knots at nearly equal distances from PV Ceph imply an E-W motion of the source, and a plasmon model fit to the knot positions gives a good fit of 22 km/s motion for the star. The second bit of damning evidence comes from a redshifted "trail" of molecular gas, pointing in the same E-W direction implied by the HH knot geometry. The third exhibit we offer in accusing PV Ceph of speeding involves the tilt apparent in the high-velocity molecular jet now emanating from the star. This tilt is best explained if the true, current, jet direction is N-S, as it is in HST images, and the star is moving at roughly 20 km/s. Tracing the motion of PV Ceph backward in time, to the nearest cluster from which it might have been ejected, we find that it is very likely to have been thrown out of the massive star-forming cluster NGC 7023 (more than 10 pc away). We propose that PV Ceph was ejected, at a speed large enough to escape NGC 7023, at least 100,000 years ago, but that it did not enter the molecular cloud in which it now finds itself until more like 35,000 years ago. Our calculations show that the currently-observable molecular outflow associated with PV Ceph is about 10,000 years old, so that the flow has had plenty of time to form while in its current molecular cloud. But, the question of what PV Ceph was doing, and what gas/disk it took along with it in the time it was traveling through the low-density region between NGC 7023 and its current home is an open question.

Magnetic fields are believed to play an important role in the evolution of molecular clouds, from their large scale structure to dense cores, protostellar envelopes, and protoplanetary disks. How important is unclear, and whether magnetic fields are the dominant force driving star formation at any scale is also unclear. In this review we examine the observational data which address these questions, with particular emphasis on high angular resolution observations. Unfortunately the data do not clarify the situation. It is clear that the fields are important, but to what degree we don't yet know. Observations to date have been limited by the sensitivity of available telescopes and instrumentation. In the future ALMA and the SKA in particular should provide great advances in observational studies of magnetic fields, and we discuss which observations are most desirable when they become available.

We compute the spectral correlation function (SCF) of 13CO J=1--0 maps of molecular cloud complexes. The SCF is a power law over approximately an order of magnitude in spatial separation in every map. The power law slope of the SCF, alpha, its normalization, S_0(1pc), and the spectral line width averaged over the whole map, sigma_v, are computed for all the observational maps. The values of alpha, S_0(1pc) and sigma_v are combined to obtain empirical correlations to be used as tests for theoretical models of molecular clouds. Synthetic spectral maps are computed from different theoretical models, including solutions of the magneto-hydrodynamic (MHD) equations with different values of the rms Mach number of the flow and stochastic models with different power spectra of the velocity field. In order to compute the radiative transfer from the MHD models it is necessary to assign the models a physical scale and a physical density. When these assignments are made according to Larson type relations the best fit to the observational correlations is obtained. Unphysical stochastic models are instead ruled out by the empirical correlations. MHD models with equipartition of magnetic and kinetic energy of turbulence do not reproduce the observational data when their average magnetic field is oriented approximately parallel to the line of sight.

We present a new set of high-resolution molecular line maps of the gas immediately surrounding various Herbig-Haro (HH) knots of the giant HH flow HH 315, from the young star PV Cephei. The observations, aimed at studying the entrainment mechanism of the 2.6 pc-long HH 315 flow, include IRAM 30m maps of the 12CO(2-1), 12CO(1-0), and 13CO(1-0) lines, with beam sizes of 11'', 21'', and 22'', respectively. We compare the morphology and the kinematics of the outflow gas, as well as the temperature and momentum distribution of the molecular outflow with those predicted by different entrainment models. With our detailed study we are able to conclude that jet bow shock entrainment by an episodic stellar wind, with a time-varying axis, produces most of the high-velocity molecular outflow observed far from the source. In addition, near PV Cephei we find evidence for a poorly collimated, wide-angle, molecular outflow and a collimated wiggling jet-like molecular outflow. We propose that the poorly collimated component is entrained by a wide-angle wind, and the collimated component is entrained by a variable jet with internal working surfaces. If this picture is true, then a stellar wind model which allows for the coexistence of a wide-angle component and a collimated (jet-like) stellar wind component is needed to explain the observed properties of the PV Ceph outflow. The wiggling axis of the redshifted molecular outflow lobe indicates that the outflow ejection axis is changing over time. We find that the time-scale of the axis variation shown by the molecular outflow lobe is about a factor of 10 less than that shown by the large-scale optical HH knots.

We present a set of detailed molecular line maps of the region associated with the giant Herbig-Haro flow HH 315, from the young star PV Cephei, aimed at studying the outflow-cloud interaction. Our study clearly shows that the HH 315 flow is effecting the kinematics of its surrounding medium, and has been able to redistribute considerable amounts of the surrounding medium-density gas in its star-forming core as well at parsec-scale distances from the source. The giant molecular outflow HH 315 is a highly asymmetric bipolar flow with a projected linear extent of about 2 pc. Our results indicate that the two outflow lobes are each interacting with the ambient medium in different ways. The southern (redshifted) lobe interacts with a dense ambient medium, very close to the young stellar outflow source, and its kinetic energy is comparable to both the turbulent and gravitational binding energy of its host cloud. In addition, we find evidence that the southern lobe is responsible for the creation of a cavity in the 13CO emission. In contrast, the northern (mainly blueshifted) outflow lobe extends farther from PV Ceph and interacts with ambient gas much less dense than the southern lobe. There is very little 13CO emission north of the outflow source, and the only prominent 13CO emission is a shell-like structure coincident with the outer edge of the northern lobe, about 1.2 pc northwest of PV Ceph. It appears that the northern lobe of the HH 315 outflow has been able to ``push'' aside a substantial fraction of the gas in the area, piling it in a dense shell-like structure at its edges. In addition, we find that the northern outflow lobe is responsible for a velocity gradient in the ambient gas.

(Abridged) We use the statistical tool known as the ``Spectral Correlation Function" [SCF] to intercompare simulations and observations of the atomic interstellar medium. The simulations considered mimic three distinct sets of physical conditions. One of them (run "ISM") is intended to represent a mixture of cool and warm atomic gas, and includes self-gravity and magnetic fields. For each simulation, H I spectral-line maps are synthesized and intercompared, both with each other, and with observations, using the SCF. We find that, when thermal broadening is large in comparison with fine-scale turbulent velocity structure, it masks sub-thermal velocity sub-structure in the synthesized spectra. The H I observations we use here for comparison are of the North Celestial Pole (NCP) Loop. None of the simulations match the NCP Loop data very well. The most realistic sets of line profiles and SCF statistics comes from artifically rescaling the velocity axis of run ISM. Without rescaling, almost all velocity structure is smeared out by thermal broadening. However, if the velocity axis is expanded by a factor of 6, the SCF distributions of run ISM an the NCP Loop match up fairly well. This means that the ratio of thermal to turbulent pressure in run ISM is much too large as it stands, and that the simulation is deficient in turbulent energy. This is a consequence of run ISM not including the effects of supernovae. We conclude that the SCF is a useful tool for understanding and fine-tuning simulations of interstellar gas, and in particular that realistic simulations of the atomic ISM need to include the effects of energetic stellar winds (e.g. supernovae) in order for the ratio of thermal-to-turbulent pressure to give spectra representative of the observed interstellar medium in our Galaxy.

We model the polarized thermal dust emission from protostellar cores that are assembled by super-sonic turbulent flows in molecular clouds. Self-gravitating cores are selected from a three dimensional simulation of super-sonic and super-Alfvenic magneto-hydrodynamic (MHD) turbulence. The polarization is computed in two ways. In model A it is assumed that dust properties and grain alignment efficiency are uniform; in model B it is assumed that grains are not aligned at visual extinction larger than 3 mag. The main results of this work are: i) Values of the degree of polarization P between 1 and 10% are typical, despite the super-Alfvenic nature of the turbulence; ii) A steep decrease of P with increasing values of the sub-mm dust continuum intensity I is always found in self--gravitating cores selected from the MHD simulations, if grains are not aligned above a certain value of visual extinction (model B); iii) The same behavior is hard to reproduce if grains are aligned independently of visual extinction (model A); iv) The Chandrasekhar-Fermi formula, corrected by a factor f=0.4, provides an approximate estimate of the average magnetic field strength in the cores. Sub-mm dust continuum polarization maps of quiescent protostellar cores and Bok globules always show a decrease in P with increasing value of I consistent with the predictions of our model B. We therefore conclude that sub-mm polarization maps of quiescent cores do not map the magnetic field inside the cores at visual extinction larger than approximately 3 mag. There is no inconsistency between the results from optical and near-IR polarized absorption of background stars, and the observed polarization of sub-mm dust continuum from quiescent cores. In both cases, grains at large visual extinction appear to be virtually unaligned.

While observational evidence for the episodic nature of young stellar outflows continues to mount, existing numerical and theoretical models of molecular outflows assume they are formed by the interaction of a non-episodic wind from a young stellar object with an ambient cloud. In this Letter we estimate and discuss the effects of episodicity on the mass-velocity and position-velocity relations observed in molecular outflows. We explain how many recent observational results disagree with the predictions of non-episodic outflow models, and we offer simple explanations for the discrepancies. In particular, we discuss how an episodic stellar outflow can steepen the power-law slope of the mass-velocity relation in a molecular outflow. And, we illustrate how an episodic outflow can produce multiple "Hubble-wedges'' in the position-velocity distribution of a molecular outflow. With a little more information than we have now, it may be possible to use the "fossil record" embedded in a molecular outflow's mass-velocity and position-velocity relations to reconstruct the history of a young stellar object's mass ejection episodes.

We present the first set of detailed molecular line maps of the region associated with the giant Herbig-Haro flow HH 300, from the young star IRAS 04239+2436. Our results indicate that the red lobe of the HH 300 flow is depositing a fair amount of momentum and kinetic energy over a notable volume of its host dark cloud. This makes HH 300 a key player in the evolution and fate of its parent cloud. The redshifted molecular outflow lobe of HH 300 is 1.1 pc long and 0.3 pc wide, and has a very clumpy structure. The density, velocity, and momentum distributions in the outflow indicate that the observed clumps arise from the prompt entrainment of ambient gas. Bow shock-shaped structures are associated with the clumps, and we hypothesize that the shocks are produced by different mass ejection episodes. Lines drawn from IRAS 04239+2436 to each of these clumps have different orientations on the plane of the sky, and we conclude that HH 300 is a precessing and episodic outflow. The observations include a map of the red lobe in the 12CO(2-1) line, with a beam size of 27'', and more extended maps of the outflow region in the 12CO(1-0) and 13CO(1-0) lines, with 45'' and 47'' beam sizes, respectively. We obtain a steep power-law mass spectrum for HH 300, which we believe is best explained by the evolution of the outflow mass kinematics. In addition, our 13CO observations show that the HH 300 flow has been able to redistribute (in space and velocity) considerable amounts of its surrounding medium-density gas.

We propose a new method to measure and map the gas scale height of nearby disk galaxies. This method is applied successfully to the Australia Telescope Compact Array interferometric HI survey of the Large Magellanic Cloud (LMC); it could also be applied to a significant number of nearby disk galaxies, thanks to the next generation of interferometric facilities, such as the extended VLA and CARMA. The method consists of computing the Spectral Correlation Function (SCF) for a spectral-line map of a face-on galaxy. The SCF quantifies the correlation between spectra at different map positions as a function of their separation, and is sensitive to the properties of both the gas mass distribution and the gas velocity field. It is likely that spatial correlation properties of the gas density and velocity fields in a galactic disk are sensitive to the value of the scale height of the gas disk. A scale-free turbulent cascade is unlikely to extend to scales much larger than the disk scale height, as the disk dynamics on those larger scales should be dominated by two dimensional motions. We find a clear feature in the SCF of the LMC HI disk, on the scale of approximately 180 pc, which we identify as the disk scale height. We are also tentatively able to map variations of the scale height over the disk.

Super-sonic turbulence fragments molecular clouds (MC) into a very complex density field with density contrasts of several orders of magnitude. A fraction of the gas is locked into dense and gravitationally bound cores, which collapse as proto-stars. This process can be studied with numerical simulations of super-sonic self-gravitating turbulence. In this work, we use numerical simulations of magneto-hydrodynamic (MHD), super-sonic, super-Alfv\'{e}nic and self-gravitating turbulence to compute the mass distribution of collapsing proto-stellar cores, which are selected as local density maxima. We find that the mass distribution of collapsing cores is consistent with the stellar initial mass function (IMF), suggesting that super-sonic turbulence may be responsible for the generation of the IMF. To support this conclusion we also show that the physical properties of the numerically selected cores are in agreement with the properties of observed NH$_3$ cores and that their magnetic field strength is consistent with Zeeman splitting measurements. In turbulent MCs, star formation occurs via the gravitational collapse of super-critical cores, formed by the turbulent flow, sub-critical cores being irrelevant for the process of star formation.

The fragmentation of molecular clouds (MC) into protostellar cores is a central aspect of the process of star formation. Because of the turbulent nature of super-sonic motions in MCs, it has been suggested that dense structures such as filaments and clumps are formed by shocks in a turbulent flow. In this work we present strong evidence in favor of the turbulent origin of the fragmentation of MCs. The most generic result of turbulent fragmentation is that dense post shock gas traces a gas component with a smaller velocity dispersion than lower density gas, since shocks correspond to regions of converging flows, where the kinetic energy of the turbulent motion is dissipated. Using synthetic maps of spectra of molecular transitions, computed from the results of numerical simulations of super--sonic turbulence, we show that the dependence of velocity dispersion on gas density generates an observable relation between the rms velocity centroid and the integrated intensity (column density), Sigma(V_0)-I, which is indeed found in the observational data. The comparison between the theoretical model (maps of synthetic 13CO spectra), with 13CO maps from the Perseus, Rosette and Taurus MC complexes, shows excellent agreement in the Sigma(V_0)-I relation. The Sigma(V_0)-I relation of different observational maps with the same total rms velocity are remarkably similar, which is a strong indication of their origin from a very general property of the fluid equations, such as the turbulent fragmentation process.

The effects of noise and sampling on the ``Spectral Correlation Function'' (SCF) introduced by Rosolowsky et al. 1999 are studied using observational data, numerical simulations of magneto-hydrodynamic turbulence, and simple models of Gaussian spectral line profiles. The most significant innovations of this paper are: i) the normalization of the SCF based on an analytic model for the effect of noise; ii) the computation of the SCF as a function of the spatial lag between spectra within a map. By computing of the dependence of the SCF on the spatial lag, S_o(dr), we have been able to conclude that: 1) S_o(dr) is a power law, with slope in the range of scales l_i-l_o. 2) The correlation outer scale, l_o, is determined by the size of the map, and no evidence for a true departure from self-similarity on large scales has been found. 3) The correlation inner scale, l_i, is a true estimate of the smallest self-similar scale in a map. 4) The spectral slope in a given region, is independent of velocity resolution (above a minimum resolution threshold), spatial resolution, and average spectrum quality. 5) Molecular transitions which trace higher gas density yield larger values of the spectral slope than transitions tracing lower gas density. 6) Nyquist sampling, bad pixels in detector arrays, and reference sharing data acquisition need to be taken into account for a correct determination of the SCF at lag dr=1. The value of the spectral slope, however, can be computed correctly without a detailed knowledge of observational procedures.

The "spectral correlation function" analysis we introduce in this paper is a new tool for analyzing spectral-line data cubes. Our initial tests, carried out on a suite of observed and simulated data cubes, indicate that the spectral correlation function [SCF] is likely to be a more discriminating statistic than other statistical methods normally applied. The SCF is a measure of similarity between neighboring spectra in the data cube. When the SCF is used to compare a data cube consisting of spectral-line observations of the ISM with a data cube derived from MHD simulations of molecular clouds, it can find differences that are not found by other analyses. The initial results presented here suggest that the inclusion of self-gravity in numerical simulations is critical for reproducing the correlation behavior of spectra in star-forming molecular clouds.

We test the recently published all-sky reddening map of Schlegel, Finkbeiner & Davis (1998 [SFD]) using the extinction study of a region in the Taurus dark cloud complex by Arce & Goodman (1999 [AG]). In their study, AG use four different techniques to measure the amount and structure of the extinction toward Taurus, and all four techniques agree very well. Thus we believe that the AG results are a truthful representation of the extinction in the region and can be used to test the reliability of the SFD reddening map. The results of our test show that the SFD all-sky reddening map, which is based on data from COBE/DIRBE and IRAS/ISSA, overestimates the reddening by a factor of 1.3 to 1.5 in regions of smooth extinction with A_V > 0.5 mag. In some regions of steep extinction gradients the SFD map underestimates the reddening value, probably due to its low spatial resolution. We expect that the astronomical community will be using the SFD reddening map extensively. We offer this Letter as a cautionary note about using the SFD map in regions of high extinction (A_V > 0.5 mag), as it might not be giving accurate reddening values there.

We present a study of the detailed distribution of extinction in a region of the Taurus dark cloud complex. Our study uses new BVR images of the region, spectral classification data for 95 stars, and IRAS Sky Survey Atlas (ISSA) 60 and 100 micron images. We study the extinction of the region in four different ways, and we present the first inter-comparison of all these methods, which are: 1) using the color excess of background stars for which spectral types are known; 2) using the ISSA 60 and 100 micron images; 3) using star counts; and 4) using an optical (V and R) version of the average color excess method used by Lada et al. (1994). We find that all four methods give generally similar results, with important exceptions. To study the structure in the dust distribution, we compare the ISSA extinction and the extinction measured for individual stars. From the comparison, we conclude that in the relatively low extinction regions studied, with 0.9 < A_V < 3.0 mag (away from filamentary dark clouds and IRAS cores), there are no fluctuations in the dust column density greater than 45% (at the 99.7% confidence level), on scales smaller than 0.2 pc. We also report the discovery of a previously unknown stellar cluster behind the Taurus dark cloud near R.A 4h19m00s, Dec. 27:30:00 (B1950)

We present a study of the polarizing power of the dust in cold dense regions (dark clouds) compared to that of dust in the general interstellar medium (ISM). Our study uses new polarimetric, optical, and spectral classification data for 36 stars to carefully study the relation between polarization percentage (p) and extinction (A_V) in the Taurus dark cloud complex. We find two trends in our p-A_V study: (1) stars background to the warm ISM show an increase in p with A_V; and (2) the percentage of polarization of stars background to cold dark clouds does not increase with extinction. We detect a break in the p-A_V relation at an extinction 1.3 +/- 0.2 mag, which we expect corresponds to a set of conditions where the polarizing power of the dust associated with the Taurus dark clouds drops precipitously. This breakpoint places important restrictions on the use of polarimetry in studying interstellar magnetic fields.

A quantitative analysis of grain alignment in the filamentary dark cloud L1755 in Ophiuchus is presented. We show, that the observed decrease of the polarization to extinction ratio for the inner parts of this quiescent dark cloud can be explained as a result of the decrease of the efficiency of grain alignment.We make quantitative estimates of grain alignment efficiency for six mechanisms, involving grains with either thermal or suprathermal rotation, interacting with either magnetic field or gaseous flow. We also make semiquantitative estimates of grain alignment by radiative torques.We show that in conditions typical of dark cloud interiors, all known major mechanisms of grain alignment fail. All the studied mechanisms predict polarization at least an order of magnitude below the currently detectable levels of order of 1%. On the contrary, in the dark cloud environments where A_v<1 the grain alignment can be much more efficient. There the alignment of suprathermally rotating grains with superparamagentic inclusions, and possibly also radiative torques account for observed polarization. These results apply to L1755, which we model in detail and probably also to B216 and other similar dark clouds. Our study suggests an explanation for the difference in results obtained through polarimetry of background starlight and polarized thermal emission from the dust itself. We conjecture that the emission polarimetry selectively reveals aligned grains in the environment far from thermodynamic equilibrium as opposed to starlight polarization studie that probe the alignment of grains all the way along the line of sight, including the interiors of dark quiescent clouds, where no alignment is possible.