Showing 1–41 of 41 results for author: Maud, L T

FAUST XVII: Super deuteration in the planet forming system IRS 63 where the streamer strikes the disk

Authors: L. Podio, C. Ceccarelli, C. Codella, G. Sabatini, D. Segura-Cox, N. Balucani, A. Rimola, P. Ugliengo, C. J. Chandler, N. Sakai, B. Svoboda, J. Pineda, M. De Simone, E. Bianchi, P. Caselli, A. Isella, Y. Aikawa, M. Bouvier, E. Caux, L. Chahine, S. B. Charnley, N. Cuello, F. Dulieu, L. Evans, D. Fedele , et al. (33 additional authors not shown)

Abstract: Recent observations suggest that planets formation starts early, in protostellar disks of $\le10^5$ yrs, which are characterized by strong interactions with the environment, e.g., through accretion streamers and molecular outflows. To investigate the impact of such phenomena on disk physical and chemical properties it is key to understand what chemistry planets inherit from their natal environment… ▽ More Recent observations suggest that planets formation starts early, in protostellar disks of $\le10^5$ yrs, which are characterized by strong interactions with the environment, e.g., through accretion streamers and molecular outflows. To investigate the impact of such phenomena on disk physical and chemical properties it is key to understand what chemistry planets inherit from their natal environment. In the context of the ALMA Large Program Fifty AU STudy of the chemistry in the disk/envelope system of Solar-like protostars (FAUST), we present observations on scales from ~1500 au to ~60 au of H$_2$CO, HDCO, and D$_2$CO towards the young planet-forming disk IRS~63. H$_2$CO probes the gas in the disk as well as in a large scale streamer (~1500 au) impacting onto the South-East (SE) disk side. We detect for the first time deuterated formaldehyde, HDCO and D$_2$CO, in a planet-forming disk, and HDCO in the streamer that is feeding it. This allows us to estimate the deuterium fractionation of H$_2$CO in the disk: [HDCO]/[H$_2$CO]$\sim0.1-0.3$ and [D$_2$CO]/[H$_2$CO]$\sim0.1$. Interestingly, while HDCO follows the H$_2$CO distribution in the disk and in the streamer, the distribution of D$_2$CO is highly asymmetric, with a peak of the emission (and [D]/[H] ratio) in the SE disk side, where the streamer crashes onto the disk. In addition, D$_2$CO is detected in two spots along the blue- and red-shifted outflow. This suggests that: (i) in the disk, HDCO formation is dominated by gas-phase reactions similarly to H$_2$CO, while (ii) D$_2$CO was mainly formed on the grain mantles during the prestellar phase and/or in the disk itself, and is at present released in the gas-phase in the shocks driven by the streamer and the outflow. These findings testify on the key role of streamers in the build-up of the disk both concerning the final mass available for planet formation and its chemical composition. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 12 pages, 10 figures, accepted for publication on A&A

Multiple chemical tracers finally unveil the intricate NGC\,1333 IRAS\,4A outflow system. FAUST XVI

Authors: Layal Chahine, Cecilia Ceccarelli, Marta De Simone, Claire J. Chandler, Claudio Codella, Linda Podio, Ana López-Sepulcre, Nami Sakai, Laurent Loinard, Mathilde Bouvier, Paola Caselli, Charlotte Vastel, Eleonora Bianchi, Nicolás Cuello, Francesco Fontani, Doug Johnstone, Giovanni Sabatini, Tomoyuki Hanawa, Ziwei E. Zhang, Yuri Aikawa, Gemma Busquet, Emmanuel Caux, Aurore Durán, Eric Herbst, François Ménard , et al. (32 additional authors not shown)

Abstract: The exploration of outflows in protobinary systems presents a challenging yet crucial endeavour, offering valuable insights into the dynamic interplay between protostars and their evolution. In this study, we examine the morphology and dynamics of jets and outflows within the IRAS\,4A protobinary system. This analysis is based on ALMA observations of SiO(5--4), H$_2$CO(3$_{0,3}$--2$_{0,3}$), and H… ▽ More The exploration of outflows in protobinary systems presents a challenging yet crucial endeavour, offering valuable insights into the dynamic interplay between protostars and their evolution. In this study, we examine the morphology and dynamics of jets and outflows within the IRAS\,4A protobinary system. This analysis is based on ALMA observations of SiO(5--4), H$_2$CO(3$_{0,3}$--2$_{0,3}$), and HDCO(4$_{1,4}$--3$_{1,3}$) with a spatial resolution of $\sim$150\,au. Leveraging an astrochemical approach involving the use of diverse tracers beyond traditional ones has enabled the identification of novel features and a comprehensive understanding of the broader outflow dynamics. Our analysis reveals the presence of two jets in the redshifted emission, emanating from IRAS\,4A1 and IRAS\,4A2, respectively. Furthermore, we identify four distinct outflows in the region for the first time, with each protostar, 4A1 and 4A2, contributing to two of them. We characterise the morphology and orientation of each outflow, challenging previous suggestions of bends in their trajectories. The outflow cavities of IRAS\,4A1 exhibit extensions of 10$''$ and 13$''$ with position angles (PA) of 0$^{\circ}$ and -12$^{\circ}$, respectively, while those of IRAS\,4A2 are more extended, spanning 18$''$ and 25$''$ with PAs of 29$^{\circ}$ and 26$^{\circ}$. We propose that the misalignment of the cavities is due to a jet precession in each protostar, a notion supported by the observation that the more extended cavities of the same source exhibit lower velocities, indicating they may stem from older ejection events. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Into the thick of it: ALMA 0.45 mm observations of HL Tau at 2 au resolution

Authors: Osmar M. Guerra-Alvarado, Carlos Carrasco-González, Enrique Macías, Nienke van der Marel, Adrien Houge, Luke T. Maud, Paola Pinilla, Marion Villenave, Yoshiharu Asaki, Elizabeth Humphreys

Abstract: Aims. To comprehend the efficiency of dust evolution within protoplanetary disks, it is crucial to conduct studies of these disks using high-resolution observations at multiple wavelengths with the Atacama Large Millimeter/submillimeter Array (ALMA). Methods. In this work, we present high-frequency ALMA observations of the HL Tau disk using its Band 9 centered at a wavelength of 0.45 mm. These obs… ▽ More Aims. To comprehend the efficiency of dust evolution within protoplanetary disks, it is crucial to conduct studies of these disks using high-resolution observations at multiple wavelengths with the Atacama Large Millimeter/submillimeter Array (ALMA). Methods. In this work, we present high-frequency ALMA observations of the HL Tau disk using its Band 9 centered at a wavelength of 0.45 mm. These observations achieve the highest angular resolution in a protoplanetary disk to date, 12 milliarcseconds (mas), allowing the study of the dust emission at scales of 2 au. We use these data to extend the previously published multi-wavelength analysis of the HL Tau disk. Results. Our new 0.45 mm data traces mainly optically thick emission, providing a tight constraint to the dust temperature profile. We derive maximum particle sizes of $\sim$1 cm from the inner disk to $\sim$60 au. Beyond this radius, we find particles between 300 $μ$m and 1 mm. Moreover, an intriguing asymmetry is observed at 32 au in the northeast inner part of the HL Tau disk at 0.45 mm. We propose that this asymmetry is the outcome of a combination of factors including the optically thick nature of the emission, the orientation of the disk, and a relatively large dust scale height of the grains. To validate this, we conducted a series of radiative transfer models using the RADMC-3D software. If this scenario is correct, our measured dust mass within 32 au would suggest a dust scale height H/R> 0.08 for the inner disk. Finally, the unprecedented resolution allowed us to probe for the first time the dust emission down to a few au scales. We observed an increase in brightness temperature inside the estimated water snowline and speculate whether it could indicate the presence of a traffic jam effect in the inner disk. Abridge △ Less

Submitted 5 April, 2024; originally announced April 2024.

Comments: Accepted for publication in A&A, 13 pages, 8 figures

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

FAUST XIII. Dusty cavity and molecular shock driven by IRS7B in the Corona Australis cluster

Authors: G. Sabatini, L. Podio, C. Codella, Y. Watanabe, M. De Simone, E. Bianchi, C. Ceccarelli, C. J. Chandler, N. Sakai, B. Svoboda, L. Testi, Y. Aikawa, N. Balucani, M. Bouvier, P. Caselli, E. Caux, L. Chahine, S. Charnley, N. Cuello, F. Dulieu, L. Evans, D. Fedele, S. Feng, F. Fontani, T. Hama , et al. (32 additional authors not shown)

