Search

We present new results of a numerical study of the pumping of 4.8 and 14.5 GHz maser of o– |$\mathrm{H_2CO}$| in the presence of a free–free radiation field. It is shown that in the presence of a free–free radiation field inversion of not only the |$1_{10}{\!-\!}1_{11}$|⁠ , but also the |$2_{11}{\!-\!}2_{12}$| and other doublet state transitions occur. Further results are presented to illustrate how, as a consequence of the pumping scheme, the inversion of the |$1_{10}{\!-\!}1_{11}$| and |$2_{11}{\!-\!}2_{12}$| transitions respond to the free–free radiation fields associated with H  ii regions with different emission measures and levels of geometric dilution with respect to the masing region. We also discuss the criticism raised in the past by various authors against the pumping of the 4.8 GHz |$\mathrm{H_2CO}$| masers by a free–free radiation field. It is argued that the rarity of the |$\mathrm{H_2CO}$| masers is not to be ascribed to the pumping scheme but to other factors such as, e.g. the evolution of the associated H  ii region or the chemical evolution of the star-forming region that determines the |$\mathrm{H_2CO}$| abundance or a combination of both. [ABSTRACT FROM AUTHOR]

Population III stars are characterized by extremely low metallicities as they are thought to be formed from a pristine gas in the early Universe. Although the existence of Population III stars is widely accepted, the lack of direct observational evidence hampers the study of the nature of the putative stars. In this article, we explore the possibilities of constraining the nature of the oldest stars by using the luminosity function of their remnants – white dwarfs. We study the formation and evolution of white dwarf populations by following star formation in a Milky Way-like galaxy using the semi-analytic model a-sloth. We derive the white dwarf luminosity function by applying a linear initial-final mass relation and Mestel's cooling model. The obtained luminosity function is generally in agreement with available observations and theoretical predictions – with an exponential increase to a maximum of |$M_{\mathrm{abs}} = 16$| and a sudden drop for |$M_{\mathrm{abs}} \gt 16$|⁠. We explore the uncertainties of our model and compare them to the observational estimates. We adopt two different models of the initial mass function of Population III stars to show that the faint end of the luminosity function imprints the signature of Population III remnants. If the feature is detected in future observations, it would provide a clue to Population III stars and would also be an indirect evidence of low- to intermediate-mass Population III stars. We discuss the challenges and prospects for detecting the signatures. [ABSTRACT FROM AUTHOR]

We explain the large variety of star formation laws in terms of one single, simple law that can be inferred from the definition of the star formation rate and basic algebra. The resulting equation, |${\rm SFR}= {\epsilon _\mathrm{ff}}\ M_\mathrm{collapse}/{\tau _{\rm ff}}$|⁠ , although it has been presented elsewhere, is interpreted in terms of clouds undergoing collapse rather than being turbulence supported, an idea that different groups have pursued this century. Under such assumption, one can explain the constancy of |${\epsilon _\mathrm{ff}}$|⁠ , the different intra-cloud correlations observed in Milky Way's molecular clouds, as well as the resolved and unresolved extragalactic relationships between SFR and a measurement of the mass in CO, HCN, and CO + HI. We also explain why the slope of the correlation changes when the orbital time |${\tau _{\rm orb}}$| is considered instead of the free-fall time, and why estimations of the free-fall time from extragalactic observations skew the correlation, providing a false sublinear correlation. We furthermore show that the apparent nearly linear correlation between the star formation rate and the dynamical equilibrium pressure in the mid-plane of the galaxies, |$P_{\rm DE}$|⁠ , is just a consequence of |$P_{\rm DE}$| values being dominated by the variation of the column density of molecular gas. All in all, we argue that the star formation law is driven by the collapse of cold, dense gas, which happens to be primarily molecular in the present Universe, and that the role of stellar feedback is just to shut down the star formation process, not to shape the star formation law. [ABSTRACT FROM AUTHOR]

Meteorites trace planet formation in the Sun's protoplanetary disk, but they also record the influence of the Sun's birth environment. Whether the Sun formed in a region like Taurus-Auriga with ∼ 10 2 stars, or a region like the Carina Nebula with ∼ 10 6 stars, matters for how large the Sun's disk was, for how long and from how far away it accreted gas from the molecular cloud, and how it acquired radionuclides like 26Al. To provide context for the interpretation of meteoritic data, we review what is known about the Sun's birth environment. Based on an inferred gas disk outer radius ≈ 50 − 90 AU, radial transport in the disk, and the abundances of noble gases in Jupiter's atmosphere, the Sun's molecular cloud and protoplanetary disk were exposed to an ultraviolet flux G 0 ∼ 30 − 3000 during its birth and first ≈ 10 Myr of evolution. Based on the orbits of Kuiper Belt objects, the Solar System was subsequently exposed to a stellar density ≈ 100 M ⊙ pc − 3 for ≈ 100 Myr, strongly implying formation in a bound cluster. These facts suggest formation in a region like the outskirts of the Orion Nebula, perhaps 2 pc from the center. The protoplanetary disk might have accreted gas for many Myr, but a few × 10 5 yr seems more likely. It probably inherited radionuclides from its molecular cloud, enriched by inputs from supernovae and especially Wolf-Rayet star winds, and acquired a typical amount of 26Al. [ABSTRACT FROM AUTHOR]

We study metal enrichment originating from stellar wind and supernovae in low-metallicity clouds by performing three-dimensional radiation hydrodynamics simulations. We find that metals ejected from stellar wind are accumulated, leading to subsequent star formation in the nitrogen-enriched gas. During this early phase, the |${\rm N/O}$| ratios are similar to observed nitrogen-enriched galaxies (⁠|${\rm [N/O]}\gtrsim 0.5$|⁠). Then, once supernovae occur, the |${\rm N/O}$| ratios decrease significantly. If the duration of star formation is comparable to the time-scale of supernovae, the mass fraction of nitrogen-enriched stars reaches half the mass of star clusters. We suggest that the mass of the star cluster needs to exceed |$\sim \!10^6$|   M |$_{\odot }$| to have multiple populations due to stellar wind, considering the condition for massive star cluster formation and the timescales of stellar evolution. [ABSTRACT FROM AUTHOR]

The evolution of giant molecular clouds (GMCs), the main sites of high-mass star formation, is an essential process to unravel the galaxy evolution. Using a GMC catalogue of M 33 from the ALMA-ACA (Atacama Large Millimeter/submillimeter Array–Atacama Compact Array) survey, we classified 848 GMCs into three types based on the association with H  ii regions and their H |$\alpha$| luminosities |$L\, (\rm{H}\alpha)$|⁠ : Type I is associated with no H  ii regions; Type II with H  ii regions of |$L\, (\rm{H}\alpha)$| |$\lt 10^{37.5}$|  erg s |$^{-1}$|⁠ ; and Type III with H  ii regions of |$L\, (\rm{H}\alpha)$|   |$\geqq$|   |$10^{37.5}$|  erg s |$^{-1}$|⁠. These criteria yield 224 Type I GMCs, 473 Type II GMCs, and 151 Type III GMCs. GMCs show changes in their physical properties according to the types; mass, radius, velocity dispersion, and |$^{13}$| CO detection rate of GMCs systematically increase from Type I to Type III, and additionally, Type III GMCs are closest to virial equilibrium. Type III GMCs show the highest spatial correlation with clusters younger than |$10\:$| Myr, Type II GMCs moderate correlation, and Type I GMCs are almost uncorrelated. We interpret that these types indicate an evolutionary sequence from Type I to Type II, and then to Type III with timescales of 4 Myr, 13 Myr, and 5 Myr, respectively, indicating a GMC lifetime of 22 Myr by assuming that a Type II GMC has the same timescale as the Large Magellanic Cloud. The evolved GMCs are concentrated on the spiral arms, while the younger GMCs are apart from the arm, both to the leading and trailing sides. This indicates that GMCs collide with each other via the spiral potential, leading to the compression of GMCs and the triggering of high-mass star formation, which may support the dynamic spiral model. Overall, we suggest that the GMC evolution concept helps illuminate the galaxy evolution, including the spiral arm formation. [ABSTRACT FROM AUTHOR]

We investigated the giant molecular clouds (GMCs) in M 74 (NGC 628), using data obtained from the PHANGS (Physics at High Angular resolution in Nearby GalaxieS) project. We applied GMC types according to the activity of star formation: Type I without star formation, Type II with H |$\alpha$| luminosity (⁠|$L_{\mathrm{H\alpha }}$|⁠) less than |$10^{37.5}\ \rm{erg\ s ^{-1}}$|⁠ , and Type III with |$L_{\mathrm{H\alpha }}$| greater than |$10^{37.5}\ \rm{erg\ s^{-1}}$|⁠. A total of 432 GMCs were identified, with 59, 201, and 172 GMCs, for Types I, II, and III, respectively. The size and mass of the GMCs range from 23 to 238 pc and |$10^{4.9}$| to |$10^{7.1}\, M_{\odot }$|⁠ , indicating that the mass and radius increase from Types I to III. Clusters younger than 4 Myr and H  ii regions are concentrated within 150 pc of a GMC, indicating a tight association between these young objects and GMCs. The virial ratio decreases from Type I to Type III, indicating that Type III GMCs are the most gravitationally relaxed among the three. We interpret that the GMCs evolve from Type I to Type III, as previously observed in the Large Magellanic Cloud. Based on a steady-state assumption, the estimated evolutionary timescales of Types I, II, and III are 1, 5, and 4 Myr, respectively. We assume that the timescale of Type III is equal to the age of the associated clusters, indicating a GMC lifetime of 10 Myr or longer. Although Chevance et al. (2020 , MNRAS, 493, 2872) investigated GMCs using the same PHANGS dataset of M 74, they did not define a GMC, reaching an evolutionary picture with a 20 Myr duration of the non-star-forming phase, which is five times longer than 4 Myr. We compare the present results with those of Chevance et al. (2020 , MNRAS, 493, 2872) and argue that defining individual GMCs is essential for understanding GMC evolution. [ABSTRACT FROM AUTHOR]

