Showing 1–36 of 36 results for author: Möller, T

High-mass star formation across the Large Magellanic Cloud I. Chemical properties and hot molecular cores observed with ALMA at 1.2 mm

Authors: Roya Hamedani Golshan, Álvaro Sánchez-Monge, Peter Schilke, Marta Sewiło, Thomas Möller, V. S. Veena, Gary A. Fuller

Abstract: To study the impact of the initial effects of metallicity (i.e., the abundance of elements heavier than helium) on star formation and the formation of different molecular species, we searched for hot molecular cores in the sub-solar metallicity environment of the Large Magellanic Cloud (LMC). We conducted an ALMA Band 6 observations of 20 fields centered on young stellar objects (YSOs) distributed… ▽ More To study the impact of the initial effects of metallicity (i.e., the abundance of elements heavier than helium) on star formation and the formation of different molecular species, we searched for hot molecular cores in the sub-solar metallicity environment of the Large Magellanic Cloud (LMC). We conducted an ALMA Band 6 observations of 20 fields centered on young stellar objects (YSOs) distributed over the LMC in order to search for hot molecular cores in this galaxy. We detected a total of 65 compact 1.2 mm continuum cores in the 20 ALMA fields and analyzed their spectra with XCLASS software. The main temperature tracers are CH3OH and SO2, with more than two transitions detected in the observed frequency ranges. Other molecular lines with high detection rates in our sample are CS , SO, H13CO+, H13CN, HC15 N, and SiO. More complex molecules, such as HNCO, HDCO, HC3N, CH3CN, and NH2CHO, and multiple transitions of SO and SO2 isotopologues showed tentative or definite detection toward a small subset of the cores. According to the chemical richness of the cores and high temperatures from the XCLASS fitting, we report the detection of four hot cores and one hot core candidate. With one new hot core detection in this study, the number of detected hot cores in the LMC increases to seven. Six out of seven hot cores detected in the LMC to date are located in the stellar bar region of this galaxy. These six hot cores show emission from complex organic molecules (COMs), such as CH3OH, CH3CN, CH3OCHO, and CH3OCH3. The only known hot core in the LMC with no detection of COMs is located outside the bar region. The metallicity in the LMC presents a shallow gradient increasing from outer regions toward the bar. We suggest that the formation of hot molecular cores containing COMs ensues from the new generation of stars forming in the more metal-rich environment of the LMC bar. △ Less

Submitted 2 May, 2024; originally announced May 2024.

PRODIGE -- Envelope to Disk with NOEMA III. The origin of complex organic molecule emission in SVS13A

Authors: T. -H. Hsieh, J. E. Pineda, D. M. Segura-Cox, P. Caselli, M. T. Valdivia-Mena, C. Gieser, M. J. Maureira, A. Lopez-Sepulcre, L. Bouscasse, R. Neri, Th. Möller, A. Dutrey, A. Fuente, D. Semenov, E. Chapillon, N. Cunningham, Th. Henning, V. Pietu, I. Jimenez-Serra, S. Marino, C. Ceccarelli

Abstract: Complex Organic Molecules (COMs) have been found toward low-mass protostars but the origins of the COM emission are still unclear. It can be associated with, for example, hot corinos, outflows, and/or accretion shock/disk atmosphere. We have conducted NOEMA observations toward SVS13A from the PROtostars & DIsks: Global Evolution (PRODIGE) program. Our previous \ce{DCN} observations reveal a possib… ▽ More Complex Organic Molecules (COMs) have been found toward low-mass protostars but the origins of the COM emission are still unclear. It can be associated with, for example, hot corinos, outflows, and/or accretion shock/disk atmosphere. We have conducted NOEMA observations toward SVS13A from the PROtostars & DIsks: Global Evolution (PRODIGE) program. Our previous \ce{DCN} observations reveal a possible infalling streamer, which may affect the chemistry of the central protobinary by inducing accretion outbursts and/or shocked gas. Here, we further analyze six O-bearing COMs: CH3OH, aGg'-(CH2OH)2, C2H5OH, CH2(OH)CHO, CH3CHO, and CH3OCHO. Although the COM emission is not spatially resolved, we constrain the source sizes to $\lesssim0.3-0.4$ arcsec (90$-$120 au) by conducting uv-domain Gaussian fitting. Interestingly, the high-spectral resolution data reveal complex line profiles with multiple peaks showing differences between these six O-bearing COMs. The LTE fitting unveils differences in excitation temperatures and emitting areas among these COMs. We further conduct multiple-velocity-component LTE fitting to decompose the line emission into different kinematic components. Up to 6 velocity components are found from the LTE modeling. The temperature, column density, and source size of these components from each COM are obtained. We find a variety in excitation temperatures ($100-500$ K) and source sizes (D$\sim10-70$ au) from these kinematic components from different COMs. The emission of each COM can trace several components and different COMs most likely trace different regions. Given this complex structure, we suggest that the central region is inhomogeneous and unlikely to be heated by only protostellar radiation. We conclude that accretion shocks induced by the large-scale infalling streamer likely exist and contribute to the complexity of the COM emission. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: 29 pages, 18 figures, accepted to A&A

Journal ref: A&A 686, A289 (2024)

PRODIGE -- Planet-forming disks in Taurus with NOEMA. I. Overview and first results for 12CO, 13CO, and C18O

Authors: D. Semenov, Th. Henning, S. Guilloteau, G. Smirnov-Pinchukov, A. Dutrey, E. Chapillon, V. Pietu, R. Franceschi, K. Schwarz, S. van Terwisga, L. Bouscasse, P. Caselli, C. Ceccarelli, N. Cunningham, A. Fuente, C. Gieser, T. -H. Hsieh, A. Lopez-Sepulcre, D. M. Segura-Cox, J. E. Pineda, M. J. Maureira, Th. Moeller, M. Tafalla, M. T. Valdivia-Mena

Abstract: We are performing a line survey of 8 planet-forming Class II disks in Taurus with the IRAM NOrthern Extended Millimeter Array (NOEMA), as a part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution; PIs: P. Caselli and Th. Henning). Compact and extended disks around T Tauri stars CI, CY, DG, DL, DM, DN, IQ Tau, and UZ Tau E are observed in ~80 lines from >20 C-, O,-… ▽ More We are performing a line survey of 8 planet-forming Class II disks in Taurus with the IRAM NOrthern Extended Millimeter Array (NOEMA), as a part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution; PIs: P. Caselli and Th. Henning). Compact and extended disks around T Tauri stars CI, CY, DG, DL, DM, DN, IQ Tau, and UZ Tau E are observed in ~80 lines from >20 C-, O,- N-, and S-bearing species. The observations in four spectral settings at 210-280 GHz with $1σ$ rms sensitivity of ~ 8-12 mJy/beam at 0.9" and 0.3 km/s resolution will be completed in 2024. The uv-visibilities are fitted with the DiskFit model to obtain key stellar and disk properties. In this paper, the combined $^{12}$CO, $^{13}$CO and C$^{18}$O $J = 2-1$ data are presented. We find that the CO fluxes and disk masses inferred from dust continuum tentatively correlate with the CO emission sizes. We constrain dynamical stellar masses, geometries, temperatures, the CO column densities and gas masses for each disk. The best-fit temperatures at 100 au are ~ 17-37 K, and decrease radially with the power-law exponent q ~ 0.05-0.76. The inferred CO column densities decrease radially with the power-law exponent p ~ 0.2-3.1. The gas masses estimated from $^{13}$CO (2-1) are ~ $0.001-0.2 M_\textrm{Sun}$. The best-fit CO column densities point to severe CO freeze-out in the disks. The DL Tau disk is an outlier, and has either stronger CO depletion or lower gas mass than the rest of the sample. The CO isotopologue ratios are roughly consistent with the observed values in disks and the low-mass star-forming regions. △ Less

Submitted 27 February, 2024; v1 submitted 22 February, 2024; originally announced February 2024.

Comments: 47 pages, 10 figures, 6 tables, accepted for publication by A&A on 22/02/2024

Journal ref: A&A 685, A126 (2024)

Density distributions, magnetic field structures and fragmentation in high-mass star formation

Authors: H. Beuther, C. Gieser, J. D. Soler, Q. Zhang, R. Rao, D. Semenov, Th. Henning, R. Pudritz, T. Peters, P. Klaassen, M. T. Beltran, A. Palau, T. Moeller, K. G. Johnston, H. Zinnecker, J. Urquhart, R. Kuiper, A. Ahmadi, A. Sanchez-Monge, S. Feng, S. Leurini, S. E. Ragan

Abstract: Methods: Observing the large pc-scale Stokes I mm dust continuum emission with the IRAM 30m telescope and the intermediate-scale (<0.1pc) polarized submm dust emission with the Submillimeter Array toward a sample of 20 high-mass star-forming regions allows us to quantify the dependence of the fragmentation behaviour of these regions depending on the density and magnetic field structures. Results… ▽ More Methods: Observing the large pc-scale Stokes I mm dust continuum emission with the IRAM 30m telescope and the intermediate-scale (<0.1pc) polarized submm dust emission with the Submillimeter Array toward a sample of 20 high-mass star-forming regions allows us to quantify the dependence of the fragmentation behaviour of these regions depending on the density and magnetic field structures. Results: We infer density distributions n~r^{-p} of the regions with typical power-law slopes p around ~1.5. There is no obvious correlation between the power-law slopes of the density structures on larger clump scales (~1pc) and the number of fragments on smaller core scales (<0.1pc). Comparing the large-scale single-dish density profiles to those derived earlier from interferometric observations at smaller spatial scales, we find that the smaller-scale power-law slopes are steeper, typically around ~2.0. The flattening toward larger scales is consistent with the star-forming regions being embedded in larger cloud structures that do not decrease in density away from a particular core. Regarding the magnetic field, for several regions it appears aligned with filamentary structures leading toward the densest central cores. Furthermore, we find different polarization structures with some regions exhibiting central polarization holes whereas other regions show polarized emission also toward the central peak positions. Nevertheless, the polarized intensities are inversely related to the Stokes I intensities. We estimate magnetic field strengths between ~0.2 and ~4.5mG, and we find no clear correlation between magnetic field strength and the fragmentation level of the regions. Comparison of the turbulent to magnetic energies shows that they are of roughly equal importance in this sample. The mass-to-flux ratios range between ~2 and ~7, consistent with collapsing star-forming regions. △ Less

Submitted 20 November, 2023; originally announced November 2023.

Comments: Accepted for Astronomy & Astrophysics, 14 pages, 14 figures plus appendices, also download option at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

The Detection of Higher-Order Millimeter Hydrogen Recombination Lines in the Large Magellanic Cloud

Authors: Marta Sewiło, Kazuki Tokuda, Stan E. Kurtz, Steven B. Charnley, Thomas Möller, Jennifer Wiseman, C. -H. Rosie Chen, Remy Indebetouw, Álvaro Sánchez-Monge, Kei E. I. Tanaka, Peter Schilke, Toshikazu Onishi, Naoto Harada