Abstract: The origin of the chemical diversity observed around low-mass protostars probably resides in the earliest history of these systems. We aim to investigate the impact of protostellar feedback on the chemistry and grain growth in the circumstellar medium of multiple stellar systems. In the context of the ALMA Large Program FAUST, we present high-resolution (50 au) observations of CH$_3$OH, H$_2$CO, a… ▽ More The origin of the chemical diversity observed around low-mass protostars probably resides in the earliest history of these systems. We aim to investigate the impact of protostellar feedback on the chemistry and grain growth in the circumstellar medium of multiple stellar systems. In the context of the ALMA Large Program FAUST, we present high-resolution (50 au) observations of CH$_3$OH, H$_2$CO, and SiO and continuum emission at 1.3 mm and 3 mm towards the Corona Australis star cluster. Methanol emission reveals an arc-like structure at $\sim$1800 au from the protostellar system IRS7B along the direction perpendicular to the major axis of the disc. The arc is located at the edge of two elongated continuum structures that define a cone emerging from IRS7B. The region inside the cone is probed by H$_2$CO, while the eastern wall of the arc shows bright emission in SiO, a typical shock tracer. Taking into account the association with a previously detected radio jet imaged with JVLA at 6 cm, the molecular arc reveals for the first time a bow shock driven by IRS7B and a two-sided dust cavity opened by the mass-loss process. For each cavity wall, we derive an average H$_2$ column density of $\sim$7$\times$10$^{21}$ cm$^{-2}$, a mass of $\sim$9$\times$10$^{-3}$ M$_\odot$, and a lower limit on the dust spectral index of $1.4$. These observations provide the first evidence of a shock and a conical dust cavity opened by the jet driven by IRS7B, with important implications for the chemical enrichment and grain growth in the envelope of Solar System analogues. △ Less

Submitted 2 April, 2024; v1 submitted 26 March, 2024; originally announced March 2024.

Comments: 12 pages, 8 figures, 3 tables. Accepted Letter in Astronomy & Astrophysics

ALMA High-frequency Long Baseline Campaign in 2021: Highest Angular Resolution Submillimeter Wave Images for the Carbon-rich Star R Lep

Authors: Yoshiharu Asaki, Luke T. Maud, Harold Francke, Hiroshi Nagai, Dirk Petry, Edward B. Fomalont, Elizabeth Humphreys, Anita M. S. Richards, Ka Tat Wong, William Dent, Akihiko Hirota, Jose Miguel Fernandez, Satoko Takahashi, Antonio S. Hales

Abstract: The Atacama Large Millimeter/submillimeter Array (ALMA) was used in 2021 to image the carbon-rich evolved star R Lep in Bands 8-10 (397-908 GHz) with baselines up to 16 km. The goal was to validate the calibration, using band-to-band (B2B) phase referencing with a close phase calibrator J0504-1512, 1.2 deg from R Lep in this case, and the imaging procedures required to obtain the maximum angular r… ▽ More The Atacama Large Millimeter/submillimeter Array (ALMA) was used in 2021 to image the carbon-rich evolved star R Lep in Bands 8-10 (397-908 GHz) with baselines up to 16 km. The goal was to validate the calibration, using band-to-band (B2B) phase referencing with a close phase calibrator J0504-1512, 1.2 deg from R Lep in this case, and the imaging procedures required to obtain the maximum angular resolution achievable with ALMA. Images of the continuum emission and the hydrogen cyanide (HCN) maser line at 890.8 GHz, from the J=10-9 transition between the (1110) and (0400) vibrationally excited states, achieved angular resolutions of 13, 6, and 5 mas in Bands 8-10, respectively. Self-calibration (self-cal) was used to produce ideal images as to compare with the B2B phase referencing technique. The continuum emission was resolved in Bands 9 and 10, leaving too little flux for self-cal of the longest baselines, so these comparisons are made at coarser resolution. Comparisons showed that B2B phase referencing provided phase corrections sufficient to recover 92%, 83%, and 77% of the ideal image continuum flux densities. The HCN maser was sufficiently compact to obtain self-cal solutions in Band 10 for all baselines (up to 16 km). In Band 10, B2B phase referencing as compared to the ideal images recovered 61% and 70% of the flux density for the HCN maser and continuum, respectively. △ Less

Submitted 18 October, 2023; v1 submitted 14 October, 2023; originally announced October 2023.

Comments: 37 pages, 12 figures, 9 tables, accepted by ApJ (Aug 30, 2023)

The ALMA Interferometric Pipeline Heuristics

Authors: Todd R. Hunter, Remy Indebetouw, Crystal L. Brogan, Kristin Berry, Chin-Shin Chang, Harold Francke, Vincent C. Geers, Laura Gómez, John E. Hibbard, Elizabeth M. Humphreys, Brian R. Kent, Amanda A. Kepley, Devaky Kunneriath, Andrew Lipnicky, Ryan A. Loomis, Brian S. Mason, Joseph S. Masters, Luke T. Maud, Dirk Muders, Jose Sabater, Kanako Sugimoto, László Szűcs, Eugene Vasiliev, Liza Videla, Eric Villard , et al. (3 additional authors not shown)

Abstract: We describe the calibration and imaging heuristics developed and deployed in the ALMA interferometric data processing pipeline, as of ALMA Cycle 9. The pipeline software framework is written in Python, with each data reduction stage layered on top of tasks and toolkit functions provided by the Common Astronomy Software Applications package. This framework supports a variety of tasks for observator… ▽ More We describe the calibration and imaging heuristics developed and deployed in the ALMA interferometric data processing pipeline, as of ALMA Cycle 9. The pipeline software framework is written in Python, with each data reduction stage layered on top of tasks and toolkit functions provided by the Common Astronomy Software Applications package. This framework supports a variety of tasks for observatory operations, including science data quality assurance, observing mode commissioning, and user reprocessing. It supports ALMA and VLA interferometric data along with ALMA and NRO45m single dish data, via different stages and heuristics. In addition to producing calibration tables, calibrated measurement sets, and cleaned images, the pipeline creates a WebLog which serves as the primary interface for verifying the data quality assurance by the observatory and for examining the contents of the data by the user. Following the adoption of the pipeline by ALMA Operations in 2014, the heuristics have been refined through annual development cycles, culminating in a new pipeline release aligned with the start of each ALMA Cycle of observations. Initial development focused on basic calibration and flagging heuristics (Cycles 2-3), followed by imaging heuristics (Cycles 4-5), refinement of the flagging and imaging heuristics with parallel processing (Cycles 6-7), addition of the moment difference analysis to improve continuum channel identification (2020 release), addition of a spectral renormalization stage (Cycle 8), and improvement in low SNR calibration heuristics (Cycle 9). In the two most recent Cycles, 97% of ALMA datasets were calibrated and imaged with the pipeline, ensuring long-term automated reproducibility. We conclude with a brief description of plans for future additions, including self-calibration, multi-configuration imaging, and calibration and imaging of full polarization data. △ Less

Submitted 25 July, 2023; v1 submitted 12 June, 2023; originally announced June 2023.

Comments: Accepted for publication by Publications of the Astronomical Society of the Pacific, 65 pages, 20 figures, 10 tables, 2 appendices. Small corrections and additions from the proof stage have been applied

Journal ref: Publications of the Astronomical Society of the Pacific, volume 135, number 1049, 074501, 2023 July 24

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)

Updates to ALMA Site Properties: using the ESO-Allegro Phase RMS database -- ALMA Memo 624

Authors: Luke T. Maud, Andrés F. Pérez-Sánchez, Yoshiharu Asaki, Felix Stoehr, Bill Dent, María Díaz Trigo

Abstract: We present a long-term overview of the atmospheric phase stability at the Atacama Large Millimeter/submillimeter Array (ALMA) site, using >5 years of data, that acts as the successor to the studies summarized two decades ago by Evans et al 2003. Importantly, we explore the atmospheric variations, the `phase RMS', and associated metadata of over 17000 accrued ALMA observations taken since Cycle 3 (… ▽ More We present a long-term overview of the atmospheric phase stability at the Atacama Large Millimeter/submillimeter Array (ALMA) site, using >5 years of data, that acts as the successor to the studies summarized two decades ago by Evans et al 2003. Importantly, we explore the atmospheric variations, the `phase RMS', and associated metadata of over 17000 accrued ALMA observations taken since Cycle 3 (2015) by using the Bandpass calibrator source scans. We indicate the temporal phase RMS trends for average baseline lengths of 500, 1000, 5000, and 10000m, in contrast to the old stability studies that used a single 300m baseline phase monitor system. At the ALMA site, on the Chajnantor plateau, we report the diurnal variations and monthly changes in the phase RMS on ALMA relevant baselines lengths, measured directly from data, and we reaffirm such trends in atmospheric transmission (via Precipitable Water Vapour - PWV). We confirm that day observations have respectively higher phase RMS and PWV in contrast to night, while the monthly variations show Chilean winter (June - August) providing the best, high-frequency and long-baseline observing conditions - low (stable) phase RMS and low PWV. Yet, not all good phase stability condition occur when the PWV is low. Measurements of the phase RMS as a function of short timescales, 30 to 240s, that tie with typical target source scan times, and as a function of baseline length indicate that phase variations are smaller for short timescales and baselines and larger for longer timescales and baselines. We illustrate that fast-switching phase-referencing techniques, that allow short target scan times, could work well in reducing the phase RMS to suitable levels specifically for high-frequencies (Band 8, 9 and 10), long-baselines, and the two combined. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 34 pages, 19 Figures, 10 Tables ALMA Memo 624: https://library.nrao.edu/public/memos/alma/main/memo624.pdf