W3 is one of the most outstanding regions of high-mass star formation in the outer solar circle, and includes two active star-forming clouds: W3 Main and W3(OH). Based on a new analysis of the |${^{12}\text{CO}(J = 2-1)}$| data obtained at |$38^{\prime \prime }$| resolution, we have found three clouds that have molecular masses from 2000 to |$8000\, {M_\odot }$| at velocities |$-50\:\rm{km\: s^{-1}}$|⁠ , |$-43\:\rm{km\:s^{-1}}$|⁠ , and |$-39\:\rm{km\:s^{-1}}$|⁠. The |$-43\:\rm{km\:s^{-1}}$| cloud is the most massive one, overlapping with the |$-39\:\rm{km\:s^{-1}}$| cloud and the |$-50\:\rm{km\:s^{-1}}$| cloud toward W3 Main and W3(OH), respectively. In W3 Main and W3(OH), we have found typical signatures of a cloud–cloud collision, i.e. the complementary distribution with/without a displacement between the two clouds and/or a V-shape in the position–velocity diagram. We frame a hypothesis that a cloud–cloud collision triggered the high-mass star formation in each region. The collision in W3 Main involves the |$-39\:\rm{km\:s^{-1}}$| cloud and the |$-43\:\rm{km\:s^{-1}}$| cloud. The collision likely produced a cavity in the |$-43\:\rm{km\:s^{-1}}$| cloud that has a size similar to the |$-39\:\rm{km\:s^{-1}}$| cloud and triggered the formation of young high-mass stars in IC 1795 |$2\:$| Myr ago. We suggest that the |$-39\:\rm{km\:s^{-1}}$| cloud is still triggering the high-mass objects younger than |$1\:$| Myr currently embedded in W3 Main. On the other hand, another collision between the |$-50\:\rm{km\:s^{-1}}$| cloud and the |$-43\:\rm{km\:s^{-1}}$| cloud likely formed the heavily embedded objects in W3(OH) within |$\sim\! 0.5\:$| Myr ago. The present results favour an idea that cloud–cloud collisions are common phenomena not only in the inner solar circle but also in the outer solar circle, where the number of reported cloud–cloud collisions is yet limited (Fukui et al. 2021, PASJ, 73, S1). [ABSTRACT FROM AUTHOR]

Molecular deuteration is a powerful diagnostic tool for probing the physical conditions and chemical processes in astrophysical environments. In this work, we focus on formaldehyde deuteration in the protobinary system NGC 1333 IRAS 4A, located in the Perseus molecular cloud. Using high-resolution (⁠|$\sim$| 100 au) ALMA (The Atacama Large Millimeter/submillimeter Array) observations, we investigate the [D |$_2$| CO]/[HDCO] ratio along the cavity walls of the outflows emanating from IRAS 4A1. Our analysis reveals a consistent decrease in the deuteration ratio (from |$\sim$| 60-20 per cent to |$\sim$| 10 per cent) with increasing distance from the protostar (from |$\sim$| 2000 to |$\sim$| 4000 au). Given the large measured [D |$_2$| CO]/[HDCO], both HDCO and D |$_2$| CO are likely injected by the shocks along the cavity walls into the gas-phase from the dust mantles, formed in the previous prestellar phase. We propose that the observed [D |$_2$| CO]/[HDCO] decrease is due to the density profile of the prestellar core from which NGC 1333 IRAS 4A was born. When considering the chemical processes at the base of formaldehyde deuteration, the IRAS 4A's prestellar precursor had a predominantly flat density profile within 3000 au and a decrease of density beyond this radius. [ABSTRACT FROM AUTHOR]

We use |$1-4$|  arcsec (⁠|$300-1200$| au) resolution |$^{12}$| CO (2 − 1) data from the MASSES (Mass Assembly of Stellar Systems and their Evolution with the Submillimeter Array) project to measure the projected opening angles of 46 protostellar outflows in the Perseus Molecular Cloud, 37 of which are measured with sufficiently high confidence to use in further analysis. We find that there is a statistically significant difference in the distributions of outflow opening angles for Classes 0 and I outflows, with a distinct lack of both wide-angle Class 0 outflows and highly collimated Class I outflows. Synthesizing our results with several previous studies, we find that outflows widen with age through the Class 0 stage but do not continue to widen in the Class I stage. The maximum projected opening angle reached is approximately 90 |$^{\circ }$| |$\pm$| 20 |$^{\circ }$|⁠ , with the transition between widening and remaining constant occurring near the boundary between the Classes 0 and I phases of evolution. While the volume fractions occupied by these outflows are no more than a few tens of per cent of the total core volume, at most, recent theoretical work suggests outflows may still be capable of playing a central role in setting the low star formation efficiencies of 25  per cent–50  per cent observed on core scales. [ABSTRACT FROM AUTHOR]

Cold (⁠|$\sim$| 10 K) and dense (⁠|$\sim 10^{5}$| cm |$^{-3}$|⁠) cores of gas and dust within molecular clouds, known as starless and dynamically evolved pre-stellar cores, are the birthplaces of low-mass (M |$\le$| few M |$_\odot$|⁠) stars. As detections of interstellar complex organic molecules, or COMs, in starless cores has increased, abundance comparisons suggest that some COMs might be seeded early in the star formation process and inherited to later stages (i.e. protostellar discs and eventually comets). To date observations of COMs in starless cores have been limited, with most detections reported solely in the Taurus molecular cloud. It is therefore still a question whether different environments affect abundances. We have surveyed 35 starless and pre-stellar cores in the Perseus molecular cloud with the Arizona Radio Observatory (ARO) 12 m telescope detecting both methanol, CH |$_3$| OH, and acetaldehyde, CH |$_3$| CHO, in 100 per cent and 49 per cent of the sample, respectively. In the sub-sample of 15 cores where CH |$_3$| CHO was detected at |$\gt 3\sigma$| (⁠|$\sim$| 18 mK) with the ARO 12 m, follow-up observations with the Yebes 40 m telescope were carried out. Detections of formic acid, t -HCOOH, ketene, H |$_2$| CCO, methyl cyanide, CH |$_3$| CN, vinyl cyanide, CH |$_2$| CHCN, methyl formate, HCOOCH |$_3$|⁠ , and dimethyl ether, CH |$_3$| OCH |$_3$|⁠ , are seen in at least 20 per cent of the cores. We discuss detection statistics, calculate column densities, and compare abundances across various stages of low-mass star formation. Our findings have more than doubled COM detection statistics in cold cores and show COMs are prevalent in the gas before star and planet formation in the Perseus molecular cloud. [ABSTRACT FROM AUTHOR]

Fuelling star formation in large, discy galaxies requires a continuous supply of gas accreting into star-forming regions. Previously, we characterized this accretion in four Milky Way mass galaxies (⁠|$M_{\rm halo}\sim 10^{12}{\rm M}_{\odot }$|⁠) in the FIRE-2 cosmological zoom-in simulations. At |$z\sim 0$|⁠ , we found that gas within the inner circumgalactic medium (iCGM) approaches the disc with comparable angular momentum (AM) to the disc edge, joining in the outer half of the gaseous disc. Within the disc, gas moves inwards at velocities of |$\sim$| 1–5 km s |$^{-1}$| while fully rotationally supported. In this study, we analyse the torques that drive these flows. In all cases studied, we find that the torques in discs enable gas accreted near the disc edge to transport inwards and fuel star formation in the central few kpc. The primary sources of torque come from gravity, hydrodynamical forces, and the sub-grid |$P \mathrm{ d}V$| work done by supernova (SN) remnants interacting with gas on |$\lesssim$| 10 pc scales. These SNe remnant interactions induce negative torques within the inner disc and positive torques in the outer disc. The gas–gas gravitational, hydro, and 'feedback' torques transfer AM outwards to where accreting gas joins the disc, playing an important role in driving inflows and regulating disc structure. Gravitational torques from stars and dark matter provide an AM sink within the innermost regions of the disc and iCGM, respectively. Feedback torques are dominant within the disc, while gravitational and hydrodynamical torques have similar significance depending on the system/region. Torques from viscous shearing, magnetic forces, stellar winds, and radiative transfer are less significant. [ABSTRACT FROM AUTHOR]

This paper analyses spectropolarimetric observations of the classical T Tauri star (CTTS) GM Aurigae collected with SPIRou, the near-infrared spectropolarimeter at the Canada–France–Hawaii Telescope, as part of the SLS and SPICE Large Programs. We report for the first time results on the large-scale magnetic field at the surface of GM Aur using Zeeman Doppler imaging. Its large-scale magnetic field energy is almost entirely stored in an axisymmetric poloidal field, which places GM Aur close to other CTTSs with similar internal structures. A dipole of about 730 G dominates the large-scale field topology, while higher order harmonics account for less than 30 per cent of the total magnetic energy. Overall, we find that the main difference between our three reconstructed maps (corresponding to sequential epochs) comes from the evolving tilt of the magnetic dipole, likely generated by non-stationary dynamo processes operating in this largely convective star rotating with a period of about 6 d. Finally, we report a |$5.5\sigma$| detection of a signal in the activity-filtered radial velocity data of semi-amplitude |$110\pm 20$|  m s |$^{-1}$|  at a period of |$8.745\pm 0.009$|  d. If attributed to a close-in planet in the inner accretion disc of GM Aur, it would imply that this planet candidate has a minimum mass of |$1.10 \pm 0.30\, M_\mathrm{Jup}$| and orbits at a distance of |$0.082 \pm 0.002$|  au. [ABSTRACT FROM AUTHOR]

Context. Star formation in the central molecular zone (CMZ) is suppressed with respect to that of the Galactic disk, and this is likely related to its high turbulent environment. This turbulence impedes the potential detection of prestellar cores. Aims. We present the temperature, density, and spatial structure of the CMZ molecular cloud G+0.693-0.027, which has been proposed to host a prestellar cluster in the Sgr B2 region. Methods. We analysed multiple HC3N rotational transitions that were observed with the IRAM 30m, APEX, Yebes 40m, and GBT radio telescopes, together with SMA+APEX spatially resolved maps. Results. The spectral shape of HC3N lines shows three distinct velocity components: a broad component with a line width of 23 km s−1 (C1), and two narrow components with line widths of 7.2 and 8.8 km s−1 (C2 and C3). This suggests that a fraction of the molecular gas in this cloud is undergoing turbulence dissipation. From an initial local thermodynamic equilibrium analysis, we found column densities of N = (6.54 ± 0.07) × 1014 cm−2, (9 ± 3) × 1014 cm−2, and (3.6 ± 0.7) × 1013 cm−2 for C1, C2, and C3, respectively. These values were used as input for a subsequent non-local thermodynamic equilibrium analysis, in which we found H2 densities of 2 × 104 cm−3, 5 × 104 cm−3, and 4 × 105 cm−3 and kinetic temperatures of 140 K, 30 K, and 80 K for C1, C2, and C3, respectively. The spatially resolved maps confirm that the colder and high-density condensations C2 and C3, which peak in the 70–85 km s−1 velocity range, have deconvolved sizes of 9″ (0.36 pc) and 7.6″ (0.3 pc), respectively, and are embedded in a more diffuse and warmer gas (C1). Conclusions. The larger-scale structure of the Sgr B2 region, consistently with previous works, shows a hole at 40–50 km s−1 that is likely due to a small cloud that shocked the Sgr B2 region and is spatially related with a massive cloud at 60–80 km s−1. We propose that the impacting small cloud sequentially triggered the formation of Sgr B2(M), (N), and (S) and the condensations in G+0.693-0.027 during its passage. The two condensations are in a post-shocked environment that has undergone internal fragmentation. Based on the analysis of their masses and the virial parameters, C2 might expand, while C3 might further fragment or collapse. [ABSTRACT FROM AUTHOR]

