Source author record

Alyssa A. Goodman

Alyssa A. Goodman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

35works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

35 published item(s)

preprint2022arXiv

An Interferometric View of H-MM1. I. Direct Observation of NH3 Depletion

Spectral lines of ammonia, NH$_3$, are useful probes of the physical conditions in dense molecular cloud cores. In addition to advantages in spectroscopy, ammonia has also been suggested to be resistant to freezing onto grain surfaces, which should make it a superior tool for studying the interior parts of cold, dense cores. Here we present high-resolution NH$_3$ observations with the Very Large Array (VLA) and Green Bank Telescope (GBT) towards a prestellar core. These observations show an outer region with a fractional NH$_3$ abundance of X(NH$_3$) = (1.975$\pm$0.005)$\times 10^{-8}$ ($\pm 10\%$ systematic), but it also reveals that after all, the X(NH$_3$) starts to decrease above a H$_2$ column density of $\approx 2.6 \times 10^{22}$ cm$^{-2}$. We derive a density model for the core and find that the break-point in the fractional abundance occurs at the density n(H$_2$) $\sim 2\times10^5$ cm$^{-3}$, and beyond this point the fractional abundance decreases with increasing density, following the power law $n^{-1.1}$. This power-law behavior is well reproduced by chemical models where adsorption onto grains dominates the removal of ammonia and related species from the gas at high densities. We suggest that the break-point density changes from core to core depending on the temperature and the grain properties, but that the depletion power law is anyway likely to be close to $n^{-1}$ owing to the dominance of accretion in the central parts of starless cores.

preprint2022arXiv

Building the molecular cloud population: the role of cloud mergers

We study the physical drivers of slow molecular cloud mergers within a simulation of a Milky Way-like galaxy in the moving-mesh code Arepo, and determine the influence of these mergers on the mass distribution and star formation efficiency of the galactic cloud population. We find that 83 per cent of these mergers occur at a relative velocity below 5 km/s, and are associated with large-scale atomic gas flows, driven primarily by (1) expanding bubbles of hot, ionised gas caused by supernova explosions and (2) galactic rotation. The major effect of these mergers is to aggregate molecular mass into higher-mass clouds: mergers account for over 50 per cent of the molecular mass contained in clouds of mass M > 2 x 10^6 Msun. These high-mass clouds have higher densities, internal velocity dispersions and instantaneous star formation efficiencies than their unmerged, lower-mass precursors. As such, the mean instantaneous star formation efficiency in our simulated galaxy, with its merger rate of just 1 per cent of clouds per Myr, is 25 per cent higher than in a similar population of clouds containing no mergers.

preprint2022arXiv

Star formation near the Sun is driven by expansion of the Local Bubble

For decades we have known that the Sun lies within the Local Bubble, a cavity of low-density, high-temperature plasma surrounded by a shell of cold, neutral gas and dust. However, the precise shape and extent of this shell, the impetus and timescale for its formation, and its relationship to nearby star formation have remained uncertain, largely due to low-resolution models of the local interstellar medium. Leveraging new spatial and dynamical constraints from the Gaia space mission, here we report an analysis of the 3D positions, shapes, and motions of dense gas and young stars within 200 pc of the Sun. We find that nearly all the star-forming complexes in the solar vicinity lie on the surface of the Local Bubble and that their young stars show outward expansion mainly perpendicular to the bubble's surface. Tracebacks of these young stars' motions support a scenario where the origin of the Local Bubble was a burst of stellar birth and then death (supernovae) taking place near the bubble's center beginning 14 Myr ago. The expansion of the Local Bubble created by the supernovae swept up the ambient interstellar medium into an extended shell that has now fragmented and collapsed into the most prominent nearby molecular clouds, in turn providing robust observational support for the theory of supernova-driven star formation.

preprint2021arXiv

Dissecting the super-critical filaments embedded in the 0.5 pc subsonic region of Barnard 5