Salt-bearing disk candidates around high-mass young stellar objects

Authors: Adam Ginsburg, Brett A. McGuire, Patricio Sanhueza, Fernando Olguin, Luke T Maud, Kei E. I. Tanaka, Yichen Zhang, Henrik Beuther, Nick Indriolo

Abstract: Molecular lines tracing the orbital motion of gas in a well-defined disk are valuable tools for inferring both the properties of the disk and the star it surrounds. Lines that arise only from a disk, and not also from the surrounding molecular cloud core that birthed the star or from the outflow it drives, are rare. Several such emission lines have recently been discovered in one example case, tho… ▽ More Molecular lines tracing the orbital motion of gas in a well-defined disk are valuable tools for inferring both the properties of the disk and the star it surrounds. Lines that arise only from a disk, and not also from the surrounding molecular cloud core that birthed the star or from the outflow it drives, are rare. Several such emission lines have recently been discovered in one example case, those from NaCl and KCl salt molecules. We studied a sample of 23 candidate high-mass young stellar objects (HMYSOs) in 17 high-mass star-forming regions to determine how frequently emission from these species is detected. We present five new detections of water, NaCl, KCl, PN, and SiS from the innermost regions around the objects, bringing the total number of known briny disk candidates to nine. Their kinematic structure is generally disk-like, though we are unable to determine whether they arise from a disk or outflow in the sources with new detections. We demonstrate that these species are spatially coincident in a few resolved cases and show that they are generally detected together, suggesting a common origin or excitation mechanism. We also show that several disks around HMYSOs clearly do not exhibit emission in these species. Salty disks are therefore neither particularly rare in high-mass disks, nor are they ubiquitous. △ Less

Submitted 4 November, 2022; originally announced November 2022.

Comments: accepted to ApJ

FAUST VI. VLA 1623--2417 B: a new laboratory for astrochemistry around protostars on 50 au scale

Authors: C. Codella, A. López-Sepulcre, S. Ohashi, C. J. Chandler, M. De Simone, L. Podio, C. Ceccarelli, N. Sakai, F. Alves, A. Durán, D. Fedele, L. Loinard, S. Mercimek, N. Murillo, E. Bianchi, M. Bouvier, G. Busquet, P. Caselli, F. Dulieu, S. Feng, T. Hanawa, D. Johnstone, B. Lefloch, L. T. Maud, G. Moellenbrock , et al. (3 additional authors not shown)

Abstract: The ALMA interferometer, with its unprecedented combination of high-sensitivity and high-angular resolution, allows for (sub-)mm wavelength mapping of protostellar systems at Solar System scales. Astrochemistry has benefited from imaging interstellar complex organic molecules in these jet-disk systems. Here we report the first detection of methanol (CH3OH) and methyl formate (HCOOCH3) emission tow… ▽ More The ALMA interferometer, with its unprecedented combination of high-sensitivity and high-angular resolution, allows for (sub-)mm wavelength mapping of protostellar systems at Solar System scales. Astrochemistry has benefited from imaging interstellar complex organic molecules in these jet-disk systems. Here we report the first detection of methanol (CH3OH) and methyl formate (HCOOCH3) emission towards the triple protostellar system VLA1623-2417 A1+A2+B, obtained in the context of the ALMA Large Program FAUST. Compact methanol emission is detected in lines from Eu = 45 K up to 61 K and 537 K towards components A1 and B, respectively. LVG analysis of the CH3OH lines towards VLA1623-2417 B indicates a size of 0.11-0.34 arcsec (14-45 au), a column density N(CH3OH) = 10^16-10^17 cm-2, kinetic temperature > 170 K, and volume density > 10^8 cm-3. An LTE approach is used for VLA1623-2417 A1, given the limited Eu range, and yields Trot < 135 K. The methanol emission around both VLA1623-2417 A1 and B shows velocity gradients along the main axis of each disk. Although the axial geometry of the two disks is similar, the observed velocity gradients are reversed. The CH3OH spectra from B shows two broad (4-5 km s-1) peaks, which are red- and blue-shifted by about 6-7 km s-1 from the systemic velocity. Assuming a chemically enriched ring within the accretion disk, close to the centrifugal barrier, its radius is calculated to be 33 au. The methanol spectra towards A1 are somewhat narrower (about 4 km s-1), implying a radius of 12-24 au. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: MNRAS, in press

FAUST III. Misaligned rotations of the envelope, outflow, and disks in the multiple protostellar system of VLA 1623$-$2417

Authors: Satoshi Ohashi, Claudio Codella, Nami Sakai, Claire J. Chandler, Cecilia Ceccarelli, Felipe Alves, Davide Fedele, Tomoyuki Hanawa, Aurora Durán, Cécile Favre, Ana López-Sepulcre, Laurent Loinard, Seyma Mercimek, Nadia M. Murillo, Linda Podio, Yichen Zhang, Yuri Aikawa, Nadia Balucani, Eleonora Bianchi, Mathilde Bouvier, Gemma Busquet, Paola Caselli, Emmanuel Caux, Steven Charnley, Spandan Choudhury , et al. (47 additional authors not shown)

Abstract: We report a study of the low-mass Class-0 multiple system VLA 1623AB in the Ophiuchus star-forming region, using H$^{13}$CO$^+$ ($J=3-2$), CS ($J=5-4$), and CCH ($N=3-2$) lines as part of the ALMA Large Program FAUST. The analysis of the velocity fields revealed the rotation motion in the envelope and the velocity gradients in the outflows (about 2000 au down to 50 au). We further investigated the… ▽ More We report a study of the low-mass Class-0 multiple system VLA 1623AB in the Ophiuchus star-forming region, using H$^{13}$CO$^+$ ($J=3-2$), CS ($J=5-4$), and CCH ($N=3-2$) lines as part of the ALMA Large Program FAUST. The analysis of the velocity fields revealed the rotation motion in the envelope and the velocity gradients in the outflows (about 2000 au down to 50 au). We further investigated the rotation of the circum-binary VLA 1623A disk as well as the VLA 1623B disk. We found that the minor axis of the circum-binary disk of VLA 1623A is misaligned by about 12 degrees with respect to the large-scale outflow and the rotation axis of the envelope. In contrast, the minor axis of the circum-binary disk is parallel to the large-scale magnetic field according to previous dust polarization observations, suggesting that the misalignment may be caused by the different directions of the envelope rotation and the magnetic field. If the velocity gradient of the outflow is caused by rotation, the outflow has a constant angular momentum and the launching radius is estimated to be $5-16$ au, although it cannot be ruled out that the velocity gradient is driven by entrainments of the two high-velocity outflows. Furthermore, we detected for the first time a velocity gradient associated with rotation toward the VLA 16293B disk. The velocity gradient is opposite to the one from the large-scale envelope, outflow, and circum-binary disk. The origin of its opposite gradient is also discussed. △ Less

Submitted 18 January, 2022; originally announced January 2022.

Comments: 27 pages, 21 figures, 2 Tables, Accepted for publication in ApJ

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)

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)

FAUST II. Discovery of a Secondary Outflow in IRAS 15398-3359: Variability in Outflow Direction during the Earliest Stage of Star Formation?

Authors: Yuki Okoda, Yoko Oya, Logan Francis, Doug Johnstone, Shu-ichiro Inutsuka, Cecilia Ceccarelli, Claudio Codella, Claire Chandler, Nami Sakai, Yuri Aikawa, Felipe Alves, Nadia Balucani, Eleonora Bianchi, Mathilde Bouvier, Paola Caselli, Emmanuel Caux, Steven Charnley, Spandan Choudhury, Marta De Simone, Francois Dulieu, Aurora Durán, Lucy Evans, Cécile Favre, Davide Fedele, Siyi Feng , et al. (44 additional authors not shown)

