Showing 1–35 of 35 results for author: Ahmadi, A

The accretion burst of the massive young stellar object G323.46 -0.08

Authors: V. Wolf, B. Stecklum, A. Caratti o Garatti, P. A. Boley, Ch. Fischer, T. Harries, J. Eislöffel, H. Linz, A. Ahmadi, J. Kobus, X. Haubois, A. Matter, P. Cruzalebes

Abstract: Accretion bursts from low-mass young stellar objects (YSOs) are known for many decades. In recent years, the first accretion bursts of massive YSOs (MYSOs) have been observed. These phases of intense protostellar growth are of particular importance for studying massive star formation. Bursts of MYSOs are accompanied by flares of Class II methanol masers (hereafter masers), caused by an increase in… ▽ More Accretion bursts from low-mass young stellar objects (YSOs) are known for many decades. In recent years, the first accretion bursts of massive YSOs (MYSOs) have been observed. These phases of intense protostellar growth are of particular importance for studying massive star formation. Bursts of MYSOs are accompanied by flares of Class II methanol masers (hereafter masers), caused by an increase in exciting mid-infrared (MIR) emission. The G323.46$-$0.08 (hereafter G323) event extends the small sample of known MYSO bursts. Maser observations of the MYSO G323 show evidence of a flare, which was presumed to be caused by an accretion burst. This should be verified with IR data. We used time-dependent radiative transfer (TDRT) to characterize the heating and cooling timescales for eruptive MYSOs and to infer the main burst parameters. The G323 accretion burst is confirmed. It reached its peak in late 2013/early 2014 with a Ks-band increase of 2.5mag. TDRT indicates that the duration of the thermal afterglow in the far-infrared (FIR) can exceed the burst duration by years. The latter was proved by SOFIA observations, which indicate a flux increase of $(14.2\pm4.6)$% at $70\, \rm μm$ and $(8.5\pm6.1)$% at $160\, μ$m in 2022 (2 years after the burst end). A one-sided light echo emerged that was propagating into the interstellar medium. The G323 burst is probably the most energetic MYSO burst observed so far. Within $8.4 \rm \, yrs$, an energy of $(0.9\pm_{0.8}^{2.5}) \times 10^{47}\,\rm erg$ was released. The short timescale points to the accretion of a compact body, while the burst energy corresponds to an accumulated mass of at least $(7\pm_{6}^{20})\,M_{Jup}$ and possibly even more if the protostar is bloated. In this case, the accretion event might have triggered protostellar pulsations, which give rise to the observed maser periodicity. △ Less

Submitted 21 May, 2024; v1 submitted 16 May, 2024; originally announced May 2024.

ALMA view of the L1448-mm protostellar system on disk scales: CH$_3$OH and H$^{13}$CN as new disk wind tracers

Authors: P. Nazari, B. Tabone, A. Ahmadi, S. Cabrit, E. F. van Dishoeck, C. Codella, J. Ferreira, L. Podio, Ł. Tychoniec, M. L. van Gelder

Abstract: Protostellar disks are known to accrete, however, the exact mechanism that extracts the angular momentum and drives accretion in the low-ionization "dead" region of the disk is under debate. In recent years, magneto-hydrodynamic (MHD) disk winds have become a popular solution. Yet, observations of these winds require both high spatial resolution (${\sim}10$s au) and high sensitivity, which has res… ▽ More Protostellar disks are known to accrete, however, the exact mechanism that extracts the angular momentum and drives accretion in the low-ionization "dead" region of the disk is under debate. In recent years, magneto-hydrodynamic (MHD) disk winds have become a popular solution. Yet, observations of these winds require both high spatial resolution (${\sim}10$s au) and high sensitivity, which has resulted in only a handful of MHD disk wind candidates so far. In this work we present high angular resolution (${\sim}30$ au) ALMA observations of the emblematic L1448-mm protostellar system and find suggestive evidence for an MHD disk wind. The disk seen in dust continuum (${\sim}0.9$ mm) has a radius of ${\sim}23$ au. Rotating infall signatures in H$^{13}$CO$^+$ indicate a central mass of $0.4\pm 0.1$ M$_\odot$ and a centrifugal radius similar to the dust disk radius. Above the disk, we unveil rotation signatures in the outflow traced by H$^{13}$CN, CH$_3$OH, and SO lines and find a kinematical structure consistent with theoretical predictions for MHD disk winds. This is the first detection of an MHD disk wind candidate in H$^{13}$CN and CH$_3$OH. The wind launching region estimated from cold MHD wind theory extends out to the disk edge. The magnetic lever arm parameter would be $λ_φ \simeq 1.7$, in line with recent non-ideal MHD disk models. The estimated mass-loss rate is ${\sim}4$ times the protostellar accretion rate ($\dot{M}_{\rm acc} \simeq 2 \times 10^{-6} M_{\odot}/yr$) and suggests that the rotating wind could carry enough angular momentum to drive disk accretion. △ Less

Submitted 28 February, 2024; originally announced February 2024.

Comments: Accepted for publication in A&A

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

Constraining the gas mass of Herbig disks using CO isotopologues

Authors: L. M. Stapper, M. R. Hogerheijde, E. F. van Dishoeck, L. Lin, A. Ahmadi, A. S. Booth, S. L. Grant, K. Immer, M. Leemker, A. F. Pérez-Sánchez

Abstract: The total disk mass sets the formation potential for exoplanets. Carbon-monoxide (CO) has been used as a gas mass tracer in T Tauri disks, but was found to be less abundant than expected due to freeze-out and chemical conversion of CO on the surfaces of cold dust grains. The disks around more massive intermediate mass pre-main sequence stars called Herbig disks are likely to be warmer, allowing fo… ▽ More The total disk mass sets the formation potential for exoplanets. Carbon-monoxide (CO) has been used as a gas mass tracer in T Tauri disks, but was found to be less abundant than expected due to freeze-out and chemical conversion of CO on the surfaces of cold dust grains. The disks around more massive intermediate mass pre-main sequence stars called Herbig disks are likely to be warmer, allowing for the possibility of using CO as a more effective total gas mass tracer. Using ALMA archival data and new NOEMA data of 12CO, 13CO, and C18O transitions of 35 Herbig disks within 450 pc, the masses are determined using the thermo-chemical code Dust And LInes (DALI). The majority of Herbig disks for which 13CO and C18O are detected are optically thick in both. Computing the gas mass using a simple optically thin relation between line flux and column density results in an underestimate of the gas mass of at least an order of magnitude compared to the masses obtained with DALI. The inferred gas masses with DALI are consistent with a gas-to-dust ratio of at least 100. These gas-to-dust ratios are two orders of magnitude higher compared to those found for T Tauri disks using similar techniques, even over multiple orders of magnitude in dust mass, illustrating the importance of chemical conversion of CO in colder T Tauri disks. Similar high gas-to-dust ratios are found for Herbig group I and II disks. Since group II disks have dust masses comparable to T Tauri disks, their higher CO gas masses illustrate the determining role of temperature. Compared to debris disks, Herbig disks have four orders of magnitude higher gas masses. At least one Herbig disk, HD 163296, has a detected molecular disk wind, but our investigation has not turned up other detections of the CO disk wind in spite of similar sensitivities. △ Less

Submitted 6 December, 2023; originally announced December 2023.

Comments: Accepted for publication in Astronomy and Astrophysics. 25 pages, 19 figures, plus appendices

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

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)

Physical and chemical complexity in high-mass star-forming regions with ALMA. I. Overview and evolutionary trends of physical properties

Authors: Caroline Gieser, Henrik Beuther, Dmitry Semenov, Aida Ahmadi, Thomas Henning, Molly R. A. Wells