Context. Understanding the launching mechanism of winds and jets remains one of the fundamental challenges in astrophysics. The Protostellar Outflows at the EarliesT Stages (POETS) survey has recently mapped the 3D velocity field of the protostellar winds in a sample (37) of luminous young stellar objects (YSOs) at scales of 10–100 au via very long baseline interferometry (VLBI) observations of the 22 GHz water masers. In most of the targets, the distribution of the 3D maser velocities can be explained in terms of a magnetohydrodynamic (MHD) disk wind (DW). Aims. Our goal is to assess the launching mechanism of the protostellar wind in the YSO IRAS 21078+5211, the most promising MHD DW candidate from the POETS survey. Methods. We have performed multi-epoch Very Long Baseline Array (VLBA) observations of the 22 GHz water masers in IRAS 21078+5211 to determine the 3D velocities of the gas flowing along several wind streamlines previously identified at a linear resolution of ~1 au. Results. Near the YSO at small separations along (xl ≤ 150 au) and across (R ≤ 40 au) the jet axis, water masers trace three individual DW streamlines. By exploiting the 3D kinematic information of the masers, we determined the launch radii of these streamlines with an accuracy of ~1 au, and they lie in the range of 10–50 au. At increasingly greater distances along the jet (110 au ≤ xl ≤ 220 au), the outflowing gas speeds up while it collimates close to the jet axis. Magneto-centrifugal launching in a radially extended MHD DW appears to be the only viable process to explain the fast (up to 60 km s−1) and collimated (down to 10°) velocities of the wind in correspondence with launch radii ranging between 10 and 50 au. At larger separations from the jet axis (R ≥ 100 au), the water masers trace a slow (≤20 km s−1), radially expanding arched shock-front with kinematics inconsistent with magneto-centrifugal launching. Our resistive-magnetohydrodynamical simulations indicate that this shock-front could be driven by magnetic pressure. Conclusions. The results obtained in IRAS 21078+5211 demonstrate that VLBI observations of the 22 GHz water masers can reliably determine the launching mechanism of protostellar winds. [ABSTRACT FROM AUTHOR]

Context. Turbulence is a key component of molecular cloud structure. It is usually described by a cascade of energy down to the dissipation scale. The power spectrum for subsonic incompressible turbulence is ∝k−5/3, while for supersonic turbulence it is ∝k−2. Aims. We determine the power spectrum in an actively star-forming molecular cloud, from parsec scales down to the expected magnetohydrodynamic (MHD) wave cutoff (dissipation scale). Methods. We analyzed observations of the nearby NGC 1333 star-forming region in three different tracers to cover the different scales from ∼10 pc down to 20 mpc. The largest scales are covered with the low-density gas tracer 13CO (1–0) obtained with a single dish, the intermediate scales are covered with single-dish observations of the C18O (3–2) line, while the smallest scales are covered in H13CO+ (1–0) and HNC (1–0) with a combination of NOEMA interferometer and IRAM 30m single-dish observations. The complementarity of these observations enables us to generate a combined power spectrum covering more than two orders of magnitude in spatial scale. Results. We derive the power spectrum in an active star-forming region spanning more than 2 decades of spatial scales. The power spectrum of the intensity maps shows a single power-law behavior, with an exponent of 2.9 ± 0.1 and no evidence of dissipation. Moreover, there is evidence that the power spectrum of the ions to have more power at smaller scales than the neutrals, which is opposite to the theoretical expectations. Conclusions. We show new possibilities for studying the dissipation of energy at small scales in star-forming regions provided by interferometric observations. [ABSTRACT FROM AUTHOR]

Context. Recent millimetre and sub-millimetre observations have unveiled elongated and asymmetric structures around protostars. These structures, referred to as streamers, often exhibit coherent velocity gradients, seemingly indicating a directed gas flow towards the protostars. However, their origin and role in star formation remain uncertain. Aims. The protostellar core Per-emb-2, located in Barnard 1, has a relatively large streamer of 104 au that is more prominent in emission from carbon-chain molecules. We aim to unveil the formation mechanism of this streamer. Methods. We conducted mapping observations towards Per-emb-2 using the Nobeyama 45 m telescope. We targeted carbon-chain molecular lines such as CCS, HC3N, and HC5N at 45 GHz. Results. Using astrodendro, we identified one protostellar and four starless cores, including three new detections, on the Herschel column density map. The starless and protostellar cores are more or less gravitationally bound. We discovered strong CCS and HC3N emissions extending from the north to the south, appearing to bridge the gap between the protostellar core and the starless core to its north. This bridge spans 3 × 104 au with velocities of 6.5–7.0 km s−1. The velocity gradient of the bridge is opposite that of the streamer. Thus, the streamer is unlikely to be connected to the bridge, suggesting that the former does not have an accretion origin. Conclusions. We propose that a collision between a spherical core and the filament has shaped the density structure in this region, consequently triggering star formation within the head-tail-shaped core. In this core-filament collision scenario, the collision appears to have fragmented the filament into two structures. The streamer is a bow structure, while the bridge is a remnant of the shock-compressed filament. Thus, we conclude that the Per-emb-2 streamer does not significantly contribute to the mass accumulation towards the protostar. [ABSTRACT FROM AUTHOR]

Context. Eclipsing compact hierarchical triples (CHTs) are systems in which a tertiary star orbits an eclipsing binary (EB) in an orbit of fewer than 1000 days. In a CHT, all three stars exist in a space that is less than 5 AU in radius. A low-mass CHT is an interesting case through which we can understand the formation of multiple stars and planets at such small scales. Aims. In this study, we combine spectroscopy and photometry to estimate the orbital, stellar, and atmospheric parameters of stars in a sample of CHTs. Using the complete set of parameters, we aim to constrain the metallicity and age of the systems. Methods. We used time-series spectroscopy to obtain radial velocities (RVs) and disentangled spectra. Using RV modelling, EB light curve modelling, and spectral analysis, we estimated the metallicities and temperatures. Using isochrone fitting, we constrained the ages of the system. We then combined observations of masses, outer eccentricities (e2), orbital periods, and age estimates of the systems from the literature. We compared the distributions of e2, and the tertiary mass ratio, q3 = M3/(M1 + M2), for three different metallicity ranges and two age ranges. Results. We have estimated the masses, radii, temperatures, metallicities, and ages of 12 stars in four CHTs. The CHT CD-32 6459 shows signs of von Zeipel-Lidov-Kozai oscillations, while CD-62 1257 can evolve to form a triple common envelope. The rest of the CHTs are old and have an M-dwarf tertiary. We find that the q3 distribution for CHTs with sub-solar metallicity has a uniform distribution but the systems with solar and above-solar metallicity peak between 0.5 and 1. When dividing them according to their ages, we find the q3 of old systems to be around 0.5. The eccentricity, e2, favours a value of around 0.3 irrespective of metallicity or age. The distributions of q3 and e2 resemble the distributions of the mass ratio and eccentricity of close field binaries. [ABSTRACT FROM AUTHOR]

Context. A rich inventory of complex organic molecules (COMs) has been observed in high abundances in the gas phase toward Class 0 protostars. It has been suggested that these molecules are formed in ices and sublimate in the warm inner envelope close to the protostar. However, only the most abundant COM, methanol (CH3OH), had been firmly detected in ices before the era of the James Webb Space Telescope (JWST). Now, it is possible to detect the interstellar ices of other COMs and constrain their ice column densities quantitatively. Aims. We aim to determine the column densities of several oxygen-bearing COMs (O-COMs) in both gas and ice for two low-mass protostellar sources, NGC 1333 IRAS 2A (hereafter IRAS 2A) and B1-c, as case studies in our JWST Observations of Young proto-Stars (JOYS+) program. By comparing the column density ratios with respect to CH3OH between both phases measured in the same sources, we can probe the evolution of COMs from ice to gas in the early stages of star formation. Methods. The column densities of COMs in gas and ice were derived by fitting the spectra observed by the Atacama Large Millimeter/submillimeter Array (ALMA) and the JWST/Mid-InfraRed Instrument-Medium Resolution Spectroscopy (MIRI-MRS), respectively. The gas-phase emission lines were fit using local thermal equilibrium models, and the ice absorption bands were fit by matching the infrared spectra measured in laboratories. The column density ratios of four O-COMs (CH3CHO, C2H5OH, CH3OCH3, and CH3OCHO) with respect to CH3OH were compared between ice and gas in IRAS 2A and B1-c. Results. We were able to fit the fingerprint range of COM ices between 6.8 and 8.8 μm in the JWST/MIRI-MRS spectra of B1-c using similar components to the ones recently used for NGC 1333 IRAS 2A. We claim detection of CH4, OCN−, HCOO−, HCOOH, CH3CHO, C2H5OH, CH3OCH3, CH3OCHO, and CH3COCH3 in B1-c, and upper limits have been estimated for SO2, CH3COOH, and CH3CN. The total abundance of O-COM ices is constrained to be 15% with respect to H2O ice, 80% of which is dominated by CH3OH. The comparison of O-COM ratios with respect to CH3OH between ice and gas shows two different cases. On the one hand, the column density ratios of CH3OCHO and CH3OCH3 match well between the two phases, which may be attributed to a direct inheritance from ice to gas or strong chemical links with CH3OH. On the other hand, the ice ratios of CH3CHO and C2H5OH with respect to CH3OH are higher than the gas ratios by 1–2 orders of magnitude. This difference can be explained by gas-phase reprocessing following sublimation, or different spatial distributions of COMs in the envelope, which is an observational effect resulting from ALMA and JWST tracing different components in a protostellar system. Conclusions. The firm detection of COM ices other than CH3OH is reported in another well-studied low-mass protostar, B1-c, following the recent detection in NGC 1333 IRAS 2A. The column density ratios of four O-COMs with respect to CH3OH show both similarities and differences between gas and ice. Although the straightforward explanations would be the direct inheritance from ice to gas and the gas-phase reprocessing, respectively, other possibilities such as different spatial distributions of molecules cannot be excluded. [ABSTRACT FROM AUTHOR]