Abstract: We report the first extragalactic detection of the higher-order millimeter hydrogen recombination lines ($Δn>2$). The $γ$-, $ε$-, and $η$-transitions have been detected toward the millimeter continuum source N105-1A in the star-forming region N105 in the Large Magellanic Cloud (LMC) with the Atacama Large Millimeter/submillimeter Array (ALMA). We use the H40$α$ line, the brightest of the detected… ▽ More We report the first extragalactic detection of the higher-order millimeter hydrogen recombination lines ($Δn>2$). The $γ$-, $ε$-, and $η$-transitions have been detected toward the millimeter continuum source N105-1A in the star-forming region N105 in the Large Magellanic Cloud (LMC) with the Atacama Large Millimeter/submillimeter Array (ALMA). We use the H40$α$ line, the brightest of the detected recombination lines (H40$α$, H36$β$, H50$β$, H41$γ$, H57$γ$, H49$ε$, H53$η$, and H54$η$), and/or the 3 mm free-free continuum emission to determine the physical parameters of N105-1A (the electron temperature, emission measure, electron density, and size) and study ionized gas kinematics. We compare the physical properties of N105-1A to a large sample of Galactic compact and ultracompact (UC) H II regions and conclude that N105-1A is similar to the most luminous ($L>10^5$ $L_{\odot}$) UC H II regions in the Galaxy. N105-1A is ionized by an O5.5 V star, it is deeply embedded in its natal molecular clump, and likely associated with a (proto)cluster. We incorporate high-resolution molecular line data including CS, SO, SO$_2$, and CH$_3$OH ($\sim$0.12 pc), and HCO$^{+}$ and CO ($\sim$0.087 pc) to explore the molecular environment of N105-1A. Based on the CO data, we find evidence for a cloud-cloud collision that likely triggered star formation in the region. We find no clear outflow signatures, but the presence of filaments and streamers indicates on-going accretion onto the clump hosting the UC H II region. Sulfur chemistry in N105-1A is consistent with the accretion shock model predictions. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Comments: 51 pages, 30 figures, 2 tables (including appendices); accepted for publication in The Astrophysical Journal (ApJ)

The physical and chemical structure of Sagittarius B2 VIIIa. Dust and ionized gas contributions to the full molecular line survey of 47 hot cores

Authors: T. Möller, P. Schilke, Á. Sánchez-Monge, A. Schmiedeke, F. Meng

Abstract: Sagittarius B2 (Sgr B2) is a giant molecular cloud complex in the central molecular zone of our Galaxy hosting several sites of high-mass star formation. The two main centers of activity are Sgr B2(M) and Sgr B2(N), which contain 27 and 20 continuum sources, respectively. Our analysis aims to be a comprehensive modeling of each core spectrum, where we take the complex interaction between molecular… ▽ More Sagittarius B2 (Sgr B2) is a giant molecular cloud complex in the central molecular zone of our Galaxy hosting several sites of high-mass star formation. The two main centers of activity are Sgr B2(M) and Sgr B2(N), which contain 27 and 20 continuum sources, respectively. Our analysis aims to be a comprehensive modeling of each core spectrum, where we take the complex interaction between molecular lines, dust attenuation, and free-free emission arising from HII regions into account. In this work, we determine the dust and, if HII regions are contained, the parameters of the free-free thermal emission of the ionized gas for each core, and derive a self-consistent description of the continuum levels of each core. Using the high sensitivity of ALMA, we characterize the physical and chemical structure of these continuum sources and gain better insight into the star formation process within the cores. We used ALMA to perform an unbiased spectral line survey of all 47 sources in ALMA band 6 with a frequency coverage from 211 GHz to 275 GHz. In order to model the free-free continuum contribution of a specific core, we fit the contained recombination lines (RRLs) to obtain the electron temperatures and the emission measures, where we use an extended XCLASS program to describe RRLs and free-free continuum simultaneously. In contrast to previous analyses, we derived the corresponding parameters here not only for each core, but also for their local surrounding envelope, and determined their physical properties. The distribution of RRLs we found in the core spectra closely fits the distribution of HII regions described in previous analyses. For the cores we determine average dust temperatures of around 236 K (Sgr B2(M)) and 225 K (Sgr B2(N)), while the electronic temperatures are located in a range between 3800 K and 23800 K. △ Less

Submitted 24 July, 2023; v1 submitted 12 July, 2023; originally announced July 2023.

Journal ref: A&A 676, A121 (2023)

Kinematics and stability of high-mass protostellar disk candidates at sub-arcsecond resolution -- Insights from the IRAM NOEMA large program CORE

Authors: Aida Ahmadi, H. Beuther, F. Bosco, C. Gieser, S. Suri, J. C. Mottram, R. Kuiper, Th. Henning, Á. Sánchez-Monge, H. Linz, R. E. Pudritz, D. Semenov, J. M. Winters, T. Möller, M. T. Beltrán, T. Csengeri, R. Galván-Madrid, K. G. Johnston, E. Keto, P. D. Klaassen, S. Leurini, S. N. Longmore, S. L. Lumsden, L. T. Maud, L. Moscadelli , et al. (6 additional authors not shown)

Abstract: The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in… ▽ More The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in the 1.37 mm wavelength regime in both continuum and line emission, reaching $\sim$0.4" resolution (800 au at 2 kpc). Using the dense gas tracer CH$_3$CN, we find velocity gradients across 13 cores perpendicular to the directions of bipolar molecular outflows, making them excellent disk candidates. Specific angular momentum ($j$) radial profiles are on average $\sim10^{-3}$ km /s pc and follow $j \propto r^{1.7}$, consistent with a poorly resolved rotating and infalling envelope/disk model. Fitting the velocity profiles with a Keplerian model, we find protostellar masses in the range of $\sim 10-25$ $M_\odot$. Modelling the level population of CH$_3$CN lines, we present temperature maps and find median gas temperatures in the range $70-210$ K. We create Toomre $Q$ maps to study the stability of the disks and find almost all (11 of 13) disk candidates to be prone to fragmentation due to gravitational instabilities at the scales probed by our observations. In particular, disks with masses greater than $\sim10-20\%$ of the mass of their host (proto)stars are Toomre unstable, and more luminous protostellar objects tend to have disks that are more massive and hence more prone to fragmentation. Our finings show that most disks around high-mass protostars are prone to disk fragmentation early in their formation due to their high disk to stellar mass ratio. This impacts the accretion evolution of high-mass protostars which will have significant implications for the formation of the most massive stars. △ Less

Submitted 3 May, 2023; v1 submitted 28 April, 2023; originally announced May 2023.

Comments: 27 pages, 12 figures, 6 appendices - accepted for publication in Astronomy and Astrophysics

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

The star formation history of the Sco-Cen association. Coherent star formation patterns in space and time

Authors: Sebastian Ratzenböck, Josefa E. Großschedl, João Alves, Núria Miret-Roig, Immanuel Bomze, John Forbes, Alyssa Goodman, Alvaro Hacar, Doug Lin, Stefan Meingast, Torsten Möller, Martin Piecka, Laura Posch, Alena Rottensteiner, Cameren Swiggum, Catherine Zucker

Abstract: We reconstruct the star formation history of the Sco-Cen OB association using a novel high-resolution age map of the region. We develop an approach to produce robust ages for Sco-Cen's recently identified 37 stellar clusters using the SigMA algorithm. The Sco-Cen star formation timeline reveals four periods of enhanced star formation activity, or bursts, remarkably separated by about 5 Myr. Of the… ▽ More We reconstruct the star formation history of the Sco-Cen OB association using a novel high-resolution age map of the region. We develop an approach to produce robust ages for Sco-Cen's recently identified 37 stellar clusters using the SigMA algorithm. The Sco-Cen star formation timeline reveals four periods of enhanced star formation activity, or bursts, remarkably separated by about 5 Myr. Of these, the second burst, which occurred about 15 million years ago, is by far the dominant, and most of Sco-Cen's stars and clusters were in place by the end of this burst. The formation of stars and clusters in Sco-Cen is correlated but not linearly, implying that more stars were formed per cluster during the peak of the star formation rate. Most of the clusters that are large enough to have supernova precursors were formed during the 15 Myr period. Star and cluster formation activity has been continuously declining since then. We have clear evidence that Sco-Cen formed from the inside out and contains 100-pc long chains of contiguous clusters exhibiting well-defined age gradients, from massive older clusters to smaller young clusters. These observables suggest an important role for feedback in forming about half of Sco-Cen stars, although follow-up work is needed to quantify this statement. Finally, we confirm that the Upper-Sco age controversy discussed in the literature during the last decades is solved: the nine clusters previously lumped together as Upper-Sco, a benchmark region for planet formation studies, exhibit a wide range of ages from 3 to 19 Myr. △ Less

Submitted 14 May, 2023; v1 submitted 15 February, 2023; originally announced February 2023.

Comments: 20 pages, 15 figures, 3 tables, submitted to A&A

Journal ref: A&A 678, A71 (2023)

HyGAL: Characterizing the Galactic ISM with observations of hydrides and other small molecules II. The absorption line survey with the IRAM 30 m telescope

Authors: W. -J. Kim, P. Schilke, D. A. Neufeld, A. M. Jacob, Á. Sánchez-Monge, D. Seifried, B. Godard, K. M. Menten, S. Walch, E. Falgarone, V. S. Veena, S. Bialy, T. Möller, F. Wyrowski

Abstract: As a complement to the HyGAL Stratospheric Observatory for Infrared Astronomy Legacy Program, we report the results of a ground-based absorption line survey of simple molecules in diffuse and translucent Galactic clouds. Using the Institut de Radioastronomie Millimétrique (IRAM) 30 m telescope, we surveyed molecular lines in the 2 mm and 3 mm wavelength ranges toward 15 millimeter continuum source… ▽ More As a complement to the HyGAL Stratospheric Observatory for Infrared Astronomy Legacy Program, we report the results of a ground-based absorption line survey of simple molecules in diffuse and translucent Galactic clouds. Using the Institut de Radioastronomie Millimétrique (IRAM) 30 m telescope, we surveyed molecular lines in the 2 mm and 3 mm wavelength ranges toward 15 millimeter continuum sources. These sources, which are all massive star-forming regions located mainly in the first and second quadrants of the Milky Way, form the subset of the HyGAL sample that can be observed by the IRAM 30 m telescope. We detected HCO$^+$ absorption lines toward 14 sightlines, toward which we identified 78 foreground cloud components, as well as lines from HCN, HNC, C$_2$H, and c-C$_3$H$_2$ toward most sightlines. In addition, CS and H$_2$S absorption lines are found toward at least half of the continuum sources. Static Meudon photodissociation region (PDR) isobaric models that consider ultraviolet-dominated chemistry were unable to reproduce the column densities of all seven molecular species by just a factor of a few, except for H$_2$S. The inclusion of other formation routes driven by turbulent dissipation could possibly explain the observed high column densities of these species in diffuse clouds. There is a tentative trend for H$_2$S and CS abundances relative to H$_2$ to be larger in diffuse clouds ($X$(H$_2$S) and $X$(CS) $\sim 10^{-8} - 10^{-7}$) than in translucent clouds ($X$(H$_2$S) and $X$(CS) $\sim 10^{-9} - 10^{-8}$) toward a small sample; however, a larger sample is required in order to confirm this trend. The derived H$_2$S column densities are higher than the values predicted from the isobaric PDR models, suggesting that chemical desorption of H$_2$S from sulfur-containing ice mantles may play a role in increasing the H$_2$S abundance. △ Less