Abstract: In this study, we investigate how physical properties, such as the density and temperature profiles, evolve on core scales through the evolutionary sequence during high-mass star formation ranging from protostars in cold infrared dark clouds to evolved UCHII regions. We observed 11 high-mass star-forming regions with ALMA at 3 mm wavelengths. Based on the 3 mm continuum morphology and recombinatio… ▽ More In this study, we investigate how physical properties, such as the density and temperature profiles, evolve on core scales through the evolutionary sequence during high-mass star formation ranging from protostars in cold infrared dark clouds to evolved UCHII regions. We observed 11 high-mass star-forming regions with ALMA at 3 mm wavelengths. Based on the 3 mm continuum morphology and recombination line emission, tracing locations with free-free (ff) emission, the fragmented cores analyzed in this study are classified into either dust or dust+ff cores. In addition, we resolve three cometary UCHII regions with extended 3 mm emission that is dominated by free-free emission. The temperature structure and radial profiles (T~r^-q ) are determined by modeling molecular emission of CH3CN and CH313CN with XCLASS and by using the HCN-to- HNC intensity ratio as probes for the gas kinetic temperature. The density profiles (n~r^-p ) are estimated from the 3 mm continuum visibility profiles. The masses M and H2 column densities N(H2) are then calculated from the 3 mm dust continuum emission. Results. We find a large spread in mass and peak H2 column density in the detected sources ranging from 0.1-150 Msun and 10^23 - 10^26 cm-2 , respectively. Including the results of the CORE and CORE-extension studies (Gieser et al. 2021, 2022) to increase the sample size, we find evolutionary trends on core scales for the temperature power-~law index q increasing from 0.1 to 0.7 from infrared dark clouds to UCHII regions, while for the the density power-law index p on core scales, we do not find strong evidence for an evolutionary trend. However, we find that on the larger clump scales throughout these evolutionary phases the density profile flattens from p = 2.2 to p = 1.2. (abridged) △ Less

Submitted 24 April, 2023; v1 submitted 14 April, 2023; originally announced April 2023.

Comments: 21 pages, 10 figures, accepted for publication in A&A

Journal ref: A&A 674, A160 (2023)

N-bearing complex organics toward high-mass protostars: Constant ratios pointing to formation in similar pre-stellar conditions across a large mass range

Authors: P. Nazari, J. D. Meijerhof, M. L. van Gelder, A. Ahmadi, E. F. van Dishoeck, B. Tabone, D. Langeroodi, N. F. W. Ligterink, J. Jaspers, M. T. Beltrán, G. A. Fuller, Á. Sánchez-Monge, P. Schilke

Abstract: No statistical study of COMs toward a large sample of high-mass protostars with ALMA has been carried out so far. We aim to study six N-bearing species: CH$_3$CN, HNCO, NH$_2$CHO, C$_2$H$_5$CN, C$_2$H$_3$CN and CH$_3$NH$_2$ in a large sample of high-mass protostars. From the ALMAGAL survey, 37 of the most line-rich hot molecular cores are selected. Next, we fit their spectra and find column densit… ▽ More No statistical study of COMs toward a large sample of high-mass protostars with ALMA has been carried out so far. We aim to study six N-bearing species: CH$_3$CN, HNCO, NH$_2$CHO, C$_2$H$_5$CN, C$_2$H$_3$CN and CH$_3$NH$_2$ in a large sample of high-mass protostars. From the ALMAGAL survey, 37 of the most line-rich hot molecular cores are selected. Next, we fit their spectra and find column densities and excitation temperatures of the above N-bearing species, in addition to CH$_3$OH. We (tentatively) detect CH$_3$NH$_2$ in $\sim32%$ of the sources. We find three groups of species when comparing their excitation temperatures: hot (NH$_2$CHO; Tex > 250 K), warm (C$_2$H$_3$CN, HN$^{13}$CO and CH$_{3}^{13}$CN; 100 K < Tex < 250 K) and cold species (CH$_3$OH and CH$_3$NH$_2$; Tex < 100 K). This temperature segregation reflects the trend seen in their sublimation temperatures and validates the idea of onion-like structure of COMs around protostars. Moreover, the molecules studied here show constant column density ratios across low- and high-mass protostars with scatter less than a factor $\sim3$ around the mean. The constant column density ratios point to a common formation environment of COMs or their precursors, most likely in the pre-stellar ices. The scatter around the mean of the ratios, although small, varies depending on the species considered. This spread can either have a physical origin (source structure, line or dust optical depth) or a chemical one. Formamide is most prone to the physical effects as it is tracing the closest regions to the protostars, whereas such effects are small for other species. Assuming that all molecules form in the pre-stellar ices, the scatter variations could be explained by differences in lifetimes or physical conditions of the pre-stellar clouds. If the pre-stellar lifetimes are the main factor, they should be similar for low- and high-mass protostars. △ Less

Submitted 23 August, 2022; originally announced August 2022.

Comments: 49 pages, 46 figures, Accepted for publication in A&A

Journal ref: A&A 668, A109 (2022)

Methanol deuteration in high-mass protostars

Authors: M. L. van Gelder, J. Jaspers, P. Nazari, A. Ahmadi, E. F. van Dishoeck, M. T. Beltrán, G. A. Fuller, Á. Sánchez-Monge, P. Schilke

Abstract: The deuteration of molecules forming in the ices such as methanol (CH$_3$OH) is sensitive to the physical conditions during their formation in dense cold clouds and can be probed through observations of deuterated methanol in hot cores. Observations with ALMA containing transitions of CH$_3$OH, CH$_2$DOH, CHD$_2$OH, $^{13}$CH$_3$OH, and CH$_3^{18}$OH are investigated. The column densities of CH… ▽ More The deuteration of molecules forming in the ices such as methanol (CH$_3$OH) is sensitive to the physical conditions during their formation in dense cold clouds and can be probed through observations of deuterated methanol in hot cores. Observations with ALMA containing transitions of CH$_3$OH, CH$_2$DOH, CHD$_2$OH, $^{13}$CH$_3$OH, and CH$_3^{18}$OH are investigated. The column densities of CH$_2$DOH, CHD$_2$OH, and CH$_3$OH are determined for all sources, where the column density of CH$_3$OH is derived from optically thin $^{13}$C and $^{18}$O isotopologues. Consequently, the D/H ratio of methanol is derived taking statistical effects into account. Singly deuterated methanol (CH$_2$DOH) is detected toward 25 of the 99 sources in our sample of the high-mass protostars. Including upper limits, the $\rm (D/H)_{CH_3OH}$ ratio inferred from $N_\mathrm{CH_2DOH}/N_\mathrm{CH_3OH}$ was derived for 38 of the 99 sources and varies between $\sim10^{-3}-10^{-2}$. Including other high-mass hot cores from the literature, the mean methanol D/H ratio is $1.1\pm0.7\times10^{-3}$. This is more than one order of magnitude lower than what is seen for low-mass protostellar systems ($2.2\pm1.2\times10^{-2}$). Doubly deuterated methanol (CHD$_2$OH) is detected toward 11 of the 99 sources. Including upper limits for 15 sources, the $\rm (D/H)_{CH_2DOH}$ ratios derived from $N_\mathrm{CHD_2OH}/N_\mathrm{CH_2DOH}$ are more than two orders of magnitude higher than $\rm (D/H)_{CH_3OH}$ with an average of $2.0\pm0.8\times10^{-1}$ which is similar to what is found for low-mass sources. Comparison with literature GRAINOBLE models suggests that the high-mass prestellar phases are either warm ($>20$ K) or live shorter than the free-fall timescale. In contrast, for low-mass protostars, both a low temperature of $<15$ K and a prestellar phase timescale longer than the free-fall timescale are necessary. △ Less

Submitted 12 August, 2022; originally announced August 2022.

Comments: 17 pages, 9 figures, accepted to A&A

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

The prebiotic molecular inventory of Serpens SMM1: II. The building blocks of peptide chains

