Massive-star formation involves significant ionization in the innermost regions near the central object, such as gravitationally trapped H ii regions, jets, ionized disks, or winds. Resolved observations of the associated continuum and recombination line emission are crucial for guiding theory. The next-generation Very Large Array (ngVLA) will enable unprecedented observations of thermal emission with 1 mas resolution, providing a new perspective on massive-star formation at scales down to a few astronomical units at kiloparsec distances. This work presents synthetic interferometric ngVLA observations of the free–free continuum (93 GHz band), H41 α , and H38 α recombination lines from ionized jets and disks around massive protostars. Using the sf3dmodels Python package, we generate gas distributions based on analytical models, which are then processed through the RADMC-3D radiative transfer code. Our results indicate that the ngVLA can easily resolve, both spatially and spectrally, the ionized jet from a 15 M _⊙ protostar at 700 pc, distinguishing between collimated jets and wide-angle winds, and resolving their launching radii, widths, and any substructure down to a few astronomical units. Detailed studies of radio jets launched by massive protostars will be feasible up to distances of ​​​​​​∼2 kpc. Furthermore, ngVLA will be able to study in detail the ionized disks around massive (>10 M _⊙ ) protostars up to distances from 4 to 12 kpc, resolving their kinematics and enabling the measurement of their central masses across the Galaxy. These observations can be conducted with on-source integrations of only a few hours.

Context. The fragmentation of high-mass star-forming regions depends on a variety of physical parameters, including density, the magnetic field, and turbulent gas properties. Aims. We evaluate the importance of the density and magnetic field structures in relation to the fragmentation properties during high-mass star formation. Methods. Observing the large parsec-scale Stokes I millimeter dust continuum emission with the IRAM 30 m telescope and the intermediate-scale (<0.1 pc) polarized submillimeter dust emission with the Submillimeter Array toward a sample of 20 high-mass star-forming regions allows us to quantify the dependence of the fragmentation behavior of these regions on the density and magnetic field structures. Results. Based on the IRAM 30 m data, we infer density distributions n ∝ r−p of the regions with typical power-law slopes p around ~1.5. There is no obvious correlation between the power-law slopes of the density structures on larger clump scales (~1 pc) and the number of fragments on smaller core scales (<0.1 pc). Comparing the large-scale single-dish density profiles to those derived earlier from interferometric observations at smaller spatial scales, we find that the smaller-scale power-law slopes are steeper, typically around ~2.0. The flattening toward larger scales is consistent with the star-forming regions being embedded in larger cloud structures that do not decrease in density away from a particular core. The magnetic fields of several regions appear to be aligned with filamentary structures that lead toward the densest central cores. Furthermore, we find different polarization structures; some regions exhibit central polarization holes, whereas other regions show polarized emission also toward the central peak positions. Nevertheless, the polarized intensities are inversely related to the Stokes I intensities, following roughly a power-law slope of ∝ SI−0.62. We estimate magnetic field strengths between ~0.2 and ~4.5 mG, and we find no clear correlation between magnetic field strength and the fragmentation level of the regions. A comparison of the turbulent to magnetic energies shows that they are of roughly equal importance in this sample. The mass-to-flux ratios range between ~2 and ~7, consistent with collapsing star-forming regions. Conclusions. Finding no clear correlations between the present-day large-scale density structure, the magnetic field strength, and the smaller-scale fragmentation properties of the regions, indicates that the fragmentation of high-mass star-forming regions may not be affected strongly by the initial density profiles and magnetic field properties. However, considering the limited evolutionary range and spatial scales of the presented CORE analysis, future research directions should include density structure analysis of younger regions that better resemble the initial conditions, as well as connecting the observed intermediate-scale magnetic field structure with the larger-scale magnetic fields of the parental molecular clouds.