We characterize in detail the two ~0.3 pc long filamentary structures found within the subsonic region of Barnard 5. We use combined GBT and VLA observations of the molecular lines NH$_3$(1,1) and (2,2) at a resolution of 1800 au, as well as JCMT continuum observations at 850 and 450 $μ$m at a resolution of 4400 au and 3000 au, respectively. We find that both filaments are highly super-critical with a mean mass per unit length, $M/L$, of ~80 M$_\odot$ pc$^{-1}$, after background subtraction, with local increases reaching values of ~150 M$_\odot$ pc$^{-1}$. This would require a magnetic field strength of ~500 $μ$G to be stable against radial collapse. We extract equidistant cuts perpendicular to the spine of the filament and fit a modified Plummer profile as well as a Gaussian to each of the cuts. The filament widths (deconvolved FWHM) range between 6500-7000 au (~0.03 pc) along the filaments. This equals ~2.0 times the radius of the flat inner region. We find an anti-correlation between the central density and this flattening radius, suggestive of contraction. Further, we also find a strong correlation between the power-law exponent at large radii and the flattening radius. We note that the measurements of these three parameters fall in a plane and derive their empirical relation. Our high-resolution observations provide direct constraints of the distribution of the dense gas within super-critical filaments showing pre- and protostellar activity.

preprint2020arXiv

A Galactic-scale gas wave in the Solar Neighborhood

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.

preprint2020arXiv

Efficient methanol production on the dark side of a prestellar core

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.

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2019arXiv

Droplets II: Internal Velocity Structures and Potential Rotational Motions in Pressure-dominated Coherent Structures

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, $β$, of $\sim 0.046$ for the droplets, and when including both droplets and larger-scale dense cores, we find $β\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.

preprint2019arXiv

Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) -- Full Data Release

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 $μ$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/.

preprint2016arXiv

A Uniform Catalog of Molecular Clouds in the Milky Way

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 $γ=-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 $γ=-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.

preprint2016arXiv

Misalignment of Outflow Axes in the Proto-Multiple Systems in Perseus

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.

preprint2015arXiv

Mass Assembly of Stellar Systems and Their Evolution with the SMA (MASSES). Multiplicity and the Physical Environment in L1448N

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.

preprint2014arXiv

The Bones of the Milky Way

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.

preprint2013arXiv

A New Approach to Developing Interactive Software Modules through Graduate Education

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

preprint2013arXiv

Quantifying Observational Projection Effects Using Molecular Cloud Simulations

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.

preprint2012arXiv

Dust SEDs in the era of Herschel and Planck: a Hierarchical Bayesian fitting technique

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.

preprint2012arXiv

Evidence for Grain Growth in Molecular Clouds: A Bayesian Examination of the Extinction Law in Perseus

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.

preprint2012arXiv

Principles of High-Dimensional Data Visualization in Astronomy

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.

preprint2011arXiv

A Bubbling Nearby Molecular Cloud: COMPLETE Shells in Perseus

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.

preprint2011arXiv

A Guide to Comparisons of Star Formation Simulations with Observations

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.

preprint2011arXiv

Classifying structures in the ISM with Support Vector Machines: the G16.05-0.57 supernova remnant

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.

preprint2011arXiv

EVLA Observations of the Barnard 5 Star-Forming Core: Embedded Filaments Revealed

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

preprint2011arXiv

Observing Turbulent Fragmentation in Simulations: Predictions for CARMA and ALMA

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.

preprint2011arXiv

The enigmatic core L1451-mm: a first hydrostatic core? or a hidden VeLLO?

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.

preprint2010arXiv

Direct observation of a sharp transition to coherence in Dense Cores

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.

preprint2010arXiv

Star formation in the Taurus filament L1495: From Dense Cores to Stars

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.

preprint2010arXiv

The Angular Momentum of Magnetized Molecular Cloud Cores: A 2D-3D Comparison

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.

preprint2010arXiv

The COMPLETE Survey of Outflows in Perseus

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.

preprint2010arXiv

The Effect of Projection on Derived Mass-Size and Linewidth-Size Relationships

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.

preprint2010arXiv

The Mass-Size Relation from Clouds to Cores. II. Solar Neighborhood Clouds

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

preprint2009arXiv

Dense Cores in Perseus: The Influence of Stellar Content and Cluster Environment

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.

preprint2009arXiv

The Effect of Line of Sight Temperature Variation and Noise on Dust Continuum Observations

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.

preprint2009arXiv

The Effect of Noise on the Dust Temperature - Spectral Index Correlation

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.

preprint2008arXiv

The "True" Column Density Distribution in Star-Forming Molecular Clouds

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]