Authors: Niels F. W. Ligterink, Aida Ahmadi, Bijaya. Luitel, Audrey Coutens, Hannah Calcutt, Łukasz Tychoniec, Harold Linnartz, Jes K. Jørgensen, Robin T. Garrod, Jordy Bouwman

Abstract: This work aims to constrain the abundances of interstellar amides, by searching for this group of prebiotic molecules in the intermediate-mass protostar Serpens SMM1-a. ALMA observations are conducted toward Serpens SMM1. A spectrum is extracted toward the SMM1-a position and analyzed with the CASSIS line analysis software for the presence of characteristic rotational lines of a number of amides a… ▽ More This work aims to constrain the abundances of interstellar amides, by searching for this group of prebiotic molecules in the intermediate-mass protostar Serpens SMM1-a. ALMA observations are conducted toward Serpens SMM1. A spectrum is extracted toward the SMM1-a position and analyzed with the CASSIS line analysis software for the presence of characteristic rotational lines of a number of amides and other molecules. NH$_{2}$CHO, NH$_{2}$CHO $ν_{12}$=1, NH$_{2}^{13}$CHO, CH$_{3}$C(O)NH$_{2}$ $ν$=0,1, CH$_{2}$DOH, CH$_{3}$CHO, and CH$_{3}$C(O)CH$_{3}$ are securely detected, while trans-NHDCHO, NH$_{2}$CDO, CH$_{3}$NHCHO $ν$=0,1, CH$_{3}$COOH, and HOCH$_{2}$CHO are tentatively identified. The results of this work are compared with detections presented in the literature. A uniform CH$_{3}$C(O)NH$_{2}$/NH$_{2}$CHO ratio is found for a group of interstellar sources with vast physical differences. A similar ratio is seen for CH$_{3}$NHCHO, based on a smaller data sample. The D/H ratio of NH$_{2}$CHO is about 1--3\% and is close to values found in the low-mass source IRAS~16293--2422B. The formation of CH$_{3}$C(O)NH$_{2}$ and NH$_{2}$CHO is likely linked. Formation of these molecules on grain surfaces during the dark cloud stage is a likely scenario. The high D/H ratio of NH$_{2}$CHO is also seen as an indication that these molecules are formed on icy dust grains. As a direct consequence, amides are expected to be present in the most pristine material from which planetary systems form, thus providing a reservoir of prebiotic material. △ Less

Submitted 19 February, 2022; originally announced February 2022.

Comments: Accepted for publication in ACS Earth and Space Chemistry; 59 pages, 10 figures

Importance of source structure on complex organics emission. I. Observations of CH$_3$OH from low-mass to high-mass protostars

Authors: M. L. van Gelder, P. Nazari, B. Tabone, A. Ahmadi, E. F. van Dishoeck, M. T. Beltrán, G. A. Fuller, N. Sakai, Á. Sánchez-Monge, P. Schilke, Y. -L. Yang, Y. Zhang

Abstract: Complex organic molecules (COMs) are often observed toward embedded Class 0 and I protostars. However, not all Class 0 and I protostars exhibit COMs emission. In this work, variations in methanol (CH$_3$OH) emission are studied to test if absence of CH$_3$OH emission can be linked to source properties. Combining both new and archival observations with ALMA and sources from the literature, a sample… ▽ More Complex organic molecules (COMs) are often observed toward embedded Class 0 and I protostars. However, not all Class 0 and I protostars exhibit COMs emission. In this work, variations in methanol (CH$_3$OH) emission are studied to test if absence of CH$_3$OH emission can be linked to source properties. Combining both new and archival observations with ALMA and sources from the literature, a sample of 184 low-mass and high-mass protostars is investigated. The warm (T > 100 K) gaseous CH$_3$OH mass, $M_{\rm CH_3OH}$, is determined for each source using primarily optically thin isotopologues. On average, Class I protostellar systems seem to have less warm $M_{\rm CH_3OH}$ ($<10^{-10}$ M$_\odot$) than younger Class 0 sources ($\sim10^{-7}$ M$_\odot$). High-mass sources in our sample show higher warm $M_{\rm CH_3OH}$ up to $10^{-7}-10^{-3}$ M$_\odot$. To take into account the effect of the source's overall mass on $M_{\rm CH_3OH}$, a normalized CH$_3$OH mass is defined as $M_{\rm CH_3OH}/M_{\rm dust,0}$, where $M_{\rm dust,0}$ is the cold + warm dust mass within a fixed radius. Excluding upper limits, a simple power-law fit to the normalized warm CH$_3$OH masses results in $M_{\rm CH_3OH}/M_{\rm dust,0}\propto L_{\rm bol}^{0.70\pm0.05}$. This is in good agreement with a simple hot core toy model which predicts that the normalized $M_{\rm CH_3OH}$ increases with $L_{\rm bol}^{0.75}$ due to the snowline moving outward. Sources for which the size of the disk is equivalent or smaller than the estimated 100 K radius agree well with the best-fit power-law model, whereas sources with significantly larger disks show up to two orders of magnitude lower normalized warm CH$_3$OH masses. Based on the latter results, we suggest that source structure such as a disk can result in colder gas and thus fewer COMs in the gas phase. Additionally, optically thick dust can hide the emission of COMs. △ Less

Submitted 25 February, 2022; v1 submitted 9 February, 2022; originally announced February 2022.

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

Journal ref: A&A 662, A67 (2022)

The sharp ALMA view of infall and outflow in the massive protocluster G31.41+0.31

Authors: M. T. Beltrán, V. M. Rivilla, R. Cesaroni, D. Galli, L. Moscadelli, A. Ahmadi, H. Beuther, S. Etoka, C. Goddi, P. D. Klaassen, R. Kuiper, M. S. N. Kumar, A. Lorenzani, T. Peters, Á. Sánchez-Monge, P. Schilke, F. van der Tak, S. Vig

Abstract: Context. To better understand the formation of high-mass stars, it is fundamental to investigate how matter accretes onto young massive stars, how it is ejected, and how all this differs from the low-mass case. The massive protocluster G31.41+0.31 is the ideal target to study all these processes because observations at millimeter and centimeter wavelengths have resolved the emission of the Main co… ▽ More Context. To better understand the formation of high-mass stars, it is fundamental to investigate how matter accretes onto young massive stars, how it is ejected, and how all this differs from the low-mass case. The massive protocluster G31.41+0.31 is the ideal target to study all these processes because observations at millimeter and centimeter wavelengths have resolved the emission of the Main core into at least four massive dust continuum sources, named A, B, C, and D, within 1" or 0.018 pc, and have identified signatures of infall and several outflows associated with the core. Aims. We study the interplay between infall and outflow in G31.41+0.31 by investigating at a spatial resolution of a few 100 au their properties and their possible impact on the core. Methods. We carried out molecular line observations of typical high-density tracers, such as CH3CN or H2CO, and shock and outflow tracers, such as SiO, with ALMA at 1.4 mm that achieved an angular resolution of 0.09" (340 au). Results. The observations have revealed inverse P-Cygni profiles in CH3CN and H2CO toward the four sources embedded in the Main core, suggesting that all of them are undergoing collapse. The infall rates, estimated from the red-shifted absorption are on the order of 1E-2 Msun/yr. The individual infall rates imply that the accretion timescale of the Main core is an order of magnitude smaller than its rotation timescale. This confirms that rotating toroids such as the G31 Main core are non-equilibrium, transient collapsing structures that need to be constantly replenished with fresh material from a large-scale reservoir. For sources B, C, and D, the infall could be accelerating inside the sources, while for source A, the presence of a second emission component complicates the interpretation. The SiO observations have revealed the presence of at least six outflows in the G31.41+0.31 star-forming region, ... △ Less

Submitted 25 January, 2022; originally announced January 2022.

Comments: 17 pages, 9 figures, 2 tables. Accepted by A&A