The formation of massive objects via gravitational collapse is relevant both for explaining the origin of the first supermassive black holes and in the context of massive star formation. Here, we analyze simulations of the formation of massive objects pursued by different groups and in various environments, concerning the formation of supermassive black holes, primordial stars, as well as present-day massive stars. We focus here particularly on the regime of small virial parameters, that is, low ratios of the initial kinetic to gravitational energy, low to moderate Mach numbers, and the phase before feedback is very efficient. We compare the outcomes of collapse under different conditions using dimensionless parameters, particularly the star formation efficiency *, the fraction * of mass in the most massive object relative to the total stellar mass, and the fraction (tot) of mass of the most massive object as a function of the total mass. We find that in all simulations analyzed here, (tot) increases as a function of *, although the steepness of the increase depends on the environment. The relation between * and * is found to be more complex and also strongly depends on the number of protostars present at the beginning of the simulations. We show that a collision parameter, estimated as the ratio of the system size divided by the typical collision length, allows us to approximately characterize whether collisions are important. A high collision parameter implies a steeper increase in the relation between (tot) and *. We analyze the statistical correlation between the dimensionless quantities using the Spearman coefficient and further confirm via a machine learning analysis that good predictions of * can be obtained from * together with a rough estimate of the collision parameter. This suggests that a good estimate of the mass of the most massive object can be obtained once the maximum efficiency for a given environment is known and an estimate for the collision parameter has been determined.

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.

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

We present Atacama Large Millimeter/submillimeter Array observations with a 800 au resolution and radiative-transfer modeling of the inner part ( r ≈ 6000 au) of the ionized accretion flow around a compact star cluster in formation at the center of the luminous ultracompact H ii region G10.6-0.4. We modeled the flow with an ionized Keplerian disk with and without radial motions in its outer part, or with an external Ulrich envelope. The Markov Chain Monte Carlo fits to the data give total stellar masses M ⋆ from 120 to 200 M ⊙ , with much smaller ionized-gas masses M ion-gas = 0.2–0.25 M ⊙ . The stellar mass is distributed within the gravitational radius R g ≈ 1000 to 1500 au, where the ionized gas is bound. The viewing inclination angle from the face-on orientation is i = 49°–56°. Radial motions at radii r > R g converge to v r ,0 ≈ 8.7 km s −1 , or about the speed of sound of ionized gas, indicating that this gas is marginally unbound at most. From additional constraints on the ionizing-photon rate and far-IR luminosity of the region, we conclude that the stellar cluster consists of a few massive stars with M star = 32–60 M ⊙ , or one star in this range of masses accompanied by a population of lower-mass stars. Any active accretion of ionized gas onto the massive (proto)stars is residual. The inferred cluster density is very large, comparable to that reported at similar scales in the Galactic center. Stellar interactions are likely to occur within the next million years.

"Statistical and thermal physics: with computer applications, 2nd edition." Contemporary Physics, ahead-of-print(ahead-of-print), p. 1

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

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

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

Context. ALMA observations at 1.4 mm and ~0.′′2 (~750 au) 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 if 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–500 au; their masses range from ~15 to ~26 M⊙; and their number densities are several 109 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 brightest in the core. The spectral energy distribution of these two sources and their masses and sizes are similar and suggest a common origin. Source C has not been detected at centimeter wavelengths, while source D has been clearly detected at 1.3 cm. Source D is likely the driving source of an E–W SiO outflow previously detected in the region, which suggests that the free-free emission might be coming from a radio jet. Conclusions. The observations have confirmed that the Main core in G31 is collapsing, that it has undergone fragmentation, and that its homogeneous appearance previously observed at short wavelengths is a consequence of both high dust opacity and insufficient angular resolution. The low level of fragmentation together with the fact that the core is moderately magnetically supercritical, suggests that G31 could have undergone a phase of magnetically regulated evolution characterized by a reduced fragmentation efficiency, eventually leading to the formation of a small number of relatively massive dense cores.

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

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