Submitted 19 December, 2022; originally announced December 2022.

Comments: 41 pages, 25 figures, 10 tables, and 5 appendices. Accepted for publication in A&A

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

Significance Mode Analysis (SigMA) for hierarchical structures. An application to the Sco-Cen OB association

Authors: Sebastian Ratzenböck, Josefa E. Großschedl, Torsten Möller, João Alves, Immanuel Bomze, Stefan Meingast

Abstract: We present a new clustering method, Significance Mode Analysis (SigMA), to extract co-spatial and co-moving stellar populations from large-scale surveys such as ESA Gaia. The method studies the topological properties of the density field in the multidimensional phase space. We validate SigMA on simulated clusters and find that it outperforms competing methods, especially in cases where many cluste… ▽ More We present a new clustering method, Significance Mode Analysis (SigMA), to extract co-spatial and co-moving stellar populations from large-scale surveys such as ESA Gaia. The method studies the topological properties of the density field in the multidimensional phase space. We validate SigMA on simulated clusters and find that it outperforms competing methods, especially in cases where many clusters are closely spaced. We apply the new method to Gaia DR3 data of the closest OB association to Earth, Scorpio-Centaurus (Sco-Cen), and find more than 13,000 co-moving young objects, with about 19% of these having a sub-stellar mass. SigMA finds 37 co-moving clusters in Sco-Cen. These clusters are independently validated by their narrow HRD sequences and, to a certain extent, by their association with massive stars too bright for Gaia, hence unknown to SigMA. We compare our results with similar recent work and find that the SigMA algorithm recovers richer populations, is able to distinguish clusters with velocity differences down to about 0.5 km s$^{-1}$, and reaches cluster volume densities as low as 0.01 sources/pc$^3$. The 3D distribution of these 37 coeval clusters implies a larger extent and volume for the Sco-Cen OB association than typically assumed in the literature. Additionally, we find the association to be more actively star-forming and dynamically more complex than previously thought. We confirm that the star-forming molecular clouds in the Sco-Cen region, namely, Ophiuchus, L134/L183, Pipe Nebula, Corona Australis, Lupus, and Chamaeleon, are part of the Sco-Cen The application of SigMA to Sco-Cen demonstrates that advanced machine learning tools applied to the superb Gaia data allows to construct an accurate census of the young populations, to quantify their dynamics, and to reconstruct the recent star formation history of the local Milky Way. △ Less

Submitted 25 November, 2022; originally announced November 2022.

Comments: Submitted after first revision to A&A

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

PRODIGE -- Envelope to Disk with NOEMA II. Small-scale temperature structure and a streamer feeding the SVS13A protobinary using CH3CN and DCN

Authors: T. -H. Hsieh, D. M. Segura-Cox, J. E. Pineda, P. Caselli, L. Bouscasse, R. Neri, A. Lopez-Sepulcre, M. T. Valdivia-Mena, M. J. Maureira, Th. Henning, G. V. Smirnov-Pinchukov, D. Semenov, Th. Möller, N. Cunningham, A. Fuente, S. Marino, A. Dutrey, M. Tafalla, E. Chapillon, C. Ceccarelli, B. Zhao

Abstract: Aims. We present high sensitivity and high-spectral resolution NOEMA observations of the Class 0/I binary system SVS13A, composed of the low-mass protostars VLA4A and VLA4B with a separation of ~90 au. VLA4A is undergoing an accretion burst that enriches the chemistry of the surrounding gas. This gives us an excellent opportunity to probe the chemical and physical conditions as well as the accreti… ▽ More Aims. We present high sensitivity and high-spectral resolution NOEMA observations of the Class 0/I binary system SVS13A, composed of the low-mass protostars VLA4A and VLA4B with a separation of ~90 au. VLA4A is undergoing an accretion burst that enriches the chemistry of the surrounding gas. This gives us an excellent opportunity to probe the chemical and physical conditions as well as the accretion process. Methods. We observe the (12K-11K) lines of CH3CN and CH313CN, the DCN (3-2) line, and the C18O (2-1) line toward SVS13A using NOEMA. Results. We find complex line profiles at disk scales which cannot be explained by a single component or pure Keplerian motion. By adopting two velocity components to model the complex line profiles, we find that the temperatures and densities are significantly different between these two components. This suggests that the physical conditions of the emitting gas traced via CH3CN can change dramatically within the circumbinary disk. In addition, combining our observations of DCN (3-2) with previous ALMA high-angular-resolution observations, we find that the binary system (or VLA4A) might be fed by an infalling streamer from envelope scales (~700 au). If this is the case, this streamer contributes to the accretion of material onto the system with a rate of at least 1.4x10-6 Msun yr-1. Conclusions. We conclude that the CH3CN emission in SVS13A traces hot gas from a complex structure. This complexity might be affected by a streamer that is possibly infalling and funneling material into the central region. △ Less

Submitted 25 December, 2022; v1 submitted 9 November, 2022; originally announced November 2022.

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

Journal ref: A&A 669, A137 (2023)

The physical and chemical structure of Sagittarius B2 -- VI. UCHII regions in Sgr B2

Authors: Fanyi Meng, Álvaro Sánchez-Monge, Peter Schilke, Adam Ginsburg, Chris De Pree, Nazar Budaiev, Desmond Jeff, Anika Schmiedeke, Andreas Schwörer, V. S. Veena, Thomas Möller

Abstract: The giant molecular cloud Sagittarius B2 (hereafter SgrB2) is the most massive region with ongoing high-mass star formation in the Galaxy. Two ultra-compact HII (UCHII) regions were identified in SgrB2's central hot cores, SgrB2(M) and SgrB2(N). Our aim is to characterize the properties of the HII regions in the entire SgrB2 cloud. Comparing the HII regions and the dust cores, we aim to depict the… ▽ More The giant molecular cloud Sagittarius B2 (hereafter SgrB2) is the most massive region with ongoing high-mass star formation in the Galaxy. Two ultra-compact HII (UCHII) regions were identified in SgrB2's central hot cores, SgrB2(M) and SgrB2(N). Our aim is to characterize the properties of the HII regions in the entire SgrB2 cloud. Comparing the HII regions and the dust cores, we aim to depict the evolutionary stages of different parts of SgrB2. We use the Very Large Array in its A, CnB, and D configurations, and in the frequency band C (~6 GHz) to observe the whole SgrB2 complex. Using ancillary VLA data at 22.4 GHz and ALMA data at 96 GHz, we calculated the physical parameters of the UCHII regions and their dense gas environment. We identify 54 UCHII regions in the 6 GHz image, 39 of which are also detected at 22.4 GHz. Eight of the 54 UCHII regions are newly discovered. The UCHII regions have radii between $0.006 {\rm pc}$ and $0.04 {\rm pc}$, and have emission measure between $10^{6} {\rm pc\,cm^{-6}}$ and $10^{9} {\rm pc\,cm^{-6}}$. The UCHII regions are ionized by stars of types from B0.5 to O6. We found a typical gas density of $\sim10^6-10^9 {\rm cm^{-3}}$ around the UCHII regions. The pressure of the UCHII regions and the dense gas surrounding them are comparable. The expansion timescale of these UCHII regions is determined to be $\sim10^4-10^5 {\rm yr}$. The percentage of the dust cores that are associated with HII regions are 33%, 73%, 4%, and 1% for SgrB2(N), SgrB2(M), SgrB2(S), and SgrB2(DS), respectively. Two-thirds of the dust cores in SgrB2(DS) are associated with outflows. The electron densities of the UCHII regions we identified are in agreement with that of typical UCHII regions, while the radii are smaller than those of the typical UCHII regions. The dust cores in SgrB2(N) are more evolved than in SgrB2(DS) but younger than in SgrB2(M). △ Less

Submitted 18 August, 2022; v1 submitted 16 August, 2022; originally announced August 2022.

Comments: 17 pages, 15 figure, accepted to A&A

Journal ref: A&A 666, A31 (2022)

Clustered star formation at early evolutionary stages. Physical and chemical analysis of the young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953

Authors: C. Gieser, H. Beuther, D. Semenov, S. Suri, J. D. Soler, H. Linz, J. Syed, Th. Henning, S. Feng, T. Möller, A. Palau, J. M. Winters, M. T. Beltrán, R. Kuiper, L. Moscadelli, P. Klaassen, J. S. Urquhart, T. Peters, S. N. Longmore, Á. Sánchez-Monge, R. Galván-Madrid, R. E. Pudritz, K. G. Johnston

Abstract: We aim to characterize the physical and chemical properties of fragmented cores during the earliest evolutionary stages in the very young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953. NOEMA 1.3 mm data are used in combination with archival mid- and far-infrared observations to construct and fit the SEDs of individual fragmented cores. The radial density profiles are inferred from t… ▽ More We aim to characterize the physical and chemical properties of fragmented cores during the earliest evolutionary stages in the very young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953. NOEMA 1.3 mm data are used in combination with archival mid- and far-infrared observations to construct and fit the SEDs of individual fragmented cores. The radial density profiles are inferred from the 1.3 mm continuum visibility profiles and the radial temperature profiles are estimated from H2CO rotation temperature maps. Molecular column densities are derived with the line fitting tool XCLASS. The physical and chemical properties are combined by applying the physical-chemical model MUSCLE in order to constrain the chemical timescales of a few line-rich cores. The morphology and spatial correlations of the molecular emission are analyzed using the HOG method. The mid-infrared data show that both regions contain a cluster of young stellar objects. Bipolar molecular outflows are observed in the CO 2-1 transition toward the strong mm cores indicating protostellar activity. We find strong molecular emission of SO, SiO, H2CO, and CH3OH in locations which are not associated with the mm cores. These shocked knots can be either associated with the bipolar outflows or, in the case of ISOSS J23053+5953, with a colliding flow that creates a large shocked region between the mm cores. The mean chemical timescale of the cores is lower (20 000 yr) compared to that of the sources of the more evolved CORE sample (60 000 yr). With the HOG method, we find that the spatial emission of species tracing the extended emission and of shock-tracing molecules are well correlated within transitions of these groups. △ Less

Submitted 14 October, 2021; v1 submitted 5 October, 2021; originally announced October 2021.

Comments: 34 pages, 17 figures, accepted for publication in A&A

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

Disk fragmentation in high-mass star formation. High-resolution observations towards AFGL 2591-VLA 3

