Showing 1–50 of 80 results for author: Kainulainen, J

Deep learning denoising by dimension reduction: Application to the ORION-B line cubes

Authors: Lucas Einig, Jérôme Pety, Antoine Roueff, Paul Vandame, Jocelyn Chanussot, Maryvonne Gerin, Jan H. Orkisz, Pierre Palud, Miriam Garcia Santa-Maria, Victor de Souza Magalhaes, Ivana Bešlić, Sébastien Bardeau, Emeric E. Bron, Pierre Chainais, Javier R Goicoechea, Pierre Gratier, Viviana Guzman Veloso, Annie Hughes, Jouni Kainulainen, David Languignon, Rosine Lallement, François Levrier, Dariuscz C. Lis, Harvey Liszt, Jacques Le Bourlot , et al. (7 additional authors not shown)

Abstract: Context. The availability of large bandwidth receivers for millimeter radio telescopes allows the acquisition of position-position-frequency data cubes over a wide field of view and a broad frequency coverage. These cubes contain much information on the physical, chemical, and kinematical properties of the emitting gas. However, their large size coupled with inhomogenous signal-to-noise ratio (SNR… ▽ More Context. The availability of large bandwidth receivers for millimeter radio telescopes allows the acquisition of position-position-frequency data cubes over a wide field of view and a broad frequency coverage. These cubes contain much information on the physical, chemical, and kinematical properties of the emitting gas. However, their large size coupled with inhomogenous signal-to-noise ratio (SNR) are major challenges for consistent analysis and interpretation.Aims. We search for a denoising method of the low SNR regions of the studied data cubes that would allow to recover the low SNR emission without distorting the signals with high SNR.Methods. We perform an in-depth data analysis of the 13 CO and C 17 O (1 -- 0) data cubes obtained as part of the ORION-B large program performed at the IRAM 30m telescope. We analyse the statistical properties of the noise and the evolution of the correlation of the signal in a given frequency channel with that of the adjacent channels. This allows us to propose significant improvements of typical autoassociative neural networks, often used to denoise hyperspectral Earth remote sensing data. Applying this method to the 13 CO (1 -- 0) cube, we compare the denoised data with those derived with the multiple Gaussian fitting algorithm ROHSA, considered as the state of the art procedure for data line cubes.Results. The nature of astronomical spectral data cubes is distinct from that of the hyperspectral data usually studied in the Earth remote sensing literature because the observed intensities become statistically independent beyond a short channel separation. This lack of redundancy in data has led us to adapt the method, notably by taking into account the sparsity of the signal along the spectral axis. The application of the proposed algorithm leads to an increase of the SNR in voxels with weak signal, while preserving the spectral shape of the data in high SNR voxels.Conclusions. The proposed algorithm that combines a detailed analysis of the noise statistics with an innovative autoencoder architecture is a promising path to denoise radio-astronomy line data cubes. In the future, exploring whether a better use of the spatial correlations of the noise may further improve the denoising performances seems a promising avenue. In addition, △ Less

Submitted 24 July, 2023; originally announced July 2023.

Journal ref: A&A 677, A158 (2023)

Gas kinematics around filamentary structures in the Orion B cloud

Authors: Mathilde Gaudel, Jan H. Orkisz, Maryvonne Gerin, Jérôme Pety, Antoine Roueff, Antoine Marchal, François Levrier, Marc-Antoine Miville-Deschênes, Javier R. Goicoechea, Evelyne Roueff, Franck Le Petit, Victor de Souza Magalhaes, Pierre Palud, Miriam G. Santa-Maria, Maxime Vono, Sébastien Bardeau, Emeric Bron, Pierre Chainais, Jocelyn Chanussot, Pierre Gratier, Viviana Guzman, Annie Hughes, Jouni Kainulainen, David Languignon, Jacques Le Bourlot , et al. (5 additional authors not shown)

Abstract: Understanding the initial properties of star-forming material and how they affect the star formation process is key. From an observational point of view, the feedback from young high-mass stars on future star formation properties is still poorly constrained. In the framework of the IRAM 30m ORION-B large program, we obtained observations of the translucent and moderately dense gas, which we used t… ▽ More Understanding the initial properties of star-forming material and how they affect the star formation process is key. From an observational point of view, the feedback from young high-mass stars on future star formation properties is still poorly constrained. In the framework of the IRAM 30m ORION-B large program, we obtained observations of the translucent and moderately dense gas, which we used to analyze the kinematics over a field of 5 deg^2 around the filamentary structures. We used the ROHSA algorithm to decompose and de-noise the C18O(1-0) and 13CO(1-0) signals by taking the spatial coherence of the emission into account. We produced gas column density and mean velocity maps to estimate the relative orientation of their spatial gradients. We identified three cloud velocity layers at different systemic velocities and extracted the filaments in each velocity layer. The filaments are preferentially located in regions of low centroid velocity gradients. By comparing the relative orientation between the column density and velocity gradients of each layer from the ORION-B observations and synthetic observations from 3D kinematic toy models, we distinguish two types of behavior in the dynamics around filaments: (i) radial flows perpendicular to the filament axis that can be either inflows (increasing the filament mass) or outflows and (ii) longitudinal flows along the filament axis. The former case is seen in the Orion B data, while the latter is not identified. We have also identified asymmetrical flow patterns, usually associated with filaments located at the edge of an HII region. This is the first observational study to highlight feedback from HII regions on filament formation and, thus, on star formation in the Orion B cloud. This simple statistical method can be used for any molecular cloud to obtain coherent information on the kinematics. △ Less

Submitted 25 November, 2022; originally announced November 2022.

Comments: 45 pages. Abridged abstract. Accepted by Astronomy and Astrophysics

Journal ref: A&A 670, A59 (2023)

Dust extinction map of the Galactic plane based on the VVV survey data

Authors: Miaomiao Zhang, Jouni Kainulainen

Abstract: Dust extinction is one of the most reliable tracers of the gas distribution in the Milky Way. The near-infrared (NIR) Vista Variables in the Via Lactea (VVV) survey enables extinction mapping based on stellar photometry over a large area in the Galactic plane. We devise a novel extinction mapping approach, XPNICER, by bringing together VVV photometric catalogs, stellar parameter data from StarHors… ▽ More Dust extinction is one of the most reliable tracers of the gas distribution in the Milky Way. The near-infrared (NIR) Vista Variables in the Via Lactea (VVV) survey enables extinction mapping based on stellar photometry over a large area in the Galactic plane. We devise a novel extinction mapping approach, XPNICER, by bringing together VVV photometric catalogs, stellar parameter data from StarHorse catalogs, and previously published Xpercentile and PNICER extinction mapping techniques. We apply the approach to the VVV survey area, resulting in an extinction map that covers the Galactic disk between 295 and 350 degrees at longitude and -2 to 2 degrees at latitude, and the Galactic bulge between -10 and 5 degrees at latitude. The map has 30 arcseconds spatial resolution and it traces extinctions typically up to about 10-20 mag of visual extinction and maximally up to Av~30 mag. We compare our map to previous dust based maps, concluding that it provides a high-fidelity extinction-based map, especially in its ability to recover both the diffuse dust component of the Galaxy and moderately extincted giant molecular cloud regions. The map is especially useful as independent, extinction-based data on the Galactic dust distribution and applicable for a wide range of studies from individual molecular clouds to the studies of the Galactic stellar populations. △ Less

Submitted 18 October, 2022; originally announced October 2022.

Comments: accepted by MNRAS, 37 pages

Bird's eye view of molecular clouds in the Milky Way: II. Cloud kinematics from sub-pc to kpc scales

Authors: Andri Spilker, Jouni Kainulainen, Jan Orkisz

Abstract: The kinematics of molecular gas are crucial for setting the stage for star formation. One key question related to the kinematic properties of gas is how they depend on the spatial scale. We aim to describe the CO spectra, velocity dispersions, and especially the linewidth-size relation, of molecular gas from cloud (parsec-) scales to kiloparsec scales in a complete region within the Milky Way disk… ▽ More The kinematics of molecular gas are crucial for setting the stage for star formation. One key question related to the kinematic properties of gas is how they depend on the spatial scale. We aim to describe the CO spectra, velocity dispersions, and especially the linewidth-size relation, of molecular gas from cloud (parsec-) scales to kiloparsec scales in a complete region within the Milky Way disk. We utilise the census of molecular clouds within 2 kpc from our earlier work, together with CO emission data for them from the literature. We study the kinematics and the Larson's relations for the sample of individual clouds. We also mimic a face-on view of the Milky Way and analyse the kinematics of the clouds within apertures of 0.25-2 kpc in size. In this way, we describe the scale-dependency of the CO gas kinematics and Larson's relations. We describe the spectra of CO gas at cloud scales and in apertures between 0.25-2 kpc in our survey area. The spectra within the apertures are relatively symmetric but show non-Gaussian high-velocity wings. At cloud-scales, our sample shows a linewidth-size relation σ_v=1.5*R^{0.3\pm0.1} with a large scatter. The mass-size relation in the sample of clouds is M_{CO}= 794*R^{1.5\pm0.5}. The relations are also present for the apertures at kpc-scales. A suggestive dependency on galactic environment is seen, with apertures closer to the Galactic centre and the Sagittarius spiral arm having slightly higher velocity dispersions. We explore the possible effect of a diffuse component in the survey area, and find that such a component would widen the CO spectra and could flatten the linewidth-size relation. Understanding the nature of the possible diffuse CO component and its effects on observations is crucial for connecting Galactic and extragalactic data. △ Less

Submitted 13 September, 2022; originally announced September 2022.

Comments: 8 pages, 10 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 667, A110 (2022)

3D shape explains star formation mystery of California and Orion A

Authors: Sara Rezaei Kh., Jouni Kainulainen

Abstract: The new Gaia data release (EDR3) with improved astrometry has opened a new era in studying our Milky Way in fine detail. We use Gaia EDR3 astrometry together with 2MASS and WISE photometry to study two of the most massive molecular clouds in the solar vicinity: Orion A and California. Despite having remarkable similarities in the plane of the sky in terms of shape, size, and extinction, California… ▽ More The new Gaia data release (EDR3) with improved astrometry has opened a new era in studying our Milky Way in fine detail. We use Gaia EDR3 astrometry together with 2MASS and WISE photometry to study two of the most massive molecular clouds in the solar vicinity: Orion A and California. Despite having remarkable similarities in the plane of the sky in terms of shape, size, and extinction, California has an order of magnitude lower star formation efficiency. We use our state-of-the-art dust mapping technique to derive the detailed three-dimensional (3D) structure of the two clouds, taking into account both distance and extinction uncertainties, and a full 3D spatial correlation between neighbouring points. We discover that, despite the apparent filamentary structure in the plane of the sky, California is a flat 120-pc-long sheet extending from 410 to 530 $pc$. We show that not only Orion A and California differ substantially in their 3D shapes, but also Orion A has considerably higher density substructures in 3D than California. This result presents a compelling reason why the two clouds have different star formation activities. We also demonstrate how the viewing angle of California can substantially change the cloud's position in the Kennicutt-Schmidt relation. This underlines the importance of 3D information in interpreting star formation relations and challenges studies that rely solely on the column density thresholds to determine star formation activities in molecular clouds. Finally, we provide accurate distance estimates to multiple lines of sight towards various parts of the two clouds. △ Less

Submitted 22 April, 2022; originally announced April 2022.

Comments: Accepted for publication in ApJ Letters

Initial Conditions for Star Formation: A Physical Description of the Filamentary ISM

Authors: A. Hacar, S Clark, F. Heitsch, J. Kainulainen, G. Panopoulou, D. Seifried, R. Smith

Abstract: The interstellar medium contains filamentary structure over a wide range of scales. Understanding the role of this structure, both as a conduit of gas across the scales and a diagnostic tool of local physics, is a major focus of star formation studies. We review recent progress in studying filamentary structure in the ISM, interpreting its properties in terms of physical processes, and exploring f… ▽ More The interstellar medium contains filamentary structure over a wide range of scales. Understanding the role of this structure, both as a conduit of gas across the scales and a diagnostic tool of local physics, is a major focus of star formation studies. We review recent progress in studying filamentary structure in the ISM, interpreting its properties in terms of physical processes, and exploring formation and evolution scenarios. We include structures from galactic-scale filaments to tenth-of-a-parsec scale filaments, comprising both molecular and atomic structures, from both observational and theoretical perspectives. In addition to the literature overview, we assemble a large amount of catalogue data from different surveys and provide the most comprehensive census of filamentary structures to date. Our census consists of 22 803 filamentary structures, facilitating a holistic perspective and new insights. We use our census to conduct a meta-analysis, leading to a description of filament properties over four orders of magnitudes in length and eight in mass. Our analysis emphasises the hierarchical and dynamical nature of filamentary structures. Filaments do not live in isolation, nor they generally resemble static structures close to equilibrium. We propose that accretion during filament formation and evolution sets some of the key scaling properties of filaments. This highlights the role of accretion during filament formation and evolution and also in setting the initial conditions for star formation. Overall, the study of filamentary structures during the past decade has been observationally driven. While great progress has been made on measuring the basic properties of filaments, our understanding of their formation and evolution is clearly lacking. In this context, we identify a number of directions and questions we consider most pressing for the field. △ Less

Submitted 17 March, 2022; originally announced March 2022.

Comments: 36 pages, 10 figures. To appear in Protostars and Planets VII, Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura

The width of Herschel filaments varies with distance

Authors: G. V. Panopoulou, S. E. Clark, A. Hacar, F. Heitsch, J. Kainulainen, E. Ntormousi, D. Seifried, R. J. Smith

Abstract: Context: Filamentary structures in nearby molecular clouds have been found to exhibit a characteristic width of 0.1 pc, as observed in dust emission. Understanding the origin of this universal width has become a topic of central importance in the study of molecular cloud structure and the early stages of star formation. Aims: We investigate how the recovered widths of filaments depend on the dista… ▽ More Context: Filamentary structures in nearby molecular clouds have been found to exhibit a characteristic width of 0.1 pc, as observed in dust emission. Understanding the origin of this universal width has become a topic of central importance in the study of molecular cloud structure and the early stages of star formation. Aims: We investigate how the recovered widths of filaments depend on the distance from the observer, by using previously published results from the Herschel Gould Belt Survey. Methods: We obtained updated estimates on the distances to nearby molecular clouds observed with Herschel by using recent results based on 3D dust extinction mapping and Gaia. We examined the widths of filaments from individual clouds separately, as opposed to treating them as a single population. We used these per-cloud filament widths to search for signs of variation amongst the clouds of the previously published study. Results: We find a significant dependence of the mean per-cloud filament width with distance. The distribution of mean filament widths for nearby clouds is incompatible with that of farther away clouds. The mean per-cloud widths scale with distance approximately as 4--5 times the beam size. We examine the effects of resolution by performing a convergence study of a filament profile in the Herschel image of the Taurus Molecular Cloud. We find that resolution can severely affect the shapes of radial profiles over the observed range of distances. Conclusions: We conclude that the data are inconsistent with 0.1 pc being the universal characteristic width of filaments. △ Less

Submitted 8 August, 2022; v1 submitted 15 November, 2021; originally announced November 2021.

Comments: An erratum to this manuscript has been published here: https://www.aanda.org/articles/aa/full_html/2022/07/aa42281e-21/aa42281e-21.html

Journal ref: A&A 657, L13 (2022)

The effect of viewing angle on the Kennicutt-Schmidt relation of the local molecular clouds

Authors: Jouni Kainulainen, Sara Rezaei Kh., Andri Spilker, Jan Orkisz

Abstract: The Gaia data give us an unprecedented view to the 3-dimensional (3D) structure of molecular clouds in the Solar neighbourhood. We study how the projected areas and masses of clouds, and consequently the Kennicutt-Schmidt relation (KS-relation), depend on the viewing angle. We derive the probability distributions of the projected areas and masses for nine clouds within 400 pc from the Sun using 3D… ▽ More The Gaia data give us an unprecedented view to the 3-dimensional (3D) structure of molecular clouds in the Solar neighbourhood. We study how the projected areas and masses of clouds, and consequently the Kennicutt-Schmidt relation (KS-relation), depend on the viewing angle. We derive the probability distributions of the projected areas and masses for nine clouds within 400 pc from the Sun using 3D dust distribution data from the literature. We find that the viewing angle can have a dramatic effect on the observed areas and masses of individual clouds. The joint probability distributions of the areas and masses are strongly correlated, relatively flat, and can show multiple peaks. The typical ranges and 50% quartiles of the distributions are roughly 100-200% and 20-80% of the median value, respectively, making viewing angle effects important for all individual clouds. The threshold value used to define the cloud areas is also important; our analysis suggests that the clouds become more anisotropic for smaller thresholds (larger scales). On average, the areas and masses of the plane-of-the-sky and face-on projections agree, albeit with a large scatter. This suggests that sample averages of areas and masses are relatively free of viewing angle effects, which is important to facilitate comparisons of extragalactic and galactic data. Ultimately, our results demonstrate that a cloud's location in the KS-relation is affected by viewing angle in a non-trivial manner. However, the KS-relation of our sample as a whole is not strongly affected by these effects, because the co-variance of the areas and masses causes the observed mean column density to remain relatively constant. △ Less

Submitted 15 November, 2021; originally announced November 2021.

Comments: Accepted to A&A Letters, 15 pages including appendices

Journal ref: A&A 659, L6 (2022)

Bird's eye view of molecular clouds in the Milky Way: I. Column density and star formation from sub-pc to kpc scales

Authors: Andri Spilker, Jouni Kainulainen, Jan Orkisz

Abstract: Describing how the properties of the interstellar medium combine across size-scales is crucial for understanding star formation scaling laws and connecting Galactic and extragalactic data of molecular clouds. We describe how the statistical structure of clouds, and its connection to star formation, changes from sub-parsec to kiloparsec scales in a complete region within the Milky Way disk. We buil… ▽ More Describing how the properties of the interstellar medium combine across size-scales is crucial for understanding star formation scaling laws and connecting Galactic and extragalactic data of molecular clouds. We describe how the statistical structure of clouds, and its connection to star formation, changes from sub-parsec to kiloparsec scales in a complete region within the Milky Way disk. We build a census of molecular clouds within 2 kpc from the Sun using literature. We examine the dust-based column density probability distributions (N-PDFs) of the clouds and their relation to star formation traced by young stellar objects (YSOs). We then examine our survey region from the outside, within apertures of varying sizes, and describe how the N-PDFs and their relation to star formation changes with the size-scale. The N-PDFs of the clouds are not well described by any single simple model; use of any single model may bias the interpretation of the N-PDFs. The top-heaviness of the N-PDFs correlates with star formation activity, and the correlation changes with Galactic environment (spiral-/inter-arm regions). We find that the density contrast of clouds may be more intimately linked to star formation than the dense gas mass fraction. The aperture-averaged N-PDFs vary with the size-scale and are more top-heavy for larger apertures. The top-heaviness of the aperture N-PDFs correlates with star formation activity up to roughly 0.5 kpc size-scale, depending on the environment. Our results suggest that the relations between cloud structure and star formation are environment specific and best captured by relative quantities (e.g., the density contrast). Finally, we show how the density structures of individual clouds give rise to a kpc-scale Kennicutt-Schmidt relationship as a combination of sampling effects and blending of different galactic environments. △ Less

Submitted 2 September, 2021; v1 submitted 10 August, 2021; originally announced August 2021.

Comments: 19 pages and 19 figures in main text. Accepted for publication in A&A

Journal ref: A&A 653, A63 (2021)

The SEDIGISM survey: first data release and overview of the Galactic structure

Authors: F. Schuller, J. S. Urquhart, T. Csengeri, D. Colombo, A. Duarte-Cabral, M. Mattern, A. Ginsburg, A. R. Pettitt, F. Wyrowski, L. Anderson, F. Azagra, P. Barnes, M. Beltran, H. Beuther, S. Billington, L. Bronfman, R. Cesaroni, C. Dobbs, D. Eden, M. -Y. Lee, S. -N. X. Medina, K. M. Menten, T. Moore, F. M. Montenegro-Montes, S. Ragan , et al. (35 additional authors not shown)

Abstract: The SEDIGISM (Structure, Excitation and Dynamics of the Inner Galactic Interstellar Medium) survey used the APEX telescope to map 84 deg^2 of the Galactic plane between l = -60 deg and l = +31 deg in several molecular transitions, including 13CO(2-1) and C18O(2-1), thus probing the moderately dense (~10^3 cm^-3) component of the interstellar medium. With an angular resolution of 30'' and a typical… ▽ More The SEDIGISM (Structure, Excitation and Dynamics of the Inner Galactic Interstellar Medium) survey used the APEX telescope to map 84 deg^2 of the Galactic plane between l = -60 deg and l = +31 deg in several molecular transitions, including 13CO(2-1) and C18O(2-1), thus probing the moderately dense (~10^3 cm^-3) component of the interstellar medium. With an angular resolution of 30'' and a typical 1-sigma sensitivity of 0.8-1.0 K at 0.25 km/s velocity resolution, it gives access to a wide range of structures, from individual star-forming clumps to giant molecular clouds and complexes. The coverage includes a good fraction of the first and fourth Galactic quadrants, allowing us to constrain the large scale distribution of cold molecular gas in the inner Galaxy. In this paper we provide an updated overview of the full survey and the data reduction procedures used. We also assess the quality of these data and describe the data products that are being made publicly available as part of this first data release (DR1). We present integrated maps and position-velocity maps of the molecular gas and use these to investigate the correlation between the molecular gas and the large scale structural features of the Milky Way such as the spiral arms, Galactic bar and Galactic centre. We find that approximately 60 per cent of the molecular gas is associated with the spiral arms and these appear as strong intensity peaks in the derived Galactocentric distribution. We also find strong peaks in intensity at specific longitudes that correspond to the Galactic centre and well known star forming complexes, revealing that the 13CO emission is concentrated in a small number of complexes rather than evenly distributed along spiral arms. △ Less

Submitted 2 December, 2020; originally announced December 2020.

Comments: MNRAS, in press

The SEDIGISM survey: Molecular clouds in the inner Galaxy

Authors: A. Duarte-Cabral, D. Colombo, J. S. Urquhart, A. Ginsburg, D. Russeil, F. Schuller, L. D. Anderson, P. J. Barnes, M. T. Beltran, H. Beuther, S. Bontemps, L. Bronfman, T. Csengeri, C. L. Dobbs, D. Eden, A. Giannetti, J. Kauffmann, M. Mattern, S. -N. X. Medina, K. M. Menten, M. -Y. Lee, A. R. Pettitt, M. Riener, A. J. Rigby, A. Trafficante , et al. (35 additional authors not shown)

Abstract: We use the 13CO(2-1) emission from the SEDIGISM high-resolution spectral-line survey of the inner Galaxy, to extract the molecular cloud population with a large dynamic range in spatial scales, using the SCIMES algorithm. This work compiles a cloud catalogue with a total of 10663 molecular clouds, 10300 of which we were able to assign distances and compute physical properties. We study some of the… ▽ More We use the 13CO(2-1) emission from the SEDIGISM high-resolution spectral-line survey of the inner Galaxy, to extract the molecular cloud population with a large dynamic range in spatial scales, using the SCIMES algorithm. This work compiles a cloud catalogue with a total of 10663 molecular clouds, 10300 of which we were able to assign distances and compute physical properties. We study some of the global properties of clouds using a science sample, consisting of 6664 well resolved sources and for which the distance estimates are reliable. In particular, we compare the scaling relations retrieved from SEDIGISM to those of other surveys, and we explore the properties of clouds with and without high-mass star formation. Our results suggest that there is no single global property of a cloud that determines its ability to form massive stars, although we find combined trends of increasing mass, size, surface density and velocity dispersion for the sub-sample of clouds with ongoing high-mass star formation. We then isolate the most extreme clouds in the SEDIGISM sample (i.e. clouds in the tails of the distributions) to look at their overall Galactic distribution, in search for hints of environmental effects. We find that, for most properties, the Galactic distribution of the most extreme clouds is only marginally different to that of the global cloud population. The Galactic distribution of the largest clouds, the turbulent clouds and the high-mass star-forming clouds are those that deviate most significantly from the global cloud population. We also find that the least dynamically active clouds (with low velocity dispersion or low virial parameter) are situated further afield, mostly in the least populated areas. However, we suspect that part of these trends may be affected by some observational biases, and thus require further follow up work in order to be confirmed. △ Less

Submitted 2 December, 2020; originally announced December 2020.

Comments: 25 pages (+ appendices, 15 pages), 26 figures, MNRAS

SEDIGISM-ATLASGAL: Dense Gas Fraction and Star Formation Efficiency Across the Galactic Disk

Authors: J. S. Urquhart, C. Figura, J. R. Cross, M. R. A. Wells, T. J. T. Moore, D. J. Eden, S. E. Ragan, A. R. Pettitt, A. Duarte-Cabral, D. Colombo, F. Schuller, T. Csengeri, M. Mattern, H. Beuther, K. M. Menten, F. Wyrowski, L. D. Anderson, P. J. Barnes, M. T. Beltrán, S. J. Billington, L. Bronfman, A. Giannetti, J. Kainulainen, J. Kauffmann, M. -Y. Lee , et al. (10 additional authors not shown)

Abstract: By combining two surveys covering a large fraction of the molecular material in the Galactic disk we investigate the role the spiral arms play in the star formation process. We have matched clumps identified by ATLASGAL with their parental GMCs as identified by SEDIGISM, and use these giant molecular cloud (GMC) masses, the bolometric luminosities, and integrated clump masses obtained in a concurr… ▽ More By combining two surveys covering a large fraction of the molecular material in the Galactic disk we investigate the role the spiral arms play in the star formation process. We have matched clumps identified by ATLASGAL with their parental GMCs as identified by SEDIGISM, and use these giant molecular cloud (GMC) masses, the bolometric luminosities, and integrated clump masses obtained in a concurrent paper to estimate the dense gas fractions (DGF$_{\rm gmc}=\sum M_{\rm clump}/M_{\rm gmc}$) and the instantaneous star forming efficiencies (i.e., SFE$_{\rm gmc} = \sum L_{\rm clump}/M_{\rm gmc}$). We find that the molecular material associated with ATLASGAL clumps is concentrated in the spiral arms ($\sim$60% found within $\pm$10 km s$^{-1}$ of an arm). We have searched for variations in the values of these physical parameters with respect to their proximity to the spiral arms, but find no evidence for any enhancement that might be attributable to the spiral arms. The combined results from a number of similar studies based on different surveys indicate that, while spiral-arm location plays a role in cloud formation and HI to H$_2$ conversion, the subsequent star formation processes appear to depend more on local environment effects. This leads us to conclude that the enhanced star formation activity seen towards the spiral arms is the result of source crowding rather than the consequence of a any physical process. △ Less

Submitted 2 December, 2020; originally announced December 2020.

Comments: 15 pages, 12 figures. Full list of affiliations can be found at the end of the paper

Quantitative inference of the $H_2$ column densities from 3 mm molecular emission: A case study towards Orion B

Authors: Pierre Gratier, Jérôme Pety, Emeric Bron, Antoine Roueff, Jan H. Orkisz, Maryvonne Gerin, Victor de Souza Magalhaes, Mathilde Gaudel, Maxime Vono, Sébastien Bardeau, Jocelyn Chanussot, Pierre Chainais, Javier R. Goicoechea, Viviana V. Guzmán, Annie Hughes, Jouni Kainulainen, David Languignon, Jacques Le Bourlot, Franck Le Petit, François Levrier, Harvey Liszt, Nicolas Peretto, Evelyne Roueff, Albrecht Sievers

Abstract: Molecular hydrogen being unobservable in cold molecular clouds, the column density measurements of molecular gas currently rely either on dust emission observation in the far-IR or on star counting. (Sub-)millimeter observations of numerous trace molecules are effective from ground based telescopes, but the relationships between the emission of one molecular line and the H2 column density (NH2) is… ▽ More Molecular hydrogen being unobservable in cold molecular clouds, the column density measurements of molecular gas currently rely either on dust emission observation in the far-IR or on star counting. (Sub-)millimeter observations of numerous trace molecules are effective from ground based telescopes, but the relationships between the emission of one molecular line and the H2 column density (NH2) is non-linear and sensitive to excitation conditions, optical depths, abundance variations due to the underlying physico-chemistry. We aim to use multi-molecule line emission to infer NH2 from radio observations. We propose a data-driven approach to determine NH2 from radio molecular line observations. We use supervised machine learning methods (Random Forests) on wide-field hyperspectral IRAM-30m observations of the Orion B molecular cloud to train a predictor of NH2, using a limited set of molecular lines as input, and the Herschel-based dust-derived NH2 as ground truth output. For conditions similar to the Orion B molecular cloud, we obtain predictions of NH2 within a typical factor of 1.2 from the Herschel-based estimates. An analysis of the contributions of the different lines to the predictions show that the most important lines are $^{13}$CO(1-0), $^{12}$CO(1-0), C$^{18}$O(1-0), and HCO$^+$(1-0). A detailed analysis distinguishing between diffuse, translucent, filamentary, and dense core conditions show that the importance of these four lines depends on the regime, and that it is recommended to add the N$_2$H$^+$(1-0) and CH$_3$OH(20-10) lines for the prediction of NH2 in dense core conditions. This article opens a promising avenue to directly infer important physical parameters from the molecular line emission in the millimeter domain. The next step will be to try to infer several parameters simultaneously (e.g., NH2 and far-UV illumination field) to further test the method. [Abridged] △ Less

Submitted 31 August, 2020; originally announced August 2020.

Comments: Accepted for publication in A&A

Journal ref: A&A 645, A27 (2021)

Tracers of the ionization fraction in dense and translucent gas: I. Automated exploitation of massive astrochemical model grids

Authors: Emeric Bron, Evelyne Roueff, Maryvonne Gerin, Jérôme Pety, Pierre Gratier, Franck Le Petit, Viviana Guzman, Jan H. Orkisz, Victor de Souza Magalhaes, Mathilde Gaudel, Maxime Vono, Sébastien Bardeau, Pierre Chainais, Javier R. Goicoechea, Annie Hughes, Jouni Kainulainen, David Languignon, Jacques Le Bourlot, François Levrier, Harvey Liszt, Karin Öberg, Nicolas Peretto, Antoine Roueff, Albrecht Sievers

Abstract: The ionization fraction plays a key role in the physics and chemistry of the neutral interstellar medium, from controlling the coupling of the gas to the magnetic field to allowing fast ion-neutral reactions that drive interstellar chemistry. Most estimations of the ionization fraction have relied on deuterated species such as DCO+, whose detection is limited to dense cores representing an extreme… ▽ More The ionization fraction plays a key role in the physics and chemistry of the neutral interstellar medium, from controlling the coupling of the gas to the magnetic field to allowing fast ion-neutral reactions that drive interstellar chemistry. Most estimations of the ionization fraction have relied on deuterated species such as DCO+, whose detection is limited to dense cores representing an extremely small fraction of the volume of the giant molecular clouds they are part of. As large field-of-view hyperspectral maps become available, new tracers may be found. We search for the best observable tracers of the ionization fraction based on a grid of astrochemical models. We build grids of models that sample randomly a large space of physical conditions (unobservable quantities such as gas density, temperature, etc.) and compute the corresponding observables (line intensities, column densities) and the ionization fraction. We estimate the predictive power of each potential tracer by training a Random Forest model to predict the ionization fraction from that tracer, based on these model grids. In both translucent medium and cold dense medium conditions, several observable tracers with very good predictive power for the ionization fraction are found. Several tracers in cold dense medium conditions are found to be better and more widely applicable than the traditional DCO+/HCO+ ratio. We also provide simpler analytical fits for estimating the ionization fraction from the best tracers, and for estimating the associated uncertainties. We discuss the limitations of the present study and select a few recommended tracers in both types of conditions. The method presented here is very general and can be applied to the measurement of any other quantity of interest (cosmic ray flux, elemental abundances, etc.) from any type of model (PDR models, time-dependent chemical models, etc.). (abridged) △ Less

Submitted 27 July, 2020; originally announced July 2020.

Comments: 29 pages, 15 figures. Accepted for publication in A&A

Journal ref: A&A 645, A28 (2021)

Ubiquitous velocity fluctuations throughout the molecular interstellar medium

Authors: J. D. Henshaw, J. M. D. Kruijssen, S. N. Longmore, M. Riener, A. K. Leroy, E. Rosolowsky, A. Ginsburg, C. Battersby, M. Chevance, S. E. Meidt, S. C. O. Glover, A. Hughes, J. Kainulainen, R. S. Klessen, E. Schinnerer, A. Schruba, H. Beuther, F. Bigiel, G. A. Blanc, E. Emsellem, T. Henning, C. N. Herrera, E. W. Koch, J. Pety, S. E. Ragan , et al. (1 additional authors not shown)

Abstract: The density structure of the interstellar medium (ISM) determines where stars form and release energy, momentum, and heavy elements, driving galaxy evolution. Density variations are seeded and amplified by gas motion, but the exact nature of this motion is unknown across spatial scale and galactic environment. Although dense star-forming gas likely emerges from a combination of instabilities, conv… ▽ More The density structure of the interstellar medium (ISM) determines where stars form and release energy, momentum, and heavy elements, driving galaxy evolution. Density variations are seeded and amplified by gas motion, but the exact nature of this motion is unknown across spatial scale and galactic environment. Although dense star-forming gas likely emerges from a combination of instabilities, convergent flows, and turbulence, establishing the precise origin is challenging because it requires quantifying gas motion over many orders of magnitude in spatial scale. Here we measure the motion of molecular gas in the Milky Way and in nearby galaxy NGC 4321, assembling observations that span an unprecedented spatial dynamic range ($10^{-1}{-}10^3$ pc). We detect ubiquitous velocity fluctuations across all spatial scales and galactic environments. Statistical analysis of these fluctuations indicates how star-forming gas is assembled. We discover oscillatory gas flows with wavelengths ranging from $0.3{-}400$ pc. These flows are coupled to regularly-spaced density enhancements that likely form via gravitational instabilities. We also identify stochastic and scale-free velocity and density fluctuations, consistent with the structure generated in turbulent flows. Our results demonstrate that ISM structure cannot be considered in isolation. Instead, its formation and evolution is controlled by nested, interdependent flows of matter covering many orders of magnitude in spatial scale. △ Less

Submitted 3 July, 2020; originally announced July 2020.

Comments: Published in Nature Astronomy on July 6th 2020. This is the authors' version before final edits. Includes methods and supplementary information. Link to the NA publication: https://www.nature.com/articles/s41550-020-1126-z

Autonomous Gaussian decomposition of the Galactic Ring Survey. II. The Galactic distribution of 13CO

Authors: M. Riener, J. Kainulainen, J. D. Henshaw, H. Beuther

Abstract: Knowledge about the distribution of CO emission in the Milky Way is essential to understand the impact of Galactic environment on the formation and evolution of structures in the interstellar medium. However, currently our insight about the fraction of CO in spiral arm and interarm regions is still limited by large uncertainties in assumed rotation curve models or distance determination techniques… ▽ More Knowledge about the distribution of CO emission in the Milky Way is essential to understand the impact of Galactic environment on the formation and evolution of structures in the interstellar medium. However, currently our insight about the fraction of CO in spiral arm and interarm regions is still limited by large uncertainties in assumed rotation curve models or distance determination techniques. In this work we use a Bayesian approach to obtain the current best assessment of the distribution of 13CO from the Galactic Ring Survey. We performed two different distance estimates that either included or excluded a model for Galactic features. We also include a prior for the solution of the kinematic distance ambiguity that was determined from a compilation of literature distances and an assumed size-linewidth relationship. We find that the fraction of 13CO emission associated with spiral arm features varies from 76% to 84% between the two distance runs. The vertical distribution of the gas is concentrated around the Galactic midplane showing FWHM values of ~75 pc. We do not find any significant difference between gas emission properties associated with spiral arm and interarm features. In particular the distribution of velocity dispersion values of gas emission in spurs and spiral arms is very similar. We detect a trend of higher velocity dispersion values with increasing heliocentric distance, which we attribute to beam averaging effects caused by differences in spatial resolution. We argue that the true distribution of the gas emission is likely more similar to a combination of our two distance results, and highlight the importance of using complementary distance estimations to safeguard against the pitfalls of any single approach. We conclude that the methodology presented in this work is a promising way to determine distances to gas emission features in Galactic plane surveys. △ Less

Submitted 5 June, 2020; originally announced June 2020.

Comments: Accepted for publication in A&A. 27 pages, 31 figures (main text: 18 pages, 21 figures). The PDF contains multi-layer figures that can be displayed with Adobe Acrobat. Table 2 with the distance results is available in its entirety at the CDS. The Python wrapper for the Bayesian Distance Calculator (BD_wrapper) is freely available at GitHub under https://github.com/mriener/BD_wrapper

Journal ref: A&A 640, A72 (2020)

C18O, 13CO, and 12CO abundances and excitation temperatures in the Orion B molecular cloud: An analysis of the precision achievable when modeling spectral line within the Local Thermodynamic Equilibrium approximation

Authors: Antoine Roueff, Maryvonne Gerin, Pierre Gratier, Francois Levrier, Jerome Pety, Mathilde Gaudel, Javier R. Goicoechea, Jan H. Orkisz, Victor de Souza Magalhaes, Maxime Vono, Sebastien Bardeau, Emeric Bron, Jocelyn Chanussot, Pierre Chainais, Viviana V. Guzman, Annie Hughes, Jouni Kainulainen, David Languignon, Jacques Le Bourlot, Franck Le Petit, Harvey S. Liszt, Antoine Marchal, Marc-Antoine Miville-Deschenes, Nicolas Peretto, Evelyne Roueff , et al. (1 additional authors not shown)

Abstract: CO isotopologue transitions are routinely observed in molecular clouds to probe the column density of the gas, the elemental ratios of carbon and oxygen, and to trace the kinematics of the environment. We aim at estimating the abundances, excitation temperatures, velocity field and velocity dispersions of the three main CO isotopologues towards a subset of the Orion B molecular cloud. We use the C… ▽ More CO isotopologue transitions are routinely observed in molecular clouds to probe the column density of the gas, the elemental ratios of carbon and oxygen, and to trace the kinematics of the environment. We aim at estimating the abundances, excitation temperatures, velocity field and velocity dispersions of the three main CO isotopologues towards a subset of the Orion B molecular cloud. We use the Cramer Rao Bound (CRB) technique to analyze and estimate the precision of the physical parameters in the framework of local-thermodynamic-equilibrium excitation and radiative transfer with an additive white Gaussian noise. We propose a maximum likelihood estimator to infer the physical conditions from the 1-0 and 2-1 transitions of CO isotopologues. Simulations show that this estimator is unbiased and efficient for a common range of excitation temperatures and column densities (Tex > 6 K, N > 1e14 - 1e15 cm-2). Contrary to the general assumptions, the different CO isotopologues have distinct excitation temperatures, and the line intensity ratios between different isotopologues do not accurately reflect the column density ratios. We find mean fractional abundances that are consistent with previous determinations towards other molecular clouds. However, significant local deviations are inferred, not only in regions exposed to UV radiation field but also in shielded regions. These deviations result from the competition between selective photodissociation, chemical fractionation, and depletion on grain surfaces. We observe that the velocity dispersion of the C18O emission is 10% smaller than that of 13CO. The substantial gain resulting from the simultaneous analysis of two different rotational transitions of the same species is rigorously quantified. The CRB technique is a promising avenue for analyzing the estimation of physical parameters from the fit of spectral lines. △ Less

Submitted 17 May, 2020; originally announced May 2020.

Comments: 27 pages, 23 PDF figures. Accepted for publication in A&A. Uses aa latex macro

Journal ref: A&A 645, A26 (2021)

Cloud formation in the atomic and molecular phase: HI self absorption (HISA) towards a Giant Molecular Filament

Authors: Y. Wang, S. Bihr, H. Beuther, M. R. Rugel, J. D. Soler, J. Ott, J. Kainulainen, N. Schneider, R. S. Klessen, S. C. O. Glover, N. M. McClure-Griffiths, P. F. Goldsmith, K. G. Johnston, K. M. Menten, S. Ragan, L. D. Anderson, J. S. Urquhart, H. Linz, N. Roy, R. J. Smith, F. Bigiel, T. Henning, S. N. Longmore

Abstract: Molecular clouds form from the atomic phase of the interstellar medium. However, characterizing the transition between the atomic and the molecular interstellar medium (ISM) is a difficult observational task. Here we address cloud formation processes by combining HSIA with molecular line data. One scenario proposed by numerical simulations is that the column density probability density functions (… ▽ More Molecular clouds form from the atomic phase of the interstellar medium. However, characterizing the transition between the atomic and the molecular interstellar medium (ISM) is a difficult observational task. Here we address cloud formation processes by combining HSIA with molecular line data. One scenario proposed by numerical simulations is that the column density probability density functions (N-PDF) evolves from a log-normal shape at early times to a power-law-like shape at later times. In this paper, we study the cold atomic component of the giant molecular filament GMF38a (d=3.4 kpc, length$\sim230$ pc). We identify an extended HISA feature, which is partly correlated with the 13CO emission. The peak velocities of the HISA and 13CO observations agree well on the eastern side of the filament, whereas a velocity offset of approximately 4 km s$^{-1}$ is found on the western side. The sonic Mach number we derive from the linewidth measurements shows that a large fraction of the HISA, which is ascribed to the cold neutral medium (CNM), is at subsonic and transonic velocities. The column density of the CNM is on the order of 10$^{20}$ to 10$^{21}$ cm$^{-2}$. The column density of molecular hydrogen is an order of magnitude higher. The N-PDFs from HISA (CNM), HI emission (WNM+CNM), and 13CO (molecular component) are well described by log-normal functions, which is in agreement with turbulent motions being the main driver of cloud dynamics. The N-PDF of the molecular component also shows a power law in the high column-density region, indicating self-gravity. We suggest that we are witnessing two different evolutionary stages within the filament. The eastern subregion seems to be forming a molecular cloud out of the atomic gas, whereas the western subregion already shows high column density peaks, active star formation and evidence of related feedback processes. △ Less

Submitted 3 January, 2020; originally announced January 2020.

Comments: 22 pages, 25 figures, accepted for publication by A&A

Journal ref: A&A 634, A139 (2020)

Autonomous Gaussian decomposition of the Galactic Ring Survey. I. Global statistics and properties of the 13CO emission data

Authors: M. Riener, J. Kainulainen, H. Beuther, J. D. Henshaw, J. H. Orkisz, Y. Wang

Abstract: The analysis of large molecular line surveys of the Galactic plane is essential for our understanding of the gas kinematics on Galactic scales, in particular its link with the formation and evolution of dense structures in the interstellar medium. An approximation of the emission peaks with Gaussian functions allows for an efficient and straightforward extraction of useful physical information con… ▽ More The analysis of large molecular line surveys of the Galactic plane is essential for our understanding of the gas kinematics on Galactic scales, in particular its link with the formation and evolution of dense structures in the interstellar medium. An approximation of the emission peaks with Gaussian functions allows for an efficient and straightforward extraction of useful physical information contained in the shape and Doppler-shifted frequency of the emission lines contained in these enormous data sets. In this work we present an overview and first results of a Gaussian decomposition of the entire Galactic Ring Survey (GRS) 13CO (1-0) data that consists of about 2.3 million spectra. We performed the decomposition with the fully automated GaussPy+ algorithm and fitted about 4.6 million Gaussian components to the GRS spectra. We discuss the statistics of the fit components and relations between the fitted intensities, velocity centroids, and velocity dispersions. We find that the magnitude of the velocity dispersion values increase toward the inner Galaxy and around the Galactic midplane, which we speculate is partly due to the influence of the Galactic bar and regions with higher non-thermal motions located in the midplane, respectively. We also use our decomposition results to infer global properties of the gas emission and find that the number of fit components used per spectrum is indicative for the amount of structure along the line of sight. We find that the emission lines from regions located on the far side of the Galaxy show increased velocity dispersion values, likely due to beam averaging effects. We demonstrate how this trend has the potential to aid in characterising Galactic structure by disentangling emission that is belonging to the nearby Aquila Rift molecular cloud from emission that is more likely associated with the Perseus and Outer spiral arms. △ Less

Submitted 15 November, 2019; originally announced November 2019.

Comments: Accepted for publication in A&A. 26 pages, 25 figures (main text: 16 pages, 16 figures; Appendix: 10 pages, 9 figures). The PDF contains multi-layer figures that can be displayed with Adobe Acrobat. Table 1 with the full decomposition results will be made publicly available soon; earlier access to the data is possible upon request from the corresponding author (M. Riener). Comments are welcome

Journal ref: A&A 633, A14 (2020)

Deep PSF photometric catalog of the VVV survey data

Authors: Miaomiao Zhang, Jouni Kainulainen

Abstract: The Vista Variables in the Via Lactea (VVV) survey has performed a multi-epoch near-infrared imaging of the inner Galactic plane. High-fidelity photometric catalogs are needed to utilize the data. We aim at producing a deep, point-spread-function (PSF) photometric catalog for the VVV survey J, H, and Ks band data. Specifically, we aim at taking advantage of all the epochs of the survey to reach hi… ▽ More The Vista Variables in the Via Lactea (VVV) survey has performed a multi-epoch near-infrared imaging of the inner Galactic plane. High-fidelity photometric catalogs are needed to utilize the data. We aim at producing a deep, point-spread-function (PSF) photometric catalog for the VVV survey J, H, and Ks band data. Specifically, we aim at taking advantage of all the epochs of the survey to reach high limiting magnitudes. We develop an automatic PSF-fitting pipeline based on the DaoPHOT algorithm and perform photometry on the stacked VVV images in J, H, and Ks bands. We present a PSF photometric catalog in the Vega system that contains about 926 million sources in the J, H, and Ks filters. About 10% of the sources are flagged as possible spurious detections. The 5 sigma limiting magnitudes of the sources with high reliability are about 20.8, 19.5, and 18.7 mag in the J, H, and Ks band, respectively, depending on the local crowding condition. Our photometric catalog reaches on average about one magnitude deeper than the previously released PSF DoPHOT photometric catalog. It also includes less spurious detections. There are significant differences in the brightnesses of faint sources between our catalog and the previously released one. The likely origin of these differences is in the different photometric algorithms that are utilized; it is not straightforward to assess which catalog is more accurate in which situations. Our new catalog is beneficial especially for science goals that require high limiting magnitudes; our catalog reaches such in fields that have a relatively uniform source number density. Overall, the limiting magnitudes and completeness are different in the fields with different crowding conditions. △ Less

Submitted 31 October, 2019; v1 submitted 18 October, 2019; originally announced October 2019.

Comments: 19 pages, 22 figures, and 2 tables. Accepted by A&A

Journal ref: A&A 632, A85 (2019)

Bow Magnetic Morphology Surrounding Filamentary Molecular Clouds? 3D Magnetic Field Structure of Orion-A

Authors: M. Tahani, R. Plume, J. C. Brown, J. D. Soler, J. Kainulainen

Abstract: Using a method based on Faraday rotation measurements, Tahani et al. 2018 find the line-of-sight component of magnetic fields in Orion-A and show that their direction changes from the eastern side of this filamentary structure to its western side. Three possible magnetic field morphologies that can explain this reversal across the Orion-A region are toroidal, helical, and bow-shaped morphologies.… ▽ More Using a method based on Faraday rotation measurements, Tahani et al. 2018 find the line-of-sight component of magnetic fields in Orion-A and show that their direction changes from the eastern side of this filamentary structure to its western side. Three possible magnetic field morphologies that can explain this reversal across the Orion-A region are toroidal, helical, and bow-shaped morphologies. In this paper we construct simple models to represent these three morphologies and compare them with the available observational data to find the most probable morphology(ies). To compare the observations with the models, we use probability values and a Monte-Carlo analysis to determine the most likely magnetic field morphology among these three morphologies. We find that the bow morphology has the highest probability values and that our Monte-Carlo analysis suggests that the bow morphology is more likely. We suggest that the bow morphology is the most likely and the most natural of the three morphologies that could explain a magnetic field reversal across the Orion-A filamentary structure (i.e., bow, helical and toroidal morphologies). △ Less

Submitted 16 October, 2019; originally announced October 2019.

Comments: 11 pages, 15 Figures

Journal ref: A&A 632, A68 (2019)

OH maser emission in the THOR survey of the northern Milky Way

Authors: H. Beuther, A. Walsh, Y. Wang, M. Rugel, J. Soler, H. Linz, R. S. Klessen, L. D. Anderson, J. S. urquhart, S. C. O. Glover, S. J. Billington, J. Kainulainen, K. M. Menten, N. Roy, S. N. Longmore, F. Bigiel

Abstract: Context: OH masers trace diverse physical processes, from the expanding envelopes around evolved stars to star-forming regions or supernovae remnants. Aims: We identify the ground-state OH masers at 18cm wavelength in the area covered by ``The HI/OH/Recombination line survey of the Milky Way (THOR)''. We present a catalogue of all OH maser features and their possible associated environments. Metho… ▽ More Context: OH masers trace diverse physical processes, from the expanding envelopes around evolved stars to star-forming regions or supernovae remnants. Aims: We identify the ground-state OH masers at 18cm wavelength in the area covered by ``The HI/OH/Recombination line survey of the Milky Way (THOR)''. We present a catalogue of all OH maser features and their possible associated environments. Methods: The THOR survey covers longitude and latitude ranges of 14.3<l<66.8 and b<1.25 deg. All OH ground state lines at 1612, 1665, 1667 and 1720MHz have been observed, employing the Very Large Array. The spatial resolution of the data varies between 12.5'' and 19'', the spectral resolution is 1.5km/s, and the rms sensitivity of the data is ~10mJy/beam per channel. Results: We identify 1585 individual maser spots distributed over 807 maser sites. Based on different criteria from spectral profiles to literature comparison, we try to associate the maser sites with astrophysical source types. Approximately 51\% of the sites exhibit the double-horned 1612MHz spectra typically emitted from the expanding shells of evolved stars. The separations of the two main velocity features of the expanding shells typically vary between 22 and 38km/s. In addition to this, at least 20% of the maser sites are associated with star-forming regions. While the largest fraction of 1720MHz maser spots (21 out of 53) is associated with supernova remnants, a significant fraction of the 1720MHz maser spots (17) are also associated with star-forming regions. We present comparisons to the thermal 13CO(1-0) emission as well as to other surveys of class II CH3OH and H2O maser emission. The catalogue attempts to present associations to astrophysical sources where available, and the full catalogue is available in electronic form. △ Less

Submitted 26 July, 2019; originally announced July 2019.

Comments: The main paper contains 13 pages and 7 pages. However, the full paper with appendix where all maser sites are shown and also the full catalogue is presented can be found at http://www.mpia.de/thor/Data_%26_Publications.html

Journal ref: A&A 628, A90 (2019)

Dynamics of cluster-forming hub-filament systems: The case of the high-mass star-forming complex Monoceros R2

Authors: S. P. Trevino-Morales, A. Fuente, A. Sanchez-Monge, J. Kainulainen, P. Didelon, S. Suri, N. Schneider, J. Ballesteros-Paredes, Y. -N. Lee, P. Hennebelle, P. Pilleri, M. Gonzalez-Garcia, C. Kramer, S. Garcia-Burillo, A. Luna, J. R. Goicoechea, P. Tremblin, S. Geen

Abstract: High-mass stars and star clusters commonly form within hub-filament systems. Monoceros R2, harbors one of the closest such systems, making it an excellent target for case studies. We investigate the morphology, stability and dynamical properties of the hub-filament system on basis of 13CO and C18O observations obtained with the IRAM-30m telescope and H2 column density maps derived from Herschel du… ▽ More High-mass stars and star clusters commonly form within hub-filament systems. Monoceros R2, harbors one of the closest such systems, making it an excellent target for case studies. We investigate the morphology, stability and dynamical properties of the hub-filament system on basis of 13CO and C18O observations obtained with the IRAM-30m telescope and H2 column density maps derived from Herschel dust emission observations. We identified the filamentary network and characterized the individual filaments as either main (converging into the hub) or secondary (converging to a main filament) filaments. The main filaments have line masses of 30-100 Msun/pc and show signs of fragmentation. The secondary filaments have line masses of 12-60 Msun/pc and show fragmentation only sporadically. In the context of Ostriker's hydrostatic filament model, the main filaments are thermally super-critical. If non-thermal motions are included, most of them are trans-critical. Most of the secondary filaments are roughly trans-critical regardless of whether non-thermal motions are included or not. From the main filaments, we estimate a mass accretion rate of 10(-4)-10(-3) Msun/pc into the hub. The secondary filaments accrete into the main filaments with a rate of 0.1-0.4x10(-4) Msun/pc. The main filaments extend into the hub. Their velocity gradients increase towards the hub, suggesting acceleration of the gas. We estimate that with the observed infall velocity, the mass-doubling time of the hub is ~2.5 Myr, ten times larger than the free-fall time, suggesting a dynamically old region. These timescales are comparable with the chemical age of the HII region. Inside the hub, the main filaments show a ring- or a spiral-like morphology that exhibits rotation and infall motions. One possible explanation for the morphology is that gas is falling into the central cluster following a spiral-like pattern. △ Less

Submitted 8 July, 2019; originally announced July 2019.

Comments: The paper contains 19 pages of main text and 27 pages of appendix. The main text includes 15 figures and 2 tables. The appendix has 14 images, 5 regular size tables and a long table

Journal ref: A&A 629, A81 (2019)

GaussPy+: A fully automated Gaussian decomposition package for emission line spectra

Authors: Manuel Riener, Jouni Kainulainen, Jonathan D. Henshaw, Jan H. Orkisz, Claire E. Murray, Henrik Beuther

Abstract: Our understanding of the dynamics of the interstellar medium is informed by the study of the detailed velocity structure of emission line observations. One approach to study the velocity structure is to decompose the spectra into individual velocity components; this leads to a description of the dataset that is significantly reduced in complexity. However, this decomposition requires full automati… ▽ More Our understanding of the dynamics of the interstellar medium is informed by the study of the detailed velocity structure of emission line observations. One approach to study the velocity structure is to decompose the spectra into individual velocity components; this leads to a description of the dataset that is significantly reduced in complexity. However, this decomposition requires full automation lest it becomes prohibitive for large datasets, such as Galactic plane surveys. We developed GaussPy+, a fully automated Gaussian decomposition package that can be applied to emission line datasets, especially large surveys of HI and isotopologues of CO. We built our package upon the existing GaussPy algorithm and significantly improved its performance for noisy data. New functionalities of GaussPy+ include: i) automated preparatory steps, such as an accurate noise estimation, which can also be used as standalone applications; ii) an improved fitting routine; iii) an automated spatial refitting routine that can add spatial coherence to the decomposition results by refitting spectra based on neighbouring fit solutions. We thoroughly tested the performance of GaussPy+ on synthetic spectra and a test field from the Galactic Ring Survey. We found that GaussPy+ can deal with cases of complex emission and even low to moderate signal-to-noise values. △ Less

