2019

Jorma Harju, Jaime E. Pineda, Anton I. Vasyunin, Paola Caselli, Stella S. R. Offner, Alyssa A. Goodman, and alia. 2019. “Efficient methanol desorption in shear instability.” arXiv, 1903, 11298. Publisher's VersionAbstract
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.
Hope How-Huan Chen, Jaime E. Pineda, Stella S. R. Offner, Alyssa A. Goodman, and alia. 2019. “Droplets II: Internal Velocity Structures and Potential Rotational Motions in Coherent Cores.” arXiv, 1908, 04367. Publisher's VersionAbstract
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 ∼0.046 for the droplets, and when including both droplets and larger-scale dense cores, we find β∼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.
Michelle Ntampaka, Camille Avestruz, Steven Boada, Joao Caldeira, Jessi Cisewski-Kehe, Rosanne Di Stefano, Cora Dvorkin, August E. Evrard, Arya Farahi, Doug Finkbeiner, Shy Genel, Alyssa Goodman, and alia. 2019. “The Role of Machine Learning in the Next Decade of Cosmology.” Bulletin of the American Astronomical Society, 51, 14. Publisher's VersionAbstract
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.
C. Zucker, J. S. Speagle, E. Schlafly, G. M. Green, D. P. Finkbeiner, A. Goodman, and J. Alves. 2019. “VizieR Online Data Catalog: Distances to molecular clouds in SFR (Zucker+, 2020).” VizieR Online Data Catalog, J/A+A/633/A51.Abstract

Distances to ~60 star-forming regions in Reipurth (2008, Star Formation Handbook, vols I and II) have been computed using stellar photometry and Gaia DR2 parallax measurements. Usually, several distance estimates are taken across each cloud.

 

For each sightline, the median distance (d50) is provided, plus the 16th and 84th percentiles on the distance probability distribution function. There is an additional systematic uncertainty, which is unknown but estimated to be ~5% in distance for clouds <1.5kpc, ~10% in distance for clouds >1.5kpc, and ~7% in distance for the southern clouds Lupus, Chamaeleon, and Corona Australis. These should be added in quadrature with the statistical uncertainties reported in the table. In addition to the distances, ancillary model parameters used in our fit are also included (e.g. the amount of foreground extinction "f"). See Section 3.2.1 and Section 3.2.2 in Zucker et al. (2019ApJ...879..125Z) for a complete description of model parameters.
H. Houghton, P. Udomprasert, S. Sunbury, E. Wright, A. Goodman, E. Johnson, and A. Bishop. 2019. “Cultivating Curiosity with Life in the Universe and WorldWide Telescope.” Astronomical Society of the Pacific, 524, 273.Abstract
When students encounter complex topics like the search for extraterrestrial life, questions abound - thoughtful, unpredictable, and often profound. Despite this thriving curiosity, the first step to be able to explore complex questions is developing the capacity to verbalize a meaningful question. The WorldWide Telescope Ambassadors team designed an out-of-school curriculum called Life in the Universe, which engages middle school-aged students in the science and scientific process of the search for distant life. Students practice generating meaningful questions, which will guide them through the science content, as groups of students build to culminating capstone projects. Results from surveys administered to participating students indicate gains in curiosity in science, as well as in seeing oneself as successful in science.
Catherine Zucker, Rowan Smith, and Alyssa Goodman. 12/2019. “Synthetic Large-Scale Galactic Filaments: On Their Formation, Physical Properties, and Resemblance to Observations.” The Astrophysical Journal , 887, Pp. 186. Publisher's VersionAbstract
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 ≳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 kpc2 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 ≈ 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.
P. Udomprasert, H. Houghton, S. Sunbury, J. Plummer, E. Wright, A. Goodman, E. Johnson, H. Zhang, A. Vaishampayan, and K. Cho. 11/2019. “Visualizing Seasons and Moon Phases with WorldWide Telescope.” Advancing Astronomy for All: ASP 2018 ASP Conference Series, 524, Pp. 125. Publisher's VersionAbstract

WorldWide Telescope (WWT) is a powerful visualization program that allows users to connect Earth-based and space-based views of the Sun- Earth-Moon system. By blending hands-on physical activities with WWT's virtual models, students can visualize spatially complex concepts like seasons, Moon phases, and eclipses. In this workshop, we will demonstrate how WWT and the physical models are used together in our WWT ThinkSpace curriculum, developed with funding from the National Science Foundation. We will also present student learning outcomes based on written assessments and student interviews.