Journal ref: A&A 659, A81 (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)

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)

Fragmentation in the massive G31.41+0.31 protocluster

Authors: M. T. Beltrán, V. M. Rivilla, R. Cesaroni, L. T. Maud, D. Galli, L. Moscadelli, A. Lorenzani, A. Ahmadi, H. Beuther, T. Csengeri, S. Etoka, C. Goddi, P. D. Klaassen, R. Kuiper, M. S. N. Kumar, T. Peters, Á. Sánchez-Monge, P. Schilke, F. van der Tak, S. Vig, H. Zinnecker

Abstract: Context. ALMA observations at 1.4 mm and 0.2'' (750au) angular resolution of the Main core in the high-mass star forming region G31.41+0.31 have revealed a puzzling scenario: on the one hand, the continuum emission looks very homogeneous and the core appears to undergo solid-body rotation, suggesting a monolithic core stabilized by the magnetic field; on the other hand, rotation and infall speed u… ▽ More Context. ALMA observations at 1.4 mm and 0.2'' (750au) angular resolution of the Main core in the high-mass star forming region G31.41+0.31 have revealed a puzzling scenario: on the one hand, the continuum emission looks very homogeneous and the core appears to undergo solid-body rotation, suggesting a monolithic core stabilized by the magnetic field; on the other hand, rotation and infall speed up toward the core center, where two massive embedded free-free continuum sources have been detected, pointing to an unstable core having undergone fragmentation. Aims. To establish whether the Main core is indeed monolithic or its homogeneous appearance is due to a combination of large dust opacity and low angular resolution, we carried out millimeter observations at higher angular resolution and different wavelengths. Methods. We carried out ALMA observations at 1.4 mm and 3.5 mm that achieved angular resolutions of 0.1''(375 au) and 0.075'' (280 au), respectively. VLA observations at 7 mm and 1.3 cm at even higher angular resolution, 0.05'' (190 au) and 0.07'' (260 au), respectively, were also carried out to better study the nature of the free-free continuum sources detected in the core. Results. The millimeter continuum emission of the Main core has been clearly resolved into at least four sources, A, B, C, and D, within 1'', indicating that the core is not monolithic. The deconvolved radii of the dust emission of the sources, estimated at 3.5 mm, are 400-500au, their masses range from 15 to 26 Msun, and their number densities are several 1E9 cm-3. Sources A and B, located closer to the center of the core and separated by 750 au, are clearly associated with two free-free continuum sources, likely thermal radio jets, and are the brightest in the core. The spectral energy distribution of these two sources and their masses and sizes are similar and suggest a common origin. △ Less

Submitted 8 March, 2021; originally announced March 2021.

Comments: 18 pages,9 figures, 5 tables. Accepted by A&A