Submitted 25 June, 2019; originally announced June 2019.

Comments: Accepted for publication in A&A. 36 pages, 33 figures (main text: 23 pages, 18 figures; Appendix: 13 pages, 15 figures). The GaussPy+ package is written in Python 3 and is freely available at GitHub under https://github.com/mriener/gausspyplus

Young massive star cluster formation in the Galactic Centre is driven by global gravitational collapse of high-mass molecular clouds

Authors: A. T. Barnes, S. N. Longmore, A. Avison, Y. Contreras, A. Ginsburg, J. D. Henshaw, J. M. Rathborne, D. L. Walker, J. Alves, J. Bally, C. Battersby, M. T. Beltrán, H. Beuther, G. Garay, L. Gomez, J. Jackson, J. Kainulainen, J. M. D. Kruijssen, X. Lu, E. A. C. Mills, J. Ott, T. Peters

Abstract: Young massive clusters (YMCs) are the most compact, high-mass stellar systems still forming at the present day. The precursor clouds to such systems are, however, rare due to their large initial gas mass reservoirs and rapid dispersal timescales due to stellar feedback. Nonetheless, unlike their high-z counterparts, these precursors are resolvable down to the sites of individually forming stars, a… ▽ More Young massive clusters (YMCs) are the most compact, high-mass stellar systems still forming at the present day. The precursor clouds to such systems are, however, rare due to their large initial gas mass reservoirs and rapid dispersal timescales due to stellar feedback. Nonetheless, unlike their high-z counterparts, these precursors are resolvable down to the sites of individually forming stars, and hence represent the ideal environments in which to test the current theories of star and cluster formation. Using high angular resolution (1$^{\prime\prime}$ / 0.05pc) and sensitivity ALMA observations of two YMC progenitor clouds in the Galactic Centre, we have identified a suite of molecular line transitions -- e.g. c-C$_{3}$H$_{2} $($7-6$) -- that are believed to be optically thin, and reliably trace the gas structure in the highest density gas on star-forming core scales. We conduct a virial analysis of the identified core and proto-cluster regions, and show that half of the cores (5/10) and both proto-clusters are unstable to gravitational collapse. This is the first kinematic evidence of global gravitational collapse in YMC precursor clouds at such an early evolutionary stage. The implications are that if these clouds are to form YMCs, then they likely do so via the "conveyor-belt" mode, whereby stars continually form within dispersed dense gas cores as the cloud undergoes global gravitational collapse. The concurrent contraction of both the cluster-scale gas and embedded (proto)stars ultimately leads to the high (proto)stellar density in YMCs. △ Less