Abstract: We have observed the very low-mass Class 0 protostar IRAS 15398-3359 at scales ranging from 50 au to 1800 au, as part of the ALMA Large Program FAUST. We uncover a linear feature, visible in H2CO, SO, and C18O line emission, which extends from the source along a direction almost perpendicular to the known active outflow. Molecular line emission from H2CO, SO, SiO, and CH3OH further reveals an arc-… ▽ More We have observed the very low-mass Class 0 protostar IRAS 15398-3359 at scales ranging from 50 au to 1800 au, as part of the ALMA Large Program FAUST. We uncover a linear feature, visible in H2CO, SO, and C18O line emission, which extends from the source along a direction almost perpendicular to the known active outflow. Molecular line emission from H2CO, SO, SiO, and CH3OH further reveals an arc-like structure connected to the outer end of the linear feature and separated from the protostar, IRAS 15398-3359, by 1200 au. The arc-like structure is blue-shifted with respect to the systemic velocity. A velocity gradient of 1.2 km/s over 1200 au along the linear feature seen in the H2CO emission connects the protostar and the arc-like structure kinematically. SO, SiO, and CH3OH are known to trace shocks, and we interpret the arc-like structure as a relic shock region produced by an outflow previously launched by IRAS 15398-3359. The velocity gradient along the linear structure can be explained as relic outflow motion. The origins of the newly observed arc-like structure and extended linear feature are discussed in relation to turbulent motions within the protostellar core and episodic accretion events during the earliest stage of protostellar evolution. △ Less

Submitted 18 January, 2021; originally announced January 2021.

ALMA High-frequency Long-baseline Campaign in 2017: A Comparison of the Band-to-band and In-band Phase Calibration Techniques and Phase-calibrator Separation Angles

Authors: Luke T. Maud, Yoshiharu Asaki, Edward B. Fomalont, William R. F. Dent, Akihiko Hirota, Satoki Matsushita, Neil M. Phillips, John M. Carpenter, Satoko Takahashi, Eric Villard, Tsuyoshi Sawada, Stuartt Corder

Abstract: The Atacama Large millimeter/submillimeter Array (ALMA) obtains spatial resolutions of 15 to 5 milli-arcsecond (mas) at 275-950GHz (0.87-0.32mm) with 16km baselines. Calibration at higher-frequencies is challenging as ALMA sensitivity and quasar density decrease. The Band-to-Band (B2B) technique observes a detectable quasar at lower frequency that is closer to the target, compared to one at the ta… ▽ More The Atacama Large millimeter/submillimeter Array (ALMA) obtains spatial resolutions of 15 to 5 milli-arcsecond (mas) at 275-950GHz (0.87-0.32mm) with 16km baselines. Calibration at higher-frequencies is challenging as ALMA sensitivity and quasar density decrease. The Band-to-Band (B2B) technique observes a detectable quasar at lower frequency that is closer to the target, compared to one at the target high-frequency. Calibration involves a nearly constant instrumental phase offset between the frequencies and the conversion of the temporal phases to the target frequency. The instrumental offsets are solved with a differential-gain-calibration (DGC) sequence, consisting of alternating low and high frequency scans of strong quasar. Here we compare B2B and in-band phase referencing for high-frequencies ($>$289GHz) using 2-15km baselines and calibrator separation angles between $\sim$0.68 and $\sim$11.65$^{\circ}$. The analysis shows that: (1) DGC for B2B produces a coherence loss $<$7% for DGC phase RMS residuals $<$30$^{\circ}$. (2) B2B images using close calibrators ( $<$1.67$^{\circ}$ ) are superior to in-band images using distant ones ( $>$2.42$^{\circ}$ ). (3) For more distant calibrators, B2B is preferred if it provides a calibrator $\sim$2$^{\circ}$ closer than the best in-band calibrator. (4) Decreasing image coherence and poorer image quality occur with increasing phase calibrator separation angle because of uncertainties in the antenna positions and sub-optimal phase referencing. (5) To achieve $>$70% coherence for long-baseline (16 km) band 7 (289GHz) observations, calibrators should be within $\sim$4$^{\circ}$ of the target. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Comments: 17 Figures, 10 Tables

Journal ref: 2020ApJS..250...18M

FAUST I. The hot corino at the heart of the prototypical Class I protostar L1551 IRS5

Authors: E. Bianchi, C. J. Chandler, C. Ceccarelli, C. Codella, N. Sakai, A. López-Sepulcre, L. T. Maud, G. Moellenbrock, B. Svoboda, Y. Watanabe, T. Sakai, F. Ménard, Y. Aikawa, F. Alves, N. Balucani, M. Bouvier, P. Caselli, E. Caux, S. Charnley, S. Choudhury, M. De Simone, F. Dulieu, A. Durán, L. Evans, C. Favre , et al. (41 additional authors not shown)

Abstract: The study of hot corinos in Solar-like protostars has been so far mostly limited to the Class 0 phase, hampering our understanding of their origin and evolution. In addition, recent evidence suggests that planet formation starts already during Class I phase, which, therefore, represents a crucial step in the future planetary system chemical composition. Hence, the study of hot corinos in Class I p… ▽ More The study of hot corinos in Solar-like protostars has been so far mostly limited to the Class 0 phase, hampering our understanding of their origin and evolution. In addition, recent evidence suggests that planet formation starts already during Class I phase, which, therefore, represents a crucial step in the future planetary system chemical composition. Hence, the study of hot corinos in Class I protostars has become of paramount importance. Here we report the discovery of a hot corino towards the prototypical Class I protostar L1551 IRS5, obtained within the ALMA Large Program FAUST. We detected several lines from methanol and its isopotologues ($^{13}$CH$_{\rm 3}$OH and CH$_{\rm 2}$DOH), methyl formate and ethanol. Lines are bright toward the north component of the IRS5 binary system, and a possible second hot corino may be associated with the south component. The methanol lines non-LTE analysis constrains the gas temperature ($\sim$100 K), density ($\geq$1.5$\times$10$^{8}$ cm$^{-3}$), and emitting size ($\sim$10 au in radius). All CH$_{\rm 3}$OH and $^{13}$CH$_{\rm 3}$OH lines are optically thick, preventing a reliable measure of the deuteration. The methyl formate and ethanol relative abundances are compatible with those measured in Class 0 hot corinos. Thus, based on the present work, little chemical evolution from Class 0 to I hot corinos occurs. △ Less

Submitted 20 July, 2020; originally announced July 2020.

Comments: 6 pages, 2 figures

ALMA Band-to-band Phase Referencing: Imaging Capabilities on Long Baselines and High Frequencies

Authors: Yoshiharu Asaki, Luke T. Maud, Edward B. Fomalont, William R. F. Dent, Loreto Barcos-Muñoz, Neil M. Phillips, Akihiko Hirota, Satoko Takahashi, Stuartt Corder, John M. Carpenter, Eric Villard

Abstract: High-frequency long-baseline experiments with the Atacama Large Millimeter/submillimeter Array were organized to test the high angular resolution imaging capabilities in the submillimeter wave regime using baselines up to 16 km. Four experiments were conducted, two Band 7 (289 GHz) and two Band 8 (405 GHz) observations. Phase correction using band-to-band (B2B) phase referencing was used with a ph… ▽ More High-frequency long-baseline experiments with the Atacama Large Millimeter/submillimeter Array were organized to test the high angular resolution imaging capabilities in the submillimeter wave regime using baselines up to 16 km. Four experiments were conducted, two Band 7 (289 GHz) and two Band 8 (405 GHz) observations. Phase correction using band-to-band (B2B) phase referencing was used with a phase calibrator only 0.7deg away observed in Band 3 (96 GHz) and Band 4 (135 GHz), respectively. In Band 8, we achieved the highest resolution of 14x11 mas. We compared the synthesis images of the target quasar using 20 and 60 s switching cycle times in the phase referencing. In Band 7, the atmosphere had good stability in phase rms (<0.5 rad over 2 minutes), and there was little difference in image coherence between the 20 and 60 s switching cycle times. One Band 8 experiment was conducted under a worse phase rms condition (>1 rad over 2 minutes), which led to a significantly reduced coherence when using the 60 s switching cycle time. One of our four experiments indicates that the residual phase rms error after phase referencing can be reduced to 0.16 rad at 289 GHz in using the 20 s switching cycle time. Such conditions would meet the phase correction requirement of image coherence of >70% in Band 10, assuming a similar phase calibrator separation angle, emphasizing the need for such B2B phase referencing observing at high frequencies. △ Less

Submitted 14 July, 2020; originally announced July 2020.

Comments: 40 pages, 10 figures, 6 tables