Authors: S. Suri, H. Beuther, C. Gieser, A. Ahmadi, Á. Sánchez-Monge, J. M. Winters, H. Linz, Th. Henning, M. T. Beltrán, F. Bosco, R. Cesaroni, T. Csengeri, S. Feng, M. G. Hoare, K. G. Johnston, P. Klaasen, R. Kuiper, S. Leurini, S. Longmore, S. Lumsden, L. Maud, L. Moscadelli, T. Möller, A. Palau, T. Peters , et al. (7 additional authors not shown)

Abstract: Increasing evidence suggests that, similar to their low-mass counterparts, high-mass stars form through a disk-mediated accretion process. At the same time, formation of high-mass stars still necessitates high accretion rates, and hence, high gas densities, which in turn can cause disks to become unstable against gravitational fragmentation. We study the kinematics and fragmentation of the disk ar… ▽ More Increasing evidence suggests that, similar to their low-mass counterparts, high-mass stars form through a disk-mediated accretion process. At the same time, formation of high-mass stars still necessitates high accretion rates, and hence, high gas densities, which in turn can cause disks to become unstable against gravitational fragmentation. We study the kinematics and fragmentation of the disk around the high-mass star forming region AFGL 2591-VLA 3 which was hypothesized to be fragmenting based on the observations that show multiple outflow directions. We use a new set of high-resolution (0.19 arcsec) IRAM/NOEMA observations at 843 micron towards VLA 3 which allow us to resolve its disk, characterize the fragmentation, and study its kinematics. In addition to the 843 micron continuum emission, our spectral setup targets warm dense gas and outflow tracers such as HCN, HC$_3$N and SO$_2$, as well as vibrationally excited HCN lines. The high resolution continuum and line emission maps reveal multiple fragments with subsolar masses within the inner 1000 AU of VLA 3. Furthermore, the velocity field of the inner disk observed at 843 micron shows a similar behavior to that of the larger scale velocity field studied in the CORE project at 1.37 mm. We present the first observational evidence for disk fragmentation towards AFGL 2591-VLA 3, a source that was thought to be a single high-mass core. While the fragments themselves are low-mass, the rotation of the disk is dominated by the protostar with a mass of 10.3$\pm 1.8~M_{\odot}$. These data also show that NOEMA Band 4 can obtain the highest currently achievable spatial resolution at (sub-)mm wavelengths in observations of strong northern sources. △ Less

Submitted 10 September, 2021; originally announced September 2021.

Comments: 16 pages, 11 figures, Accepted for publication in A&A

Journal ref: A&A 655, A84 (2021)

Herschel observations of extraordinary sources: Full Herschel/HIFI molecular line survey of Sagittarius B2(M)

Authors: T. Möller, P. Schilke, A. Schmiedeke, E. A. Bergin, D. C. Lis, Á. Sánchez-Monge, A. Schwörer, C. Comito