Submitted 14 March, 2019; originally announced March 2019.

Comments: Accepted for publication in MNRAS. 22 pages (+22 pages of appendix), 10 (+21 in appendix) figures, 5 (+1 in appendix) tables

Star-forming content of the giant molecular filaments in the Milky Way

Authors: Miaomiao Zhang, Jouni Kainulainen, Michael Mattern, Min Fang, Thomas Henning

Abstract: Observations have discovered numerous giant molecular filaments (GMFs) in the Milky Way. However, their role in the Galactic star formation and Galaxy-scale evolution of dense gas is still unknown. We investigate systematically the star-forming content of all currently known GMFs. This allows us to estimate the star formation rates (SFRs) of the GMFs and to establish relationships between the SFRs… ▽ More Observations have discovered numerous giant molecular filaments (GMFs) in the Milky Way. However, their role in the Galactic star formation and Galaxy-scale evolution of dense gas is still unknown. We investigate systematically the star-forming content of all currently known GMFs. This allows us to estimate the star formation rates (SFRs) of the GMFs and to establish relationships between the SFRs and the GMF properties. We identify and classify the young stellar object (YSO) population of each GMF using multi-wavelength photometry from near- to far-infrared. We estimate the total SFRs assuming a universal and fully sampled initial mass function and luminosity function. We uniformly estimate the physical properties of 57 GMFs. The GMFs show correlations between the 13CO line width, mass, and size, similar to Larson's relations. We identify 36 394 infrared excess sources in 57 GMFs and obtain SFRs for 46 GMFs. The median SFR surface density (Sigma_SFR) and star formation efficiency (SFE) of GMFs are 0.62 Msun Myr-1 pc-2 and 1%, similar to the nearby star-forming clouds. The star formation rate per free-fall time of GMFs is between 0.002-0.05 with the median value of 0.02. We also find a strong correlation between SFR and dense gas mass that is defined as gas mass above a visual extinction of 7 mag, which suggests that the SFRs of the GMFs scale similarly with dense gas as those of nearby molecular clouds. We also find a strong correlation between the mean SFR per unit length and dense gas mass per unit length. The origin of this scaling remains unknown, calling for further studies that can link the structure of GMFs to their SF activity and explore the differences between GMFs and other molecular clouds △ Less