Journal ref: 2020 AJ 160 59

ALMA High-frequency Long Baseline Campaign in 2017: Band-to-band Phase Referencing in Submillimeter Waves

Authors: Yoshiharu Asaki, Luke T. Maud, Edward B. Fomalont, Neil M. Phillips, Akihiko Hirota, Tsuyoshi Sawada, Loreto Barcos-Muñoz, Anita M. S. Richards, William R. F. Dent, Satoko Takahashi, Stuartt Corder, John M. Carpenter, Eric Villard, Elizabeth M. Humphreys

Abstract: In 2017, an Atacama Large Millimeter/submillimeter Array (ALMA) high-frequency long baseline campaign was organized to test image capabilities with baselines up to 16 km at submillimeter (submm) wavelengths. We investigated image qualities using ALMA receiver Bands 7, 8, 9, and 10 (285-875 GHz) by adopting band-to-band (B2B) phase referencing in which a phase calibrator is tracked at a lower frequ… ▽ More In 2017, an Atacama Large Millimeter/submillimeter Array (ALMA) high-frequency long baseline campaign was organized to test image capabilities with baselines up to 16 km at submillimeter (submm) wavelengths. We investigated image qualities using ALMA receiver Bands 7, 8, 9, and 10 (285-875 GHz) by adopting band-to-band (B2B) phase referencing in which a phase calibrator is tracked at a lower frequency. For B2B phase referencing, it is expected that a closer phase calibrator to a target can be used, comparing to standard in-band phase referencing. In the first step, it is ensured that an instrumental phase offset difference between low- and high-frequency Bands can be removed using a differential gain calibration in which a phase calibrator is certainly detected while frequency switching. In the next step, comparative experiments are arranged to investigate the image quality between B2B and in-band phase referencing with phase calibrators at various separation angles. In the final step, we conducted long baseline imaging tests for a quasar at 289 GHz in Band 7 and 405 GHz in Band 8 and complex structure sources of HL Tau and VY CMa at ~670 GHz in Band 9. The B2B phase referencing was successfully applied, allowing us to achieve an angular resolution of 14x11 and 10x8 mas for HL Tau and VY CMa, respectively. There is a high probability of finding a low-frequency calibrator within 5.4 deg in B2B phase referencing, bright enough to use an 8 s scan length combined with a 7.5 GHz bandwidth. △ Less

Submitted 23 March, 2020; v1 submitted 16 March, 2020; originally announced March 2020.

Comments: 61 pages, 17 figures, 8 tables

Journal ref: 2020ApJS..247...23A

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)

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)

High-mass star formation at sub-50AU scales

Authors: H. Beuther, A. Ahmadi. J. C. Mottram, H. Linz, L. T. Maud, Th. Henning, R. Kuiper, A. J. Walsh, K. G. Johnston, S. N. Longmore

Abstract: Methods: We observed the high-mass hot core region G351.77-0.54 with ALMA and more than 16km baselines. Results: At a spatial resolution of 18/40au (depending on the distance), we identify twelve sub-structures within the inner few thousand au of the region. The brightness temperatures are high, reaching values greater 1000K, signposting high optical depth toward the peak positions. Core separat… ▽ More Methods: We observed the high-mass hot core region G351.77-0.54 with ALMA and more than 16km baselines. Results: At a spatial resolution of 18/40au (depending on the distance), we identify twelve sub-structures within the inner few thousand au of the region. The brightness temperatures are high, reaching values greater 1000K, signposting high optical depth toward the peak positions. Core separations vary between sub-100au to several 100 and 1000au. The core separations and approximate masses are largely consistent with thermal Jeans fragmentation of a dense gas core. Due to the high continuum optical depth, most spectral lines are seen in absorption. However, a few exceptional emission lines are found that most likely stem from transitions with excitation conditions above1000K. Toward the main continuum source, these emission lines exhibit a velocity gradient across scales of 100-200au aligned with the molecular outflow and perpendicular to the previously inferred disk orientation. While we cannot exclude that these observational features stem from an inner hot accretion disk, the alignment with the outflow rather suggests that it stems from the inner jet and outflow region. The highest-velocity features are found toward the peak position, and no Hubble-like velocity structure can be identified. Therefore, these data are consistent with steady-state turbulent entrainment of the hot molecular gas via Kelvin-Helmholtz instabilities at the interface between the jet and the outflow. Conclusions: Resolving this high-mass star-forming region at sub-50au scales indicates that the hierarchical fragmentation process in the framework of thermal Jeans fragmentation can continue down to the smallest accessible spatial scales. Velocity gradients on these small scales have to be treated cautiously and do not necessarily stem from disks, but may be better explained with outflow emission. △ Less

Submitted 26 November, 2018; originally announced November 2018.

Comments: Accepted by Astronomy & Astrophysics, high-resolution version available at http://www.mpia.de/homes/beuther/papers.html

Journal ref: A&A 621, A122 (2019)

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)

Evidence for the start of planet formation in a young circumstellar disk

Authors: Daniel Harsono, Per Bjerkeli, Matthijs H. D. van der Wiel, Jon P. Ramsey, Luke T. Maud, Lars E. Kristensen, Jes K. Jørgensen

Abstract: The growth of dust grains in protoplanetary disks is a necessary first step towards planet formation. This growth has been inferred via observations of thermal dust emission towards mature protoplanetary systems (age >2 million years) with masses that are, on average, similar to Neptune3. In contrast, the majority of confirmed exoplanets are heavier than Neptune. Given that young protoplanetary di… ▽ More The growth of dust grains in protoplanetary disks is a necessary first step towards planet formation. This growth has been inferred via observations of thermal dust emission towards mature protoplanetary systems (age >2 million years) with masses that are, on average, similar to Neptune3. In contrast, the majority of confirmed exoplanets are heavier than Neptune. Given that young protoplanetary disks are more massive than their mature counterparts, this suggests that planet formation starts early, but evidence for grain growth that is spatially and temporally coincident with a massive reservoir in young disks remains scarce. Here, we report observations on a lack of emission of carbon monoxide isotopologues within the inner ~15 au of a very young (age ~100,000 years) disk around the Solar-type protostar TMC1A. By using the absence of spatially resolved molecular line emission to infer the gas and dust content of the disk, we conclude that shielding by millimeter-size grains is responsible for the lack of emission. This suggests that grain growth and millimeter-size dust grains can be spatially and temporally coincident with a mass reservoir sufficient for giant planet formation. Hence, planet formation starts during the earliest, embedded phases in the life of young stars. △ Less

Submitted 25 June, 2018; originally announced June 2018.

Comments: Accepted for publication in Nature Astronomy, 3 figures, 3 extended figures

V1094 Sco: a rare giant multi-ringed disk around a T Tauri star

Authors: S. E. van Terwisga, E. F. van Dishoeck, M. Ansdell, N. van der Marel, L. Testi, J. P. Williams, S. Facchini, M. Tazzari, M. R. Hogerheijde, L. Trapman, C. F. Manara, A. Miotello, L. T. Maud, D. Harsono

Abstract: A wide variety of ring-like dust structures has been detected in protoplanetary disks, but their origin and frequency are still unclear. We characterize the structure of an extended, multi-ringed disk discovered serendipitously in the ALMA Lupus disk survey and put it in the context of the Lupus disk population. ALMA observations in Band 6 at 234 GHz and Band 7 at 328 GHz at 0.3" resolution toward… ▽ More A wide variety of ring-like dust structures has been detected in protoplanetary disks, but their origin and frequency are still unclear. We characterize the structure of an extended, multi-ringed disk discovered serendipitously in the ALMA Lupus disk survey and put it in the context of the Lupus disk population. ALMA observations in Band 6 at 234 GHz and Band 7 at 328 GHz at 0.3" resolution toward the K6 star V1094 Sco in Lupus III are presented, and its disk structure is analyzed. The spectral index $α_{mm}$ is determined in the inner 150 AU of the disk. The ALMA continuum data show a very extended disk with two gap/ring pairs. The gaps are located at 100 AU and 170 AU, the bright rings at 130 AU and 220 AU. Continuum emission is detected out to a 300 AU distance, similar to IM Lup but a factor of 5 larger than typically found for Lupus disks at this sensitivity and resolution. The bright central region of the disk (within 35 AU) is possibly optically thick at 1 mm wavelengths, and has a brightness temperature of only 13 K. The spectral index increases between the inner disk and the first ring, at the location of the first gap. Due to the low temperature of the disk midplane, snow lines can be excluded as the drivers behind the ring and gap formation in this disk. Disks the size of V1094 Sco are rare, and only 2.1+-1.5% of disks in Lupus show continuum emission beyond 200 AU. Possible connections between the large primordial disk population, transition disks, and exoplanets are discussed. △ Less