Abstract: We present a full analysis of a broadband spectral line survey of Sagittarius B2 (Main), one of the most chemically rich regions in the Galaxy located within the giant molecular cloud complex Sgr B2 in the Central Molecular Zone. Our goal is to derive the molecular abundances and temperatures of the high-mass star-forming region Sgr B2(M) and thus its physical and astrochemical conditions. Sgr B2(… ▽ More We present a full analysis of a broadband spectral line survey of Sagittarius B2 (Main), one of the most chemically rich regions in the Galaxy located within the giant molecular cloud complex Sgr B2 in the Central Molecular Zone. Our goal is to derive the molecular abundances and temperatures of the high-mass star-forming region Sgr B2(M) and thus its physical and astrochemical conditions. Sgr B2(M) was observed using the Heterodyne Instrument for the Far-Infrared (HIFI) on board the Herschel Space Observatory in a spectral line survey from 480 to 1907 GHz at a spectral resolution of 1.1 MHz, which provides one of the largest spectral coverages ever obtained toward this high-mass star-forming region in the submillimeter with high spectral resolution and includes frequencies > 1 THz unobservable from the ground. We model the molecular emission from the submillimeter to the far-IR using the XCLASS program. For each molecule, a quantitative description was determined taking all emission and absorption features of that species across the entire spectral range into account. Additionally, we derive velocity resolved ortho / para ratios for those molecules for which ortho and para resolved molecular parameters are available. Finally, the temperature and velocity distributions are analyzed and the derived abundances are compared with those obtained for Sgr B2(N) from a similar HIFI survey. A total of 92 isotopologues were identified, arising from 49 different molecules, ranging from free ions to complex organic compounds and originating from a variety of environments from the cold envelope to hot and dense gas within the cores. Sulfur dioxide, methanol, and water are the dominant contributors. For the ortho / para ratios we find deviations from the high temperature values between 13 and 27 %. In total 14 % of all lines remain unidentified. △ Less

Submitted 27 April, 2021; originally announced April 2021.

Comments: 67 pages, 102 figures, submitted to A&A

Journal ref: A&A 651, A9 (2021)

Fragmentation and kinematics in high-mass star formation: CORE-extension targeting two very young high-mass star-forming regions

Authors: H. Beuther, C. Gieser, S. Suri, H. Linz, P. Klaassen, D. Semenov, J. M. Winters, Th. Henning, J. D. Soler, J. S. Urquhart, J. Syed, S . Feng, T. Moeller, M. T. Beltran, A. Sanchez-Monge, S. N. Longmore, T. Peters, J. Ballesteros-Paredes, P. Schilke, L. Moscadelli, A. Palau, R. Cesaroni, S. Lumsden, R. Pudritz, F. Wyrowski , et al. (2 additional authors not shown)

Abstract: Context: The formation of high-mass star-forming regions from their parental gas cloud and the subsequent fragmentation processes lie at the heart of star formation research. Aims: We aim to study the dynamical and fragmentation properties at very early evolutionary stages of high-mass star formation. Methods: Employing the NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30m telescope, we… ▽ More Context: The formation of high-mass star-forming regions from their parental gas cloud and the subsequent fragmentation processes lie at the heart of star formation research. Aims: We aim to study the dynamical and fragmentation properties at very early evolutionary stages of high-mass star formation. Methods: Employing the NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30m telescope, we observed two young high-mass star-forming regions, ISOSS22478 and ISOSS23053, in the 1.3mm continuum and spectral line emission at a high angular resolution (~0.8''). Results: We resolved 29 cores that are mostly located along filament-like structures. Depending on the temperature assumption, these cores follow a mass-size relation of approximately M~r^2.0, corresponding to constant mean column densities. However, with different temperature assumptions, a steeper mass-size relation up to M~r^3.0, which would be more likely to correspond to constant mean volume densities, cannot be ruled out. The correlation of the core masses with their nearest neighbor separations is consistent with thermal Jeans fragmentation. We found hardly any core separations at the spatial resolution limit, indicating that the data resolve the large-scale fragmentation well. Although the kinematics of the two regions appear very different at first sight - multiple velocity components along filaments in ISOSS22478 versus a steep velocity gradient of more than 50km/s/pc in ISOSS23053 - the findings can be explained within the framework of a dynamical cloud collapse scenario. Conclusions: While our data are consistent with a dynamical cloud collapse scenario and subsequent thermal Jeans fragmentation, the importance of additional environmental properties, such as the magnetization of the gas or external shocks triggering converging gas flows, is nonetheless not as well constrained and would require future investigation. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: 16 pages, 16 figures, accepted for Astronomy and Astrophysics, a higher-resolution version can also be found at https://www.mpia.de/homes/beuther/papers.html

Journal ref: A&A 649, A113 (2021)

The physical and chemical structure of high-mass star-forming regions. Unraveling chemical complexity with the NOEMA large program "CORE"

Authors: C. Gieser, H. Beuther, D. Semenov, A. Ahmadi, S. Suri, T. Möller, M. T. Beltran, P. Klaassen, Q. Zhang, J. S. Urquhart, Th. Henning, S. Feng, R. Galván-Madrid, V. de Souza Magalhães, L. Moscadelli, S. Longmore, S. Leurini, R. Kuiper, T. Peters, K. M. Menten, T. Csengeri, G. Fuller, F. Wyrowski, S. Lumsden, Á. Sánchez-Monge , et al. (8 additional authors not shown)

Abstract: We use sub-arcsecond resolution ($\sim$0.4$''$) observations with NOEMA at 1.37 mm to study the dust emission and molecular gas of 18 high-mass star-forming regions. We combine the derived physical and chemical properties of individual cores in these regions to estimate their ages. The temperature structure of these regions are determined by fitting H2CO and CH3CN line emission. The density profil… ▽ More We use sub-arcsecond resolution ($\sim$0.4$''$) observations with NOEMA at 1.37 mm to study the dust emission and molecular gas of 18 high-mass star-forming regions. We combine the derived physical and chemical properties of individual cores in these regions to estimate their ages. The temperature structure of these regions are determined by fitting H2CO and CH3CN line emission. The density profiles are inferred from the 1.37 mm continuum visibilities. The column densities of 11 different species are determined by fitting the emission lines with XCLASS. Within the 18 observed regions, we identify 22 individual cores with associated 1.37 mm continuum emission and with a radially decreasing temperature profile. We find an average temperature power-law index of q = 0.4$\pm$0.1 and an average density power-law index of p = 2.0$\pm$0.2 on scales on the order of several 1 000 au. Comparing these results with values of p derived in the literature suggest that the density profiles remain unchanged from clump to core scales. The column densities relative to N(C18O) between pairs of dense gas tracers show tight correlations. We apply the physical-chemical model MUSCLE to the derived column densities of each core and find a mean chemical age of $\sim$60 000 yrs and an age spread of 20 000-100 000 yrs. With this paper we release all data products of the CORE project available at https://www.mpia.de/core. The CORE sample reveals well constrained density and temperature power-law distributions. Furthermore, we characterize a large variety in molecular richness that can be explained by an age spread confirmed by our physical-chemical modeling. The hot molecular cores show the most emission lines, but we also find evolved cores at an evolutionary stage, in which most molecules are destroyed and thus the spectra appear line-poor again. △ Less

Submitted 23 February, 2021; originally announced February 2021.

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

Journal ref: A&A 648, A66 (2021)

Multi-scale view of star formation in IRAS 21078+5211: From clump fragmentation to disk wind

Authors: L. Moscadelli, H. Beuther, A. Ahmadi, C. Gieser, F. Massi, R. Cesaroni, Á. Sánchez-Monge, F. Bacciotti, M. T. Beltrán, T. Csengeri, R. Galván-Madrid, Th. Henning, P. D. Klaassen, R. Kuiper, S. Leurini, S. N. Longmore, L. T. Maud, T. Möller, A. Palau, T. Peters, R. E. Pudritz, A. Sanna, D. Semenov, J. S. Urquhart, J. M. Winters , et al. (1 additional authors not shown)

Abstract: In the massive star-forming region IRAS 21078+5211, a highly fragmented cluster (0.1~pc in size) of molecular cores is observed, located at the density peak of an elongated (1~pc in size) molecular cloud. A small (1~km/s per 0.1~pc) LSR velocity (Vlsr) gradient is detected across the axis of the molecular cloud. Assuming we are observing a mass flow from the harboring cloud to the cluster, we deri… ▽ More In the massive star-forming region IRAS 21078+5211, a highly fragmented cluster (0.1~pc in size) of molecular cores is observed, located at the density peak of an elongated (1~pc in size) molecular cloud. A small (1~km/s per 0.1~pc) LSR velocity (Vlsr) gradient is detected across the axis of the molecular cloud. Assuming we are observing a mass flow from the harboring cloud to the cluster, we derive a mass infall rate of about 10^{-4}~M_{sun}~yr^{-1}. The most massive cores (labeled 1, 2, and 3) are found at the center of the cluster, and these are the only ones that present a signature of protostellar activity in terms of emission from high-excitation molecular lines or a molecular outflow. We reveal an extended (size about 0.1~pc), bipolar collimated molecular outflow emerging from core 1. We believe this is powered by a (previously discovered) compact (size <= 1000~au) radio jet, ejected by a YSO embedded in core 1 (named YSO-1), since the molecular outflow and the radio jet are almost parallel and have a comparable momentum rate. By means of high-excitation lines, we find a large (14~km/s over 500~au) Vlsr gradient at the position of YSO-1, oriented approximately perpendicular to the radio jet. Assuming this is an edge-on, rotating disk and fitting a Keplerian rotation pattern, we determine the YSO-1 mass to be 5.6+/-2.0~M_{sun}. The water masers (previously observed with VLBI) emerge within 100-300~au from YSO-1 and are unique tracers of the jet kinematics. Their three-dimensional (3D) velocity pattern reveals that the gas flows along, and rotates about, the jet axis. We show that the 3D maser velocities are fully consistent with the magneto-centrifugal disk-wind models predicting a cylindrical rotating jet. Under this hypothesis, we determine the jet radius to be about 16~au and the corresponding launching radius and terminal velocity to be about 2.2~au and 200~km/s, respectively. △ Less

Submitted 9 February, 2021; originally announced February 2021.

Comments: 25 pages, 16 figures

Journal ref: A&A 647, A114 (2021)

The $ρ$ Oph region revisited with Gaia EDR3

Authors: Natalie Grasser, Sebastian Ratzenböck, João Alves, Josefa Großschedl, Stefan Meingast, Catherine Zucker, Alvaro Hacar, Charles Lada, Alyssa Goodman, Marco Lombardi, John C. Forbes, Immanuel M. Bomze, Torsten Möller

Abstract: Context. Young and embedded stellar populations are important probes of the star formation process. Paradoxically, we have a better census of nearby embedded young populations than the slightly more evolved optically visible young populations. The high accuracy measurements and all-sky coverage of Gaia data are about to change this situation. Aims. This work aims to construct the most complete sam… ▽ More Context. Young and embedded stellar populations are important probes of the star formation process. Paradoxically, we have a better census of nearby embedded young populations than the slightly more evolved optically visible young populations. The high accuracy measurements and all-sky coverage of Gaia data are about to change this situation. Aims. This work aims to construct the most complete sample to date of YSOs in the $ρ$ Oph region. Methods. We compile a catalog of 1114 Ophiuchus YSOs from the literature and crossmatch it with the Gaia EDR3, Gaia-ESO and APOGEE-2 surveys. We apply a multivariate classification algorithm to this catalog to identify new, co-moving population candidates. Results. We find 191 new high-fidelity YSO candidates in the Gaia EDR3 catalog belonging to the $ρ$ Oph region. The new sources appear to be mainly Class III M-stars and substellar objects and are less extincted than the known members. We find 28 previously unknown sources with disks. The analysis of the proper motion distribution of the entire sample reveals a well-defined bimodality, implying two distinct populations sharing a similar 3D volume. The first population comprises young stars' clusters around the $ρ$ Ophiuchi star and the main Ophiuchus clouds (L1688, L1689, L1709). In contrast, the second population is older ($\sim$ 10 Myr), dispersed, has a distinct proper motion, and is possibly from the Upper Sco group. The two populations are moving away from each other at about 4.1 km/s, and will no longer overlap in about 4 Myr. Finally, we flag 17 sources in the literature as impostors, which are sources that exhibit large deviations from the average distance and proper motion properties of the $ρ$ Oph population. Our results show the importance of accurate 3D space and motion information for improved stellar population analysis. (Abridged) △ Less

Submitted 8 June, 2021; v1 submitted 28 January, 2021; originally announced January 2021.

Comments: Submitted to A&A on January 28th 2021. This is the second (revised) version of the paper. All comments welcome

Journal ref: A&A 652, A2 (2021)

The GUAPOS project: G31.41+0.31 Unbiased ALMA sPectral Observational Survey -- I. Isomers of C$_{2}$H$_{4}$O$_{2}$

Authors: C. Mininni, M. T. Beltrán, V. M. Rivilla, A. Sánchez-Monge, F. Fontani, T. Möller, R. Cesaroni, P. Schilke, S. Viti, I. Jiménez-Serra, L. Colzi, A. Lorenzani, L. Testi

Abstract: Understanding the degree of chemical complexity that can be reached in star-forming regions, together with the identification of precursors of the building blocks of life in the interstellar medium, is one of the goals of astrochemistry. Unbiased spectral surveys with large bandwidth and high spectral resolution are thus needed, to resolve line blending in chemically rich sources and identify comp… ▽ More Understanding the degree of chemical complexity that can be reached in star-forming regions, together with the identification of precursors of the building blocks of life in the interstellar medium, is one of the goals of astrochemistry. Unbiased spectral surveys with large bandwidth and high spectral resolution are thus needed, to resolve line blending in chemically rich sources and identify complex organic molecules. This kind of observations has been successfully carried out, mainly towards the Galactic Center, a region that shows peculiar environmental conditions. We present an unbiased spectral survey at 3mm of one of the most chemically rich hot molecular cores located outside the Galactic Center, in the high-mass star-forming region G31.41+0.31. In this first paper, we present the survey and discuss the detection of the 3 isomers of C$_{2}$H$_{4}$O$_{2}$: methyl formate, glycolaldehyde and acetic acid. Observations were carried out with ALMA and cover the entire Band 3 from 84 to 116 GHz with an angular resolution of $1.2^{''}$x$1.2^{''}$ and a spectral resolution of $\sim0.488$ MHz. The transitions of the 3 molecules have been analyzed with the software XCLASS. All three isomers were detected and methyl formate and acetic acid abundances in G31 are the highest detected up to now, if compared to sources in literature. The size of the emission varies among the three isomers with acetic acid showing the most compact emission while methyl formate the most extended. The comparison with chemical models suggests the necessity of grain-surface routes for the formation of methyl formate in G31, while for glycolaldehyde both scenarios could be feasible. Proposed grain-surface reaction for acetic acid is not able to reproduce the observed abundance in this work, while gas-phase scenario should be further tested due to large uncertainties. △ Less

Submitted 28 September, 2020; originally announced September 2020.

Comments: 49 pages, 29 figures, accepted for publication in A&A

Journal ref: A&A 644, A84 (2020)

Extended stellar systems in the solar neighborhood. IV. Meingast 1: the most massive stellar stream in the solar neighborhood

Authors: S. Ratzenböck, S. Meingast, J. Alves, T. Möller, I. Bomze

Abstract: Nearby stellar streams carry unique information on the dynamical evolution and disruption of stellar systems in the Galaxy, the mass distribution in the disk, and provide unique targets for planet formation and evolution studies. We revisit the stream discovered in Meingast et al (2019) to search for new members, using Gaia DR2 data and a machine learning approach. We use a bagging classifier of o… ▽ More Nearby stellar streams carry unique information on the dynamical evolution and disruption of stellar systems in the Galaxy, the mass distribution in the disk, and provide unique targets for planet formation and evolution studies. We revisit the stream discovered in Meingast et al (2019) to search for new members, using Gaia DR2 data and a machine learning approach. We use a bagging classifier of one-class Support Vector Machines to perform a search in positions and proper motions for new stream members. We use the variable prediction frequency resulting from the multitude of classifiers to estimate a stream membership criterion which we use to select high fidelity sources. We use the HR diagram and the Cartesian velocity distribution as test and validation tools. We find about 2000 stream members with high-fidelity, or about an order of magnitude more than previously known, unveiling the stream's population across the entire stellar mass spectrum, from B-stars to M-stars, including white dwarfs. We find that, apart from being slightly more metal-poor, the HRD of the stream is indistinguishable from that of the Pleiades cluster. For the mass range at which we are mostly complete, $\sim$0.2 M$_\odot$ $ < $ M $ < $ $\sim$4 M$_\odot$, we find a normal IMF, allowing us to estimate the total mass of stream to be about 2000 M$_\odot$, making this relatively young stream by far the most massive known. In addition, we identify several white dwarfs as potential stream members. The nearby Meingast 1 stream, due to its richness, age, and distance, is a new fundamental laboratory for star and planet formation and evolution studies for the poorly studied gravitationally unbound star-formation mode. We also demonstrate that One-Class Support Vector Machines can be effectively used to unveil the full stellar populations of nearby stellar systems with Gaia data. △ Less

Submitted 13 February, 2020; originally announced February 2020.

Comments: submitted to A&A on 28.1.2020

Journal ref: A&A 639, A64 (2020)

Chemical complexity in high-mass star formation: An observational and modeling case study of the AFGL 2591 VLA 3 hot core

Authors: C. Gieser, D. Semenov, H. Beuther, A. Ahmadi, J. C. Mottram, Th. Henning, M. Beltran, L. T. Maud, F. Bosco, S. Leurini, T. Peters, P. Klaassen, R. Kuiper, S. Feng, J. S. Urquhart, L. Moscadelli, T. Csengeri, S. Lumsden, J. M. Winters, S. Suri, Q. Zhang, R. Pudritz, A. Palau, K. M. Menten, R. Galvan-Madrid , et al. (8 additional authors not shown)

Abstract: We present a detailed observational and modeling study of the hot core VLA 3 in the high-mass star-forming region AFGL 2591, which is a target region of the NOrthern Extended Millimeter Array (NOEMA) large program CORE. Using NOEMA observations at 1.37 mm with an angular resolution of ~0."42 (1 400 au at 3.33 kpc), we derived the physical and chemical structure of the source. We modeled the observ… ▽ More We present a detailed observational and modeling study of the hot core VLA 3 in the high-mass star-forming region AFGL 2591, which is a target region of the NOrthern Extended Millimeter Array (NOEMA) large program CORE. Using NOEMA observations at 1.37 mm with an angular resolution of ~0."42 (1 400 au at 3.33 kpc), we derived the physical and chemical structure of the source. We modeled the observed molecular abundances with the chemical evolution code MUSCLE (MUlti Stage ChemicaL codE). Results. With the kinetic temperature tracers CH3CN and H2CO we observe a temperature distribution with a power-law index of q = 0.41+-0.08. Using the visibilities of the continuum emission we derive a density structure with a power-law index of p = 1.7+-0.1. The hot core spectra reveal high molecular abundances and a rich diversity in complex molecules. The majority of the molecules have an asymmetric spatial distribution around the forming protostar(s), which indicates a complex physical structure on scales < 1 400 au. Using MUSCLE, we are able to explain the observed molecular abundance of 10 out of 14 modeled species at an estimated hot core chemical age of ~21 100 years. In contrast to the observational analysis, our chemical modeling predicts a lower density power-law index of p < 1.4. Reasons for this discrepancy are discussed. Conclusions. Combining high spatial resolution observations with detailed chemical modeling allows us to derive a concise picture of the physical and chemical structure of the famous AFGL 2591 hot core. The next steps are to conduct a similar analysis for the whole CORE sample, and then use this analysis to constrain the chemical diversity in high-mass star formation to a much greater depth. △ Less

Submitted 11 October, 2019; originally announced October 2019.

Comments: 22 pages, 13 figures, accepted for publication in A&A

Journal ref: A&A 631, A142 (2019)

The physical and chemical structure of Sagittarius B2 -- V. Non-thermal emission in the envelope of Sgr B2

Authors: F. Meng, Á. Sánchez-Monge, P. Schilke, M. Padovani, A. Marcowith, A. Ginsburg, A. Schmiedeke, A. Schwörer, C. DePree, V. S. Veena, Th. Möller

Abstract: The giant molecular cloud Sagittarius B2 (hereafter SgrB2) is the most massive region with ongoing high-mass star formation in the Galaxy. In the southern region of the 40-pc large envelope of SgrB2, we encounter the SgrB2(DS) region which hosts more than 60 high-mass protostellar cores distributed in an arc shape around an extended HII region. We use the Very Large Array in its CnB and D configur… ▽ More The giant molecular cloud Sagittarius B2 (hereafter SgrB2) is the most massive region with ongoing high-mass star formation in the Galaxy. In the southern region of the 40-pc large envelope of SgrB2, we encounter the SgrB2(DS) region which hosts more than 60 high-mass protostellar cores distributed in an arc shape around an extended HII region. We use the Very Large Array in its CnB and D configurations, and in the frequency bands C (4--8 GHz) and X (8--12 GHz) to observe the whole SgrB2 complex. Continuum and radio recombination line maps are obtained. We detect radio continuum emission in SgrB2(DS) in a bubble-shaped structure. From 4 to 12 GHz, we derive a spectral index between -1.2 and -0.4, indicating the presence of non-thermal emission. We decompose the contribution from thermal and non-thermal emission, and find that the thermal component is clumpy and more concentrated, while the non-thermal component is more extended and diffuse. The radio recombination lines in the region are found to be not in local thermodynamic equilibrium (LTE) but stimulated by the non-thermal emission. The thermal free-free emission is likely tracing an HII region ionized by an O7 star, while the non-thermal emission can be generated by relativistic electrons created through first-order Fermi acceleration. We have developed a simple model of the SgrB2(DS) region and found that first-order Fermi acceleration can reproduce the observed flux density and spectral index. △ Less

Submitted 19 September, 2019; v1 submitted 20 August, 2019; originally announced August 2019.

Comments: Accepted for publication in Astronomy & Astrophysics. 15 pages with 12 figures and 1 table, plus Appendices with extra figures. (Abstract modified)

Journal ref: A&A 630, A73 (2019)

The Physical and chemical structure of Sagittarius B2 -- IV. Converging filaments in the high-mass cluster forming region Sgr B2(N)

Authors: A. Schwörer, Á. Sánchez-Monge, P. Schilke, T. Möller, A. Ginsburg, F. Meng, A. Schmiedeke, H. S. P. Müller, D. Lis, S. -L. Qin

Abstract: We have used an unbiased, spectral line-survey that covers the frequency range from 211 to 275 GHz and was obtained with ALMA (angular resolution of 0.4 arcsec) to study the small-scale structure of the dense gas in Sagittarius B2 (north). Eight filaments are found converging to the central hub and extending for about 0.1 pc. The spatial structure, together with the presence of the massive central… ▽ More We have used an unbiased, spectral line-survey that covers the frequency range from 211 to 275 GHz and was obtained with ALMA (angular resolution of 0.4 arcsec) to study the small-scale structure of the dense gas in Sagittarius B2 (north). Eight filaments are found converging to the central hub and extending for about 0.1 pc. The spatial structure, together with the presence of the massive central region, suggest that these filaments may be associated with accretion processes. In order to derive the kinematic properties of the gas in a chemically line-rich source like Sgr B2(N), we have developed a new tool that stacks all the detected transition lines of any molecular species. This permits to increase the signal-to-noise ratio of our observations and average out line blending effects, which are a common problem in line-rich regions. We derive velocity gradients along the filaments of about 20-100 km s$^{-1}$ pc$^{-1}$, which are 10-100 times larger than those typically found on larger scales (1 pc) in other star-forming regions. The mass accretion rates of individual filaments are about 0.05 M$_\odot$ yr$^{-1}$, which result in a total accretion rate of 0.16 M$_\odot$ yr$^{-1}$. Some filaments harbor dense cores that are likely forming stars and stellar clusters. The stellar content of these dense cores is on the order of 50% of the total mass. We conclude that the cores may merge in the center when already forming stellar clusters but still containing a significant amount of gas, resulting in a "damp" merger. The high density and mass of the central region, combined with the presence of converging filaments with high mass, high accretion rates and embedded dense cores already forming stars, suggest that Sgr B2(N) may have the potential to evolve into a super stellar cluster. △ Less

Submitted 26 June, 2019; originally announced June 2019.

Comments: Accepted for publication in Astronomy & Astrophysics. 10 pages with 8 figures and 2 tables, plus Appendices with extra figures and tables. (Abstract modified)

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

VISION - Vienna survey in Orion. III. Young stellar objects in Orion A

Authors: Josefa E. Großschedl, João Alves, Paula S. Teixeira, Hervé Bouy, Jan Forbrich, Charles J. Lada, Stefan Meingast, Álvaro Hacar, Joana Ascenso, Christine Ackerl, Birgit Hasenberger, Rainer Köhler, Karolina Kubiak, Irati Larreina, Lorenz Linhardt, Marco Lombardi, Torsten Möller

Abstract: We have extended and refined the existing young stellar object (YSO) catalogs for the Orion A molecular cloud, the closest massive star-forming region to Earth. This updated catalog is driven by the large spatial coverage (18.3 deg$^2$, $\sim$950 pc$^2$), seeing limited resolution ($\sim$0.7$"$), and sensitivity ($K_s<19$ mag) of the ESO-VISTA near-infrared survey of the Orion A cloud (VISION). Co… ▽ More We have extended and refined the existing young stellar object (YSO) catalogs for the Orion A molecular cloud, the closest massive star-forming region to Earth. This updated catalog is driven by the large spatial coverage (18.3 deg$^2$, $\sim$950 pc$^2$), seeing limited resolution ($\sim$0.7$"$), and sensitivity ($K_s<19$ mag) of the ESO-VISTA near-infrared survey of the Orion A cloud (VISION). Combined with archival mid- to far-infrared data, the VISTA data allow for a refined and more robust source selection. We estimate that among previously known protostars and pre-main-sequence stars with disks, source contamination levels (false positives) are at least $\sim$6.4% and $\sim$2.3%, respectively, mostly due to background galaxies and nebulosities. We identify 274 new YSO candidates using VISTA/Spitzer based selections within previously analyzed regions, and VISTA/WISE based selections to add sources in the surroundings, beyond previously analyzed regions. The WISE selection method recovers about 59% of the known YSOs in Orion A's low-mass star-forming part L1641, which shows what can be achieved by the all-sky WISE survey in combination with deep near-infrared data in regions without the influence of massive stars. The new catalog contains 2980 YSOs, which were classified based on the de-reddened mid-infrared spectral index into 188 protostars, 185 flat-spectrum sources, and 2607 pre-main-sequence stars with circumstellar disks. We find a statistically significant difference in the spatial distribution of the three evolutionary classes with respect to regions of high dust column-density, confirming that flat-spectrum sources are at a younger evolutionary phase compared to Class IIs, and are not a sub-sample seen at particular viewing angles. △ Less

Submitted 19 December, 2018; v1 submitted 1 October, 2018; originally announced October 2018.

Comments: 39 pages, 25 figures, Accepted for publication by A&A

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

The physical and chemical structure of Sagittarius B2. III. Radiative transfer simulations of the hot core SgrB2(M) for methyl cyanide

Authors: S. Pols, A. Schwoerer, P. Schilke, A. Schmiedeke, A. Sanchez-Monge, Th. Moeller

Abstract: We model the emission of methyl cyanide (CH3CN) lines towards the massive hot molecular core SgrB2(M). We aim at reconstructing the CH3CN abundance field and investigating the gas temperature distribution as well as the velocity field. SgrB2(M) was observed with the ALMA in a spectral line survey from 211 to 275 GHz. This frequency range includes several transitions of CH3CN (including isotopologu… ▽ More We model the emission of methyl cyanide (CH3CN) lines towards the massive hot molecular core SgrB2(M). We aim at reconstructing the CH3CN abundance field and investigating the gas temperature distribution as well as the velocity field. SgrB2(M) was observed with the ALMA in a spectral line survey from 211 to 275 GHz. This frequency range includes several transitions of CH3CN (including isotopologues and vibrationally excited states). We employ the three-dimensional radiative transfer toolbox Pandora in order to retrieve the velocity and abundance field by modeling different CH3CN lines. For this purpose, we base our model on the results of a previous study that determined the physical structure of SgrB2(M), i.e.\ the distribution of dust dense cores, ionized regions and heating sources. The morphology of the CH3CN emission can be reproduced by a molecular density field that consists of a superposition of cores with modified Plummer-like density profiles. The averaged relative abundance of CH3CN with respect to H2 ranges from 4x10^{-11} to 2x10^{-8} in the northern part of SgrB2(M) and from 2x10^{-10} to 5x10^{-7} in the southern part. In general, we find that the relative abundance of CH3CN is lower at the center of the very dense and hot cores, causing the general morphology of the CH3CN emission to be shifted with respect to the dust continuum emission. The dust temperature calculated by the radiative transfer simulation based on the available luminosity reaches values up to 900 K. However, in some regions vibrationally excited transitions of CH3CN are underestimated by the model, indicating that the predicted gas temperature, which is assumed to be equal to the dust temperature, is partly underestimated. The determination of the velocity component along the line of sight reveals that a velocity gradient from the north to the south exists in SgrB2(M). △ Less

Submitted 26 February, 2018; originally announced February 2018.

Comments: Accepted for publication in Astronomy & Astrophysics. 18 pages with 18 figures and 1 table, plus Appendices with extra figures and tables

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

The physical and chemical structure of Sagittarius B2. II. Continuum millimeter emission of SgrB2(M) and SgrB2(N) with ALMA

Authors: A. Sanchez-Monge, P. Schilke, A. Schmiedeke, A. Ginsburg, R. Cesaroni, D. C. Lis, S. -L. Qin, H. S. P. Mueller, E. Bergin, C. Comito, Th. Moeller

Abstract: The high-mass star forming sites SgrB2(M) and SgrB2(N) have been the target of numerous studies, revealing e.g. a rich chemistry. We want to characterize their physical and chemical structure using ALMA high-angular resolution observations at mm wavelengths, reaching spatial scales of about 4000 au, and covering the whole band 6 (from 211 to 275 GHz). In order to determine the continuum emission i… ▽ More The high-mass star forming sites SgrB2(M) and SgrB2(N) have been the target of numerous studies, revealing e.g. a rich chemistry. We want to characterize their physical and chemical structure using ALMA high-angular resolution observations at mm wavelengths, reaching spatial scales of about 4000 au, and covering the whole band 6 (from 211 to 275 GHz). In order to determine the continuum emission in line-rich sources, we use a new statistical method: STATCONT. We detect 27 continuum sources in SgrB2(M) and 20 in SgrB2(N). We study the continuum emission across the ALMA band 6, and compare it with previous SMA 345 GHz and VLA 40 GHz observations, to study the nature of the sources detected. The brightest sources are dominated by (partially optically thick) dust emission, while there is an important degree of contamination from ionized gas free-free emission in weaker sources. While the total mass in SgrB2(M) is distributed in many fragments, most of the mass in SgrB2(N) arises from a single object, with filamentary-like structures converging towards the center. There seems to be a lack of low-mass dense cores in both regions. We determine H2 volume densities for the cores of about 10^5-10^7 Msun pc^-3, one to two orders of magnitude higher than the stellar densities of super star clusters. In general, SgrB2(N) is chemically richer than SgrB2(M). There seems to be a correlation between the chemical richness and the mass of the fragments, with more massive clumps being more chemically rich. Both SgrB2(N) and SgrB2(M) harbour a cluster of hot molecular cores. We compare the continuum images with predictions from a detailed 3D radiative transfer model that reproduces the structure of SgrB2 from 45 pc down to 100 au. This dataset, together with ongoing projects in the range 5 to 200 GHz, better constrain the 3D structure of SgrB2, and allow us to understand its physical and chemical structure. △ Less

Submitted 6 April, 2017; originally announced April 2017.

Comments: Accepted for publication in Astronomy & Astrophysics. 20 pages with 13 figures and 4 tables, plus Appendices with extra figures and tables. (Abstract modified)

Journal ref: A&A 604, A6 (2017)

The CO-H2 van der Waals complex and complex organic molecules in cold molecular clouds: a TMC-1C survey

Authors: A. Potapov, A. Sanchez-Monge, P. Schilke, U. U. Graf, Th. Moeller, S. Schlemmer

Abstract: Almost 200 different species have been detected in the interstellar medium (ISM) during the last decades, revealing not only simple species but complex molecules with more than 6 atoms. Other exotic compounds, like the weakly-bound dimer (H2)2, have also been detected in astronomical sources like Jupiter. We aim at detecting for the first time the CO-H2 van der Waals complex in the ISM, which if d… ▽ More Almost 200 different species have been detected in the interstellar medium (ISM) during the last decades, revealing not only simple species but complex molecules with more than 6 atoms. Other exotic compounds, like the weakly-bound dimer (H2)2, have also been detected in astronomical sources like Jupiter. We aim at detecting for the first time the CO-H2 van der Waals complex in the ISM, which if detected can be a sensitive indicator for low temperatures. We use the IRAM30m telescope, located in Pico Veleta (Spain), to search for the CO-H2 complex in a cold, dense core in TMC-1C (with a temperature of 10 K). All the brightest CO-H2 transitions in the 3 mm (80-110 GHz) band have been observed with a spectral resolution of 0.5-0.7 km/s, reaching a rms noise level of 2 mK. The simultaneous observation of a broad frequency band, 16 GHz, has allowed us to conduct a serendipitous spectral line survey. No lines belonging to the CO-H2 complex have been detected. We have set up a new, more stringent upper limit for its abundance to be [CO-H2]/[CO] = 5x10^{-6}, while we expect the abundance of the complex to be in the range 10^{-8}-10^{-3}. The spectral line survey has allowed us to detect 75 lines associated with 41 different species (including isotopologues). We detect a number of complex organic species, e.g. methyl cyanide (CH3CN), methanol (CH3OH), propyne (CH3CCH) and ketene (CH2CO), associated with cold gas (excitation temperatures about 7 K), confirming the presence of these complex species not only in warm objects but also in cold regimes. △ Less

Submitted 10 July, 2016; originally announced July 2016.

Comments: 9 pages, 3 figures, 2 tables. Accepted for publication in Astronomy & Astrophysics

The physical and chemical structure of Sagittarius B2, I. Three-dimensional thermal dust and free-free continuum modeling on 100 au to 45 pc scales

Authors: A. Schmiedeke, P. Schilke, Th. Möller, Á. Sánchez-Monge, E. Bergin, C. Comito, T. Csengeri, D. C. Lis, S. Molinari, S. -L. Qin, R. Rolffs

Abstract: We model the dust and free-free continuum emission in the high-mass star-forming region Sagittarius B2 in order to reconstruct the three-dimensional density and dust temperature distribution, as a crucial input to follow-up studies of the gas velocity field and molecular abundances. We employ the three-dimensional radiative transfer program RADMC-3D to calculate the dust temperature self-consisten… ▽ More We model the dust and free-free continuum emission in the high-mass star-forming region Sagittarius B2 in order to reconstruct the three-dimensional density and dust temperature distribution, as a crucial input to follow-up studies of the gas velocity field and molecular abundances. We employ the three-dimensional radiative transfer program RADMC-3D to calculate the dust temperature self-consistently, provided a given initial density distribution. This density distribution of the entire cloud complex is then recursively reconstructed based on available continuum maps, including both single-dish and high-resolution interferometric maps covering a wide frequency range (40 GHz - 4 THz). The model covers spatial scales from 45 pc down to 100 au, i.e. a spatial dynamic range of 10^5. We find that the density distribution of Sagittarius B2 can be reasonably well fitted by applying a superposition of spherical cores with Plummer-like density profiles. In order to reproduce the spectral energy distribution, we position Sgr B2(N) along the line of sight behind the plane containing Sgr B2(M). We find that the entire cloud complex comprises a total gas mass of 8.0 x 10^6 Msun within a diameter of 45 pc, corresponding to an averaged gas density of 170 Msun/pc^3. We estimate stellar masses of 2400 Msun and 20700 Msun and luminosities of 1.8 x 10^6 Lsun and 1.2 x 10^7 Lsun for Sgr B2(N) and Sgr B2(M), respectively. We report H_2 column densities of 2.9 x 10^24 cm^-2 for Sgr B2(N) and 2.5 x 10^24 cm^-2 for Sgr B2(M) in a 40" beam. For Sgr B2(S), we derive a stellar mass of 1100 Msun, a luminosity of 6.6 x 10^5 Lsun and a H_2 column density of 2.2 x 10^24 cm^-2 in a 40" beam. We calculate a star formation efficiency of 5% for Sgr B2(N) and 50% for Sgr B2(M), indicating that most of the gas content in Sgr B2(M) has already been converted to stars or dispersed. △ Less

Submitted 6 February, 2016; originally announced February 2016.

Comments: Accepted for publication in A&A. 30 pages, 24 figures, 6 tables

Journal ref: A&A 588, A143 (2016)

eXtended CASA Line Analysis Software Suite (XCLASS)

Authors: T. Möller, C. Endres, P. Schilke

Abstract: The eXtended CASA Line Analysis Software Suite (XCLASS) is a toolbox for the Common Astronomy Software Applications package (CASA) containing new functions for modeling interferometric and single dish data. Among the tools is the myXCLASS program which calculates synthetic spectra by solving the radiative transfer equation for an isothermal object in one dimension, whereas the finite source size a… ▽ More The eXtended CASA Line Analysis Software Suite (XCLASS) is a toolbox for the Common Astronomy Software Applications package (CASA) containing new functions for modeling interferometric and single dish data. Among the tools is the myXCLASS program which calculates synthetic spectra by solving the radiative transfer equation for an isothermal object in one dimension, whereas the finite source size and dust attenuation are considered as well. Molecular data required by the myXCLASS program are taken from an embedded SQLite3 database containing entries from the Cologne Database for Molecular Spectroscopy CDMS) and JPL using the Virtual Atomic and Molecular Data Center (VAMDC) portal. Additionally, the toolbox provides an interface for the model optimizer package Modeling and Analysis Generic Interface for eXternal numerical codes (MAGIX), which helps to find the best description of observational data using myXCLASS (or another external model program), i.e., finding the parameter set that most closely reproduces the data. △ Less

Submitted 21 October, 2016; v1 submitted 17 August, 2015; originally announced August 2015.

Comments: 18 pages, 10 figures

Journal ref: A&A 598, A7 (2017)

Evolution of complex organic molecules in hot molecular cores: Synthetic spectra at (sub-)mm wavebands

Authors: Rumpa Choudhury, Peter Schilke, Gwendoline Stéphan, Edwin A. Bergin, Thomas Möller, Anika Schmiedeke, Alexander Zernickel

Abstract: Hot molecular cores (HMCs) are intermediate stages of high-mass star formation and are also known for their rich emission line spectra at (sub-)mm wavebands. The observed spectral feature of HMCs such as total number of emission lines and associated line intensities are also found to vary with evolutionary stages. We developed various 3D models for HMCs guided by the evolutionary scenarios propose… ▽ More Hot molecular cores (HMCs) are intermediate stages of high-mass star formation and are also known for their rich emission line spectra at (sub-)mm wavebands. The observed spectral feature of HMCs such as total number of emission lines and associated line intensities are also found to vary with evolutionary stages. We developed various 3D models for HMCs guided by the evolutionary scenarios proposed by recent empirical and modeling studies. We then investigated the spatio-temporal variation of temperature and molecular abundances in HMCs by consistently coupling gas-grain chemical evolution with radiative transfer calculations. We explored the effects of varying physical conditions on molecular abundances including density distribution and luminosity evolution of the central protostar(s). The time-dependent temperature structure of the hot core models provides a realistic framework for investigating the spatial variation of ice mantle evaporation as a function of evolutionary timescales. With increasing protostellar luminosity, the water ice evaporation font ($\sim$100K) expands and the spatial distribution of gas phase abundances of these COMs also spreads out. We simulated the synthetic spectra for these models at different evolutionary timescales to compare with observations. A qualitative comparison of the simulated and observed spectra suggests that these self-consistent hot core models can reproduce the notable trends in hot core spectral variation within the typical hot core timescales of 10$^{5}$ year. These models predict that the spatial distribution of various emission line maps will also expand with evolutionary time. The model predictions can be compared with high resolution observation that can probe scales of a few thousand AU in high-mass star forming regions such as from ALMA.[Abridged] △ Less

Submitted 11 January, 2015; originally announced January 2015.

Comments: accepted for publication in A&A

Journal ref: A&A 575, A68 (2015)

Ubiquitous argonium (ArH$^+$) in the diffuse interstellar medium -- a molecular tracer of almost purely atomic gas

Authors: Peter Schilke, David A. Neufeld, Holger S. P. Müller, Claudia Comito, Edwin A. Bergin, Dariusz C. Lis, Maryvonne Gerin, John H. Black, Mark Wolfire, Nick Indriolo, John. C. Pearson, Karl M. Menten, Benjamin Winkel, Alvaro Sanchez-Monge, Thomas Moeller, Benjamin Godard, Edith Falgarone

Abstract: We describe the assignment of a previously unidentified interstellar absorption line to ArH$^+$ and discuss its relevance in the context of hydride absorption in diffuse gas with a low H$_2$ fraction. The column densities along several lines of sight are determined and discussd in the framework of chemical models. The column densities of ArH$^+$ are compared to those of other species, tracing inte… ▽ More We describe the assignment of a previously unidentified interstellar absorption line to ArH$^+$ and discuss its relevance in the context of hydride absorption in diffuse gas with a low H$_2$ fraction. The column densities along several lines of sight are determined and discussd in the framework of chemical models. The column densities of ArH$^+$ are compared to those of other species, tracing interstellar medium (ISM) components with different H$_2$ abundances. Chemical models are constructed, taking UV radiation and cosmic ray ionization into account. Due to the detection of two isotopologues, $^{36}$ArH$^+$ and $^{38}$ArH$^+$, we are confident about the carrier assignment to ArH$^+$. NeH$^+$ is not detected with a limit of [NeH$^+$]/[ArH$^+$] $\le$ 0.1. The derived column densities agree well with the predictions of chemical models. ArH$^+$ is a unique tracer of gas with a fractional H$_2$ abundance of $10^{-4}- 10^{-3}$ and shows little correlation with H$_2$O$^+$, which traces gas with a fractional H$_2$ abundance of $\approx $0.1. A careful analysis of variations in the ArH$^+$, OH$^+$, H$_2$O$^+$ and HF column densities promises to be a faithful tracer of the distribution of the H$_2$ fractional abundance, providing unique information on a poorly known phase in the cycle of interstellar matter, its transition from atomic diffuse gas to dense molecular gas traced by CO emission. Abundances of these species put strong observational constraints upon magnetohydrodynamical (MHD) simulations of the interstellar medium, and potentially could evolve into a tool to characterize the ISM. Paradoxically, the ArH$^+$ molecule is a better tracer of \new{almost} purely atomic hydrogen gas than H{\sc i} itself, since H{\sc i} can also be present in gas with a significant molecular content, but ArH$^+$ singles out gas that is $>99.9$\% atomic. △ Less

Submitted 14 April, 2014; v1 submitted 31 March, 2014; originally announced March 2014.

Comments: 12 pages, 11 figures, Astronomy and Astrophysics, accepted revised version: previous version had column density of SgrB2 sources wrong

Journal ref: Astron. Astrophys. 566 (2014) Art. No. A29

Modeling and Analysis Generic Interface for eXternal numerical codes (MAGIX)

Authors: T. Möller, I. Bernst, D. Panoglou, D. Muders, V. Ossenkopf, M. Röllig, P. Schilke

Abstract: The modeling and analysis generic interface for external numerical codes (MAGIX) is a model optimizer developed under the framework of the coherent set of astrophysical tools for spectroscopy (CATS) project. The MAGIX package provides a framework of an easy interface between existing codes and an iterating engine that attempts to minimize deviations of the model results from available observationa… ▽ More The modeling and analysis generic interface for external numerical codes (MAGIX) is a model optimizer developed under the framework of the coherent set of astrophysical tools for spectroscopy (CATS) project. The MAGIX package provides a framework of an easy interface between existing codes and an iterating engine that attempts to minimize deviations of the model results from available observational data, constraining the values of the model parameters and providing corresponding error estimates. Many models (and, in principle, not only astrophysical models) can be plugged into MAGIX to explore their parameter space and find the set of parameter values that best fits observational/experimental data. MAGIX complies with the data structures and reduction tools of ALMA (Atacama Large Millimeter Array), but can be used with other astronomical and with non-astronomical data. △ Less

Submitted 24 October, 2012; originally announced October 2012.

Comments: 12 pages, 15 figures, 2 tables, paper is also available at http://www.aanda.org/articles/aa/pdf/forth/aa20063-12.pdf

Molecular line survey of the high-mass star-forming region NGC 6334I with Herschel/HIFI and the SMA

Authors: A. Zernickel, P. Schilke, A. Schmiedeke, D. C. Lis, C. L. Brogan, C. Ceccarelli, C. Comito, M. Emprechtinger, T. R. Hunter, T. Möller

Abstract: We aim at deriving the molecular abundances and temperatures of the hot molecular cores in the high-mass star-forming region NGC 6334I and consequently deriving their physical and astrochemical conditions. In the framework of the Herschel guaranteed time key program CHESS, NGC 6334I is investigated by using HIFI aboard the Herschel Space Observatory. A spectral line survey is carried out in the fr… ▽ More We aim at deriving the molecular abundances and temperatures of the hot molecular cores in the high-mass star-forming region NGC 6334I and consequently deriving their physical and astrochemical conditions. In the framework of the Herschel guaranteed time key program CHESS, NGC 6334I is investigated by using HIFI aboard the Herschel Space Observatory. A spectral line survey is carried out in the frequency range 480-1907 GHz, and auxiliary interferometric data from the SMA in the 230 GHz band provide spatial information for disentangling the different physical components contributing to the HIFI spectrum. The spectral lines are identified with the aid of former surveys and spectral line catalogs. The observed spectrum is then compared to a simulated synthetic spectrum with XCLASS, assuming local thermal equilibrium, and best fit parameters are derived using the model optimization package MAGIX. A total of 46 molecules are identified, with 31 isotopologues, resulting in about 4300 emission and absorption lines. High- energy levels of the dominant emitter methanol and vibrationally excited HCN are detected. The number of unidentified lines remains low with 75, or less than 2 percent of the lines detected. The modeling suggests that several spectral features need two or more components to be fitted properly. Other components could be assigned to cold foreground clouds or to outflows, most visible in the SiO emission. A chemical variation between the two embedded hot cores is found, with more N-bearing molecules identified in SMA1 and O-bearing molecules in SMA2. Spectral line surveys give powerful insights into the study of the interstellar medium. Different molecules trace different physical conditions like the inner hot core, the envelope, the outflows or the cold foreground clouds. The derived molecular abundances provide further constraints for astrochemical models. △ Less

Submitted 27 August, 2012; originally announced August 2012.

Comments: 30 pages including appendix, 49 figures, accepted for publication in Astronomy and Astrophysics

The quest for companions to post-common envelope binaries: I. Searching a sample of stars from the CSS and SDSS

Authors: U. Backhaus, S. Bauer, K. Beuermann, J. Diese, S. Dreizler, F. V. Hessman, T. -O. Husser, K. -H. Klapdohr, J. Moellmanns, R. Schuenecke, J. Dette, J. Dubbert, T. Miosga, A. L. Rochus Vogel, S. Simons, S. Biriuk, M. Debrah, M. Griemens, A. Hahn, T. Moeller, M. Pawlowski, M. Schweizer, A. -L. Speck, C. Zapros, T. Bollmann , et al. (17 additional authors not shown)

Abstract: As part of an ongoing collaboration between student groups at high schools and professional astronomers, we have searched for the presence of circum-binary planets in a bona-fide unbiased sample of twelve post-common envelope binaries (PCEBs) from the Catalina Sky Survey (CSS) and the Sloan Digital Sky Survey (SDSS). Although the present ephemerides are significantly more accurate than previous on… ▽ More As part of an ongoing collaboration between student groups at high schools and professional astronomers, we have searched for the presence of circum-binary planets in a bona-fide unbiased sample of twelve post-common envelope binaries (PCEBs) from the Catalina Sky Survey (CSS) and the Sloan Digital Sky Survey (SDSS). Although the present ephemerides are significantly more accurate than previous ones, we find no clear evidence for orbital period variations between 2005 and 2011 or during the 2011 observing season. The sparse long-term coverage still permits O-C variations with a period of years and an amplitude of tens of seconds, as found in other systems. Our observations provide the basis for future inferences about the frequency with which planet-sized or brown-dwarf companions have either formed in these evolved systems or survived the common envelope (CE) phase. △ Less

Submitted 10 January, 2012; originally announced January 2012.

Comments: accepted by A&A

Herschel observations of extra-ordinary sources: Detecting spiral arm clouds by CH absorption lines

Authors: S. -L. Qin, P. Schilke, C. Comito, T. Möller, R. Rolffs, H. S. P. Müller, A. Belloche, K. M. Menten, D. C. Lis, T. G. Phillips, E. A. Bergin, T. A. Bell, N. R. Crockett, G. A. Blake, S. Cabrit, E. Caux, C. Ceccarelli, J. Cernicharo, F. Daniel, M. -L. Dubernet, M. Emprechtinger, P. Encrenaz, E. Falgarone, M. Gerin, T. F. Giesen , et al. (45 additional authors not shown)

Abstract: We have observed CH absorption lines ($J=3/2, N=1 \leftarrow J=1/2, N=1$) against the continuum source Sgr~B2(M) using the \textit{Herschel}/HIFI instrument. With the high spectral resolution and wide velocity coverage provided by HIFI, 31 CH absorption features with different radial velocities and line widths are detected and identified. The narrower line width and lower column density clouds sho… ▽ More We have observed CH absorption lines ($J=3/2, N=1 \leftarrow J=1/2, N=1$) against the continuum source Sgr~B2(M) using the \textit{Herschel}/HIFI instrument. With the high spectral resolution and wide velocity coverage provided by HIFI, 31 CH absorption features with different radial velocities and line widths are detected and identified. The narrower line width and lower column density clouds show `spiral arm' cloud characteristics, while the absorption component with the broadest line width and highest column density corresponds to the gas from the Sgr~B2 envelope. The observations show that each `spiral arm' harbors multiple velocity components, indicating that the clouds are not uniform and that they have internal structure. This line-of-sight through almost the entire Galaxy offers unique possibilities to study the basic chemistry of simple molecules in diffuse clouds, as a variety of different cloud classes are sampled simultaneously. We find that the linear relationship between CH and H$_2$ column densities found at lower $A_V$ by UV observations does not continue into the range of higher visual extinction. There, the curve flattens, which probably means that CH is depleted in the denser cores of these clouds. △ Less

Submitted 12 July, 2010; originally announced July 2010.

Comments: Accepted for publication in A&A, HIFI Special Issue