Submitted 6 November, 2018; originally announced November 2018.

Comments: 74 pages, 91 figures, 7 appendix, accepted for publication in A&A

Journal ref: A&A 622, A52 (2019)

Histogram of oriented gradients: a technique for the study of molecular cloud formation

Authors: J. D. Soler, H. Beuther, M. Rugel, Y. Wang, P. C. Clark, S. C. O. Glover, P. F. Goldsmith, M. Heyer, L. D. Anderson, A. Goodman, Th. Henning, J. Kainulainen, R. S. Klessen, S. N. Longmore, N. M. McClure-Griffiths, K. M. Menten, J. C. Mottram, J. Ott, S. E. Ragan, R. J. Smith, J. S. Urquhart, F. Bigiel, P. Hennebelle, N. Roy, P. Schilke

Abstract: 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 corre… ▽ More 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. △ Less

Submitted 2 January, 2019; v1 submitted 21 September, 2018; originally announced September 2018.

Comments: 32 pages, 36 figures. Accepted for publication at A&A (28DEC2018)

Report number: AA/2018/34300

Journal ref: A&A 622, A166 (2019)

SEDIGISM: The kinematics of ATLASGAL filaments

Authors: M. Mattern, J. Kauffmann, T. Csengeri, J. S. Urquhart, S. Leurini, F. Wyrowski, A. Giannetti, P. J. Barnes, H. Beuther, L. Bronfman, A. Duarte-Cabral, T. Henning, J. Kainulainen, K. M. Menten, E. Schisano, F. Schuller