Aims. Current star formation research centers the characterization of the physical and chemical properties of massive stars, which are in the process of formation, at the spatial resolution of individual high-mass cores. Methods. We use sub-arcsecond resolution (~0.′′4) observations with the NOrthern Extended Millimeter Array at 1.37 mm to study the dust emission and molecular gas of 18 high-mass star-forming regions. With distances in the range of 0.7−5.5 kpc, this corresponds to spatial scales down to 300−2300 au that are resolved by our observations. We combined the derived physical and chemical properties of individual cores in these regions to estimate their ages. The temperature structures of these regions are determined by fitting the H 2 CO and CH 3 CN line emission. The density profiles are inferred from the 1.37 mm continuum visibilities. The column densities of 11 different species are determined by fitting the emission lines with XCLASS . Results. Within the 18 observed regions, we identified 22 individual cores with associated 1.37 mm continuum emission and with a radially decreasing temperature profile. We find an average temperature power-law index of q = 0.4 ± 0.1 and an average density power-law index of p = 2.0 ± 0.2 on scales that are on the order of several 1000 au. Comparing these results with values of p derived from the literature presumes that the density profiles remain unchanged from clump to core scales. The column densities relative to N (C 18 O) between pairs of dense gas tracers show tight correlations. We applied the physical-chemical model MUlti Stage ChemicaL codE to the derived column densities of each core and find a mean chemical age of ~60 000 yr and an age spread of 20 000−100 000 yr. With this paper, we release all data products of the CORE project. Conclusions. The CORE sample reveals well-constrained density and temperature power-law distributions. Furthermore, we characterized a large variety in molecular richness that can be explained by an age spread that is then confirmed by our physical-chemical modeling. The hot molecular cores show the greatest number of emission lines, but we also find evolved cores at an evolutionary stage in which most molecules are destroyed and, thus, the spectra appear line-poor once again.

Context. Star formation (SF) is a multi-scale process in which the mode of fragmentation of the collapsing clump on scales of 0.1–1 pc determines the mass reservoir and affects the accretion process of the individual protostars on scales of 10–100 au. Aims. We want to investigate the nearby (located at 1.63 ± 0.05 kpc) high-mass star-forming region IRAS 21078+5211 at linear scales from ~1 pc down to ~10 au. Methods. We combine the data of two recent programs: the NOrthern Extended Millimeter Array large project CORE and the Protostellar Outflows at the EarliesT Stages (POETS) survey. The former provides images of the 1 mm dust continuum and molecular line emissions with a linear resolution of ≈600 au covering a field of view up to ≈0.5 pc. The latter targets the ionized gas and 22 GHz water masers, mapping linear scales from a few 103 au down to a few astronomical units. Results. In IRAS 21078+5211, a highly fragmented cluster (size ~0.1 pc) of molecular cores is observed, located at the density peak of an elongated (size ~1 pc) molecular cloud. A small (≈1 km s−1 per 0.1 pc) LSR velocity (VLSR) gradient is detected across the major 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 ≈10−4 M⊙ 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. The masses of the young stellar objects (YSOs) inside these three cores are estimated in the range 1–6 M⊙. We reveal an extended (size ~0.1 pc), bipolar collimated molecular outflow emerging from core 1. We believe this is powered by the compact (size ≲1000 au) radio jet discovered in the POETS survey, 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−1 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⊙. The water masers observed in the POETS survey 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 ≈ 16 au and the corresponding launching radius and terminal velocity to be ≈ 2.2 au and ≈ 200 km s−1, respectively. Conclusions. Complementing high-angular resolution, centimeter and millimeter interferometric observations in thermal tracers with Very Long Baseline Interferometry of molecular masers, is invaluable in studying high-mass SF. The combination of these twodatasets allows us to connect the events that we see at large scales, as clump fragmentation and mass flows, with the physical processes identified at small scales, specifically, accretion and ejection in disk-jet systems.

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