Context. Many fast-rotating stars (rotation periods of < 2 days) are found to be unresolved binaries with separations of tens of AU. This correlation between fast rotators and binarity leads to the question of whether the formation of binary stars inherently produces fast rotators. Aims. Our goal is to understand the spin evolution of protostars and whether the formation of companions plays a role in spinning up stars. Methods. We used magneto-hydrodynamic simulations to study the formation of multiple star systems from turbulent and non-turbulent protostellar cores. We tracked the angular momentum accreted by individual star and inner disc systems by using a sink (star) particle technique. We ran a resolution study to extrapolate protostellar properties. Results. We find in all simulations that the primary star can experience a spin-up event correlated with the formation of companions, namely fragmentation into binaries or higher-order systems. The primary star can spin up by up to 84% of its pre-fragmentation angular momentum and by up to 18% of its pre-fragmentation mass-specific angular momentum. The mechanism for the spin-up is gravitational disc instabilities in the circumstellar disc around the primary star, which leads to the accretion of material with high specific angular momentum. The simulations that experience the strongest disc instabilities fragment to form companions. Simulations with weaker spin-up events experience disc instabilities triggered by a companion flyby, and the disc instability in these cases typically does not produce further fragments (i.e. they remain binary systems). Conclusions. The primary star in multiple star systems can end up with a higher spin than single stars. This is because gravitational instabilities in the circumstellar disc around the primary star can trigger a spin-up event. In the strongest spin-up events, the instability is likely to cause disc fragmentation and the formation of companions. This spin-up mechanism, coupled with shorter disc lifetimes due to truncated circumstellar discs (and thus short spin-down times), may help produce fast rotators. [ABSTRACT FROM AUTHOR]

Context. The star formation rate (SFR), the number of stars formed per unit of time, is a fundamental quantity in the evolution of the Universe. Aims. While turbulence is believed to play a crucial role in setting the SFR, the exact mechanism remains unclear. Turbulence promotes star formation by compressing the gas, but also slows it down by stabilizing the gas against gravity. Most widely used analytical models rely on questionable assumptions, including: i) integrating over the density PDF, a one-point statistical description that ignores spatial correlation, ii) selecting self-gravitating gas based on a density threshold that often ignores turbulent dispersion, iii) assuming the freefall time as the timescale for estimating SFR without considering the need to rejuvenate the density PDF, iv) assuming the density probability distribution function (PDF) to be log-normal. This leads to the reliance on fudge factors for rough agreement with simulations. Even more seriously, when a more accurate density PDF is being used, the classical theory predicts a SFR that is essentially 0. Methods. Improving upon the only existing model that incorporates the spatial correlation of the density field, we present a new analytical model that, in a companion paper, is rigorously compared against a large series of numerical simulations. We calculate the time needed to rejuvenate density fluctuations of a given density and spatial scale, revealing that it is generally much longer than the freefall time, rendering the latter inappropriate for use. Results. We make specific predictions regarding the role of the Mach number, ℳ, and the driving scale of turbulence divided by the mean Jeans length. At low to moderate Mach numbers, turbulence does not reduce and may even slightly promote star formation by broadening the PDF. However, at higher Mach numbers, most density fluctuations are stabilized by turbulent dispersion, leading to a steep drop in the SFR as the Mach number increases. A fundamental parameter is the exponent of the power spectrum of the natural logarithm of the density, ln ρ, characterizing the spatial distribution of the density field. In the high Mach regime, the SFR strongly depends on it, as lower values imply a paucity of massive, gravitationally unstable clumps. Conclusions. We provide a revised analytical model to calculate the SFR of a system, considering not only the mean density and Mach number but also the spatial distribution of the gas through the power spectrum of ln ρ, as well as the injection scale of turbulence. At low Mach numbers, the model predicts a relatively high SFR nearly independent of ℳ, whereas for high Mach, the SFR is a steeply decreasing function of ℳ. [ABSTRACT FROM AUTHOR]

A consensus prevails with regard to star-disk systems accreting most of their mass and angular momentum during the collapse of a prestellar core. However, recent results have indicated that stars experience post-collapse or late infall, during which the star and its disk are refreshed with material from the protostellar environment through accretion streamers. Apart from adding mass to the star-disk system, infall potentially supplies a substantial amount of angular momentum, as the infalling material is initially not bound to the collapsing prestellar core. We investigate the orientation of infall on star-disk systems by analyzing the properties of accreting tracer particles in three-dimensional magnetohydrodynamical (3D MHD) simulations of a molecular cloud that is (4 pc)3 in volume. In contrast to the traditional picture, where the rotational axis is inherited from the collapse of a coherent pre-stellar core, the orientation of star-disk systems changes substantially throughout the accretion process, thereby extending the possibility of primordial misalignment as the source of large obliquities. In agreement with previous results that show larger contributions of late infall for increasing stellar masses, a misaligned infall is more likely to lead to a prolonged change in orientation for stars of higher final mass. On average, brown dwarfs and very low mass stars are more likely to form and accrete all of their mass as part of a multiple system, while stars with final masses above a few 0.1 M⊙ are more likely to accrete part of their mass as single stars. Finally, we find an overall trend among our sample: the post-collapse accretion phase is more anisotropic than the early collapse phase. This result is consistent with a scenario of Bondi-Hoyle-Littletlon accretion during the post-collapse phase, while the initial collapse is less anisotropic – despite the fact that material is funneled through accretion channels. [ABSTRACT FROM AUTHOR]

Context. The cosmic-ray ionization rate (ζ2) is one of the key parameters in star formation, since it regulates the chemical and dynamical evolution of molecular clouds by ionizing molecules and determining the coupling between the magnetic field and gas. Aims. However, measurements of ζ2 in dense clouds (e.g., nH ≥ 104 cm−3) are difficult and sensitive to the model assumptions. The aim is to find a convenient analytic approach that can be used in high-mass star-forming regions (HMSFRs), especially for warm gas environments such as hot molecular cores (HMCs). Methods. We propose a new analytic approach to calculate ζ2 through HCO+, N2H+, and CO measurements. By comparing our method with various astrochemical models and with observations found in the literature, we identify the parameter space for which the analytic approach is applicable. Results. Our method gives a good approximation, to within 50%, of ζ2 in dense and warm gas (e.g., nH ≥ 104 cm−3, T = 50, 100 K) for AV ≥ 4 mag and t ≥ 2 × 104 yr at Solar metallicity. The analytic approach gives better results for higher densities. However, it starts to underestimate ζ2 at low metallicity (Z = 0.1 Z⊙) when the value is too high (ζ2 ≥ 3 × 10−15 s−1). By applying our method to the OMC-2 FIR4 envelope and the L1157-B1 shock region, we find ζ2 values of (1.0 ± 0.3) × 10−14 s−1 and (2.2 ± 0.4) × 10−16 s−1, consistent with those previously reported. Conclusions. We calculate ζ2 toward a total of 82 samples in HMSFRs, finding that the average value of ζ2 toward all HMC samples (ζ2 = (7.4±5.0)×10−16 s−1) is more than an order of magnitude higher than the theoretical prediction of cosmic-ray attenuation models, favoring the scenario that locally accelerated cosmic rays in embedded protostars should be responsible for the observed high ζ2. [ABSTRACT FROM AUTHOR]

Context. Significant progress has been made toward understanding the formation of massive (M > 8 M⊙) binaries in close orbits (with periods of less than a month). Some of the observational studies leading to this progress are the detection of a very low velocity dispersion among the massive stars in the young region M17 and the measurement of a positive trend of velocity dispersion with age in Galactic clusters. The velocity dispersion observed in M17 could be explained either by the lack of binaries among the stars in this region, which implies the highly unlikely scenario of a different formation mechanism for M17 than for other Galactic regions, or by larger binary separations than typically observed, but with a binary fraction similar to other young Galactic clusters. The latter implies that, over time, the binary components migrate toward each other. This is in agreement with the finding that the radial velocity dispersion of young Galactic clusters correlates positively with their age. Aims. We aim to determine the origin of the strikingly low velocity dispersion by determining the observed and intrinsic binary fraction of massive stars in M17 through multi-epoch spectroscopy. Methods. We performed a multi-epoch spectroscopic survey consisting of three epochs separated by days and months, respectively. We complemented this survey with existing data covering timescales of years. We determined the radial velocity of each star at each epoch by fitting the stellar absorption profiles. The velocity shifts between epochs were used to determine whether a close companion is present. Results. We determined an observed binary fraction of 27% and an intrinsic binary fraction of 87%, consistent with that of other Galactic clusters. We conclude that the low velocity dispersion is due to a large separation among the young massive binaries in M17. Our result is in agreement with a migration scenario in which massive stars are born in binaries or higher-order systems at large separation and harden within the first million years of evolution. Such an inward migration may either be driven by interaction with a remnant accretion disk or with other young stellar objects present in the system, or by dynamical interactions within the cluster. Our results imply that possibly both dynamical interactions and binary evolution are key processes in the formation of gravitational wave sources. [ABSTRACT FROM AUTHOR]

Context. The inner parts of the hot discs surrounding massive young stellar objects (MYSOs) are still barely explored due to observational limitations in terms of angular resolution, scarcity of diagnostic lines, and the embedded and rare nature of these targets. Aims. We present the first K-band spectro-interferometric observations towards the MYSO G033.3891, which based on former kinematic evidence via the CO bandhead emission is known to host an accreting disc. Methods. Using the high spectral resolution mode (R∼4000) of the GRAVITY/VLTI, we spatially resolved the emission of the inner dusty disc and the crucial gaseous interface between the star and the dusty disc. Using detailed modelling on the K-band dust continuum and tracers known to be associated with the ionised and molecular gaseous interface (Brγ, CO), we report on the smallest scales of accretion and ejection. Results. The new observations in combination with our geometric and kinematic models employed to fit former high spectral resolution observations on the source (R∼30 000; CRIRES/VLTI) allowed us to constrain the size of the inner gaseous disc both spatially and kinematically via the CO overtone emission at only 2 au. Our models reveal that both Brγ and CO emissions are located well within the dust sublimation radius (5 au) as traced by the hot 2.2 µm dust continuum. Conclusions. Our paper provides the first case study where the tiniest scales of gaseous accretion around the MYSO G033.3891 are probed both kinematically and spatially via the CO bandhead emission. This analysis of G033.3891 stands as only the second instance of such an investigation within MYSOs, underscoring the gradual accumulation of knowledge regarding how massive young stars gain their mass while further solidifying the disc nature of accretion at the smallest scales of MYSOs. [ABSTRACT FROM AUTHOR]