Abstract: Analysing the kinematics of filamentary molecular clouds is a crucial step towards understanding their role in the star formation process. Therefore, we study the kinematics of 283 filament candidates in the inner Galaxy, that were previously identified in the ATLASGAL dust continuum data. The $^{13}$CO(2 - 1) and C$^{18}$O(2 - 1) data of the SEDIGISM survey (Structure, Excitation, and Dynamics of… ▽ More Analysing the kinematics of filamentary molecular clouds is a crucial step towards understanding their role in the star formation process. Therefore, we study the kinematics of 283 filament candidates in the inner Galaxy, that were previously identified in the ATLASGAL dust continuum data. The $^{13}$CO(2 - 1) and C$^{18}$O(2 - 1) data of the SEDIGISM survey (Structure, Excitation, and Dynamics of the Inner Galactic Inter Stellar Medium) allows us to analyse the kinematics of these targets and to determine their physical properties at a resolution of 30 arcsec and 0.25 km/s. To do so, we developed an automated algorithm to identify all velocity components along the line-of-sight correlated with the ATLASGAL dust emission, and derive size, mass, and kinematic properties for all velocity components. We find two-third of the filament candidates are coherent structures in position-position-velocity space. The remaining candidates appear to be the result of a superposition of two or three filamentary structures along the line-of-sight. At the resolution of the data, on average the filaments are in agreement with Plummer-like radial density profiles with a power-law exponent of p = 1.5 +- 0.5, indicating that they are typically embedded in a molecular cloud and do not have a well-defined outer radius. Also, we find a correlation between the observed mass per unit length and the velocity dispersion of the filament of $m \sim σ_v^2$. We show that this relation can be explained by a virial balance between self-gravity and pressure. Another possible explanation could be radial collapse of the filament, where we can exclude infall motions close to the free-fall velocity. △ Less

Submitted 22 August, 2018; originally announced August 2018.

Comments: 45 pages, 20 pages for Table of derived parameters

Journal ref: A&A 619, A166 (2018)

Structure and Fragmentation of a high line-mass filament: Nessie

Authors: Michael Mattern, Jouni Kainulainen, Miaomiao Zhang, Henrik Beuther

Abstract: An increasing number of hundred-parsec scale, high line-mass filaments have been detected in the Galaxy. Their evolutionary path, including fragmentation towards star formation, is virtually unknown. We characterize the fragmentation within the Nessie filament, covering size-scales between $\sim$ 0.1-100 pc. We also connect the small-scale fragments to the star-forming potential of the cloud. We c… ▽ More An increasing number of hundred-parsec scale, high line-mass filaments have been detected in the Galaxy. Their evolutionary path, including fragmentation towards star formation, is virtually unknown. We characterize the fragmentation within the Nessie filament, covering size-scales between $\sim$ 0.1-100 pc. We also connect the small-scale fragments to the star-forming potential of the cloud. We combine near-infrared data from the VVV survey with mid-infrared GLIMPSE data to derive a high-resolution dust extinction map and apply a wavelet decomposition technique on it to analyze the fragmentation characteristics of the cloud, which are compared with predictions from fragmentation models. We compare the detected objects to those identified in $\sim$ 10 times coarser resolution from ATLASGAL data. We present a high-resolution extinction map of Nessie. We estimate the mean line-mass of Nessie to be $\sim$ 627 M$_\odot$/pc and the distance to be $\sim$ 3.5 kpc. We find that Nessie shows fragmentation at multiple size scales. The nearest-neighbour separations of the fragments at all scales are within a factor of 2 of the Jeans' length at that scale. However, the relationship between the mean densities of the fragments and their separations is significantly shallower than expected for Jeans' fragmentation. The relationship is similar to the one predicted for a filament that exhibits a Larson-like scaling between size-scale and velocity dispersion; such a scaling may result from turbulent support. Based on the number of YSOs in Nessie, we estimate that the star formation rate is $\sim$ 371 M$_\odot$/Myr; similar values result if using the number of dense cores, or the amount of dense gas, as the proxy of star formation. The star formation efficiency is 0.017. These numbers indicate that Nessie's star-forming content is comparable to the Solar neighborhood giant molecular clouds like Orion A. △ Less

Submitted 6 April, 2018; originally announced April 2018.

Journal ref: A&A 616, A78 (2018)

OH absorption in the first quadrant of the Milky Way as seen by THOR

Authors: M. R. Rugel, H. Beuther, S. Bihr, Y. Wang, J. Ott, A. Brunthaler, A. Walsh, S. C. O. Glover, P. F. Goldsmith, L. D. Anderson, N. Schneider, K. M. Menten, S. E. Ragan, J. S. Urquhart, R. S. Klessen, J. D. Soler, N. Roy, J. Kainulainen, T. Henning, F. Bigiel, R. J. Smith, F. Wyrowski, S. N. Longmore

Abstract: The hydroxyl radical (OH) is present in the diffuse molecular and partially atomic phases of the interstellar medium (ISM), but its abundance relative to hydrogen is not clear. We aim to evaluate the abundance of OH with respect to molecular hydrogen using OH absorption against cm-continuum sources over the first Galactic quadrant. This OH study is part of the HI/OH/Recombination line survey (THOR… ▽ More The hydroxyl radical (OH) is present in the diffuse molecular and partially atomic phases of the interstellar medium (ISM), but its abundance relative to hydrogen is not clear. We aim to evaluate the abundance of OH with respect to molecular hydrogen using OH absorption against cm-continuum sources over the first Galactic quadrant. This OH study is part of the HI/OH/Recombination line survey (THOR). THOR is a Karl G. Jansky Very Large Array large program of atomic, molecular and ionized gas in the range 15°$\leq$l$\leq$67° and |b|$\leq$1°. It is the highest-resolution unbiased OH absorption survey to date towards this region. We combine the derived optical depths with literature 13CO(1-0) and HI observations to determine the OH abundance. We detect absorption in the 1665 and 1667 MHz transitions for continuum sources stronger than $F_{\rm cont}\geq$0.1 Jy/beam. OH absorption is found against $\sim$15% of these continuum sources with increasing fractions for stronger sources. Most of the absorption is associated with Galactic HII regions. We find OH and 13CO gas to have similar kinematic properties. The OH abundance decreases with increasing hydrogen column density. The OH abundance with respect to the total hydrogen nuclei column density (atomic and molecular phase) is in agreement with a constant abundance for $A_V$ < 10-20. Towards the lowest column densities, we find sources that exhibit OH absorption but no 13CO emission, indicating that OH is a well suited tracer of the low column density molecular gas. We present spatially resolved OH absorption towards W43. The unbiased nature of the THOR survey opens a new window onto the gas properties of the ISM. The characterization of the OH abundance over a large range of hydrogen gas column densities contributes to the understanding of OH as a molecular gas tracer and provides a starting point for future investigations. △ Less

Submitted 16 September, 2018; v1 submitted 13 March, 2018; originally announced March 2018.

Comments: accepted for publication in Astronomy & Astrophysics; corrections before publication

Journal ref: A&A 618, A159 (2018)

Helical Magnetic Fields in Molecular Clouds? A New Method to Determine the Line-of-Sight Magnetic Field Structure in Molecular Clouds

Authors: M. Tahani, R. Plume, J. C. Brown, J. Kainulainen

Abstract: Magnetic fields pervade in the interstellar medium (ISM) and are believed to be important in the process of star formation, yet probing magnetic fields in star formation regions is challenging. We propose a new method to use Faraday rotation measurements in small scale star forming regions to find the direction and magnitude of the component of magnetic field along the line-of-sight. We test the p… ▽ More Magnetic fields pervade in the interstellar medium (ISM) and are believed to be important in the process of star formation, yet probing magnetic fields in star formation regions is challenging. We propose a new method to use Faraday rotation measurements in small scale star forming regions to find the direction and magnitude of the component of magnetic field along the line-of-sight. We test the proposed method in four relatively nearby regions of Orion A, Orion B, Perseus, and California. We use rotation measure data from the literature. We adopt a simple approach based on relative measurements to estimate the rotation measure due to the molecular clouds over the Galactic contribution. We then use a chemical evolution code along with extinction maps of each cloud to find the electron column density of the molecular cloud at the position of each rotation measure data point. Combining the rotation measures produced by the molecular clouds and the electron column density, we calculate the line-of-sight magnetic field strength and direction. In California and Orion A, we find clear evidence that the magnetic fields at one side of these filamentary structures are pointing towards us and are pointing away from us at the other side. Even though the magnetic fields in Perseus might seem to suggest the same behavior, not enough data points are available to draw such conclusions. In Orion B, as well, there are not enough data points available to detect such behavior. This behavior is consistent with a helical magnetic field morphology. In the vicinity of available Zeeman measurements in OMC-1, OMC-B, and the dark cloud Barnard 1, we find magnetic field values of $-23\pm38~μ$G, $-129\pm28~μ$G, and $32\pm101~μ$G, respectively, which are in agreement with the Zeeman Measurements. △ Less

Submitted 21 February, 2018; originally announced February 2018.

Journal ref: A&A 614, A100 (2018)

Relationship between turbulence energy and density variance in the Solar neighbourhood molecular clouds

Authors: Jouni Kainulainen, Christoph Federrath

Abstract: [abridged] We determine the relationship between turbulence energy and gas density variance for 15 molecular clouds in the Solar neighbourhood. We use the linewidths of the CO molecule as the probe of the turbulence energy (sonic Mach number, $\cal{M}_\mathrm{s}$) and three-dimensional models to reconstruct the density probability distribution function ($ρ$-PDF) of the clouds, derived using near-i… ▽ More [abridged] We determine the relationship between turbulence energy and gas density variance for 15 molecular clouds in the Solar neighbourhood. We use the linewidths of the CO molecule as the probe of the turbulence energy (sonic Mach number, $\cal{M}_\mathrm{s}$) and three-dimensional models to reconstruct the density probability distribution function ($ρ$-PDF) of the clouds, derived using near-infrared extinction and Herschel dust emission data, as the probe of the density variance ($σ_\mathrm{s}$). We find no significant correlation between $\cal{M}_\mathrm{s}$ and $σ_\mathrm{s}$ among the studied clouds, however, we also cannot rule out a weak correlation. In the context of turbulence-dominated gas, the range of the $\cal{M}_\mathrm{s}$ and $σ_\mathrm{s}$ values corresponds with the model predictions. The data cannot constrain whether or not the turbulence driving parameter, $b$, and/or thermal-to-magnetic pressure ratio, $β$, vary among the sample clouds. Most clouds are not in agreement with field strengths stronger than given by $β\lesssim 0.05$. A model with $b^2 β/ (β+1) = 0.30 \pm 0.06$ provides an adequate fit to the cloud sample as a whole. When considering the average behaviour of the sample, we can rule out three regimes: (i) strong compression combined with a weak magnetic field ($b \gtrsim 0.7$ and $β\gtrsim 3$), (ii) weak compression ($b \lesssim 0.35$), and (iii) strong magnetic field ($β\lesssim 0.1$). Including independent magnetic field strength estimates to the analysis, the data rule out solenoidal driving ($b < 0.4$) for the majority of the Solar neighbourhood clouds. However, most clouds have $b$ parameters larger than unity, which indicates a discrepancy with the turbulence-dominated picture; we discuss the possible reasons for this. △ Less

Submitted 5 November, 2017; originally announced November 2017.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 608, L3 (2017)

On the fragmentation of filaments in a molecular cloud simulation

Authors: R. -A. Chira, J. Kainulainen, J. C. Ibàñez-Mejía, Th. Henning, M. -M. Mac Low

Abstract: The fragmentation of filaments in molecular clouds has attracted a lot of attention as there seems to be a relation between the evolution of filaments and star formation. The study of the fragmentation process has been motivated by simple analytical models. However, only a few comprehensive studies have analysed the evolution of filaments using numerical simulations where the filaments form self-c… ▽ More The fragmentation of filaments in molecular clouds has attracted a lot of attention as there seems to be a relation between the evolution of filaments and star formation. The study of the fragmentation process has been motivated by simple analytical models. However, only a few comprehensive studies have analysed the evolution of filaments using numerical simulations where the filaments form self-consistently as part of molecular clouds. We address the early evolution of pc-scale filaments that form within individual clouds. We focus on three questions: How do the line masses of filaments evolve? How and when do the filaments fragment? How does the fragmentation relate to the line masses of the filaments? We examine three simulated molecular clouds formed in kpc-scale numerical simulations performed with the FLASH code. We compare the properties of the identified filaments with the predictions of analytic filament stability models. The line masses and mass fraction enclosed in the identified filaments increase continuously after the onset of self-gravity. The first fragments appear early when the line masses lie well below the critical line mass of Ostriker's hydrostatic equilibrium solution. The average line masses of filaments identified in 3D density cubes increases far more quickly than those identified in 2D column density maps. Our results suggest that hydrostatic or dynamic compression from the surrounding cloud has a significant impact on the early dynamical evolution of filaments. A simple model of an isolated, isothermal cylinder may not provide a good approach for fragmentation analysis. Caution must be exercised in interpreting distributions of properties of filaments identified in column density maps, especially in the case of low-mass filaments. Comparing or combining results from studies that use different filament finding techniques is strongly discouraged. △ Less

Submitted 28 November, 2017; v1 submitted 4 November, 2017; originally announced November 2017.

Comments: 16 pages, 17 figures, A&A, final version

Journal ref: A&A 610, A62 (2018)

The stellar content of the infalling molecular clump G286.21+0.17

Authors: M. Andersen, P. J. Barnes, J. C. Tan, J. Kainulainen, G. de Marchi

Abstract: The early evolution during massive star cluster formation is still uncertain. Observing embedded clusters at their earliest stages of formation can provide insight into the spatial and temporal distribution of the stars and thus probe different star cluster formation models. We present near-infrared imaging of an 8'*13'(5.4pc*8.7pc) region around the massive infalling clump G286.21+0.17(also known… ▽ More The early evolution during massive star cluster formation is still uncertain. Observing embedded clusters at their earliest stages of formation can provide insight into the spatial and temporal distribution of the stars and thus probe different star cluster formation models. We present near-infrared imaging of an 8'*13'(5.4pc*8.7pc) region around the massive infalling clump G286.21+0.17(also known as BYF73). The stellar content across the field is determined and photometry is derived in order to { obtain} stellar parameters for the cluster members. We find evidence for some sub-structure (on scales less than a pc diameter) within the region with apparently at least three different sub-clusters associated with the molecular clump based on differences in extinction and disk fractions. At the center of the clump we identify a deeply embedded sub-cluster. Near-infrared excess is detected for 39-44% in the two sub-clusters associated with molecular material and 27% for the exposed cluster. Using the disk excess as a proxy for age this suggests the clusters are very young. The current total stellar mass is estimated to be at least 200 Msun. The molecular core hosts a rich population of pre-main sequence stars. There is evidence for multiple events of star formation both in terms of the spatial distribution within the star forming region and possibly from the disk frequency. △ Less

Submitted 29 September, 2017; originally announced September 2017.

Comments: Submitted to AJ

Constraining the Dust Opacity Law in Three Small and Isolated Molecular Clouds

Authors: K. Webb, J. Di Francesco, S. Sadavoy, K. Thanjavur, R. Launhardt, Y. Shirley, A. Stutz, J. Abreu Vicente, J. Kainulainen

Abstract: Density profiles of isolated cores derived from thermal dust continuum emission rely on models of dust properties, such as mass opacity, which are poorly constrained. With complementary measures from near-infrared extinction maps, we can assess the reliability of commonly-used dust models. In this work, we compare Herschel-derived maps of the optical depth with equivalent maps derived from CFHT WI… ▽ More Density profiles of isolated cores derived from thermal dust continuum emission rely on models of dust properties, such as mass opacity, which are poorly constrained. With complementary measures from near-infrared extinction maps, we can assess the reliability of commonly-used dust models. In this work, we compare Herschel-derived maps of the optical depth with equivalent maps derived from CFHT WIRCAM near-infrared observations for three isolated cores: CB68, L429, and L1552. We assess the dust opacities provided from four models: OH1a, OH5a, Orm1, and Orm4. Although the consistency of the models differs between the three sources, the results suggest that the optical properties of dust in the envelopes of the cores are best described by either silicate and bare graphite grains (e.g., Orm1) or carbonaceous grains with some coagulation and either thin or no ice mantles (e.g., OH5a). None of the models, however, individually produced the most consistent optical depth maps for every source. The results suggest that either the dust in the cores is not well described by any one dust property model, the application of the dust models cannot be extended beyond the very center of the cores, or more complex SED fitting functions are necessary. △ Less

Submitted 28 September, 2017; originally announced September 2017.

Comments: Accepted to ApJ, September 27, 2017

On the effective turbulence driving mode of molecular clouds formed in disc galaxies

Authors: Keitaro Jin, Diane M. Salim, Christoph Federrath, Elizabeth J. Tasker, Asao Habe, Jouni T. Kainulainen

Abstract: We determine the physical properties and turbulence driving mode of molecular clouds formed in numerical simulations of a Milky Way-type disc galaxy with parsec-scale resolution. The clouds form through gravitational fragmentation of the gas, leading to average values for mass, radii and velocity dispersion in good agreement with observations of Milky Way clouds. The driving parameter (b) for the… ▽ More We determine the physical properties and turbulence driving mode of molecular clouds formed in numerical simulations of a Milky Way-type disc galaxy with parsec-scale resolution. The clouds form through gravitational fragmentation of the gas, leading to average values for mass, radii and velocity dispersion in good agreement with observations of Milky Way clouds. The driving parameter (b) for the turbulence within each cloud is characterised by the ratio of the density contrast (sigma_rho) to the average Mach number (Mach) within the cloud, b = sigma_rho/Mach. As shown in previous works, b ~ 1/3 indicates solenoidal (divergence-free) driving and b ~ 1 indicates compressive (curl-free) driving. We find that the average b value of all the clouds formed in the simulations has a lower limit of b > 0.2. Importantly, we find that b has a broad distribution, covering values from purely solenoidal to purely compressive driving. Tracking the evolution of individual clouds reveals that the b value for each cloud does not vary significantly over their lifetime. Finally, we perform a resolution study with minimum cell sizes of 8, 4, 2 and 1 pc and find that the average b value increases with increasing resolution. Therefore, we conclude that our measured b values are strictly lower limits and that a resolution better than 1 pc is required for convergence. However, regardless of the resolution, we find that b varies by factors of a few in all cases, which means that the effective driving mode alters significantly from cloud to cloud. △ Less

Submitted 28 March, 2017; originally announced March 2017.

Comments: 12 pages, 11 figures, accepted for publication in MNRAS, more info: https://www.mso.anu.edu.au/~chfeder/pubs/turb_driv_gal/turb_driv_gal.html

The HI/OH/Recombination line survey of the inner Milky Way (THOR)

Authors: H. Beuther, S. Bihr, M. Rugel, K. Johnston, Y. Wang, F. Walter, A. Brunthaler, A. J. Walsh, J. Ott, J. Stil, Th. Henning, T. Schierhuber, J. Kainulainen, M. Heyer, P. F. Goldsmith, L. D. Anderson, S. N. Longmore, R. S. Klessen, S. C. O. Glover, J. S. Urquhart, R. Plume, S. E. Ragan, N. Schneider, N. M. McClure-Griffiths, K. M. Menten , et al. (5 additional authors not shown)

Abstract: Context: The past decade has witnessed a large number of Galactic plane surveys at angular resolutions below 20". However, no comparable high-resolution survey exists at long radio wavelengths around 21cm in line and continuum emission. Methods: Employing the Very Large Array (VLA) in the C-array configuration and a large program, we observe the HI 21cm line, four OH lines, nineteen Halpha radio r… ▽ More Context: The past decade has witnessed a large number of Galactic plane surveys at angular resolutions below 20". However, no comparable high-resolution survey exists at long radio wavelengths around 21cm in line and continuum emission. Methods: Employing the Very Large Array (VLA) in the C-array configuration and a large program, we observe the HI 21cm line, four OH lines, nineteen Halpha radio recombination lines as well as the continuum emission from 1 to 2GHz in full polarization over a large part of the first Galactic quadrant. Results: Covering Galactic longitudes from 14.5 to 67.4deg and latitudes between +-1.25deg, we image all of these lines and the continuum at ~20" resolution. These data allow us to study the various components of the interstellar medium (ISM): from the atomic phase, traced by the HI line, to the molecular phase, observed by the OH transitions, to the ionized medium, revealed by the cm continuum and the Halpha radio recombination lines. Furthermore, the polarized continuum emission enables magnetic field studies. In this overview paper, we discuss the survey outline and present the first data release as well as early results from the different datasets. We now release the first half of the survey; the second half will follow later after the ongoing data processing has been completed. The data in fits format (continuum images and line data cubes) can be accessed through the project web-page http://www.mpia.de/thor. Conclusions: The HI/OH/Recombination line survey of the Milky Way (THOR) opens a new window to the different parts of the ISM. It enables detailed studies of molecular cloud formation, conversion of atomic to molecular gas, and feedback from HII regions as well as the magnetic field in the Milky Way. It is highly complementary to other surveys of our Galaxy, and comparing different datasets allows us to address many open questions. △ Less

Submitted 12 September, 2016; originally announced September 2016.

Comments: Accepted for Astronomy & Astrophysics, 23 pages, 17 figures, a high-resolution version can be found at http://www.mpia.de/thor

Journal ref: 2016, A&A, 595, A32

Earliest phases of star formation (EPoS): Dust temperature distributions in isolated starless cores

Authors: N. Lippok, R. Launhardt, Th. Henning, Z. Balog. H. Beuther, J. Kainulainen, O. Krause, H. Linz, M. Nielbock, S. E. Ragan, T. P. Robitaille, S. I. Sadavoy, A. Schmiedeke

Abstract: Constraining the temperature and density structure of dense molecular cloud cores is fundamental for understanding the initial conditions of star formation. We use Herschel observations of the thermal FIR dust emission from nearby isolated molecular cloud cores and combine them with ground-based submillimeter continuum data to derive observational constraints on their temperature and density struc… ▽ More Constraining the temperature and density structure of dense molecular cloud cores is fundamental for understanding the initial conditions of star formation. We use Herschel observations of the thermal FIR dust emission from nearby isolated molecular cloud cores and combine them with ground-based submillimeter continuum data to derive observational constraints on their temperature and density structure. The aim of this study is to verify the validity of a ray-tracing inversion technique developed to derive the dust temperature and density structure of isolated starless cores directly from the dust emission maps and to test if the resulting temperature and density profiles are consistent with physical models. Using this ray-tracing inversion technique, we derive the dust temperature and density structure of six isolated starless cloud cores. We employ self-consistent radiative transfer modeling to the derived density profiles, treating the ISRF as the only heating source. The best-fit values of local strength of the ISRF and the extinction by the outer envelope are derived by comparing the self-consistently calculated temperature profiles with those derived by the ray-tracing method. We find that all starless cores are significantly colder inside than outside, with the core temperatures showing a strong negative correlation with peak column density. This suggests that their thermal structure is dominated by external heating from the ISRF and shielding by dusty envelopes. The temperature profiles derived with the ray-tracing inversion method can be well-reproduced with self-consistent radiative transfer models. △ Less

Submitted 14 June, 2016; originally announced June 2016.

Comments: 16 pages, 12 figures, accepted for publication in A&A

Journal ref: A&A 592, A61 (2016)

The Distribution of Mass Surface Densities in a High-Mass Protocluster

Authors: Wanggi Lim, Jonathan C. Tan, Jouni Kainulainen, Bo Ma, Michael J. Butler

Abstract: We study the probability distribution function (PDF) of mass surface densities, $Σ$, of infrared dark cloud (IRDC) G028.37+00.07 and its surrounding giant molecular cloud. This PDF constrains the physical processes, such as turbulence, magnetic fields and self-gravity, that are expected to be controlling cloud structure and star formation activity. The chosen IRDC is of particular interest since i… ▽ More We study the probability distribution function (PDF) of mass surface densities, $Σ$, of infrared dark cloud (IRDC) G028.37+00.07 and its surrounding giant molecular cloud. This PDF constrains the physical processes, such as turbulence, magnetic fields and self-gravity, that are expected to be controlling cloud structure and star formation activity. The chosen IRDC is of particular interest since it has almost 100,000 solar masses within a radius of 8 parsecs, making it one of the most massive, dense molecular structures known and is thus a potential site for the formation of a "super star cluster." We study $Σ$ in two ways. First, we use a combination of NIR and MIR extinction maps that are able to probe the bulk of the cloud structure up to $Σ\sim1\:{\rm{g\:cm}^{-2}}\:$($A_V\simeq200$ mag). Second, we study the FIR and sub-mm dust continuum emission from the cloud utilizing Herschel PACS and SPIRE images and paying careful attention to the effects of foreground and background contamination. We find that the PDFs from both methods, applied over a $\sim20^\prime$(30 pc)-sized region that contains $\simeq1.5\times10^5\:M_\odot$ and encloses a minimum closed contour with $Σ\simeq0.013\:{\rm{g\:cm}^{-2}}\:$($A_V\simeq3$ mag), shows a log-normal shape with the peak measured at $Σ\simeq0.021\:{\rm{g\:cm}^{-2}}\:$($A_V\simeq4.7$ mag). There is tentative evidence for the presence of a high-$Σ$ power law tail that contains from $\sim3\%$ to 8\% of the mass of the cloud material. We discuss the implications of these results for the physical processes occurring in this cloud. △ Less

Submitted 7 September, 2016; v1 submitted 30 May, 2016; originally announced May 2016.

Comments: Accepted for publication in the ApJL. 9 pages, 5 figures

Appearance of Dusty Filaments at Different Viewing Angles

Authors: R. -A. Chira, R. Siebenmorgen, Th. Henning, J. Kainulainen

Abstract: Context: In the last years, there have been many studies on the omnipresence and structures of filaments in star-forming regions, as well as their role in the process of star formation. Those filaments are normally identified as elongated fibres across the plane of the sky. But how would we detect filaments that are inclined? Aims: We aim to learn more about whether, and how, total column density… ▽ More Context: In the last years, there have been many studies on the omnipresence and structures of filaments in star-forming regions, as well as their role in the process of star formation. Those filaments are normally identified as elongated fibres across the plane of the sky. But how would we detect filaments that are inclined? Aims: We aim to learn more about whether, and how, total column density or dust temperature change with respect to the line of sight. Such variations would enable observers to use dust observations to identify and study filaments at any inclination and gain more insight on the distribution and orientations of filaments within the Galactic plane. Methods: As a first step, we perform numerical calculations on simple cylindrical models to evaluate the influence of filament geometry on the average flux density. After that, we apply our three-dimensional Monte Carlo dust radiative transfer code on two models of star-forming regions and derive maps of effective total column density and dust temperature at different viewing angles. Results: We see only slight changes of average flux density for all cylinders we study. For our more complex models, we find that the effective dust temperature is not sensitive to viewing angle, while the total column density is strongly influenced, with differences exceeding an order of magnitude. The variations are not injective with the viewing angle and depend on the structure of the object. Conclusions: We conclude that there is no single quantity in our analysis that can uniquely trace the inclination and three-dimensional structure of a filament based on dust observations alone. However, observing wide variations in total column density at a given effective dust temperature is indicative of inclined filaments. △ Less

Submitted 31 May, 2016; v1 submitted 10 May, 2016; originally announced May 2016.

Comments: 19 pages, 20 figures, accepted by A&A

Journal ref: A&A 592, A90 (2016)

Resolving the fragmentation of high line-mass filaments with ALMA: the integral shaped filament in Orion A

Authors: Jouni Kainulainen, Amelia M. Stutz, Thomas Stanke, Jorge Abreu-Vicente, Henrik Beuther, Thomas Henning, Katharine G. Johnston, Tom Megeath

Abstract: We study the fragmentation of the nearest high line-mass filament, the integral shaped filament (ISF, line-mass $\sim$ 400 M$_\odot$ pc$^{-1}$) in the Orion A molecular cloud. We have observed a 1.6 pc long section of the ISF with the Atacama Large Millimetre/submillimeter Array (ALMA) at 3 mm continuum emission, at a resolution of $\sim$3" (1 200 AU). We identify from the region 43 dense cores wi… ▽ More We study the fragmentation of the nearest high line-mass filament, the integral shaped filament (ISF, line-mass $\sim$ 400 M$_\odot$ pc$^{-1}$) in the Orion A molecular cloud. We have observed a 1.6 pc long section of the ISF with the Atacama Large Millimetre/submillimeter Array (ALMA) at 3 mm continuum emission, at a resolution of $\sim$3" (1 200 AU). We identify from the region 43 dense cores with masses about a solar mass. 60% of the ALMA cores are protostellar and 40\% are starless. The nearest neighbour separations of the cores do not show a preferred fragmentation scale; the frequency of short separations increases down to 1 200 AU. We apply a two-point correlation analysis on the dense core separations and show that the ALMA cores are significantly grouped at separations below $\sim$17 000 AU and strongly grouped below $\sim$6 000 AU. The protostellar and starless cores are grouped differently: only the starless cores group strongly below $\sim$6 000 AU. In addition, the spatial distribution of the cores indicates periodic grouping of the cores into groups of $\sim$30 000 AU in size, separated by $\sim$50 000 AU. The groups coincide with dust column density peaks detected by Herschel. These results show hierarchical, two-mode fragmentation in which the maternal filament periodically fragments into groups of dense cores. Critically, our results indicate that the fragmentation models for lower line-mass filaments ($\sim$ 16 M$_\odot$ pc$^{-1}$) fail to capture the observed properties of the ISF. We also find that the protostars identified with Spitzer and Herschel in the ISF are grouped at separations below $\sim$17 000 AU. In contrast, young stars with disks do not show significant grouping. This suggests that the grouping of dense cores is partially retained over the protostar lifetime, but not over the lifetime of stars with disks. △ Less

Submitted 26 January, 2017; v1 submitted 17 March, 2016; originally announced March 2016.

Comments: Accepted to A&A

Journal ref: A&A 600, A141 (2017)

Giant molecular filaments in the Milky Way II: The fourth Galactic quadrant

Authors: J. Abreu-Vicente, S. Ragan, J. Kainulainen, Th. Henning, H. Beuther, K. Johnston

Abstract: Filamentary structures are common morphological features of the cold, molecular interstellar medium (ISM). Recent studies have discovered massive, hundred-parsec-scale filaments that may be connected to the large-scale, Galactic spiral arm structure. Addressing the nature of these Giant Molecular Filaments (GMFs) requires a census of their occurrence and properties. We perform a systematic search… ▽ More Filamentary structures are common morphological features of the cold, molecular interstellar medium (ISM). Recent studies have discovered massive, hundred-parsec-scale filaments that may be connected to the large-scale, Galactic spiral arm structure. Addressing the nature of these Giant Molecular Filaments (GMFs) requires a census of their occurrence and properties. We perform a systematic search of GMFs in the fourth Galactic quadrant and determine their basic physical properties. We identify GMFs based on their dust extinction signatures in near- and mid-infrared and velocity structure probed by ^{13}CO line emission. We use the ^{13}CO line emission and ATLASGAL dust emission data to estimate the total and dense gas masses of the GMFs. We combine our sample with an earlier sample from literature and study the Galactic environment of the GMFs. We identify nine GMFs in the fourth Galactic quadrant; six are located in the Centaurus spiral arm and three in inter-arm regions. Combining this sample with an earlier study using the same identification criteria in the first Galactic quadrant results in 16 GMFs, nine of which are located within spiral arms. The GMFs have sizes of 80-160 pc and ^{13}CO-derived masses between 5-90 x 10^{4} Msun. Their dense gas mass fractions are between 1.5-37%, being higher in the GMFs connected to spiral arms. We also compare the different GMF-identification methods and find that emission and extinction based techniques overlap only partially, highlighting the need to use both to achieve a complete census. △ Less

Submitted 25 April, 2016; v1 submitted 17 March, 2016; originally announced March 2016.

Comments: Accepted for publication in A&A

Journal ref: A&A 590, A131 (2016)

The Optical-Infrared Extinction Curve and its Variation in the Milky Way

Authors: E. F. Schlafly, A. M. Meisner, A. M. Stutz, J. Kainulainen, J. E. G. Peek, K. Tchernyshyov, H. -W. Rix, D. P. Finkbeiner, K. R. Covey, G. M. Green, E. F. Bell, W. S. Burgett, K. C. Chambers, P. W. Draper, H. Flewelling, K. W. Hodapp, N. Kaiser, E. A. Magnier, N. F. Martin, N. Metcalfe, R. J. Wainscoat, C. Waters

Abstract: The dust extinction curve is a critical component of many observational programs and an important diagnostic of the physics of the interstellar medium. Here we present new measurements of the dust extinction curve and its variation towards tens of thousands of stars, a hundred-fold larger sample than in existing detailed studies. We use data from the APOGEE spectroscopic survey in combination with… ▽ More The dust extinction curve is a critical component of many observational programs and an important diagnostic of the physics of the interstellar medium. Here we present new measurements of the dust extinction curve and its variation towards tens of thousands of stars, a hundred-fold larger sample than in existing detailed studies. We use data from the APOGEE spectroscopic survey in combination with ten-band photometry from Pan-STARRS1, 2MASS, and WISE. We find that the extinction curve in the optical through infrared is well characterized by a one-parameter family of curves described by R(V). The extinction curve is more uniform than suggested in past works, with sigma(R(V)) = 0.18, and with less than one percent of sight lines having R(V) > 4. Our data and analysis have revealed two new aspects of Galactic extinction: first, we find significant, wide-area variations in R(V) throughout the Galactic plane. These variations are on scales much larger than individual molecular clouds, indicating that R(V) variations must trace much more than just grain growth in dense molecular environments. Indeed, we find no correlation between R(V) and dust column density up to E(B-V) ~ 2. Second, we discover a strong relationship between R(V) and the far-infrared dust emissivity. △ Less

Submitted 11 February, 2016; originally announced February 2016.

Comments: 26 pages, 20 figures, comments welcome

VISION - Vienna survey in Orion I. VISTA Orion A Survey

Authors: Stefan Meingast, João Alves, Diego Mardones, Paula Teixeira, Marco Lombardi, Josefa Großschedl, Joana Ascenso, Herve Bouy, Jan Forbrich, Alyssa Goodman, Alvaro Hacar, Birgit Hasenberger, Jouni Kainulainen, Karolina Kubiak, Charles Lada, Elizabeth Lada, André Moitinho, Monika Petr-Gotzens, Lara Rodrigues, Carlos G. Román-Zúñiga

Abstract: 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… ▽ More 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. △ Less

Submitted 7 January, 2016; originally announced January 2016.

Comments: Accepted for publication in A&A, high-quality version available at http://homepage.univie.ac.at/stefan.meingast/vision.pdf

Journal ref: A&A 587, A153 (2016)

The Musca cloud: A 6 pc-long velocity-coherent, sonic filament

Authors: A. Hacar, J. Kainulainen, M. Tafalla, H. Beuther, J. Alves

Abstract: Filaments play a central role in the molecular clouds' evolution, but their internal dynamical properties remain poorly characterized. To further explore the physical state of these structures, we have investigated the kinematic properties of the Musca cloud. We have sampled the main axis of this filamentary cloud in $^{13}$CO and C$^{18}$O (2--1) lines using APEX observations. The different line… ▽ More Filaments play a central role in the molecular clouds' evolution, but their internal dynamical properties remain poorly characterized. To further explore the physical state of these structures, we have investigated the kinematic properties of the Musca cloud. We have sampled the main axis of this filamentary cloud in $^{13}$CO and C$^{18}$O (2--1) lines using APEX observations. The different line profiles in Musca shows that this cloud presents a continuous and quiescent velocity field along its $\sim$6.5 pc of length. With an internal gas kinematics dominated by thermal motions (i.e., $σ_{NT}/c_s\lesssim1$) and large-scale velocity gradients, these results reveal Musca as the longest velocity-coherent, sonic-like object identified so far in the ISM. The transonic properties of Musca present a clear departure from the predicted supersonic velocity dispersions expected in the Larson's velocity dispersion-size relationship, and constitute the first observational evidence of a filament fully decoupled from the turbulent regime over multi-parsec scales. △ Less

Submitted 8 December, 2015; v1 submitted 19 November, 2015; originally announced November 2015.

Comments: 12 pages, 6 figures; Accepted for publication in A&A

Journal ref: A&A 587, A97 (2016)

Filament Fragmentation in High-Mass Star Formation

Authors: H. Beuther, S. E. Ragan, K. Johnston, Th. Henning, A. Hacar, J. T. Kainulainen

Abstract: Aims: We resolve the length-scales for filament formation and fragmentation (res. <=0.1pc), in particular the Jeans length and cylinder fragmentation scale. Methods: We observed the prototypical high-mass star-forming filament IRDC18223 with the Plateau de Bure Interferometer (PdBI) in the 3.2mm continuum and N2H+(1-0) line emission in a ten field mosaic at a spatial resolution of ~4'' (~14000AU… ▽ More Aims: We resolve the length-scales for filament formation and fragmentation (res. <=0.1pc), in particular the Jeans length and cylinder fragmentation scale. Methods: We observed the prototypical high-mass star-forming filament IRDC18223 with the Plateau de Bure Interferometer (PdBI) in the 3.2mm continuum and N2H+(1-0) line emission in a ten field mosaic at a spatial resolution of ~4'' (~14000AU). Results: The dust continuum emission resolves the filament into a chain of at least 12 relatively regularly spaced cores. The mean separation between cores is ~0.40(+-0.18)pc. While this is approximately consistent with the fragmentation of an infinite, isothermal, gravitationally bound gas cylinder, a high mass-to-length ratio of M/l~1000M_sun/pc requires additional turbulent and/or magnetic support against radial collapse of the filament. The N2H+(1-0) data reveal a velocity gradient perpendicular to the main filament. Although rotation of the filament cannot be excluded, the data are also consistent with the main filament being comprised of several velocity-coherent sub-filaments. Furthermore, this velocity gradient perpendicular to the filament resembles recent results toward Serpens south that are interpreted as signatures of filament formation within magnetized and turbulent sheet-like structures. Lower-density gas tracers ([CI] and C18O) reveal a similar red/blueshifted velocity structure on scales around 60'' east and west of the IRDC18223 filament. This may tentatively be interpreted as a signature of the large-scale cloud and the smaller-scale filament being kinematically coupled. We do not identify a velocity gradient along the axis of the filament. This may either be due to no significant gas flows along the filamentary axis, but it may partly also be caused by a low inclination angle of the filament with respect to the plane of the sky that could minimize such signature. △ Less

Submitted 9 October, 2015; originally announced October 2015.

Comments: 12 pages, 9 figures, accepted for A&A, a higher-resolution version can be found at http://www.mpia.de/homes/beuther/papers.html

Journal ref: A&A 584, A67 (2015)

Gravitational fragmentation caught in the act: the filamentary Musca molecular cloud

Authors: Jouni Kainulainen, Alvaro Hacar, João Alves, Henrik Beuther, Hervé Bouy, Mario Tafalla

Abstract: Filamentary structures are common in molecular clouds. Explaining how they fragment to dense cores is a missing step in understanding their role in star formation. We perform a case study of whether low-mass filaments are close-to hydrostatic prior to their fragmentation, and whether their fragmentation agrees with gravitational fragmentation models. For this, we study the 6.5 pc long Musca molecu… ▽ More Filamentary structures are common in molecular clouds. Explaining how they fragment to dense cores is a missing step in understanding their role in star formation. We perform a case study of whether low-mass filaments are close-to hydrostatic prior to their fragmentation, and whether their fragmentation agrees with gravitational fragmentation models. For this, we study the 6.5 pc long Musca molecular cloud that is an ideal candidate for a filament at an early stage of fragmentation. We employ dust extinction mapping in conjunction with near-infrared data from the NEWFIRM instrument, and 870 um dust continuum emission data from the LABOCA instrument, to estimate column densities. We use the data to identify fragments from the cloud and to determine the radial density distribution of its filamentary part. We compare the cloud's morphology with 13CO and C18O line emission observed with the APEX/SHeFI instrument. The Musca cloud is pronouncedly fragmented at its ends, but harbours a remarkably well-defined, 1.6 pc long filament in its Center region. The line mass of the filament is 21-31 Ms pc^-1 and FWHM 0.07 pc. Its radial profile can be fitted with a Plummer profile that has the power-index of 2.6 \pm 11%, flatter than that of an infinite hydrostatic filament. The profile can also be fitted with a hydrostatic cylinder truncated by external pressure. These models imply a central density of 5-10 x 10^4 cm^-3. The fragments in the cloud have a mean separation of 0.4 pc, in agreement with gravitational fragmentation. These properties, together with the subsonic and velocity-coherent nature of the cloud, suggest a scenario in which an initially hydrostatic cloud is currently gravitationally fragmenting. The fragmentation has started a few tenths of a Myr ago from the cloud ends, leaving its center yet relatively non-fragmented, possibly because of gravitational focusing in a finite geometry. △ Less

Submitted 14 July, 2015; originally announced July 2015.

Comments: 14 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 586, A27 (2016)

Relationship between the column density distribution and evolutionary class of molecular clouds as viewed by ATLASGAL

Authors: J. Abreu-Vicente, J. Kainulainen, A. Stutz, T. Henning, H. Beuther

Abstract: We present the first study of the relationship between the column density distribution of molecular clouds within nearby Galactic spiral arms and their evolutionary status as measured from their stellar content. We analyze a sample of 195 molecular clouds located at distances below 5.5 kpc, identified from the ATLASGAL 870 micron data. We define three evolutionary classes within this sample: starl… ▽ More We present the first study of the relationship between the column density distribution of molecular clouds within nearby Galactic spiral arms and their evolutionary status as measured from their stellar content. We analyze a sample of 195 molecular clouds located at distances below 5.5 kpc, identified from the ATLASGAL 870 micron data. We define three evolutionary classes within this sample: starless clumps, star-forming clouds with associated young stellar objects, and clouds associated with HII regions. We find that the N(H2) probability density functions (N-PDFs) of these three classes of objects are clearly different: the N-PDFs of starless clumps are narrowest and close to log-normal in shape, while star-forming clouds and HII regions exhibit a power-law shape over a wide range of column densities and log-normal-like components only at low column densities. We use the N-PDFs to estimate the evolutionary time-scales of the three classes of objects based on a simple analytic model from literature. Finally, we show that the integral of the N-PDFs, the dense gas mass fraction, depends on the total mass of the regions as measured by ATLASGAL: more massive clouds contain greater relative amounts of dense gas across all evolutionary classes. △ Less

Submitted 2 July, 2015; originally announced July 2015.

Comments: Accepted for publication in A&A (25th June 15) 23 pages, 12 figures. Additional appendix figures will appear in the journal version of this paper

Journal ref: A&A 581, A74 (2015)

Evolution of column density distributions within Orion~A

Authors: A. M. Stutz, J. Kainulainen

Abstract: We compare the structure of star-forming molecular clouds in different regions of Orion A to determine how the column density probability distribution function (N-PDF) varies with environmental conditions such as the fraction of young protostars. A correlation between the N-PDF slope and Class 0 protostar fraction has been previously observed in a low-mass star-formation region (Perseus) by Sadavo… ▽ More We compare the structure of star-forming molecular clouds in different regions of Orion A to determine how the column density probability distribution function (N-PDF) varies with environmental conditions such as the fraction of young protostars. A correlation between the N-PDF slope and Class 0 protostar fraction has been previously observed in a low-mass star-formation region (Perseus) by Sadavoy; here we test if a similar correlation is observed in a high-mass star-forming region. We use Herschel data to derive a column density map of Orion A. We use the Herschel Orion Protostar Survey catalog for accurate identification and classification of the Orion A young stellar object (YSO) content, including the short-lived Class 0 protostars (with a $\sim$ 0.14 Myr lifetime). We divide Orion A into eight independent 13.5 pc$^2$ regions; in each region we fit the N-PDF distribution with a power-law, and we measure the fraction of Class 0 protostars. We use a maximum likelihood method to measure the N-PDF power-law index without binning. We find that the Class 0 fraction is higher in regions with flatter column density distributions. We test the effects of incompleteness, YSO misclassification, resolution, and pixel-scale. We show that these effects cannot account for the observed trend. Our observations demonstrate an association between the slope of the power-law N-PDF and the Class 0 fractions within Orion A. Various interpretations are discussed including timescales based on the Class 0 protostar fraction assuming a constant star-formation rate. The observed relation suggests that the N-PDF can be related to an "evolutionary state" of the gas. If universal, such a relation permits an evaluation of the evolutionary state from the N-PDF power-law index at much greater distances than those accesible with protostar counts. (abridged) △ Less

Submitted 20 April, 2015; originally announced April 2015.

Comments: A&A Letter, accepted

Journal ref: A&A 577, L6 (2015)

Unfolding the Laws of Star Formation: The Density Distribution of Molecular Clouds

Authors: Jouni Kainulainen, Christoph Federrath, Thomas Henning

Abstract: The formation of stars shapes the structure and evolution of entire galaxies. The rate and efficiency of this process are affected substantially by the density structure of the individual molecular clouds in which stars form. The most fundamental measure of this structure is the probability density function of volume densities (rho-PDF), which determines the star formation rates predicted with ana… ▽ More The formation of stars shapes the structure and evolution of entire galaxies. The rate and efficiency of this process are affected substantially by the density structure of the individual molecular clouds in which stars form. The most fundamental measure of this structure is the probability density function of volume densities (rho-PDF), which determines the star formation rates predicted with analytical models. This function has remained unconstrained by observations. We have developed an approach to quantify rho-PDFs and establish their relation to star formation. The rho-PDFs instigate a density threshold of star formation and allow us to quantify the star formation efficiency above it. The rho-PDFs provide new constraints for star formation theories and correctly predict several key properties of the star-forming interstellar medium. △ Less

Submitted 10 April, 2014; originally announced April 2014.

Comments: Science (2014) vol. 345. 34 pages