Ian W. Stephens, Tyler L. Bourke, Michael M. Dunham, Philip C. Myers, Riwaj Pokhrel, John J. Tobin, Héctor G. Arce, and et al. 12/2019. “Mass Assembly of Stellar Systems and Their Evolution with the SMA (MASSES)—Full Data Release.” The Astrophysical Journal. Supplement Series, 245, Pp. 21. Publisher's VersionAbstract
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 (̃4″ or ̃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), 13CO(2-1), C18O(2-1), and N2D+(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 C18O(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 C18O(2-1) line widths are 1.45 km s-1, which is higher than the C18O line widths detected toward Perseus filaments and cores. We find that N2D+(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 C18O(2-1) and N2D+(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://da taverse.harvard.edu/dataverse/MASSES.
J.D. Soler, H. Beuther, M. Rugel, Y. Wang, P. C. Clark, S. C. O. Glover, P. F. Goldsmith, and et al. 2/2019. “Histogram of Oriented Gradients: A Technique for the Study of Molecular Cloud Formation.” Astronomy and Astrophysics, 622, Pp. A166. Publisher's VersionAbstract

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 spectral line 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 characterized HOG using synthetic observations of HI and 13CO (J = 1 → 0) emission from numerical simulations of magnetohydrodynamic (MHD) turbulence leading to the formation of molecular gas after the collision of two atomic clouds. We found a significant spatial correlation between the two tracers in velocity channels where vHI ≈ v13CO, almost independent of the orientation of the collision with respect to the line of sight. Subsequently, we used HOG to investigate the spatial correlation of the HI, from The HI/OH/recombination line survey of the inner Milky Way (THOR), and the 13CO (J = 1 → 0) emission from the Galactic Ring Survey (GRS), toward the portion of the Galactic plane 33°.75 ≤l ≤ 35°.25 and |b| ≤ 1°.25. We found a significant spatial correlation between the two tracers in extended portions of the studied region. Although some of the regions with high spatial correlation are associated with HI self-absorption (HISA) features, suggesting that it is produced by the cold atomic gas, the correlation is not exclusive to this kind of region. The HOG results derived for the observational data indicate significant differences between individual regions: some show spatial correlation in channels around vHI ≈ v13CO while others present spatial correlations in velocity channels separated by a few kilometers per second. 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.

Aneta Siemiginowska, Gwendolyn Eadie, Ian Czekala, Eric Feigelson, Eric B. Ford, Vinay Kashyap, Michael Kuhn, and et al. 5/2019. “The Next Decade of Astroinformatics and Astrostatistics.” Bulletin of the American Astronomical Society, 51, Pp. 355. Publisher's VersionAbstract

Over the past century, major advances in astronomy and astrophysics have been driven by improvements in instrumentation. With the amassing of high quality data from new telescopes it is becoming clear that research in astrostatistics and astroinformatics will be necessary to develop new methodology needed in astronomy.

Laura M. Fissel, Peter A. R. Ade, Francesco E. Angilè, Peter Ashton, Steven J. Benton, Che-Yu Che, Maria Cunningham, and et al. 6/2019. “Relative Alignment between the Magnetic Field and Molecular Gas Structure in the Vela C Giant Molecular Cloud Using Low- and High-Density Tracers.” The Astrophysical Journal , 878, Pp. 110. Publisher's VersionAbstract

We compare the magnetic field orientation for the young giant molecular cloud Vela C inferred from 500 μ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 12CO and 13CO J → 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 occur between the densities traced by 13CO and by C18O J → 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 103 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.

Hope How-Huan Chen, Jaime E. Pineda, Stella S. R. Offner, Alyssa A. Goodman, Andreas Burkert, Rachel K. Friesen, Erik Rosolowsky, Samantha Scibelli, and Yancy Shirley. 12/2019. “Droplets. II. Internal Velocity Structures and Potential Rotational Motions in Pressure-Dominated Coherent Structures.” The Astrophysical Journal, 886, Pp. 119. Publisher's VersionAbstract
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 ∼0.046 for the droplets, and when including both droplets and larger-scale dense cores, we find β ∼ 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 introduced by Soler et al. 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.
Hope How-Huan Chen, Jaime E. Pineda, Alyssa A. Goodman, Andreas Burkert, Stella S. R. Offner, Rachel K. Friesen, Philip C. Myers, and et al. 6/2019. “Droplets. I. Pressure-Dominated Coherent Structures in L1688 and B18.” The Astrophysical Journal, 877, Pp. 93. Publisher's VersionAbstract
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, we identify regions of high density and near-constant, almost-thermal velocity dispersion. We reveal 18 coherent structures are revealed, 12 in L1688 and 6 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 statistically study the coherent structures in molecular clouds. The identified coherent structures have a typical radius of 0.04 pc and a typical mass of 0.4 M ☉, generally smaller than previously known coherent cores identified by Goodman et al., Caselli et al., and Pineda et al. 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 (V LSR) 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.

Pages