Context. Together with gravity, turbulence, and stellar feedback, magnetic fields (B-fields) are thought to play a critical role in the evolution of molecular clouds and star formation processes. The polarization of thermal dust emission is a popular tracer of B-fields in star-forming regions. Aims. We aim to map the morphology and measure the strength of B-fields of the nearby molecular cloud, rho Ophiuchus-A (ρ Oph-A), to understand the role of B-fields in regulating star formation and in shaping the cloud. Methods. We analyzed the far-infrared (FIR) polarization of thermal dust emission observed by SOFIA/HAWC+ at 89 and 154 μm toward the densest part of ρ Oph-A, which is irradiated by the nearby B3/4 star, Oph-S1. These FIR polarimetric maps cover an area of ~4.5′ × 4.5′ (corresponding to 0″.18 × 0″.18 pc2) with an angular resolution of 7.8″ and 13.6″ respectively. Results. The ρ Oph-A cloud exhibits well-ordered B-fields with magnetic orientations that are mainly perpendicular to the ridge of the cloud toward the densest region. We obtained a map of B-field strengths in the range of 0.2–2.5 mG, using the Davis-Chandrasekhar-Fermi (DCF) method. The B-fields are strongest at the densest part of the cloud, which is associated with the starless core SM1, and then decrease toward the outskirts of the cloud. By calculating the map of the mass-to-flux ratio, Alfvén Mach number, and plasma β parameter in ρ Oph-A, we find that the cloud is predominantly magnetically sub-critical, sub-Alfvénic, which implies that the cloud is supported by strong B-fields that dominate over gravity, turbulence, and thermal gas energy. The measured B-field strengths at the two densest subsregions using other methods that account for the compressible mode are relatively lower than that measured with the DCF method. However, these results do not significantly change our conclusions on the roles of B-fields relative to gravity and turbulence on star formation. Our virial analysis suggests that the cloud is gravitationally unbound, which is consistent with the previous detection of numerous starless cores in the cloud. By comparing the magnetic pressure with the radiation pressure from the Oph-S1 star, we find that B-fields are sufficiently strong to support the cloud against radiative feedback and to regulate the shape of the cloud. [ABSTRACT FROM AUTHOR]

The stellar initial mass function (IMF) is critical to our understanding of star formation and the effects of young stars on their environment. On large scales, it enables us to use tracers such as UV or Hα emission to estimate the star formation rate of a system and interpret unresolved star clusters across the Universe. So far, there is little firm evidence of large-scale variations of the IMF, which is thus generally considered "universal". Stars form from cores, and it is now possible to estimate core masses and compare the core mass function (CMF) with the IMF, which it presumably produces. The goal of the ALMA-IMF large programme is to measure the core mass function at high linear resolution (2700 au) in 15 typical Milky Way protoclusters spanning a mass range of 2.5 × 103 to 32.7 × 103M⊙. In this work, we used two different core extraction algorithms to extract ≈680 gravitationally bound cores from these 15 protoclusters. We adopted a per core temperature using the temperature estimate from the point-process mapping Bayesian method (PPMAP). A power-law fit to the CMF of the sub-sample of cores above the 1.64 M⊙ completeness limit (330 cores) through the maximum likelihood estimate technique yields a slope of 1.97 ± 0.06, which is significantly flatter than the 2.35 Salpeter slope. Assuming a self-similar mapping between the CMF and the IMF, this result implies that these 15 high-mass protoclusters will generate atypical IMFs. This sample currently is the largest sample that was produced and analysed self-consistently, derived at matched physical resolution, with per core temperature estimates, and cores as massive as 150 M⊙. We provide both the raw source extraction catalogues and the catalogues listing the source size, temperature, mass, spectral indices, and so on in the 15 protoclusters. [ABSTRACT FROM AUTHOR]

In this study we present 18 |$-$| 24 GHz and high-angular-resolution (0.5 arcsec) radio wavelength Australia Telescope Compact Array follow-up observations towards a sample of 39 HC H  ii region candidates. These objects, taken from a sample hosting 6.7 GHz methanol masers, were chosen due to the compact and optically thick nature of their continuum emission. We have detected 27 compact radio sources and constructed their spectral energy distributions over the 5–24 GHz range to determine the young H  ii region's physical properties, i.e. diameter, electron density |${\it n}_{\mbox{e}}$|⁠ , emission measure, Lyman continuum flux |${\it N}_{\mbox{Ly}}$|⁠ , and turnover frequency |${\nu }_{\mbox{t}}$|⁠. The flux measurements are fitted for 20 objects assuming an ionization-bounded H  ii region with uniform density model. For the remaining seven objects that lack constraints spanning both their optically thick and thin regimes, we utilize relations from the literature to determine their physical properties. Comparing these determined parameters with those of known hypercompact (HC) and ultracompact (UC) H  ii regions, we have identified 13 HC H  ii regions, six intermediate objects that fall between HC H  ii and UC H  ii regions, six UC H  ii regions and one radio jet candidate which increases the known population of HC H  ii regions by |$\sim$| 50 per cent. All the young and compact H  ii regions are embedded in dusty and dense clumps and |$\sim$| 80 per cent of the HC H  ii regions identified in this work are associated with various maser species (CH |$_3$| OH, H |$_2$| O, and OH). Four of our radio sources remain optically thick at 24 GHz; we consider these to be among the youngest HC H  ii regions. [ABSTRACT FROM AUTHOR]

We present deep, sub-arcsecond (⁠|$\sim$| 2000 au) resolution ALMA 0.82-mm observations of the former high-mass prestellar core candidate G11.92–0.61 MM2, recently shown to be an |$\sim$| 500 au-separation protobinary. Our observations show that G11.92–0.61 MM2, located in the G11.92–0.61 protocluster, lies on a filamentary structure traced by 0.82-mm continuum and N |$_2$| H |$^+$| (4-3) emission. The N |$_2$| H |$^+$| (4-3) spectra are multipeaked, indicative of multiple velocity components along the line of sight. To analyse the gas kinematics, we performed pixel-by-pixel Gaussian decomposition of the N |$_2$| H |$^+$| spectra using scousepy and hierarchical clustering of the extracted velocity components using acorns. Seventy velocity- and position-coherent clusters (called 'trees') are identified in the N |$_2$| H |$^+$| -emitting gas, with the eight largest trees accounting for |$\gt $| 60 per cent of the fitted velocity components. The primary tree, with |$\sim$| 20 per cent of the fitted velocity components, displays a roughly north–south velocity gradient along the filamentary structure traced by the 0.82-mm continuum. Analysing an |$\sim$| 0.17 pc-long substructure, we interpret its velocity gradient of |$\sim$| 10.5 km s |$^{-1}$|  pc |$^{-1}$| as tracing filamentary accretion towards MM2 and estimate a mass inflow rate of |$\sim 1.8\times 10^{-4}$| to 1.2 |$\times 10^{-3}$|  M |$_\odot$| yr |$^{-1}$|⁠. Based on the recent detection of a bipolar molecular outflow associated with MM2, accretion on to the protobinary is ongoing, likely fed by the larger scale filamentary accretion flows. If 50 per cent of the filamentary inflow reaches the protostars, each member of the protobinary would attain a mass of 8 M |$_\odot$| within |$\sim 1.6\times 10^5$| yr, comparable to the combined time-scale of the 70- μ m- and mid-infrared-weak phases derived for ATLASGAL-TOP100 massive clumps using chemical clocks. [ABSTRACT FROM AUTHOR]

The identification of complex pre-biotic molecules using millimetre and submillimetre telescopes allows us to understand how the basic building blocks of life are formed in the universe. In the interstellar medium, ethylene glycol (⁠|$({\rm CH}_2{\rm OH})_{2}$|⁠) is the simplest sugar alcohol molecule, and it is the reduced alcohol of the simplest sugar-like molecule, glycolaldehyde (⁠|${\rm CH}_{2}{\rm OHCHO}$|⁠). We present the detection of the rotational emission lines of |$aGg^{\prime }$| conformer of ethylene glycol (⁠|$({\rm CH}_2{\rm OH})_{2}$|⁠) towards the hot molecular core G358.93–0.03 MM1 using the Atacama Large Millimeter/Submillimeter Array. The estimated column density of |$aGg^{\prime }$| - |$({\rm CH}_2{\rm OH})_{2}$| towards the G358.93–0.03 MM1 is (4.5 |$\pm 0.1)\times 10^{16}$| cm |$^{-2}$| with an excitation temperature of 155 |$\pm$| 35 K. The abundance of |$aGg^{\prime }$| - |$({\rm CH}_2{\rm OH})_{2}$| with respect to |${\rm H}_{2}$| is (1.4 |$\pm 0.5)\times 10^{-8}$|⁠. Similarly, the abundances of |$aGg^{\prime }$| - |$({\rm CH}_2{\rm OH})_{2}$| with respect to |${\rm CH}_{2}{\rm OHCHO}$| and |${\rm CH}_{3}{\rm OH}$| are 3.1 |$\pm$| 0.5 and (6.1 |$\pm 0.3)\times 10^{-3}$|⁠. We compare the estimated abundance of |$aGg^{\prime }$| - |$({\rm CH}_2{\rm OH})_{2}$| with the existing three-phase warm-up chemical model abundance of |$({\rm CH}_2{\rm OH})_{2}$|⁠ , and we notice the observed abundance and modelled abundance are nearly similar. We discuss the possible formation pathways of |$aGg^{\prime }$| - |$({\rm CH}_2{\rm OH})_{2}$| towards the hot molecular cores, and we find that |$aGg^{\prime }$| - |$({\rm CH}_2{\rm OH})_{2}$| is probably created via the recombination of two |${\rm CH}_{2}{\rm OH}$| radicals on the grain surface of G358.93–0.03 MM1. [ABSTRACT FROM AUTHOR]

NGC 6357, a star-forming complex at |$\sim 1.7$|  kpc from the Sun, contains giant molecular clouds and three prominent star clusters alongside H  ii regions, very massive stars and thousands of young stellar objects in different evolutionary stages. We present a combined infrared kinematic and time domain study of the line of sight towards this region enabled by the VVVX survey. In terms of kinematics, a novel discovery emerges an asymmetrical distribution in the vector point diagram. Some stars in the sample exhibit spatial proximity to dusty regions, with their proper motions aligned with filament projections, hinting at a younger population linked to triggered star formation. However, this distribution could also stem from an asymmetric stellar expansion event within NGC 6357, warranting further investigation. Comparing these data with Gaia revealed inconsistencies likely due to high-extinction levels in the region. Additionally, owing to accretion episodes and surface cool spots, young stars display high variability. Using the |$K_{\rm s}$| -band time series data, we overcome the extreme levels of extinction towards the region, and compile a catalogue of 774 infrared light curves of young stars. Each light curve has been characterized in terms of asymmetry and periodicity, to infer the dominant underlying physical mechanism. These findings are then correlated with evolutionary stages, aiming to uncover potential age disparities among the observed stars. This study contributes to our understanding the intricate dynamics and evolutionary processes within NGC 6357, offering valuable insights into the formation and development of stellar populations within such complex environments. [ABSTRACT FROM AUTHOR]