Submitted 8 May, 2018; originally announced May 2018.

Comments: 12 pages, 7 figures. Accepted for publication in Astronomy & Astrophysics

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

Radiative transfer modelling of W33A MM1: 3-D structure and dynamics of a complex massive star forming region

Authors: Andrés F. Izquierdo, Roberto Galván-Madrid, Luke T. Maud, Melvin G. Hoare, Katharine G. Johnston, Eric R. Keto, Qizhou Zhang, Willem-Jan de Wit

Abstract: We present a composite model and radiative transfer simulations of the massive star forming core W33A MM1. The model was tailored to reproduce the complex features observed with ALMA at $\approx 0.2$ arcsec resolution in CH$_3$CN and dust emission. The MM1 core is fragmented into six compact sources coexisting within $\sim 1000$ au. In our models, three of these compact sources are better represen… ▽ More We present a composite model and radiative transfer simulations of the massive star forming core W33A MM1. The model was tailored to reproduce the complex features observed with ALMA at $\approx 0.2$ arcsec resolution in CH$_3$CN and dust emission. The MM1 core is fragmented into six compact sources coexisting within $\sim 1000$ au. In our models, three of these compact sources are better represented as disc-envelope systems around a central (proto)star, two as envelopes with a central object, and one as a pure envelope. The model of the most prominent object (Main) contains the most massive (proto)star ($M_\star\approx7~M_\odot$) and disc+envelope ($M_\mathrm{gas}\approx0.4~M_\odot$), and is the most luminous ($L_\mathrm{Main} \sim 10^4~L_\odot$). The model discs are small (a few hundred au) for all sources. The composite model shows that the elongated spiral-like feature converging to the MM1 core can be convincingly interpreted as a filamentary accretion flow that feeds the rising stellar system. The kinematics of this filament is reproduced by a parabolic trajectory with focus at the center of mass of the region. Radial collapse and fragmentation within this filament, as well as smaller filamentary flows between pairs of sources are proposed to exist. Our modelling supports an interpretation where what was once considered as a single massive star with a $\sim 10^3$ au disc and envelope, is instead a forming stellar association which appears to be virialized and to form several low-mass stars per high-mass object. △ Less

Submitted 24 April, 2018; originally announced April 2018.

Comments: 22 pages, 13 figures, 3 tables. Accepted for publication in MNRAS

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

Sulphur monoxide exposes a potential molecular disk wind from the planet-hosting disk around HD100546

Authors: Alice Booth, Catherine Walsh, Mihkel Kama, Ryan A. Loomis, Luke T. Maud, Attila Juhász

Abstract: Sulphur-bearing volatiles are observed to be significantly depleted in interstellar and circumstellar regions. This missing sulphur is postulated to be mostly locked up in refractory form. With ALMA we have detected sulphur monoxide (SO), a known shock tracer, in the HD 100546 protoplanetary disk. Two rotational transitions: $J=7_{7}-6_{6}$ (301.286 GHz) and $J=7_8-6_7$ (304.078 GHz) are detected… ▽ More Sulphur-bearing volatiles are observed to be significantly depleted in interstellar and circumstellar regions. This missing sulphur is postulated to be mostly locked up in refractory form. With ALMA we have detected sulphur monoxide (SO), a known shock tracer, in the HD 100546 protoplanetary disk. Two rotational transitions: $J=7_{7}-6_{6}$ (301.286 GHz) and $J=7_8-6_7$ (304.078 GHz) are detected in their respective integrated intensity maps. The stacking of these transitions results in a clear 5$σ$ detection in the stacked line profile. The emission is compact but is spectrally resolved and the line profile has two components. One component peaks at the source velocity and the other is blue-shifted by 5 km s$^{-1}$. The kinematics and spatial distribution of the SO emission are not consistent with that expected from a purely Keplerian disk. We detect additional blue-shifted emission that we attribute to a disk wind. The disk component was simulated using LIME and a physical disk structure. The disk emission is asymmetric and best fit by a wedge of emission in the north east region of the disk coincident with a `hot-spot' observed in the CO $J=3-2$ line. The favoured hypothesis is that a possible inner disk warp (seen in CO emission) directly exposes the north-east side of the disk to heating by the central star, creating locally the conditions to launch a disk wind. Chemical models of a disk wind will help to elucidate why the wind is particularly highlighted in SO emission and whether a refractory source of sulphur is needed. An alternative explanation is that the SO is tracing an accretion shock from a circumplanetary disk associated with the proposed protoplanet embedded in the disk at 50 au. We also report a non-detection of SO in the protoplanetary disk around HD 97048. △ Less

Submitted 19 December, 2017; v1 submitted 16 December, 2017; originally announced December 2017.

Comments: Accepted for publication in A&A on 12/12/17

Journal ref: A&A 611, A16 (2018)

Phase correction for ALMA - Investigating water vapour radiometer scaling:The long-baseline science verification data case study

Authors: L. T. Maud, R. P. J. Tilanus, T. A. van Kempen, M. R. Hogerheijde, M. Schmalzl, I. Yoon, Y. Contreras, M. C. Toribio, Y. Asaki, W. R. F. Dent, E. Fomalont, S. Matsushita

Abstract: The Atacama Large millimetre/submillimetre Array (ALMA) makes use of water vapour radiometers (WVR), which monitor the atmospheric water vapour line at 183 GHz along the line of sight above each antenna to correct for phase delays introduced by the wet component of the troposphere. The application of WVR derived phase corrections improve the image quality and facilitate successful observations in… ▽ More The Atacama Large millimetre/submillimetre Array (ALMA) makes use of water vapour radiometers (WVR), which monitor the atmospheric water vapour line at 183 GHz along the line of sight above each antenna to correct for phase delays introduced by the wet component of the troposphere. The application of WVR derived phase corrections improve the image quality and facilitate successful observations in weather conditions that were classically marginal or poor. We present work to indicate that a scaling factor applied to the WVR solutions can act to further improve the phase stability and image quality of ALMA data. We find reduced phase noise statistics for 62 out of 75 datasets from the long-baseline science verification campaign after a WVR scaling factor is applied. The improvement of phase noise translates to an expected coherence improvement in 39 datasets. When imaging the bandpass source, we find 33 of the 39 datasets show an improvement in the signal-to-noise ratio (S/N) between a few to ~30 percent. There are 23 datasets where the S/N of the science image is improved: 6 by <1%, 11 between 1 and 5%, and 6 above 5%. The higher frequencies studied (band 6 and band 7) are those most improved, specifically datasets with low precipitable water vapour (PWV), <1mm, where the dominance of the wet component is reduced. Although these improvements are not profound, phase stability improvements via the WVR scaling factor come into play for the higher frequency (>450 GHz) and long-baseline (>5km) observations. These inherently have poorer phase stability and are taken in low PWV (<1mm) conditions for which we find the scaling to be most effective. A promising explanation for the scaling factor is the mixing of dry and wet air components, although other origins are discussed. We have produced a python code to allow ALMA users to undertake WVR scaling tests and make improvements to their data. △ Less

Submitted 11 July, 2017; originally announced July 2017.

Comments: 27 pages, 11 figures, 5 tables - Accepted to A&A

Journal ref: A&A 605, A121 (2017)

The ALMA view of W33A: A Spiral Filament Feeding the Candidate Disc in MM1-Main

Authors: L. T. Maud, M. G. Hoare, R. Galván-Madrid, Q. Zhang, E. Keto, K. G. Johnston, J. E. Pineda

Abstract: We targeted the massive star forming region W33A using the Atacama Large Sub/Millimeter Array (ALMA) in band 6 (230 GHz) and 7 (345 GHz) to search for a sub-1000au disc around the central O-type massive young stellar object (MYSO) W33A MM1-Main. Our data achieves a resolution of ~0.2" (~500au) and resolves the central core, MM1, into multiple components and reveals complex and filamentary structur… ▽ More We targeted the massive star forming region W33A using the Atacama Large Sub/Millimeter Array (ALMA) in band 6 (230 GHz) and 7 (345 GHz) to search for a sub-1000au disc around the central O-type massive young stellar object (MYSO) W33A MM1-Main. Our data achieves a resolution of ~0.2" (~500au) and resolves the central core, MM1, into multiple components and reveals complex and filamentary structures. There is strong molecular line emission covering the entire MM1 region. The kinematic signatures are inconsistent with only Keplerian rotation although we propose that the shift in the emission line centroids within ~1000au of MM1-Main could hint at an underlying compact disc with Keplerian rotation. We cannot however rule out the possibility of an unresolved binary or multiple system. A putative smaller disc could be fed by the large scale spiral `feeding filament' we detect in both gas and dust emission. We also discuss the nature of the now-resolved continuum sources. △ Less