Journal ref: A&A 648, A100 (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 prebiotic molecular inventory of Serpens SMM1 I. An investigation of the isomers CH$_{3}$NCO and HOCH$_{2}$CN

Authors: N. F. W. Ligterink, A. Ahmadi, A. Coutens, Ł. Tychoniec H. Calcutt, E. F. van Dishoeck, H. Linnartz, J. K. Jørgensen, R. T. Garrod, J. Bouwman

Abstract: Methyl isocyanate (CH$_{3}$NCO) and glycolonitrile (HOCH$_{2}$CN) are isomers and prebiotic molecules that are involved in the formation of peptide structures and the nucleobase adenine, respectively. ALMA observations of the intermediate-mass Class 0 protostar Serpens SMM1-a and ALMA-PILS data of the low-mass Class 0 protostar IRAS~16293B are used. Spectra are analysed with the CASSIS line analys… ▽ More Methyl isocyanate (CH$_{3}$NCO) and glycolonitrile (HOCH$_{2}$CN) are isomers and prebiotic molecules that are involved in the formation of peptide structures and the nucleobase adenine, respectively. ALMA observations of the intermediate-mass Class 0 protostar Serpens SMM1-a and ALMA-PILS data of the low-mass Class 0 protostar IRAS~16293B are used. Spectra are analysed with the CASSIS line analysis software package in order to identify and characterise molecules. CH$_{3}$NCO, HOCH$_{2}$CN, and various other molecules are detected towards SMM1-a. HOCH$_{2}$CN is identified in the PILS data towards IRAS~16293B in a spectrum extracted at a half-beam offset position from the peak continuum. CH$_{3}$NCO and HOCH$_{2}$CN are equally abundant in SMM1-a at [X]/[CH$_{3}$OH] of 5.3$\times$10$^{-4}$ and 6.2$\times$10$^{-4}$, respectively. A comparison between SMM1-a and IRAS~16293B shows that HOCH$_{2}$CN and HNCO are more abundant in the former source, but CH$_{3}$NCO abundances do not differ significantly. Data from other sources are used to show that the [CH$_{3}$NCO]/[HNCO] ratio is similar in all these sources within $\sim$10\%. The new detections of CH$_{3}$NCO and HOCH$_{2}$CN are additional evidence for a large interstellar reservoir of prebiotic molecules that can contribute to the formation of biomolecules on terrestrial planets. A plausible formation pathway for HOCH$_{2}$CN is the thermal Strecker-like reaction of CN$^{-}$ with H$_{2}$CO. The similar [CH$_{3}$NCO]/[HNCO] ratios indicate that these two species either are chemically related or their formation is affected by physical conditions in the same way. The relatively high abundances of HOCH$_{2}$CN and HNCO in SMM1-a may be explained by a prolonged stage of relatively warm ice mantles, where thermal and energetic processing of HCN in the ice results in the efficient formation of both species. △ Less

Submitted 31 December, 2020; originally announced December 2020.

Comments: Accepted for publication in Astronomy and Astrophysics. 32 pages, 26 figures

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

Does the magnetic field suppress fragmentation in massive dense cores?

Authors: Aina Palau, Qizhou Zhang, Josep M. Girart, Junhao Liu, Ram Rao, Patrick M. Koch, Robert Estalella, Huei-Ru Vivien Chen, Hauyu Baobab Liu, Keping Qiu, Zhi-Yun Li, Luis A. Zapata, Sylvain Bontemps, Paul T. P. Ho, Henrik Beuther, Tao-Chung Ching, Hiroko Shinnaga, Aida Ahmadi

Abstract: Theoretical and numerical works indicate that a strong magnetic field should suppress fragmentation in dense cores. However, this has never been tested observationally in a relatively large sample of fragmenting massive dense cores. Here we use the polarization data obtained in the Submillimeter Array Legacy Survey of Zhang et al. to build a sample of 18 massive dense cores where both fragmentatio… ▽ More Theoretical and numerical works indicate that a strong magnetic field should suppress fragmentation in dense cores. However, this has never been tested observationally in a relatively large sample of fragmenting massive dense cores. Here we use the polarization data obtained in the Submillimeter Array Legacy Survey of Zhang et al. to build a sample of 18 massive dense cores where both fragmentation and magnetic field properties are studied in a uniform way. We measured the fragmentation level, Nmm, within the field of view common to all regions, of 0.15 pc, with a mass sensitivity of about 0.5 Msun, and a spatial resolution of about 1000 au. In order to obtain the magnetic field strength using the Davis-Chandrasekhar-Fermi method, we estimated the dispersion of the polarization position angles, the velocity dispersion of the H13CO+(4-3) gas, and the density of each core, all averaged within 0.15 pc. A strong correlation is found between Nmm and the average density of the parental core, although with significant scatter. When large-scale systematic motions are separated from the velocity dispersion and only the small-scale (turbulent) contribution is taken into account, a tentative correlation is found between Nmm and the mass-to-flux ratio, as suggested by numerical and theoretical works. △ Less

Submitted 9 May, 2021; v1 submitted 22 October, 2020; originally announced October 2020.

Comments: accepted in ApJ

Gravity and rotation drag the magnetic field in high-mass star formation

Authors: H. Beuther, J. D. Soler, H. Linz, Th. Henning, C. Gieser, R. Kuiper, W. Vlemmings, P. Hennebelle, S. Feng, R. Smith, A. Ahmadi

Abstract: The formation of hot stars out of the cold interstellar medium lies at the heart of astrophysical research. Understanding the importance of magnetic fields during star formation remains a major challenge. With the advent of the Atacama Large Millimeter Array, the potential to study magnetic fields by polarization observations has tremendously progressed. However, the major question remains how muc… ▽ More The formation of hot stars out of the cold interstellar medium lies at the heart of astrophysical research. Understanding the importance of magnetic fields during star formation remains a major challenge. With the advent of the Atacama Large Millimeter Array, the potential to study magnetic fields by polarization observations has tremendously progressed. However, the major question remains how much magnetic fields shape the star formation process or whether gravity is largely dominating. Here, we show that for the high-mass star-forming region G327.3 the magnetic field morphology appears to be dominantly shaped by the gravitational contraction of the central massive gas core where the star formation proceeds. We find that in the outer parts of the region, the magnetic field is directed toward the gravitational center of the region. Filamentary structures feeding the central core exhibit U-shaped magnetic field morphologies directed toward the gravitational center as well, again showing the gravitational drag toward the center. The inner part then shows rotational signatures, potentially associated with an embedded disk, and there the magnetic field morphology appears to be rotationally dominated. Hence, our results demonstrate that for this region gravity and rotation are dominating the dynamics and shaping the magnetic field morphology. △ Less

Submitted 12 October, 2020; originally announced October 2020.

Comments: 10 pages, 4 figures, accepted for the Astrophysical Journal, also available at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

Multi-wavelength modelling of the circumstellar environment of the massive proto-star AFGL 2591 VLA 3

Authors: F. A. Olguin, M. G. Hoare, K. G. Johnston, F. Motte, H. -R. V. Chen, H. Beuther, J. C. Mottram, A. Ahmadi, C. Gieser, D. Semenov, T. Peters, A. Palau, P. D. Klaassen, R. Kuiper, Á. Sánchez-Monge, Th. Henning

Abstract: We have studied the dust density, temperature and velocity distributions of the archetypal massive young stellar object (MYSO) AFGL 2591. Given its high luminosity ($L=2 \times 10^5$ L$_\odot$) and distance ($d=3.3$ kpc), AFGL 2591 has one of the highest $\sqrt{L}/d$ ratio, giving better resolved dust emission than any other MYSO. As such, this paper provides a template on how to use resolved mult… ▽ More We have studied the dust density, temperature and velocity distributions of the archetypal massive young stellar object (MYSO) AFGL 2591. Given its high luminosity ($L=2 \times 10^5$ L$_\odot$) and distance ($d=3.3$ kpc), AFGL 2591 has one of the highest $\sqrt{L}/d$ ratio, giving better resolved dust emission than any other MYSO. As such, this paper provides a template on how to use resolved multi-wavelength data and radiative transfer to obtain a well-constrained 2-D axi-symmetric analytic rotating infall model. We show for the first time that the resolved dust continuum emission from Herschel 70 $μ$m observations is extended along the outflow direction, whose origin is explained in part from warm dust in the outflow cavity walls. However, the model can only explain the kinematic features from CH$_3$CN observations with unrealistically low stellar masses ($<15$ M$_\odot$), indicating that additional physical processes may be playing a role in slowing down the envelope rotation. As part of our 3-step continuum and line fitting, we have identified model parameters that can be further constrained by specific observations. High-resolution mm visibilities were fitted to obtain the disc mass (6 M$_\odot$) and radius (2200 au). A combination of SED and near-IR observations were used to estimate the luminosity and envelope mass together with the outflow cavity inclination and opening angles. △ Less

Submitted 12 August, 2020; originally announced August 2020.

Comments: 34 pages, 22 figures, 15 tables, 4 appendices. Accepted for publication in MNRAS

Spiral arms and instability within the AFGL 4176 mm1 disc

Authors: K. G. Johnston, M. G. Hoare, H. Beuther, R. Kuiper, N. D. Kee, H. Linz, P. Boley, L. T. Maud, A. Ahmadi, T. P. Robitaille

Abstract: We present high-resolution (30 mas or 130 au at 4.2 kpc) Atacama Large Millimeter/submillimeter Array observations at 1.2 mm of the disc around the forming O-type star AFGL 4176 mm1. The disc (AFGL 4176 mm1-main) has a radius of ~1000 au and contains significant structure, most notably a spiral arm on its redshifted side. We fitted the observed spiral with logarithmic and Archimedean spiral models… ▽ More We present high-resolution (30 mas or 130 au at 4.2 kpc) Atacama Large Millimeter/submillimeter Array observations at 1.2 mm of the disc around the forming O-type star AFGL 4176 mm1. The disc (AFGL 4176 mm1-main) has a radius of ~1000 au and contains significant structure, most notably a spiral arm on its redshifted side. We fitted the observed spiral with logarithmic and Archimedean spiral models. We find that both models can describe its structure, but the Archimedean spiral with a varying pitch angle fits its morphology marginally better. As well as signatures of rotation across the disc, we observe gas arcs in CH$_3$CN that connect to other millimetre continuum sources in the field, supporting the picture of interactions within a small cluster around AFGL 4176 mm1-main. Using local thermodynamic equilibrium modelling of the CH$_3$CN K-ladder, we determine the temperature and velocity field across the disc, and thus produce a map of the Toomre stability parameter. Our results indicate that the outer disc is gravitationally unstable and has already fragmented or is likely to fragment in the future, possibly producing further companions. These observations provide evidence that disc fragmentation is one possible pathway towards explaining the high fraction of multiple systems around high-mass stars. △ Less

Submitted 18 February, 2020; v1 submitted 21 November, 2019; originally announced November 2019.

Comments: Accepted by Astronomy & Astrophysics Letters, 8 pages (2 pages appendices), 11 figures

Journal ref: A&A 634, L11 (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)

Disk Kinematics and Stability in High-Mass Star Formation: Linking Simulations and Observations

Authors: Aida Ahmadi, Rolf Kuiper, Henrik Beuther

Abstract: In the disk-mediated accretion scenario for the formation of the most massive stars, gravitational instabilities in the disk can force it to fragment. We investigate the effects of inclination and spatial resolution on observable kinematics and stability of disks in high-mass star formation. We study a high-resolution 3D radiation-hydrodynamic simulation that leads to the fragmentation of a massiv… ▽ More In the disk-mediated accretion scenario for the formation of the most massive stars, gravitational instabilities in the disk can force it to fragment. We investigate the effects of inclination and spatial resolution on observable kinematics and stability of disks in high-mass star formation. We study a high-resolution 3D radiation-hydrodynamic simulation that leads to the fragmentation of a massive disk. Using RADMC-3D we produce 1.3 mm continuum and CH3CN line cubes at different inclinations. The model is set to different distances and synthetic observations are created for ALMA at ~80 mas resolution and NOEMA at ~0.3''. The synthetic ALMA observations resolve all fragments and their kinematics well. The synthetic NOEMA observations at 800 pc (~300 au resolution) are able to resolve the fragments, while at 2000 pc (~800 au resolution) only a single slightly elongated structure is observed. The position-velocity (PV) plots show the differential rotation of material best in the edge-on views. As the observations become less resolved, the inner high-velocity components of the disk become blended with the envelope and the PV plots resemble rigid-body-like rotation. Protostellar mass estimates from PV plots of poorly resolved observations are therefore overestimated. We fit the emission of CH3CN lines and produce maps of gas temperature with values in the range of 100-300 K. Studying the Toomre stability of the disks in the resolved observations, we find Q values below the critical value for stability against gravitational collapse at the positions of the fragments and the arms connecting the fragments. For the poorly resolved observations we find low Q values in the outskirts of the disk. Therefore we are able to predict that the disk is unstable and fragmenting even in poorly resolved observations. This conclusion is true regardless of knowledge about the inclination of the disk. △ Less

Submitted 15 October, 2019; v1 submitted 9 September, 2019; originally announced September 2019.

Comments: 21 pages, 15 figures, accepted for publication in Astronomy and Astrophysics - corrected captions of Fig. 13 & 15

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

Fragmentation, rotation and outflows in the high-mass star-forming region IRAS 23033+5951. A case study of the IRAM NOEMA large program CORE

Authors: F. Bosco, H. Beuther, A. Ahmadi, J. C. Mottram, R. Kuiper, H. Linz, L. Maud, J. M. Winters, T. Henning, S. Feng, T. Peters, D. Semenov, P. D. Klaassen, P. Schilke, J. S. Urquhart, M. T. Beltrán, S. L. Lumsden, S. Leurini, L. Moscadelli, R. Cesaroni, Á. Sánchez-Monge, A. Palau, R. Pudritz, F. Wyrowski, S. Longmore , et al. (1 additional authors not shown)

Abstract: The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. U… ▽ More The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. Using the Northern Extended Millimeter Array (NOEMA) in three configurations and the IRAM 30-m single-dish telescope at 220GHz, we probe the gas and dust emission at an angular resolution of ~0.45arcsec, corresponding to 1900au. In the mm continuum emission, we identify a protostellar cluster with at least four mm-sources, where three of them show a significantly higher peak intensity well above a signal-to-noise ratio of 100. Hierarchical fragmentation from large to small spatial scales is discussed. Two fragments are embedded in rotating structures and drive molecular outflows, traced by $^{13}$CO (2-1) emission. The velocity profiles across two of the cores are similar to Keplerian but are missing the highest velocity components close to the center of rotation, which is a common phenomena from observations like these, and other rotation scenarios are not excluded entirely. Position-velocity diagrams suggest protostellar masses of ~6 and 19M$_\sun$. Rotational temperatures from fitting CH$_3$CN ($12_K-11_K$) spectra are used for estimating the gas temperature and by that the disk stability against gravitational fragmentation, utilizing Toomre's $Q$ parameter. [We] identify only one candidate disk to be unstable against gravitational instability caused by axisymmetric perturbations. The dominant sources cover different evolutionary stages within the same maternal gas clump. The appearance of rotation and outflows of the cores are similar to those found in low-mass star-forming regions. △ Less

Submitted 9 July, 2019; originally announced July 2019.

Comments: 23 pages, 16 figures, 5 tables, accepted for publication in A&A

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

Substructures in the Keplerian disc around the O-type (proto)star G17.64+0.16

Authors: L. T. Maud, R. Cesaroni, M. S. N. Kumar, V. M. Rivilla, A. Ginsburg, P. D. Klaassen, D. Harsono, A. Sanchez-Monge, A. Ahmadi, V. Allen, M. T. Beltran, H. Beuther, R. Galvan-Madrid, C. Goddi, M. G. Hoare, M. R. Hogerheijde, K. G. Johnston, R. Kuiper, L. Moscadelli, T. Peters, L. Testi, F. F. S. van der Tak, W. J. de Wit

Abstract: We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission… ▽ More We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission, especially visible as arc structures to the north and south. The Keplerian kinematics are most prominently seen in the vibrationally excited water line, H2O (Eu=3461.9K). The mass of the central source found by modelling the Keplerian rotation is consistent with 45+/-10Mo. The H30alpha (231.9GHz) radio-recombination line and the SiO (5-4) molecular line were detected at up to the 10 sigma$ level. The estimated disc mass is 0.6-2.6Mo under the optically thin assumption. Analysis of the Toomre Q parameter, in the optically thin regime, indicates that the disc stability is highly dependent on temperature. The disc currently appears stable for temperatures >150K, this does not preclude that the substructures formed earlier through disc fragmentation. △ Less

Submitted 15 June, 2019; originally announced June 2019.

Comments: 9 Total, inc references, appendix, long author list, 4 figures (inc appendix)

Journal ref: A&A 627, L6 (2019)

IRAS23385+6053: An embedded massive cluster in the making

Authors: R. Cesaroni, H. Beuther, A. Ahmadi, M. T. Beltran, T. Csengeri, R. Galvan-Madrid, C. Gieser, T. Henning, K. G. Johnston, P. D. Klaassen, R. Kuiper, S. Leurini, H. Linz, S. Longmore, S. L. Lumsden, L. T. Maud, L. Moscadelli, J. C. Mottram, A. Palau, T. Peters, R. E. Pudritz, A. Sanchez-Monge, P. Schilke, D. Semenov, S. Suri , et al. (4 additional authors not shown)

Abstract: This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additi… ▽ More This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additional single-dish maps, merged with the interferometric data to recover the extended emission. Our correlator setup covered a number of lines from well-known hot core tracers and a few outflow tracers. The angular (~0.45"$-$0.9") and spectral (0.5 km/s) resolutions were sufficient to resolve the clump in IRAS23385+6053 and investigate the existence of large-scale motions due to rotation, infall, or expansion. We find that the clump splits into six distinct cores when observed at sub-arcsecond resolution. These are identified through their 1.4 mm continuum and molecular line emission. We produce maps of the velocity, line width, and rotational temperature from the methanol and methyl cyanide lines, which allow us to investigate the cores and reveal a velocity and temperature gradient in the most massive core. We also find evidence of a bipolar outflow, possibly powered by a low-mass star. We present the tentative detection of a circumstellar self-gravitating disk lying in the most massive core and powering a large-scale outflow previously known in the literature. In our scenario, the star powering the flow is responsible for most of the luminosity of IRAS23385+6053 (~$3000~L_\odot$). The other cores, albeit with masses below the corresponding virial masses, appear to be accreting material from their molecular surroundings and are possibly collapsing or on the verge of collapse. We conclude that we are observing a sample of star-forming cores that is bound to turn into a cluster of massive stars. △ Less

Submitted 27 May, 2019; originally announced May 2019.

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

Information gains from Monte Carlo Markov Chains

Authors: Ahmad Mehrabi, A. Ahmadi

Abstract: In this paper, we present a novel method for computing the relative entropy as well as the expected relative entropy using an MCMC chain. The relative entropy from information theory can be used to quantify differences in posterior distributions of a pair of experiments. In cosmology, the relative entropy has been proposed as an interesting tool for model selection, experiment design, forecasting… ▽ More In this paper, we present a novel method for computing the relative entropy as well as the expected relative entropy using an MCMC chain. The relative entropy from information theory can be used to quantify differences in posterior distributions of a pair of experiments. In cosmology, the relative entropy has been proposed as an interesting tool for model selection, experiment design, forecasting and measuring information gain from subsequent experiments. In contrast to Gaussian distributions, these quantities are not generally available analytically and one needs to use numerical methods to estimate them which are certainly computationally expensive. We propose a method and provide its python package to estimate the relative entropy as well as expected relative entropy from a posterior sample. We consider the linear Gaussian model to check the accuracy of our code. Our results indicate that the relative error is below $0.2\%$ for sample size larger than $10^5$ in the linear Gaussian model. In addition, we study the robustness of our code in estimating the expected relative entropy in this model. △ Less

Submitted 26 April, 2019; originally announced April 2019.

Comments: The code is available in https://github.com/ahmadiphy/MCKLdivergence

Chasing discs around O-type (proto)stars - ALMA evidence for an SiO disc and disc wind from G17.64+0.16

Authors: L. T. Maud, R. Cesaroni, M. S. N. Kumar, F. F. S. van der Tak, V. Allen, M. G. Hoare, P. D. Klaassen, D. Harsono, M. R. Hogerheijde, Á. Sánchez-Monge, P. Schilke, A. Ahmadi, M. T. Beltrán, H. Beuther, T. Csengeri, S. Etoka, G. Fuller, R. Galván-Madrid, C. Goddi, Th. Henning, K. G. Johnston, R. Kuiper, S. Lumsden, L. Moscadelli, J. C. Mottram , et al. (6 additional authors not shown)

Abstract: We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peak… ▽ More We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. We detect SiO (5-4) emission that is marginally resolved and elongated in a direction perpendicular to the large-scale outflow seen in the 13CO (2-1) line using the main ALMA array in conjunction with the Atacama Compact Array (ACA). Morphologically, the SiO appears to represent a disc-like structure. Using parametric models we show that the position-velocity profile of the SiO is consistent with the Keplerian rotation of a disc around an object between 10-30Mo in mass, only if there is also radial expansion from a separate structure. The radial motion component can be interpreted as a disc wind from the disc surface. Models with a central stellar object mass between 20 and 30Mo are the most consistent with the stellar luminosity (100000 Lo) and indicative of an O-type star. The H30a millimetre recombination line (231.9GHz) is also detected, but spatially unresolved, and is indicative of a very compact, hot, ionised region co-spatial with the dust continuum core. Accounting for all observables, we suggest that G17.64 is consistent with a O-type young stellar object in the final stages of protostellar assembly, driving a wind, but that has not yet developed into a compact HII region. The existance and detection of the disc in G17.64 is likely related to its isolated and possibly more evolved nature, traits which may underpin discs in similar sources. △ Less

Submitted 9 October, 2018; originally announced October 2018.

Comments: 23 pages, 13 fig (inc 3 Appendix figures)

Journal ref: A&A 620, A31 (2018)

Core fragmentation and Toomre stability analysis of W3(H2O): A case study of the IRAM NOEMA large program CORE

Authors: A. Ahmadi, H. Beuther, J. C. Mottram, F. Bosco, H. Linz, Th. Henning, J. M. Winters, R. Kuiper, R. Pudritz, Á. Sánchez-Monge, E. Keto, M. Beltran, S. Bontemps, R. Cesaroni, T. Csengeri, S. Feng, R. Galvan-Madrid, K. G. Johnston, P. Klaassen, S. Leurini, S. N. Longmore, S. Lumsden, L. T. Maud, K. M. Menten, L. Moscadelli , et al. (8 additional authors not shown)

Abstract: The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.3… ▽ More The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.37 mm wavelength regime in both line and dust continuum emission. We present the spectral line setup of the CORE survey and a case study for W3(H2O). At ~0.35" (700 AU at 2 kpc) resolution, the W3(H2O) clump fragments into two cores (West and East), separated by ~2300 AU. Velocity shifts of a few km/s are observed in the dense-gas tracer, CH3CN, across both cores, consistent with rotation and perpendicular to the directions of two bipolar outflows, one emanating from each core. The kinematics of the rotating structure about W3(H2O) W shows signs of differential rotation of material, possibly in a disk-like object. The observed rotational signature around W3(H2O) E may be due to a disk-like object, an unresolved binary (or multiple) system, or a combination of both. We fit the emission of CH3CN (12-11) K = 4-6 and derive a gas temperature map with a median temperature of ~165 K across W3(H2O). We create a Toomre Q map to study the stability of the rotating structures against gravitational instability. The rotating structures appear to be Toomre unstable close to their outer boundaries, with a possibility of further fragmentation in the differentially-rotating core W3(H2O) W. Rapid cooling in the Toomre-unstable regions supports the fragmentation scenario. Combining millimeter dust continuum and spectral line data toward the famous high-mass star-forming region W3(H2O), we identify core fragmentation on large scales, and indications for possible disk fragmentation on smaller spatial scales. △ Less

Submitted 1 August, 2018; originally announced August 2018.

Comments: 23 pages, 26 figures, accepted for publication in Astronomy and Astrophysics

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

Fragmentation and disk formation during high-mass star formation: The IRAM NOEMA (Northern Extended Millimeter Array) large program CORE

Authors: H. Beuther, J. C. Mottram, A. Ahmadi, F. Bosco, H. Linz, Th. Henning, P. Klaassen, J. M. Winters, L. T. Maud, R. Kuiper, D. Semenov, C. Gieser, T. Peters, J. S. Urquhart, R. Pudritz, S. E. Ragan, S. Feng, E. Keto, S. Leurini, R. Cesaroni, M. Beltran, A. Palau, A. Sanchez-Monge, R. Galvan-Madrid, Q. Zhang , et al. (8 additional authors not shown)

Abstract: Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high… ▽ More Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high angular resolution (~0.4'') for a sample of 20 well-known high-mass star-forming regions with distances below 5.5kpc and luminosities larger than 10^4Lsun. Results: We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000AU. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions: The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation. △ Less

Submitted 3 May, 2018; originally announced May 2018.

Comments: 29 pages, 17 figures, Astronomy & Astrophysics in press, for a higher-resolution version of the paper see http://www.mpia.de/homes/beuther/papers.html

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

A New Look at an Old Cluster: The Membership, Rotation, and Magnetic Activity of Low-Mass Stars in the 1.3-Gyr-Old Open Cluster NGC 752

Authors: Marcel Agüeros, Emily Bowsher, John Bochanski, Phill Cargile, Kevin Covey, Stephanie Douglas, Adam Kraus, Alisha Kundert, Nicholas Law, Aida Ahmadi, Héctor Arce

Abstract: The nearby open cluster NGC 752 presents a rare opportunity to study stellar properties at ages >1 Gyr. However, constructing a membership catalog for it is challenging; most surveys have been limited to identifying its giants and dwarf members earlier than mid-K. We supplement past membership catalogs with candidates selected with updated photometric and proper-motion criteria, generating a list… ▽ More The nearby open cluster NGC 752 presents a rare opportunity to study stellar properties at ages >1 Gyr. However, constructing a membership catalog for it is challenging; most surveys have been limited to identifying its giants and dwarf members earlier than mid-K. We supplement past membership catalogs with candidates selected with updated photometric and proper-motion criteria, generating a list of 258 members, a >50% increase over previous catalogs. Using a Bayesian framework to fit MESA Isochrones & Stellar Tracks evolutionary models to literature photometry and the Tycho-Gaia Astrometric Solution data available for 59 cluster members, we infer the age of, and distance to, NGC 752: 1.34$\pm$0.06 Gyr and 438$_{-6}^{+8}$ pc. We also report the results of our optical monitoring of the cluster using the Palomar Transient Factory. We obtain rotation periods for 12 K and M cluster members, the first periods measured for such low-mass stars with a well-constrained age >1 Gyr. We compare these new periods to data from the younger clusters Praesepe and NGC 6811, and to a theoretical model for angular-momentum loss, to examine stellar spin down for low-mass stars over their first 1.3 Gyr. While on average NGC 752 stars are rotating more slowly than their younger counterparts, the difference is not significant. Finally, we use our spectroscopic observations to measure Halpha for cluster stars, finding that members earlier than $\approx$M2 are magnetically inactive, as expected at this age. Forthcoming Gaia data should solidify and extend the membership of NGC 752 to lower masses, thereby increasing its importance for studies of low-mass stars. △ Less

Submitted 5 April, 2018; originally announced April 2018.

Comments: Accepted to ApJ. This version 23 pages, 15 figures; on-line version will include 11 more figures as well as light curve data for the 12 rotators we identify

Journal ref: ApJ 2018 862 33

Accelerating infall and rotational spin-up in the hot molecular core G31.41+0.31

Authors: M. T. Beltrán, R. Cesaroni, V. M. Rivilla, Á. Sánchez-Monge, L. Moscadelli, A. Ahmadi, V. Allen, H. Beuther, S. Etoka, D. Galli, R. Galván-Madrid, C. Goddi, K. G. Johnston, A. Kölligan, R. Kuiper, M. S. N. Kumar, L. T. Maud, J. C. Mottram, T. Peters, P. Schilke, L. Testi, F. van der Tak, C. M. Walmsley

Abstract: As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41+0.31 with ALMA at 1.4 mm with an angular resolution of~0.22" (~1700 au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is as… ▽ More As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41+0.31 with ALMA at 1.4 mm with an angular resolution of~0.22" (~1700 au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is associated with two free-free continuum sources. The Main core looks featureless and homogeneous in dust continuum emission and does not present any hint of fragmentation. Each transition of CH3CN and CH3OCHO, both ground and vibrationally excited, as well as those of CH3CN isotopologues, shows a clear velocity gradient along the NE-SW direction, with velocity linearly increasing with distance from the center, consistent with solid-body rotation. However, when comparing the velocity field of transitions with different upper level energies, the rotation velocity increases with increasing energy of the transition, which suggests that the rotation speeds up towards the center. Spectral lines towards the dust continuum peak show an inverse P-Cygni profile that supports the existence of infall in the core. The infall velocity increases with the energy of the transition suggesting that the infall is accelerating towards the center of the core, consistent with gravitational collapse. Despite the monolithic appearance of the Main core, the presence of red-shifted absorption, the existence of two embedded free-free sources at the center, and the rotational spin-up are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum. Therefore, the most likely explanation for the monolithic morphology is that the large opacity of the dust emission prevents the detection of any inhomogeneity in the core. △ Less

Submitted 21 June, 2018; v1 submitted 14 March, 2018; originally announced March 2018.

Comments: 28 pages, 16 figures, and 7 tables. Accepted for publication in A&A

Why are rapidly rotating M dwarfs in the Pleiades so (infra)red? New period measurements confirm rotation-dependent color offsets from the cluster sequence

Authors: Kevin R. Covey, Marcel A. Agüeros, Nicholas M. Law, Jiyu Liu, Aida Ahmadi, Russ Laher, David Levitan, Branimir Sesar, Jason Surace

Abstract: Stellar rotation periods measured in open clusters have proved to be extremely useful for studying stars' angular momentum content and rotationally driven magnetic activity, which are both age- and mass-dependent processes. While period measurements have been obtained for hundreds of solar-mass members of the Pleiades, period measurements exist for only a few low-mass ($<$0.5 M$_{\odot}$) members… ▽ More Stellar rotation periods measured in open clusters have proved to be extremely useful for studying stars' angular momentum content and rotationally driven magnetic activity, which are both age- and mass-dependent processes. While period measurements have been obtained for hundreds of solar-mass members of the Pleiades, period measurements exist for only a few low-mass ($<$0.5 M$_{\odot}$) members of this key laboratory for stellar evolution theory. To fill this gap, we report rotation periods for 132 low-mass Pleiades members (including nearly 100 with M $\leq$ 0.45 M$_{\odot}$), measured from photometric monitoring of the cluster conducted by the Palomar Transient Factory in late 2011 and early 2012. These periods extend the portrait of stellar rotation at 125 Myr to the lowest-mass stars and re-establish the Pleiades as a key benchmark for models of the transport and evolution of stellar angular momentum. Combining our new rotation periods with precise BVIJHK photometry reported by Stauffer et al. and Kamai et al., we investigate known anomalies in the photometric properties of K and M Pleiades members. We confirm the correlation detected by Kamai et al. between a star's rotation period and position relative to the main sequence in the cluster's color-magnitude diagram. We find that rapid rotators have redder (V-K) colors than slower rotators at the same V, indicating that rapid and slow rotators have different binary frequencies and/or photospheric properties. We find no difference in the photometric amplitudes of rapid and slow rotators, indicating that asymmetries in the longitudinal distribution of starspots do not scale grossly with rotation rate. △ Less

Submitted 26 February, 2016; v1 submitted 26 January, 2016; originally announced January 2016.

Comments: Accepted for publication in the Astrophysical Journal; 27 pages, 24 figures, 5 tables

Exoplanets from the Arctic: The First Wide-Field Survey at 80 Degrees North

Authors: Nicholas M. Law, Raymond Carlberg, Pegah Salbi, Wai-Hin Wayne Ngan, Aida Ahmadi, Eric Steinbring, Richard Murowinski, Suresh Sivanandam, Wolfgang Kerzendorf

Abstract: Located within 10 degrees of the North Pole, northern Ellesmere Island offers continuous darkness in the winter months. This capability can greatly enhance the detection efficiency of planetary transit surveys and other time domain astronomy programs. We deployed two wide-field cameras at 80 degrees North, near Eureka, Nunavut, for a 152-hour observing campaign in February 2012. The 16-megapixel-c… ▽ More Located within 10 degrees of the North Pole, northern Ellesmere Island offers continuous darkness in the winter months. This capability can greatly enhance the detection efficiency of planetary transit surveys and other time domain astronomy programs. We deployed two wide-field cameras at 80 degrees North, near Eureka, Nunavut, for a 152-hour observing campaign in February 2012. The 16-megapixel-camera systems were based on commercial f/1.2 lenses with 70mm and 42mm apertures, and they continuously imaged 504 and 1,295 square degrees respectively. In total, the cameras took over 44,000 images and produced better-than-1% precision light curves for approximately 10,000 stars. We describe a new high-speed astrometric and photometric data reduction pipeline designed for the systems, test several methods for the precision flat-fielding of images from very-wide-angle cameras, and evaluate the cameras' image qualities. We achieved a scintillation-limited photometric precision of 1-2% in each 10s exposure. Binning the short exposures into 10-min chunks provided a photometric stability of 2-3 millimagnitudes, sufficient for the detection of transiting exoplanets around the bright stars targeted by our survey. We estimate that the cameras, when operated over the full arctic winter, will be capable of discovering several transiting exoplanets around bright (V<9) stars. △ Less

Submitted 8 November, 2012; originally announced November 2012.

Comments: 10 pages, 11 figures, AJ in press

New Exoplanet Surveys in the Canadian High Arctic at 80 Degrees North

Authors: Nicholas M. Law, Suresh Sivanandam, Richard Murowinski, Raymond Carlberg, Wayne Ngan, Pegah Salbi, Aida Ahmadi, Eric Steinbring, Mark Halman, James Graham

Abstract: Observations from near the Eureka station on Ellesmere Island, in the Canadian High Arctic at 80 degrees North, benefit from 24-hour darkness combined with dark skies and long cloud-free periods during the winter. Our first astronomical surveys conducted at the site are aimed at transiting exoplanets; compared to mid-latitude sites, the continuous darkness during the Arctic winter greatly improves… ▽ More Observations from near the Eureka station on Ellesmere Island, in the Canadian High Arctic at 80 degrees North, benefit from 24-hour darkness combined with dark skies and long cloud-free periods during the winter. Our first astronomical surveys conducted at the site are aimed at transiting exoplanets; compared to mid-latitude sites, the continuous darkness during the Arctic winter greatly improves the survey's detection efficiency for longer-period transiting planets. We detail the design, construction, and testing of the first two instruments: a robotic telescope, and a set of very wide-field imaging cameras. The 0.5m Dunlap Institute Arctic Telescope has a 0.8-square-degree field of view and is designed to search for potentially habitable exoplanets around low-mass stars. The very wide field cameras have several-hundred-square-degree fields of view pointed at Polaris, are designed to search for transiting planets around bright stars, and were tested at the site in February 2012. Finally, we present a conceptual design for the Compound Arctic Telescope Survey (CATS), a multiplexed transient and transit search system which can produce a 10,000-square-degree snapshot image every few minutes throughout the Arctic winter. △ Less

Submitted 28 August, 2012; originally announced August 2012.

Comments: 11 pages, 6 figures, SPIE vol 8444, 2012