We present the initial data for the |$(J = 3 \rightarrow 2)$| transition of |$^{13}\text{CO}$| obtained from the central molecular zone (CMZ) of the Milky Way as part of the CO Heterodyne Inner Milky Way Plane Survey 2 (CHIMPS2). Covering |$359^{\circ } \le l \le 1^{\circ }$| and |$|b| \le 0.5^{\circ }$| with an angular resolution of 19 arcsec, velocity resolution of 1 km s |$^{-1}$|⁠ , and rms |$\Delta {T_{\rm A}^{*}} = 0.59\, \mathrm{K}$| at these resolutions, our observations unveil the complex structure of the CMZ molecular gas in improved detail. Complemented by the |$\rm {^{12}CO}$| CHIMPS2 data, we estimate a median optical depth of |$\tau _{13} = 0.087$|⁠. The preliminary analysis yields a median |$^{13}\text{CO}$| column-density range equal to |$N(^{13}{\rm CO}) = 2{\!-\!}5 \times 10^{18}\, \mathrm{cm}^{-2}$|⁠ , median H |$_{2}$| column density equal to |$N(\mathrm{H_{2}}) = 4 \times 10^{22}\, \mathrm{cm}^{-2}$| to |$1 \times 10^{23}\, \mathrm{cm}^{-2}$|⁠. We derive |$N({\rm H_{2}})$| -based total mass estimates of |$M({\rm H}_{2}) = 2{\!-\!}6 \times 10^{7}\, \mathrm{M}_{\odot }$|⁠ , in agreement with previous studies. We analyse the relationship between the integrated intensity of |$^{13}\text{CO}$| and the surface density of compact sources identified by Herschel Hi-GAL, and find that younger Hi-GAL sources detected at 500  |$\rm{\mu m}$| but not at 70  |$\rm{\mu m}$| follow the dense gas of the CMZ more closely than those that are bright at 70  |$\rm{\mu m}$|⁠. The latter, actively star-forming sources appear to be more associated with material in the foreground spiral arms. [ABSTRACT FROM AUTHOR]

We present the first large-scale 3D kinematic study of |$\sim$| 2700 spectroscopically confirmed young stars (⁠|$\lt $| 20 Myr) in 18 star clusters and OB associations (hereafter groups) from the combination of Gaia astrometry and Gaia –ESO Survey spectroscopy. We measure 3D velocity dispersions for all groups, which range from 0.61 to 7.4 km s |$^{-1}$| (1D velocity dispersions of 0.35–4.3 km s |$^{-1}$| ). We find the majority of groups have anisotropic velocity dispersions, suggesting they are not dynamically relaxed. From the 3D velocity dispersions, measured radii, and estimates of total mass, we estimate the virial state and find that all systems are super-virial when only the stellar mass is considered, but that some systems are sub-virial when the mass of the molecular cloud is taken into account. We observe an approximately linear correlation between the 3D velocity dispersion and the group mass, which would imply that the virial state of groups scales as the square root of the group mass. However, we do not observe a strong correlation between virial state and group mass. In agreement with their virial state, we find that nearly all of the groups studied are in the process of expanding and that the expansion is anisotropic, implying that groups were not spherical prior to expansion. One group, Rho Oph, is found to be contracting and in a sub-virial state (when the mass of the surrounding molecular cloud is considered). This work provides a glimpse of the potential of the combination of Gaia and data from the next generation of spectroscopic surveys. [ABSTRACT FROM AUTHOR]

Three-dimensional (3D) bubble structure of the Sgr-B molecular-cloud complex is derived by a kinematical analysis of CO-line archival cube data of the Galactic Centre (GC) observed with the Nobeyama 45-m telescope. The line-of-sight depth is estimated by applying the face-on transformation method of radial velocity to the projected distance on the Galactic plane considering the Galactic rotation of the central molecular zone (CMZ). The 3D complex exhibits a conical-horn structure with the Sgr-B2 cloud located in the farthest end on the line of sight at radial velocity |$v_{\rm lsr} \sim 70$|  km s |$^{-1}$|⁠ , and the entire complex composes a lopsided bubble opening toward the Sun at |$v_{\rm lsr}\sim 50$| to 30 km s |$^{-1}$|⁠. The line-of-sight extent of the complex is |$\sim 100$|  pc according to the large velocity extent for several tens of km s |$^{-1}$| from Sgr-B2 to the outskirts. The entire complex exhibits a flattened conical bubble with full sizes |$\sim 40 \ {\rm pc} \times 20 \ {\rm pc} \times 100 \ {\rm pc}$| in the l, b and line-of-sight directions, respectively. Based on the 3D analysis, we propose a formation scenario of the giant molecular bubble structure due to a galactic bow shock, and suggest that the star formation in Sgr-B2 was enhanced by dual-side compression (DSC) of the B2 cloud by the Galactic shock wave from up-stream and expanding H  ii region from the down-stream side of the GC Arm I in Galactic rotation. [ABSTRACT FROM AUTHOR]

We present a catalogue of dense cores and filaments in a |$3.8^\circ \times 2.4^\circ$| field around the TMC1 region of the Taurus molecular cloud. The catalogue was created using photometric data from the Herschel SPIRE and PACS instruments in the 70, 160, 250, 350, and 500 μ m continuum bands. Extended structure in the region was reconstructed from a Herschel column density map. Power spectra and probability density functions (PDFs) of this structure are presented. The PDF splits into lognormal and power-law forms, with the high-density power-law component associated primarily with the central part of TMC1. The total mass in the mapped region is 2000 M |$_\odot$|⁠ , of which 34 per cent is above an extinction of |$A_V\sim 3$| mag – a level that appears as a break in the PDF and as the minimum column density at which dense cores are found. A total of 35 dense filaments were extracted from the column density map. These have a typical full width at half-maximum (FWHM) width of 0.07 pc, but the TMC1 filament itself has a mean FWHM of |$\sim 0.13$| pc. The thermally supercritical filaments in the region are aligned orthogonal to the prevailing magnetic field direction. Derived properties for the supercritical TMC1 filament support the scenario of it being relatively young. A catalogue of 44 robust and candidate prestellar cores is created and is assessed to be complete down to 0.1 M |$_\odot$|⁠. The combined prestellar core mass function for the TMC1 and L1495 regions is well fit by a single lognormal distribution and is comparable to the standard initial mass function. [ABSTRACT FROM AUTHOR]

The majority of stars are in binary/multiple systems. How such systems form in turbulent, magnetized cores of molecular clouds in the presence of non-ideal magnetohydrodynamic (MHD) effects remains relatively underexplored. Through athena ++-based non-ideal MHD adaptive mesh refinement simulations with ambipolar diffusion, we show that the collapsing protostellar envelope is dominated by dense gravo-magneto-sheetlets, a turbulence-warped version of the classic pseudodisc produced by anisotropic magnetic resistance to the gravitational collapse, in agreement with previous simulations of turbulent, magnetized single-star formation. The sheetlets feed mass, magnetic fields, and angular momentum to a Dense ROtation-Dominated (DROD) structure, which fragments into binary/multiple systems. This DROD fragmentation scenario is a more dynamic variant of the traditional disc fragmentation scenario for binary/multiple formation, with dense spiral filaments created by inhomogeneous feeding from the highly structured larger-scale sheetlets rather than the need for angular momentum transport, which is dominated by magnetic braking. Provided that the local material is sufficiently demagnetized, with a plasma- |$\beta$| of 10 or more, collisions between the dense spiralling filaments play a key role in facilitating gravitational collapse and stellar companion formation by pushing the local magnetic Toomre parameter |$Q_\mathrm{m}$| below unity. This mechanism can naturally produce in situ misaligned systems on the 100-au scale, often detected with high-resolution Atacama Large Millimeter Array (ALMA) observations. Our simulations also highlight the importance of non-ideal MHD effects, which affect whether fragmentation occurs and, if so, the masses and orbital parameters of the stellar companions formed. [ABSTRACT FROM AUTHOR]

We study the wide-binary eccentricity (e) distribution in young star clusters and the role of turbulence in setting the form of the e distribution using magnetohydrodynamical simulations of star cluster formation. The simulations incorporate gravity, turbulence, magnetic fields, protostellar heating, and jets/outflows. We find that (1) simulations that employ purely compressive turbulence driving produce binaries with a superthermal e distribution [ |$\alpha \gt 1$| in |$p(e) \propto e^\alpha$| ], while simulations with purely solenoidal driving or natural mixture of driving modes produce subthermal/thermal distributions (⁠|$\alpha \le$| 1), (2) the e distribution over the full range of binary separations in our simulations is set at the early stages of the star cluster formation process, (3) while binaries (separation of |$r_{\mathrm{pair}} \le 1000\, \mathrm{AU}$|⁠) have subthermal to thermal e distributions (⁠|$\alpha \sim 0.8$|⁠), wide binaries (⁠|$r_{\mathrm{pair}} \gt 1000\, \mathrm{AU}$|⁠) have a superthermal distribution (⁠|$\alpha \sim 1.8$|⁠), and (4) low-mass binary systems (system masses of |$M_{\mathrm{sys}} \le 0.8\, \mathrm{M_\odot }$|⁠) have a highly superthermal distribution (⁠|$\alpha \sim 2.4$|⁠), whereas high-mass systems (⁠|$M_{\mathrm{sys}} \gt 0.8\, \mathrm{M_\odot }$|⁠) exhibit a subthermal/thermal distribution (⁠|$\alpha \sim 0.8$|⁠). The binary eccentricity distribution is often modelled as a thermal distribution. However, our results suggest that the e distribution depends on the range of separation of the sampled binaries, which agrees with the findings from recent Gaia observations. We conclude that the dependence of the e distribution on the binary separation and mass is linked to the binary formation mechanism governed by the turbulent properties of the parent cloud. [ABSTRACT FROM AUTHOR]

We present the overview and first results from the North-PHASE Legacy Survey, which follows six young clusters for five years, using the 2 deg |$^2$| FoV of the JAST80 telescope from the Javalambre Observatory (Spain). North-PHASE investigates stellar variability on time-scales from days to years for thousands of young stars distributed over entire clusters. This allows us to find new YSO, characterize accretion, and study inner disc evolution within the cluster context. Each region (Tr 37, Cep OB3, IC 5070, IC 348, NGC 2264, and NGC 1333) is observed in six filters (SDSS griz, u band, and J0660, which covers H |$\alpha$|⁠), detecting cluster members as well as field variable stars. Tr 37 is used to prove feasibility and optimize the variability analysis techniques. In Tr 37, variability reveals 50 new YSO, most of them proper motion outliers. North-PHASE independently confirms the youth of astrometric members, efficiently distinguishes accreting and non-accreting stars, reveals the extent of the cluster populations along Tr37/IC 1396 bright rims, and detects variability resulting from rotation, dips, and irregular bursts. The proper motion outliers unveil a more complex star formation history than inferred from Gaia alone, and variability highlights previously hidden proper motion deviations in the surrounding clouds. We also find that non-YSO variables identified by North-PHASE cover a different variability parameter space and include long-period variables, eclipsing binaries, RR Lyr, and |$\delta$| Scuti stars. These early results also emphasize the power of variability to complete the picture of star formation where it is missed by astrometry. [ABSTRACT FROM AUTHOR]

