On the ALMA observability of nascent massive multiple systems formed by gravitational instability

This is the final version. Available from Oxford University Press via the DOI in this record.Massive young stellar object (MYSOs) form during the collapse of high-mass pre-stellar cores, where infalling molecular material is accreted through a centrifugally-balanced accretion disc that is subject to efficient gravitational instabilities. In the resulting fragmented accretion disc of the MYSO, gaseous clumps and low-mass stellar companions can form, which will influence the future evolution of massive protostars in the Hertzsprung-Russell diagram. We perform dust continuum radiative transfer calculations and compute synthetic images of disc structures modelled by the gravito-radiation-hydrodynamics simulation of a forming MYSO, in order to investigate the Atacama Large Millimeter/submillimeter Array (ALMA) observability of circumstellar gaseous clumps and forming multiple systems. Both spiral arms and gaseous clumps located at ~a few 100 au from the protostar can be resolved by interferometric ALMA Cycle 7 C43-8 and C43-10 observations at band 6 (1.2 mm), using a maximal 0.015" beam angular resolution and at least 10-30 min exposure time for sources at distances of 1-2 kpc. Our study shows that substructures are observable regardless of their viewing geometry or can be inferred in the case of an edge-viewed disc. The observation probability of the clumps increases with the gradually increasing efficiency of gravitational instability at work as the disc evolves. As a consequence, large discs around MYSOs close to the zero-age-main-sequence line exhibit more substructures than at the end of the gravitational collapse. Our results motivate further observational campaigns devoted to the close surroundings of the massive protostars S255IR-NIRS3 and NGC 6334I-MM1, whose recent outbursts are a probable signature of disc fragmentation and accretion variability.European Research Council (ERC)Science and Technology Facilities Council (STFC)Russian Science FoundationSwiss National Science Foundatio

The burst mode of accretion in massive star formation with stellar inertia