Submitted 24 January, 2017; originally announced January 2017.

ALMA Long Baseline Campaigns: Phase Characteristics of Atmosphere at Long Baselines in the Millimeter and Submillimeter Wavelengths

Authors: Satoki Matsushita, Yoshiharu Asaki, Edward B. Fomalont, Koh-Ichiro Morita, Denis Barkats, Richard E. Hills, Ryohei Kawabe, Luke T. Maud, Bojan Nikolic, Remo P. J. Tilanus, Catherine Vlahakis, Nicholas D. Whyborn

Abstract: This paper presents the first detailed investigation of the characteristics of mm/submm phase fluctuation and phase correction methods obtained using ALMA with baseline lengths up to ~15 km. Most of the spatial structure functions (SSFs) show that the phase fluctuation increases as a function of baseline length, with a power-law slope of ~0.6. In many cases, we find that the slope becomes shallowe… ▽ More This paper presents the first detailed investigation of the characteristics of mm/submm phase fluctuation and phase correction methods obtained using ALMA with baseline lengths up to ~15 km. Most of the spatial structure functions (SSFs) show that the phase fluctuation increases as a function of baseline length, with a power-law slope of ~0.6. In many cases, we find that the slope becomes shallower (average of ~0.2-0.3) at baseline lengths longer than ~1 km, namely showing a turn-over in SSF. The phase correction method using water vapor radiometers (WVRs) works well, especially for the cases where PWV >1 mm, which reduces the degree of phase fluctuations by a factor of two in many cases. However, phase fluctuations still remain after the WVR phase correction, suggesting the existence of other turbulent constituent that cause the phase fluctuation. This is supported by occasional SSFs that do not exhibit any turn-over; these are only seen when the PWV is low or after WVR phase correction. This means that the phase fluctuation caused by this turbulent constituent is inherently smaller than that caused by water vapor. Since there is no turn-over in the SSF up to the maximum baseline length of ~15 km, this turbulent constituent must have scale height of 10 km or more, and thus cannot be water vapor, whose scale height is around 1 km. This large scale height turbulent constituent is likely to be water ice or a dry component. Excess path length fluctuation after the WVR phase correction at a baseline length of 10 km is large (>200 micron), which is significant for high frequency (>450 GHz or <700 micron) observations. These results suggest the need for an additional phase correction method, such as fast switching, in addition to the WVR phase correction. We simulated the fast switching, and the result suggests that it works well, with shorter cycle times linearly improving the coherence. △ Less

Submitted 11 January, 2017; originally announced January 2017.

Comments: 12 pages, 7 figures. PASP, in press

The Winds from HL Tau

Authors: Pamela D. Klaassen, Joseph C. Mottram, Luke T. Maud, Attila Juhasz