This paper presents a comprehensive multiwavelength investigation of the intermediate-age open cluster NGC 1912. We identified 401 member stars for NGC 1912 using K-nearest neighbour and Gaussian Mixture Model machine learning techniques. The fundamental parameters determined for the cluster are metallicity (z) = 0.0141 |$\pm$| 0.0006, log(age) = 8.519 |$\pm$| 0.040, binary fraction (⁠|$b_{\text{frac}}$|⁠) = 0.445 |$\pm$| 0.035, visual extinction (⁠|$A_V$|⁠) = 0.864 |$\pm$| 0.033 mag, total-to-selective extinction ratio (⁠|$R_V$|⁠) = 3.007 |$\pm$| 0.096, distance = 1096 |$\pm$| 15 pc. The central density (⁠|$\rho _o$|⁠), core radius (⁠|$r_c$|⁠), tidal radius (⁠|$r_t$|⁠), and background density (⁠|$\rho _{\mathrm{bg}}$|⁠) for the cluster are obtained as 0.643 |$\pm$| 0.067 stars arcmin |$^{-2}$| , 7.743 |$\pm$| 0.784 arcmin, 84.633 |$\pm$| 7.188 arcmin, and 0.001 |$\pm$| 0.000 stars arcmin |$^{-2}$| , respectively. Dynamical mass segregation study revealed a concentration of some apparently massive stars towards the cluster centre. Polarimetric analysis of 90 stars in |$BVR_cI_c$| wavelength bands revealed that the observed polarization in NGC 1912 is not solely attributed to interstellar dust. The Galactic magnetic field is dominant in this region. Dust grain sizes within the intracluster region are found to be relatively smaller than that of the general ISM. The extinction map revealed the existence of a potential 'interstellar bubble' located close to the cluster region. Towards the periphery of the probable bubble region, which displays an ovoid-like geometry, 16 probable young stellar objects have been detected, suggesting the ongoing triggering of secondary star formation events in this region. [ABSTRACT FROM AUTHOR]

First stars play crucial roles in the development of the Universe, influencing events like cosmic reionization and the chemcal enrichment of the intergalactic medium. While first stars are conventionally thought to form at around |$z \sim 20$| –30 in the standard |$\Lambda$| cold dark matter (⁠|$\Lambda$| CDM) cosmology, observational constraints on small-scale (⁠|$\\lt $| Mpc) density fluctuations remain limited, possibly differing significantly from the scale-invariant fluctuations assumed in the |$\Lambda $| CDM model. Should this be the case, the formation of first stars could occur much earlier than typically predicted. In this study, we investigate the formation process of first stars in the extremely early epochs of |$z \gtrsim 100$| in the post-recombination Universe. At such early times, the effects of the warm cosmic microwave background (CMB) become significant. We calculate the collapse of primordial star-forming clouds using a one-zone thermo-chemical model that accounts for CMB influences on radiative heating, Compton cooling, and photodissociation reactions. We found that the impact of the CMB on the evolution is limited at |$z \lesssim 100 $|⁠ , with the temperature evolution closely resembling the conventional model. However, within the range |$100 \lesssim z \lesssim 400$|⁠ , the formation of H |$_{2}$| via the H |$^{-}$| channel is impeded by H |$^{-}$| photo-detachment induced by the CMB, leading to higher temperatures compared to standard thermal evolution. Consequently, first stars with masses exceeding |$1000\, {M}_{\odot }$| can emerge at |$z \gtrsim 100$|⁠. Furthermore, at |$z \gtrsim 500$|⁠ , the temperature evolution becomes nearly isothermal at several thousand Kelvin solely due to atomic cooling, as H |$_{2}$| formation is entirely suppressed, including the less-efficient H |$_2^{+}$| channel, which is blocked by H |$_2^{+}$| photodissociation. In such cases, supermassive stars with masses around |$\sim 10^{5}\, {M}_{\odot }$| are expected to form solely via atomic cooling. These findings emphasize the significant variation in the typical mass of the first stars depending on the epoch of formation. [ABSTRACT FROM AUTHOR]

Recent cosmological hydrodynamic simulations have suggested that the first stars in the Universe often form as binary or multiple systems. However, previous studies typically overlooked the potential influence of magnetic fields during this process, assuming them to be weak and minimally impactful. Emerging theoretical investigations, however, propose an alternative perspective, suggesting that turbulent dynamo effects within first-star forming clouds can generate strong magnetic fields. In this study, we perform three-dimensional ideal magnetohydrodynamics simulations, starting from the gravitational collapse of a turbulent cloud core to the early accretion phase, where disk fragmentation frequently occurs. Our findings reveal that turbulent magnetic fields, if they reach an equipartition level with turbulence energy across all scales during the collapse phase, can significantly affect the properties of the multiple systems. Specifically, both magnetic pressure and torques contribute to disk stabilization, leading to a reduction in the number of fragments, particularly for low-mass stars. Additionally, our observations indicate the launching of protostellar jets driven by magnetic pressure of toroidal fields, although their overall impact on star formation dynamics appears to be minor. Given the case with which seed magnetic fields amplify to the full equipartition level, our results suggest that magnetic fields likely play a significant role in shaping the initial mass function of the first stars, highlighting the importance of magnetic effects on star formation in the early Universe. [ABSTRACT FROM AUTHOR]

In our previous study (Tsukamoto et al. 2023b, PASJ, 75, 835), we investigated the formation and early evolution of protoplanetary disks with 3D non-ideal magnetohydrodynamics simulations considering dust growth, and found that the modified equations of the conventional steady accretion disk model that consider magnetic braking, dust growth, and ambipolar diffusion reproduce the disk structure (such as density and vertical magnetic field) obtained from simulations very well. In this paper, as a sequel to our previous study, we analytically investigate the structure and evolution of protoplanetary disks corresponding to Class 0/I young stellar objects using the modified steady accretion disk model combining an analytical model of envelope accretion. We estimate that the disk radius is several astronomical units at the disk formation epoch and increases to several hundred astronomical units at the end of the accretion phase. The disk mass is estimated to be |$0.01 \lesssim M_{\rm disk} \lesssim 0.1 \, M_\odot$| for a disk with a radius of several tens of astronomical units and a mass accretion rate of |$\dot{M}_{\rm disk} \sim 10^{-6} \, M_\odot \,\, {\rm yr^{-1}}$|⁠. These estimates seems to be consistent with recent observations. We also found that, with typical disk ionization rates (ζ ≳ 10−19 s−1) and a moderate mass accretion rate (⁠|$\dot{M}_{\rm disk}\gtrsim 10^{-8} \, M_\odot \,\, {\rm yr^{-1}}$|⁠), magnetorotational instability is suppressed in the disk because of low plasma β and efficient ambipolar diffusion. We argue that the radial profile of specific angular momentum (or rotational velocity) at the disk outer edge should be continuously connected to that of the envelope if the disk evolves by magnetic braking, and should be discontinuous if the disk evolves by an internal angular momentum transport process such as gravitational instability or magnetorotational instability. Future detailed observations of the specific angular momentum profile around the disk outer edge are important for understanding the angular momentum transport mechanism of protoplanetary disks. [ABSTRACT FROM AUTHOR]

In order to develop a complete theory of star formation, one essentially needs to know two things: what collapses and how long it takes. This is the second paper in a series, where we query how long a parcel of gas takes to collapse and the process it undergoes. We embed pseudo-Lagrangian tracer particles in simulations of collapsing molecular clouds, identify the particles that end in dense knots, and then examine the collapse history of the gas. We find a nearly universal behaviour of cruise-then-collapse, wherein a core stays at intermediate densities for a significant fraction of its life before finally collapsing. We identify time immediately before each core collapses, |$t_{\rm {sing}}$|⁠ , and examine how it transitions to high density. We find that the time to collapse is uniformly distributed between |$0.25 t_{\rm {ff}}$| and the end of the simulation at |$\sim\!\! 1 t_{\rm {ff}}$|⁠ , and that the duration of collapse is universally short, |$\Delta t \sim 0.1 t_{\rm {ff}}$|⁠ , where |$t_{\rm {ff}}$| is the free-fall time at the mean density. We describe the collapse in three stages: collection, hardening, and singularity. Collection sweeps low-density gas into moderate density. Hardening brings kinetic and gravitational energies into quasi-equipartition. Singularity is the free-fall collapse, forming an envelope in rough energy balance and central overdensity in |$\sim\!\! 0.1 t_{\rm {ff}}$|⁠. [ABSTRACT FROM AUTHOR]

Context. Interstellar surface chemistry is a complex process that occurs in icy layers that have accumulated onto grains of different sizes. The efficiency of the surface processes often depends on the immediate environment of the adsorbed molecules. Aims. We investigated how gas-grain chemistry changes when the surface molecule binding energy is modified, depending on the properties of the surface. Methods. In a gas-grain astrochemical model, molecular binding energy gradually changes for three different environments –(1) the bare grain surface, (2) polar water-dominated ices, and (3) weakly polar carbon monoxide-dominated ices. In addition to diffusion, evaporation, and chemical desorption, photodesorption was also made binding energy-dependent, in line with experimental results. These phenomena occur in a collapsing prestellar core model that considers five grain sizes with ices arranged into four layers. Results. Variable desorption energy moderately affects gas-grain chemistry. Bare-grain effects slow down ice accumulation, while easier diffusion of molecules on weakly polar ices promotes the production of carbon dioxide. Efficient chemical desorption from bare grains significantly delays the appearance of the first ice monolayer. Conclusions. The combination of multiple aspects of grain surface chemistry creates a gas-ice balance that is different from simpler models. The composition of the interstellar ices is regulated by several binding-energy dependent desorption mechanisms. Their actions overlap in time and space, explaining the similar proportions of major ice components (water and carbon oxides) observed in all directions. [ABSTRACT FROM AUTHOR]