This is the author accepted manuscript. The final version is available from Oxford University Press via the DOI in this recordData availability: This research made use of the PLUTO code developed at the University of Torino by A. Mignone (http://plutocode.ph.unito.it/). The figures have been produced using the Matplotlib plotting library for the Python programming language (https://matplotlib.org/). The data underlying this article will be shared on reasonable request to the corresponding author.The burst mode of accretion in massive star formation is a scenario linking the initial gravitational collapse of parent pre-stellar cores to the properties of their gravitationally unstable discs and of their accretion-driven bursts. In this study, we present a series of high-resolution 3D radiation-hydrodynamics numerical simulations for young massive stars formed out of collapsing 100 Mo molecular cores spinning with several values of the ratio of rotational-to-gravitational energies beta=5%-9%. The models include the indirect gravitational potential caused by disc asymmetries. We find that this modifies the barycenter of the disc, causing significant excursions of the central star position, which we term stellar wobbling. The stellar wobbling slows down and protracts the development of gravitational instability in the disc, reducing the number and magnitude of the accretion-driven bursts undergone by the young massive stars, whose properties are in good agreement with that of the burst monitored from the massive protostar M17 MIR. Including stellar wobbling is therefore important for accurate modeling disc structures. Synthetic ALMA interferometric images in the millimeter waveband show that the outcomes of efficient gravitational instability such as spiral arms and gaseous clumps can be detected for as long as the disc is old enough and has already entered the burst mode of accretion.Science and Technology Facilities Council (STFC)Ministry of Science and Higher Education of the Russian FederationEuropean Research Council (ERC)MOS

Are There Age Spreads in Star Forming Regions?

A luminosity spread at a given effective temperature is ubiquitously seen in the Hertzsprung-Russell (HR) diagrams of young star forming regions and often interpreted in terms of a prolonged period (>=10 Myr) of star formation. I review the evidence that the observed luminosity spreads are genuine and not caused by astrophysical sources of scatter. I then address whether the luminosity spreads necessarily imply large age spreads, by comparing HR diagram ages with ages from independent clocks such as stellar rotation rate, the presence of circumstellar material and lithium depletion. I argue that whilst there probably is a true luminosity dispersion, there is little evidence to support age spreads larger than a few Myr. This paradox could be resolved by brief periods of rapid accretion during the class I pre main-sequence phase.Comment: To appear in the proceedings of JENAM10: Star Clusters in the Era of Large Surveys, 8 page

A 1.3 mm SMA Survey of 29 Variable Young Stellar Objects

Young stellar objects (YSOs) may undergo periods of active accretion (outbursts), during which the protostellar accretion rate is temporarily enhanced by a few orders of magnitude. Whether or not these accretion outburst YSOs possess similar dust/gas reservoirs to each other, and whether or not their dust/gas reservoirs are similar as quiescent YSOs, are issues not yet clarified. The aim of this work is to characterize the millimeter thermal dust emission properties of a statistically significant sample of long and short duration accretion outburst YSOs (i.e., FUors and EXors) and the spectroscopically identified candidates of accretion outbursting YSOs (i.e., FUor-like objects). We have carried out extensive Submillimeter Array (SMA) observations mostly at $\sim$225 GHz (1.33 mm) and $\sim$272 GHz (1.10 mm), from 2008 to 2017. We covered accretion outburst YSOs located at $$3-$\sigma$ significance. Detected sources except for the two cases of V883 Ori and NGC 2071 MM3 were observed with $\sim$1$"$ angular resolution. Overall our observed targets show a systematically higher millimeter luminosity distribution than those of the $M_{*}>$0.3 $M_{\odot}$ Class II YSOs in the nearby ($\lesssim$400 pc) low-mass star-forming molecular clouds (e.g., Taurus, Lupus, Upp Scorpio, and Chameleon I). In addition, at 1 mm our observed confirmed binaries or triple-system sources are systematically fainter than the rest of the sources even though their 1 mm fluxes are broadly distributed. We may have detected $\sim$30-60\% millimeter flux variability from V2494 Cyg and V2495 Cyg, from the observations separated by $\sim$1 year

A concordant scenario to explain FU Orionis from deep centimeter and millimeter interferometric observations

© ESO, 2017. Aims. The aim of this work is to constrain properties of the disk around the archetype FU Orionis object, FU Ori, with as good as ~25 au resolution. Methods. We resolved FU Ori at 29-37 GHz using the Karl G. Jansky Very Large Array (JVLA) in the A-Array configuration, which provided the highest possible angular resolution to date at this frequency band (~0\hbox{\farcs}07). We also performed complementary JVLA 8-10 GHz observations, Submillimeter Array (SMA) 224 GHz and 272 GHz observations, and compared these with archival Atacama Large Millimeter Array (ALMA) 346 GHz observations to obtain the spectral energy distributions (SEDs). Results. Our 8-10 GHz observations do not find evidence for the presence of thermal radio jets, and constrain the radio jet/wind flux to at least 90 times lower than the expected value from the previously reported bolometric luminosity-radio luminosity correlation. The emission at frequencies higher than 29 GHz may be dominated by the two spatially unresolved sources, which are located immediately around FU Ori and its companion FU Ori S, respectively. Their deconvolved radii at 33 GHz are only a few au, which is two orders of magnitude smaller in linear scale than the gaseous disk revealed by the previous Subaru-HiCIAO 1.6 μm coronagraphic polarization imaging observations. We are struck by the fact that these two spatially compact sources contribute to over 50% of the observed fluxes at 224 GHz, 272 GHz, and 346 GHz. The 8-346 GHz SEDs of FU Ori and FU Ori S cannot be fit by constant spectral indices (over frequency), although we cannot rule out that it is due to the time variability of their (sub)millimeter fluxes. Conclusions. The more sophisticated models for SEDs considering the details of the observed spectral indices in the millimeter bands suggest that the >29 GHz emission is contributed by a combination of free-free emission from ionized gas and thermal emission from optically thick and optically thin dust components. We hypothesize that dust in the innermost parts of the disks (â0.1 au) has been sublimated, and thus the disks are no longer well shielded against the ionizing photons. The estimated overall gas and dust mass based on SED modeling, can be as high as a fraction of a solar mass, which is adequate for developing disk gravitational instability. Our present explanation for the observational data is that the massive inflow of gas and dust due to disk gravitational instability or interaction with a companion/intruder, was piled up at the few-Au scale due to the development of a deadzone with negligible ionization. The piled up material subsequently triggered the thermal instability and the magnetorotational instability when the ionization fraction in the inner sub-Au scale region exceeded a threshold value, leading to the high protostellar accretion rate

A concordant scenario to explain FU Orionis from deep centimeter and millimeter interferometric observations

© ESO, 2017. Aims. The aim of this work is to constrain properties of the disk around the archetype FU Orionis object, FU Ori, with as good as ~25 au resolution. Methods. We resolved FU Ori at 29-37 GHz using the Karl G. Jansky Very Large Array (JVLA) in the A-Array configuration, which provided the highest possible angular resolution to date at this frequency band (~0\hbox{\farcs}07). We also performed complementary JVLA 8-10 GHz observations, Submillimeter Array (SMA) 224 GHz and 272 GHz observations, and compared these with archival Atacama Large Millimeter Array (ALMA) 346 GHz observations to obtain the spectral energy distributions (SEDs). Results. Our 8-10 GHz observations do not find evidence for the presence of thermal radio jets, and constrain the radio jet/wind flux to at least 90 times lower than the expected value from the previously reported bolometric luminosity-radio luminosity correlation. The emission at frequencies higher than 29 GHz may be dominated by the two spatially unresolved sources, which are located immediately around FU Ori and its companion FU Ori S, respectively. Their deconvolved radii at 33 GHz are only a few au, which is two orders of magnitude smaller in linear scale than the gaseous disk revealed by the previous Subaru-HiCIAO 1.6 μm coronagraphic polarization imaging observations. We are struck by the fact that these two spatially compact sources contribute to over 50% of the observed fluxes at 224 GHz, 272 GHz, and 346 GHz. The 8-346 GHz SEDs of FU Ori and FU Ori S cannot be fit by constant spectral indices (over frequency), although we cannot rule out that it is due to the time variability of their (sub)millimeter fluxes. Conclusions. The more sophisticated models for SEDs considering the details of the observed spectral indices in the millimeter bands suggest that the >29 GHz emission is contributed by a combination of free-free emission from ionized gas and thermal emission from optically thick and optically thin dust components. We hypothesize that dust in the innermost parts of the disks (â0.1 au) has been sublimated, and thus the disks are no longer well shielded against the ionizing photons. The estimated overall gas and dust mass based on SED modeling, can be as high as a fraction of a solar mass, which is adequate for developing disk gravitational instability. Our present explanation for the observational data is that the massive inflow of gas and dust due to disk gravitational instability or interaction with a companion/intruder, was piled up at the few-Au scale due to the development of a deadzone with negligible ionization. The piled up material subsequently triggered the thermal instability and the magnetorotational instability when the ionization fraction in the inner sub-Au scale region exceeded a threshold value, leading to the high protostellar accretion rate