Abstract: Outflowing motions, whether a wind launched from the disk, a jet launched from the protostar, or the entrained molecular outflow, appear to be an ubiquitous feature of star formation. These outwards motions have a number of root causes, and how they manifest is intricately linked to their environment as well as the process of star formation itself. Using the ALMA Science Verification data of HL… ▽ More Outflowing motions, whether a wind launched from the disk, a jet launched from the protostar, or the entrained molecular outflow, appear to be an ubiquitous feature of star formation. These outwards motions have a number of root causes, and how they manifest is intricately linked to their environment as well as the process of star formation itself. Using the ALMA Science Verification data of HL Tau, we investigate the high velocity molecular gas being removed from the system as a result of the star formation process. We aim to place these motions in context with the optically detected jet, and the disk. With these high resolution ($\sim 1"$) ALMA observations of CO (J=1-0), we quantify the outwards motions of the molecular gas. We find evidence for a bipolar outwards flow, with an opening angle, as measured in the red-shifted lobe, starting off at 90$^\circ$, and narrowing to 60$^\circ$ further from the disk, likely because of magnetic collimation. Its outwards velocity, corrected for inclination angle is of order 2.4 km s$^{-1}$. △ Less

Submitted 25 April, 2016; originally announced April 2016.

Comments: 8 Pages, 8 figures, 3 tables. To appear in MNRAS

Submillimeter Array Observations of NGC 2264-C: Molecular Outflows and Driving Sources

Authors: Nichol Cunningham, Stuart L. Lumsden, Claudia J. Cyganowski, Luke T. Maud, Cormac Purcell

Abstract: We present 1.3mm Submillimeter Array (SMA) observations at $\sim$3$^{\prime\prime}$ resolution towards the brightest section of the intermediate/massive star forming cluster NGC 2264-C. The millimetre continuum emission reveals ten 1.3mm continuum peaks, of which four are new detections. The observed frequency range includes the known molecular jet/outflow tracer SiO (5-4), thus providing the firs… ▽ More We present 1.3mm Submillimeter Array (SMA) observations at $\sim$3$^{\prime\prime}$ resolution towards the brightest section of the intermediate/massive star forming cluster NGC 2264-C. The millimetre continuum emission reveals ten 1.3mm continuum peaks, of which four are new detections. The observed frequency range includes the known molecular jet/outflow tracer SiO (5-4), thus providing the first high resolution observations of SiO towards NGC 2264-C. We also detect molecular lines of twelve additional species towards this region, including CH$_3$CN, CH$_3$OH, SO, H$_2$CO, DCN, HC$_3$N, and $^{12}$CO. The SiO (5-4) emission reveals the presence of two collimated, high velocity (up to 30kms$^{-1}$ with respect to the systemic velocity) bi-polar outflows in NGC 2264-C. In addition, the outflows are traced by emission from $^{12}$CO, SO, H$_2$CO, and CH$_3$OH. We find an evolutionary spread between cores residing in the same parent cloud. The two unambiguous outflows are driven by the brightest mm continuum cores, which are IR-dark, molecular line weak, and likely the youngest cores in the region. Furthermore, towards the RMS source AFGL 989-IRS1, the IR-bright and most evolved source in NGC 2264-C, we observe no molecular outflow emission. A molecular line rich ridge feature, with no obvious directly associated continuum source, lies on the edge of a low density cavity and may be formed from a wind driven by AFGL 989-IRS1. In addition, 229GHz class I maser emission is detected towards this feature. △ Less

Submitted 18 February, 2016; originally announced February 2016.

Comments: 29 pages, 12 figures

A distance limited sample of massive molecular outflows

Authors: L. T. Maud, T. J. T. Moore, S. L. Lumsden, J. C. Mottram, J. S. Urquhart, M. G. Hoare

Abstract: We have observed 99 mid-infrared-bright, massive young stellar objects and compact HII regions drawn from the Red MSX source (RMS) survey in the J=3$-$2 transition of $^{12}$CO and $^{13}$CO, using the James Clerk Maxwell Telescope. 89 targets are within 6 kpc of the Sun, covering a representative range of luminosities and core masses. These constitute a relatively unbiased sample of bipolar molec… ▽ More We have observed 99 mid-infrared-bright, massive young stellar objects and compact HII regions drawn from the Red MSX source (RMS) survey in the J=3$-$2 transition of $^{12}$CO and $^{13}$CO, using the James Clerk Maxwell Telescope. 89 targets are within 6 kpc of the Sun, covering a representative range of luminosities and core masses. These constitute a relatively unbiased sample of bipolar molecular outflows associated with massive star formation. Of these, 59, 17 and 13 sources (66, 19 and 15 percent) are found to have outflows, show some evidence of outflow, and have no evidence of outflow, respectively. The time-dependent parameters of the high-velocity molecular flows are calculated using a spatially variable dynamic timescale. The canonical correlations between the outflow parameters and source luminosity are recovered and shown to scale with those of low-mass sources. For coeval star formation we find the scaling is consistent with all the protostars in an embedded cluster providing the outflow force, with massive stars up to $\sim$30 M$_{\odot}$ generating outflows. Taken at face value, the results support the model of a scaled-up version of the accretion-related outflow-generation mechanism associated with discs and jets in low-mass objects with time-averaged accretion rates of $\sim$10$^{-3}$ M$_{\odot}$ yr$^{-1}$ onto the cores. However, we also suggest an alternative model, in which the molecular outflow dynamics are dominated by the entrained mass and are unrelated to the details of the acceleration mechanism. We find no evidence that outflows contribute significantly to the turbulent kinetic energy of the surrounding dense cores. △ Less

Submitted 1 September, 2015; originally announced September 2015.

Journal ref: 2015MNRAS.453..645M

A distance limited sample of massive star forming cores from the RMS survey

Authors: L. T. Maud, S. L. Lumsden, T. J. T. Moore, J. C. Mottram, J. S. Urquhart, A. Cicchini

Abstract: We analyse C$^{18}$O ($J=3-$2) data from a sample of 99 infrared-bright massive young stellar objects (MYSOs) and compact HII regions that were identified as potential molecular-outflow sources in the Red MSX source (RMS) survey. We extract a distance limited (D $<$ 6 kpc) sample shown to be representative of star formation covering the transition between the source types. At the spatial resolutio… ▽ More We analyse C$^{18}$O ($J=3-$2) data from a sample of 99 infrared-bright massive young stellar objects (MYSOs) and compact HII regions that were identified as potential molecular-outflow sources in the Red MSX source (RMS) survey. We extract a distance limited (D $<$ 6 kpc) sample shown to be representative of star formation covering the transition between the source types. At the spatial resolution probed, Larson-like relationships are found for these cores, though the alternative explanation, that Larson's relations arise where surface-density-limited samples are considered, is also consistent with our data. There are no significant differences found between source properties for the MYSOs and HII regions, suggesting that the core properties are established prior to the formation of massive stars, which subsequently have little impact at the later evolutionary stages investigated. There is a strong correlation between dust-continuum and C$^{18}$O-gas masses, supporting the interpretation that both trace the same material in these IR-bright sources. A clear linear relationship is seen between the independently established core masses and luminosities. The position of MYSOs and compact HII regions in the mass-luminosity plane is consistent with the luminosity expected a cluster of protostars when using a $\sim$40 percent star-formation efficiency and indicates that they are at a similar evolutionary stage, near the end of the accretion phase. △ Less

Submitted 1 September, 2015; originally announced September 2015.

Journal ref: 2015MNRAS.452..637M

Disk evolution in the solar neighborhood. I Disk frequencies from 1 to 100 Myr

Authors: Álvaro Ribas, Bruno Merín, Hervé Bouy, Luke T. Maud

Abstract: We study the evolution of circumstellar disks in 22 young (1 to 100 Myr) nearby (within 500 pc) associations over the entire mass spectrum using photometry covering from the optical to the mid-infrared. We compiled a catalog of 2340 spectroscopically-confirmed members of these nearby associations. We analyzed their spectral energy distributions and searched for excess related to the presence of pr… ▽ More We study the evolution of circumstellar disks in 22 young (1 to 100 Myr) nearby (within 500 pc) associations over the entire mass spectrum using photometry covering from the optical to the mid-infrared. We compiled a catalog of 2340 spectroscopically-confirmed members of these nearby associations. We analyzed their spectral energy distributions and searched for excess related to the presence of protoplanetary disks in a homogeneous way. Sensitivity limits and spatial completeness were also considered. We derive disk fractions as probed by mid-infrared excess in these regions. The unprecedented size of our sample allows us to confirm the timescale of disk decay reported in the literature and to find new trends. The fraction of excess sources increases systematically if measured at longer wavelengths. Disk percentages derived using different wavelength ranges should therefore be compared with caution. The dust probed at 22-24 um evolves slower than that probed at shorter wavelengths (3.4-12 um). Assuming an exponential decay, we derive a timescale tau=4.2-5.8 Myr at 22-24 um for primordial disks, compared to 2-3 Myr at shorter wavelength (3.4-12 um). Primordial disks disappear around 10 Myr, matching in time a brief increase of the number of 'evolved' disks. The increase in timescale of excess decay at longer wavelength is compatible with inside-out disk clearing scenarios. The increased timescale of decay and larger dispersion in the distribution of disk fractions at 22-24 um suggest that the inner and outer zones evolve differently, the latter potentially following a variety of evolutionary paths. The drop of primordial disks and the coincident rise of evolved disks at 10 Myr are compatible with planet formation theories suggesting that the disappearance of the gas is immediately followed by the dynamical stirring of the disk. △ Less

Submitted 11 December, 2013; v1 submitted 2 December, 2013; originally announced December 2013.

Comments: Accepted for publication in A&A: 13 November 2013. 12 pages, 7 figures

Investigation of dust properties of the proto-planetary nebula IRAS 18276-1431

Authors: K. Murakawa, H. Izumiura, R. D. Oudmaijer, L. T. Maud

Abstract: We investigate the circumstellar dust properties of the oxygen-rich bipolar proto-planetary nebula IRAS 18276-1431 by means of two-dimensional radiative transfer simulations of the circumstellar dust shell. The model geometry is assumed to have a torus and an envelope. The parameters of the dust and the dust shell are constrained by comparing the SED and NIR intensity and polarisation data with th… ▽ More We investigate the circumstellar dust properties of the oxygen-rich bipolar proto-planetary nebula IRAS 18276-1431 by means of two-dimensional radiative transfer simulations of the circumstellar dust shell. The model geometry is assumed to have a torus and an envelope. The parameters of the dust and the dust shell are constrained by comparing the SED and NIR intensity and polarisation data with the models. The polarisation in the envelope reaches 50 -- 60 % and is nearly constant in the H and K_S bands in the observations. This weak wavelength dependence of the polarisation can be reproduced with a grain size distribution function for the torus: 0.05 micron <= a with n(a)=a^{-(p=5.5)}exp(-a/{a_c=0.3 micron}). The power index p is significantly steeper than that for interstellar dust. Similar results have also been found in some other PPNs and suggest that mechanisms that grind down large particles may also have acted when the dust particles formed. The spectral opacity index beta is found to be 0.6+/-0.5 from the millimeter fluxes. This low value indicates the presence of large dust grains in the torus. We discuss two possible dust models for the torus. One has a size distribution function of 1.0 micron <= a <= a_max=5,000.0 micron with n(a)=a^{-(p=2.5)} and the other is 1.0 micron <= a <= a_max=10,000.0 micron with n(a)=a^{-(p=3.5)}. The former has beta of 0.633, but we are not able to find reasonable geometry parameters to fit the SED in the infrared. The latter has beta of 1.12, but reproduces the SED better over a wide wavelength range. With this dust model, the geometric parameters are estimated as follows: the inner and outer radii are 30 AU and 1000 AU and the torus mass is 3.0 M_sun. Assuming an expansion velocity of 15 kms^{-1}, the torus formation time and mass-loss rate are found to be \sim300 yrs and \sim10^{-2}M_sun yr^{-1} respectively. △ Less

Submitted 28 January, 2013; originally announced January 2013.

Comments: 9 pages, 5 figures, accepted for publication in MNRAS

CO bandhead emission of massive young stellar objects: determining disc properties

Authors: J. D. Ilee, H. E. Wheelwright, R. D. Oudmaijer, W. J. de Wit, L. T. Maud, M. G. Hoare, S. L. Lumsden, T. J. T. Moore, J. S. Urquhart, J. C. Mottram

Abstract: Massive stars play an important role in many areas of astrophysics, but numerous details regarding their formation remain unclear. In this paper we present and analyse high resolution (R ~ 30,000) near-infrared 2.3 micron spectra of 20 massive young stellar objects from the RMS database, in the largest such study of CO first overtone bandhead emission to date. We fit the emission under the assumpt… ▽ More Massive stars play an important role in many areas of astrophysics, but numerous details regarding their formation remain unclear. In this paper we present and analyse high resolution (R ~ 30,000) near-infrared 2.3 micron spectra of 20 massive young stellar objects from the RMS database, in the largest such study of CO first overtone bandhead emission to date. We fit the emission under the assumption it originates from a circumstellar disc in Keplerian rotation. We explore three approaches to modelling the physical conditions within the disc - a disc heated mainly via irradiation from the central star, a disc heated mainly via viscosity, and a disc in which the temperature and density are described analytically. We find that the models described by heating mechanisms are inappropriate because they do not provide good fits to the CO emission spectra. We therefore restrict our analysis to the analytic model, and obtain good fits to all objects that possess sufficiently strong CO emission, suggesting circumstellar discs are the source of this emission. On average, the temperature and density structure of the discs correspond to geometrically thin discs, spread across a wide range of inclinations. Essentially all the discs are located within the dust sublimation radius, providing strong evidence that the CO emission originates close to the central protostar, on astronomical unit scales. In addition, we show that the objects in our sample appear no different to the general population of MYSOs in the RMS database, based on their near- and mid-infrared colours. The combination of observations of a large sample of MYSOs with CO bandhead emission and our detailed modelling provide compelling evidence of the presence of small scale gaseous discs around such objects, supporting the scenario in which massive stars form via disc accretion. △ Less

Submitted 3 December, 2012; originally announced December 2012.

Comments: 14 pages, 8 figures, accepted for publication in MNRAS