Episodic accretion in protostars leads to luminosity outbursts that end up heating their surroundings. This rise in temperature pushes the snow lines back, enabling the desorption of chemical species from dust grain surfaces, which may significantly alter the chemical history of the accreting envelope. However, a limited number of extensive chemical surveys of eruptive young stars have been performed thus far. In the present study, we carry out a large spectral survey of the binary Class I protostar L1551 IRS 5, known to be a FUor-like object, in the 3 mm and 2mm bands with the IRAM-30m telescope. As a result, we detected more than 400 molecular lines. The source displays a great chemical richness with the detection of 75 species, including isotopologues. Among these species, there are 13 hydrocarbons, 25 N-bearing species, 30 O-bearing species, 15 S-bearing species, 12 deuterated molecules, and a total of 10 complex organic molecules (l-C4H2, CH3CCH, CH2DCCH, CH3CHO, CH3CN, CH3OCH3, CH3OCHO, CH3OH, CH2DOH, and HC5N). With the help of local thermodynamic equilibrium (LTE) and non-LTE models, we determined the column densities of most molecules as well as excitation and kinetic temperatures. While most of those molecules trace the cold envelope (≲20 K), the OCS and CH3OH emission arise from the warm (>100 K) innermost (<2″) regions. We compared the chemical inventory of L1551 IRS 5 and its column density ratios, including isotopic ratios, with other protostellar sources. A broad chemical diversity is seen among Class I objects. More observations with both single-dish telescopes and interferometers are needed to characterize the diversity in a larger sample of protostars, while more astrochemical models would help explain this diversity, in addition to the impact of luminosity outbursts on the chemistry of protostellar envelopes. [ABSTRACT FROM AUTHOR]

Context. Signposts of early planet formation are ubiquitous in substructured young discs. Dense, hot, and high-pressure regions that formed during the gravitational collapse process, integral to star formation, facilitate dynamical mixing of dust within the protostellar disc. This provides an incentive to constrain the role of gas and dust interaction and resolve potential zones of dust concentration during star and disc formation stages. Aims. We explore whether the thermal and dynamical conditions that developed during protostellar disc formation can generate gas flows that efficiently mix and transport the well-coupled gas and dust components. Methods. We simulated the collapse of dusty molecular cloud cores with the hydrodynamics code PLUTO augmented with radiation transport and self-gravity. We used a two-dimensional axisymmetric geometry and followed the azimuthal component of the velocity. The dust was treated as Lagrangian particles that are subject to drag from the gas, whose motion is computed on a Eulerian grid. We considered 1, 10, and 100 µm-sized neutral, spherical dust grains. Importantly, the equation of state accurately includes molecular hydrogen dissociation. We focus on molecular cloud core masses of 1 and 3 M⊙ and explore the effects of different initial rotation rates and cloud core sizes. Results. Our study underlines mechanisms for the early transport of dust from the inner hot disc regions via the occurrence of two transient gas motions, namely meridional flow and outflow. The vortical flow fosters dynamical mixing and retention of dust, while the thermal pressure driven outflow replenishes dust in the outer disc. Notably, these phenomena occur regardless of the initial cloud core mass, size, and rotation rate. Conclusions. Young dynamical precursors to planet-forming discs exhibit regions with complex hydrodynamical gas features and high-temperature structures. These can play a crucial role in concentrating dust for subsequent growth into protoplanets. Dust transport, especially, from sub-au scales surrounding the protostar to the outer relatively cooler parts, offers an efficient pathway for thermal reprocessing during pre-stellar core collapse. [ABSTRACT FROM AUTHOR]

Context. The detections of narrow channels of accretion toward protostellar disks, known as streamers, have increased in number in the last few years. However, it is unclear whether streamers are a common feature around protostars that were previously missed, or if they are a rare phenomenon. Aims. Our goals are to obtain the incidence of streamers toward a region of clustered star formation and to trace the origins of their gas to determine whether they originate within the filamentary structure of molecular clouds or from beyond. Methods. We used combined observations of the nearby NGC 1333 star-forming region, carried out with the NOEMA interferometer and the IRAM 30m single dish. Our observations cover the area between the systems IRAS 4 and SVS 13. We traced the chemically fresh gas within NGC 1333 with HC3N molecular gas emission and the structure of the fibers in this region with N2H+ emission. We fit multiple velocity components in both maps and used clustering algorithms to recover velocity-coherent structures. Results. We find streamer candidates toward 7 out of 16 young stellar objects within our field of view. This represents an incidence of approximately 40% of young stellar objects with streamer candidates in a clustered star-forming region. The incidence increases to about 60% when we only considered embedded protostars. All streamers are found in HC3N emission. Conclusions. Given the different velocities between HC3N and N2H+ emission, and because by construction, N2H+ traces the fiber structure, we suggest that the gas that forms the streamers comes from outside the fibers. This implies that streamers can connect cloud material that falls onto the filaments with protostellar disk scales. [ABSTRACT FROM AUTHOR]

Context. The presence of grains is key to the synthesis of molecules in the interstellar medium that cannot form in the gas phase due to its low density and temperature conditions. In these reactions, the role of the grains is to enhance the encounter rate of the reactive species on their surfaces and to dissipate the energy excess of largely exothermic reactions, but less is known about their role as chemical catalysts; namely, bodies that provide low activation energy pathways with enhanced reaction rates. Different refractory materials with catalytic properties, such as those containing space-abundant d-block transition metals like iron (Fe), are present in astrophysical environments. Aims. Here, we report for first time mechanistic insights into the Fischer-Tropsch-type (FTT) synthesis of ethanol (CH3CH2OH), through ketene (CH2CO) and acetaldehyde (CH3CHO) intermediates, and methane (CH4) via a chain growing mechanism using a single-Fe atom supported on silica (SiO2) surfaces as a heterogeneous astrocatalyst. Methods. Quantum chemical simulations based on extended periodic surfaces were carried out to characterize the potential energy surfaces of the FTT chain growing mechanism. Calculations of the binding energies of reaction intermediates and products and Rice–Ramsperger–Kassel–Marcus kinetic calculations were performed to evaluate catalytic efficiencies and determine the feasibility of the reactions in different astrophysical environments. Results. Mechanistic studies demonstrate that the FTT chain growing mechanism enters into direct competition with FTT methanol formation, since formation of the CH2 chain growth initiator is feasible. The coupling of the CH2 with CO (forming ketene) and subsequent H2 additions yield acetaldehyde and finally ethanol, while direct H2 addition to CH2 produces methane. Thermodynamically, both processes are largely exergonic, but they present energy barriers that require external energy inputs to be overcome. Kinetic calculations demonstrate the strong temperature dependency of the FTT processes as tunneling does not dominate. Conclusions. The results could explain the presence of CH3CH2OH and CH4 in diverse astrophysical regions where current models fail to reproduce their observational quantities. The evidence that the chain growing mechanism is operating opens a new reactivity paradigm toward the formation of complex organic molecules, which is constrained by the temperature-dependent behaviour of the FTT reactions and by making their energy features a crucial aspect. [ABSTRACT FROM AUTHOR]

The kinetic temperature structure of the massive filament DR21 within the Cygnus X molecular cloud complex has been mapped using the IRAM 30 m telescope. This mapping employed the para-H2CO triplet (JKaKc = 303−202, 322−221, and 321–220) on a scale of ~0.1 pc. By modeling the averaged line ratios of para-H2CO 322–221/303–202 and 321–220/303 –202 with RADEX under non local thermodynamic equilibrium (LTE) assumptions, the kinetic temperature of the dense gas was derived, which ranges from 24 to 114 K, with an average temperature of 48.3 ± 0.5 K at a density of n(H2)= 105 cm−3. In comparison to temperature measurements using NH3 (1, 1)/(2,2) and far-infrared (FIR) wavelengths, the para-H2CO(3–2) lines reveal significantly higher temperatures. The dense clumps in various regions appear to correlate with the notable kinetic temperature (Tkin ≳ 50 K) of the dense gas traced by H2CO. Conversely, the outskirts of the DR21 filament display lower temperature distributions (Tkin < 50 K). Among the four dense cores (N44, N46, N48, and N54), temperature gradients are observed on a scale of ~0.1–0.3 pc. This suggests that the warm dense gas traced by H2CO is influenced by internal star formation activity. With the exception of the dense core N54, the temperature profiles of these cores were fitted with power-law indices ranging from −0.3 to −0.5, with a mean value of approximately −0.4. This indicates that the warm dense gas probed by H2CO is heated by radiation emitted from internally embedded protostar(s) and/or clusters. While there is no direct evidence supporting the idea that the dense gas is heated by shocks resulting from a past explosive event in the DR21 region on a scale of ~0.1 pc, our measurements of H2CO toward the DR21W1 region provide compelling evidence that the dense gas in this specific area is indeed heated by shocks originating from the western DR21 flow. Higher temperatures as traced by H2CO appear to be associated with turbulence on a scale of ~0.1 pc. The physical parameters of the dense gas as determined from H2CO lines in the DR21 filament exhibit aremarkable similarity to the results obtained in OMC-1 and N113, albeit on a scale of approximately 0.1–0.4 pc. This may imply that the physical mechanisms governing the dynamics and thermodynamics of dense gas traced by H2CO in diverse star formation regions may be dominated by common underlying principles despite variations in specific environmental conditions. [ABSTRACT FROM AUTHOR]

Context. Understanding circumstellar disks is of prime importance in astrophysics; however, their birth process remains poorly constrained due to observational and numerical challenges. Recent numerical works have shown that the small-scale physics, often wrapped into a sub-grid model, play a crucial role in disk formation and evolution. This calls for a combined approach in which both the protostar and circumstellar disk are studied in concert. Aims. We aim to elucidate the small-scale physics and constrain sub-grid parameters commonly chosen in the literature by resolving the star-disk interaction. Methods. We carried out a set of very high resolution 3D radiative-hydrodynamics simulations that self-consistently describe the collapse of a turbulent, dense molecular cloud core to stellar densities. We studied the birth of the protostar, the circumstellar disk, and its early evolution (< 6 yr after protostellar formation). Results. Following the second gravitational collapse, the nascent protostar quickly reaches breakup velocity and sheds its surface material, thus forming a hot (∼103 K), dense, and highly flared circumstellar disk. The protostar is embedded within the disk such that material can flow without crossing any shock fronts. The circumstellar disk mass quickly exceeds that of the protostar, and its kinematics are dominated by self-gravity. Accretion onto the disk is highly anisotropic, and accretion onto the protostar mainly occurs through material that slides on the disk surface. The polar mass flux is negligible in comparison. The radiative behavior also displays a strong anisotropy, as the polar accretion shock was shown to be supercritical, whereas its equatorial counterpart is subcritical. We also find a remarkable convergence of our results with respect to initial conditions. Conclusions. These results reveal the structure and kinematics in the smallest spatial scales relevant to protostellar and circumstellar disk evolution. They can be used to describe accretion onto regions commonly described by sub-grid models in simulations studying larger-scale physics. [ABSTRACT FROM AUTHOR]