C. Gieser1,??, H. Beuther1, D. Semenov1,2, A. Ahmadi3, S. Suri1, T. Möller4, M. T. Beltrán5, P. Klaassen6, Q. Zhang7, J. S. Urquhart8, Th. Henning1, S. Feng9,10,11, R. Galván-Madrid12, V. de Souza Magalhães13, L. Moscadelli5, S. Longmore14, S. Leurini15, R. Kuiper16, T. Peters17, K. M. Menten18, T. Csengeri19, G. Fuller20, F. Wyrowski18, S. Lumsden21, Á. Sánchez-Monge4, L. Maud22, H. Linz1, A. Palau12, P. Schilke4, J. Pety13,23, R. Pudritz24, J. M. Winters13, and V. Piétu13

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

We present new Jansky Very Large Array (VLA) images of the central region of the W49A star-forming region at 3.6 cm and at 7 mm at resolutions of 015 (1650 au) and 004 (440 au), respectively. The 3.6 cm data reveal new morphological detail in the ultracompact H ii region population, as well as several previously unknown and unresolved sources. In particular, source A shows elongated, edge-brightened bipolar lobes, indicative of a collimated outflow, and source E is resolved into three spherical components. We also present VLA observations of radio recombination lines at 3.6 cm and 7 mm, and IRAM Northern Extended Millimeter Array (NOEMA) observations at 1.2 mm. Three of the smallest ultracompact H ii regions (sources A, B2, and G2) all show broad kinematic linewidths, with Delta V-FWHM greater than or similar to 40 km s(-1). A multi-line analysis indicates that broad linewidths remain after correcting for pressure broadening effects, suggesting the presence of supersonic flows. Substantial changes in linewidth over the 21 yr time baseline at both 3.6 cm and 7 mm are found for source G2. At 3.6 cm, the linewidth of G2 changed from 31.7 1.8 km s(-1) to 55.6 2.7 km s(-1), an increase of +23.9 3.4 km s(-1). The G2 source was previously reported to have shown a 3.6 cm continuum flux density decrease of 40% between 1994 and 2015. This source sits near the center of a very young bipolar outflow whose variability may have produced these changes.

Context.High-mass star formation typically takes place in a crowded environment, with a higher likelihood of young forming stars affecting and being affected by their surroundings and neighbours, as well as links between different physical scales affecting the outcome. However, observational studies are often focused on either clump or disc scales exclusively.Aims.We explore the physical and chemical links between clump and disc scales in the high-mass star formation region W3 IRS4, a region that contains a number of different evolutionary phases in the high-mass star formation process, as a case-study for what can be achieved as part of the IRAM NOrthern Extended Millimeter Array (NOEMA) large programme named CORE: “Fragmentation and disc formation in high-mass star formation”.Methods.We present 1.4 mm continuum and molecular line observations with the IRAM NOEMA interferometer and 30 m telescope, which together probe spatial scales from ~0.3−20′′ (600−40 000 AU or 0.003−0.2 pc at 2 kpc, the distance to W3). As part of our analysis, we usedXCLASSto constrain the temperature, column density, velocity, and line-width of the molecular emission lines.Results.The W3 IRS4 region includes a cold filament and cold cores, a massive young stellar object (MYSO) embedded in a hot core, and a more evolved ultra-compact (UC)H IIregion, with some degree of interaction between all components of the region that affects their evolution. A large velocity gradient is seen in the filament, suggesting infall of material towards the hot core at a rate of 10−3−10−4M⊙yr−1, while the swept up gas ring in the photodissociation region around the UCH IIregion may be squeezing the hot core from the other side. There are no clear indications of a disc around the MYSO down to the resolution of the observations (600 AU). A total of 21 molecules are detected, with the abundances and abundance ratios indicating that many molecules were formed in the ice mantles of dust grains at cooler temperatures, below the freeze-out temperature of CO (≲35 K). This contrasts with the current bulk temperature of ~50 K, which was obtained from H2CO.Conclusions.CORE observations allow us to comprehensively link the different structures in the W3 IRS4 region for the first time. Our results argue that the dynamics and environment around the MYSO W3 IRS4 have a significant impact on its evolution. This context would be missing if only high resolution or continuum observations were available.