Search

Recent high-resolution observations at millimeter (mm) and sub-mm reveal a diverse spatial distribution for sub-pc scale dense cores within star-forming regions, ranging from clustered to aligned arrangements. To address the increasing volume of observational and simulation data, we introduce "alignment parameters" as a quantitative and reproducible method to automatically assess core alignment. We first demonstrate the effectiveness of these parameters by applying them to artificial test clumps and comparing the results with labels from visual inspection. A threshold value is then proposed to differentiate between "clustered" and "aligned" categories. Subsequently, we apply these parameters to dense cores identified from a sample of ALMA 1.3 mm dust continuum images in high-mass star-forming regions. Analysis exploring correlations between alignment parameters and clump properties rules out the presence of moderate or strong correlation, indicating that clump properties do not appear to strongly influence the outcome of fragmentation. One possible explanation for this is that the fragmentation process is chaotic, meaning that small variations in initial conditions can lead to significant differences in fragmentation outcomes, thus obscuring any direct link between clump properties and core alignment/distribution., Comment: 19 pages, 16 figures, 2 tables. Accepted for publication in the Astrophysical Journal. GitHub repository for the used code: https://github.com/Anchewei/ALPara/releases/tag/v1.1

Deuteration is sensitive to environmental conditions in star-forming regions. To investigate NH$_2$D chemistry, we compared the spatial distribution of ortho-NH$_2$D $1_{11}^s-1_{01}^a$, NH$_3$(1,1) and NH$_3$(2,2) in 12 late-stage massive star-forming regions. By averaging several pixels along the spatial slices of ortho-NH$_2$D $1_{11}^s-1_{01}^a$, we obtained the deuterium fractionation of NH$_3$. In seven targets, the deuterium fractionation of NH$_3$ shows a decreasing trend with increasing rotational temperature. This trend is less clear in the remaining five sources, likely due to limited spatial resolution. However, when considering all 12 sources together, the anticorrelation between NH$_3$ deuterium fractionation and rotational temperature becomes less significant, suggesting that other physical parameters may influence the fractionation. Additionally, we found that the region of highest deuterium fractionation of NH$_3$ is offset from the NH$_3$ peak in each source, likely because the temperature is higher near the NH$_3$ peaks and NH$_2$D may be depleted from the gas phase as the molecular cloud core evolves, as well as the increased release of CO from grains into the gas phase., Comment: Accepted for publication in ApJ. 34 pages, 25 figures, 7 tables

(Abridged) JWST continues to deliver incredibly detailed infrared (IR) images of star forming regions in the Milky Way and beyond. IR emission from star-forming regions is very spectrally rich due to emission from gas-phase atoms, ions, and polycyclic aromatic hydrocarbons (PAHs). Physically interpreting IR images of these regions relies on assumptions about the underlying spectral energy distribution in the imaging bandpasses. We aim to provide empirical prescriptions linking line, PAH, and continuum intensities from JWST images, to facilitate the interpretation of JWST images in a wide variety of contexts. We use JWST PDRs4All Near-Infrared Camera (NIRCam) and Mid-Infrared Instrument (MIRI) imaging and Near-Infrared Spectrograph (NIRSpec) integral field unit (IFU) and MIRI Medium Resolution Spectrograph (MRS) spectroscopic observations of the Orion Bar photodissociation region (PDR), to directly compare and cross-calibrate imaging and IFU data at ~100 AU resolution over a region where the radiation field and ISM environment evolves from the hot ionized gas to the cold molecular gas. We measure the relative contributions of line, PAH, and continuum emission to the NIRCam and MIRI filters as functions of local physical conditions. We provide empirical prescriptions based on NIRCam and MIRI images to derive intensities of emission lines and PAH features. Within the range of the environments probed in this study, these prescriptions accurately predict Pa-alpha, Br-alpha, PAH 3.3 um and 11.2 um intensities, while those for FeII 1.644 um, H_2 1--0 S(1) 2.12 um and 1--0 S(9) 4.96 um, and PAH 7.7 um show more complicated environmental dependencies. Linear combinations of JWST NIRCam and MIRI images provide effective tracers of ionized gas, H_2, and PAH emission in PDRs. We expect these recipes to be useful for both the Galactic and extragalactic communities., Comment: 31 pages, 24 figures, submitted to A&A

Understanding the interstellar chemistry in low-metallicity environments is crucial to unveil physical and chemical processes in the past Galaxy or those in high-redshift galaxies, where the metallicity was significantly lower compared to the present-day solar neighborhood. This is also important for the understanding of the diversity of the chemical evolution in various regions of our Galaxy. Nearby low-metallicity laboratories, such as the outskirts of our Galaxy, the Large and Small Magellanic Clouds, and the other gas-rich dwarf galaxies in the Local Group, will provide important insights for this purpose. In the last decade, there has been great progress in astrochemical studies of interstellar molecules in low-metallicity star-forming regions. Do molecular abundances simply scale with the metallicity? If not, which processes govern the chemistry in the low-metallicity interstellar medium? In this proceeding, I will discuss the role of metallicity in the chemical evolution of star-forming regions based on recent observations of interstellar molecules in low-metallicity environments., Comment: 17 pages, 7 figures, 4 tables. To be published in the proceedings of the IAU Symposium 383: Astrochemistry VIII - From the First Galaxies to the Formation of Habitable Worlds (Kavli-IAU)

Massive star-forming regions exhibit a rich chemistry with complex gas distributions, especially on small scales. While surveys have yielded constraints on typical gas conditions, they often have coarse spatial resolution and limited bandwidths. Thus, to establish an interpretative framework for these efforts, detailed observations that simultaneously provide high sensitivity, spatial resolution, and large bandwidths for a subset of diverse sources are needed. Here, we present wideband (32 GHz) Submillimeter Array observations of four high-mass star-forming regions (G28.20-0.05, G20.08-0.14 N, G35.58-0.03, W33 Main) at subarcsecond resolution, where we detect and spatially-resolve 100s of lines from over 60 molecules, including many complex organic molecules (COMs). The chemical richness of our sample is consistent with an evolutionary sequence from the line-rich hot cores and HC HII regions of G28.20-0.05 and G20.08-0.14 N to the more chemically-modest UC HII regions in G35.58-0.03, followed by the molecule-poor HII region W33 Main. We detect lines across a range of excitation conditions (Eu=20 to >800 K) and from numerous isotopologues, which enables robust estimates of gas properties. We derive nearly constant COM column density ratios that agree with literature values in other low- and high-mass protostellar cores, supporting the idea that COM abundances are set during the pre-stellar phase. In all regions, we identify spatial offsets among different molecular families, due to a combination of source physical structure and chemistry. In particular, we find potential evidence of carbon grain sublimation in G28.20-0.05 and identify an elemental oxygen gradient and rich sulfur-chemistry in G35.58-0.03. Overall, these results demonstrate that the SMA's wide bandwidth is a powerful tool to untangle the complex molecular gas structures associated with massive star formation., Comment: 47 pages, 36 figures, accepted for publication in ApJS. Data and image cubes available at https://zenodo.org/records/13342640

Recent research provides compelling evidence that the decay of short-lived radioisotopes (SLRs), such as 26Al, provided the bulk of energy for heating and desiccation of volatile-rich planetesimals in the early Solar System. However, it remains unclear whether the early Solar System was highly enriched relative to other planetary systems with similar formation characteristics. While the Solar System possesses an elevated level of SLR enrichment compared to the interstellar medium, determining SLR enrichment of individual protoplanetary disks observationally has not been performed and is markedly more difficult. We use N-body simulations to estimate enrichment of SLRs in star-forming regions through two likely important SLR sources, stellar winds from massive stars and supernovae. We vary the number of stars and the radii of the star-forming regions and implement two models of stellar wind SLR propagation for the radioisotopes 26Al and 60Fe. We find that for 26Al enrichment the Solar System is at the upper end of the expected distribution, while for the more supernovae dependent isotope 60Fe we find that the Solar System is comparatively very highly enriched. Furthermore, combined with our previous research, these results suggest that the statistical role of 26Al-driven desiccation on exoplanet bulk composition may be underestimated in typical interpretations of the low-mass exoplanet census, and that 60Fe is even less influential as a source of heating than previously assumed., Comment: 15 pages, 6 figures, 1 table

Current galaxy formation simulations often approximate star-formation, necessitating models of star-forming regions to produce observables. In the first paper of the series, we introduced TODDLERS, a time-resolved model of UV-mm emission from star-forming regions implemented in the radiative transfer code SKIRT. This work uses SKIRT-TODDLERS to produce synthetic observations, demonstrating its potential through observables related to star-formation and comparing results with existing models in SKIRT. We calculate broadband and line emission maps for 30 Milky Way-like galaxies from the Auriga simulation at z=0. Analyzing FUV and IR data, we calculate kpc-resolved IR correction factors (k_IR), quantifying the ratio of FUV luminosity absorbed by dust to reprocessed IR luminosity. We use IR maps to calculate kpc-scale MIR (8 micron / 24 micron) and FIR (70 micron / 500 micron) colors. H alpha and H beta line maps are used to study the Balmer decrement and dust correction. We also verify the fidelity of our model FIR fine structure lines as SFR indicators. We find that the Auriga integrated UV-mm SEDs show higher FUV/NUV attenuation and lower 24 micron emission when using TODDLERS instead of the existing models in SKIRT, alleviating tensions with observations. The light-weighted mean k_IR increases with aperture and inclination, while its correlation with kpc-resolved sSFR is weaker than literature values. Kpc-scale MIR-FIR colors agree excellently with local observations, with anti-correlation degree varying by galaxy morphology. The Balmer decrement effectively corrects for dust, with attenuation law varying with dust amount. H alpha emission attenuation levels are comparable to high-density regions in state-of-the-art simulations. FIR fine-structure line emission-based luminosity-SFR relations align with global observations, [CII] showing best agreement., Comment: 22 pages, 20 figures. Accepted for publication in Astronomy and Astrophysics on October 21, 2024

The production of silicon monoxide (SiO) can be considered as a fingerprint of shock interaction. In this work, we use high-sensitivity observations of the SiO (2-1) and H$^{13}$CO$^{+}$ (1-0) emission to investigate the broad and narrow SiO emission toward 146 massive star-forming regions in the ATOMS survey. We detected SiO emission in 136 regions and distinguished broad and narrow components across the extension of 118 sources (including 58 UC $H_{II}$ regions) with an average angular resolution of 2.5$^{\prime}$$^{\prime}$. The derived SiO luminosity ($L_{SiO}$) across the whole sample shows that the majority of $L_{SiO}$ (above 66$\%$) can be attributed to broad SiO, indicating its association with strong outflows. The comparison of the ALMA SiO images with the filamentary skeletons identified from H$^{13}$CO$^{+}$ and in the infrared data (at 4.5, 8, and 24 $mu$m), further confirms that most SiO emission originates from outflows. However, note that for nine sources in our sample, the observed SiO emission may be generated by expanding UC $H_{II}$ regions. There is a moderate positive correlation between the bolometric luminosity ($L_{bol}$) and $L_{SiO}$ for both components (narrow and broad). The UC $H_{II}$ sources show a weaker positive correlation between $L_{bol}$ and $L_{SiO}$ and higher $L_{SiO}$ compared to the sources without UC $H_{II}$ regions. These results imply that the SiO emission from UC $H_{II}$ sources might be affected by UV-photochemistry induced by UC $H_{II}$ regions., Comment: 23 pages, 14 figures

Context: Star-forming regions, stellar associations, and open clusters are fundamental stellar systems where predictions from star-formation theories can be robustly contrasted with observations. Aims: We aim to provide the astrophysical community with a free and open-source code to infer the phase-space (i.e. positions and velocities) parameters of stellar systems with $\lesssim$1000 stars based on \textit{Gaia} astrometry and possibly observed radial velocities. Methods: We upgrade an existing Bayesian hierarchical model and extend it to model 3D (positions) and 6D (positions and velocities) stellar coordinates and system parameters with a flexible variety of statistical models, including a linear velocity field. This velocity field allows for the inference of internal kinematics, including expansion, contraction, and rotation. Results: We extensively validated our statistical models using realistic simulations that mimic the properties of the \textit{Gaia} Data Release 3. We applied \textit{Kalkayotl} to $\beta$-Pictoris, the Hyades, and Praesepe, recovering parameter values compatible with those from the literature. In particular, we found an expansion age of $19.1\pm1.0$ Myr for $\beta$-Pictoris and rotational signal of $32\!\pm\!11\,\rm{m\,s^{-1}\,pc^{-1}}$ for the Hyades and that Praesepe's rotation reported in the literature comes from its periphery. Conclusions: The robust and flexible Bayesian hierarchical model that we make publicly available here represents a step forward in the statistical modelling of stellar systems. The products it delivers, such as expansion, contraction, rotation, and velocity dispersions, can be directly contrasted with predictions from star-formation theories., Comment: 20 pages, 8 figures, 7 tables, associated code at https://github.com/olivares-j/Kalkayotl

Fragmentation and evolution for the molecular shells of the compact HII regions are less explored compared to their evolved counterparts. We map nine compact HII regions with a typical diameter of 0.4 pc that are surrounded by molecular shells traced by CCH. Several to a dozen dense gas fragments probed by H13CO+ are embedded in these molecular shells. These gas fragments, strongly affected by the HII region, have a higher surface density, mass, and turbulence than those outside the shells but within the same pc-scale natal clump. These features suggest that the shells swept up by the early HII regions can enhance the formation of massive dense structures that may host the birth of higher-mass stars. We examine the formation of fragments and find that fragmentation of the swept-up shell is unlikely to occur in these early HII regions, by comparing the expected time scale of shell fragmentation with the age of HII region. We propose that the appearance of gas fragments in these shells is probably the result of sweeping up pre-existing fragments into the molecular shell that has not yet fragmented. Taken together, this work provides a basis for understanding the interplay of star-forming sites with an intricate environment containing ionization feedback such as those observed in starburst regions., Comment: Accepted by MNRAS, 24 pages, 14 figures

Understanding the fragmentation of the gas cloud and the formation of massive stars remains one of the most challenging questions of modern astrophysical research. Either the gas fragmentation in a Jeans-like fashion, after which the fragments grow through accretion, or the fragmentation length is larger than the Jeans length from the start. Despite significant observational efforts, a consensus has not been reached. The key is to infer the initial density distribution upon which gravitational fragmentation occurs. Since cores are the products of the fragmentation process, the distances between adjacent cores serve as a scale indicator. Based on this observation, we propose a Delaunay triangulation-based approach to infer the density structure before the fragmentation and establish the link between density distribution and gas fragmentation length. We find that at low density, the fragmenting is Jeans-like, and at high densities, the core separations are larger than the prediction of the Jeans fragmentation. This super-Jeans fragmentation, which often occurs in groups, is responsible for the clustered formation of massive stars.

We present high-spatial-resolution (0.7 to 1.0 arcsec) submillimeter observations of continuum and molecular lines of CH3OCHO, CH3OCH3, and H2CCO toward 11 high-mass star-forming regions using the Atacama Large Millimetre/submillimetre Array (ALMA). A total of 19 separate cores from 9 high-mass star-forming regions are found to be line-rich, including high-, intermediate-, and low-mass line-rich cores. The three molecules are detected in these line-rich cores. We map the emission of CH3OCHO, CH3OCH3, and H2CCO in 9 high-mass star-forming regions. The spatial distribution of the three molecules is very similar and concentrated in the areas of intense continuum emission. We also calculate the rotation temperatures, column densities, and abundances of CH3OCHO, CH3OCH3, and H2CCO under the local thermodynamic equilibrium (LTE) assumption. The abundances relative to H2 and CH3OH, and line widths of the three molecules are significantly correlated. The abundances relative to H2, temperatures and line widths of the three molecules tend to be higher in cores with higher mass and outflows detected. The possible chemical links of the three molecules are discussed., Comment: 35 pages, 9 figures

Single-dish observations suggest that the abundances of organic species in star-forming regions of the outer Galaxy, characterised by sub-Solar metallicities, are comparable to those found in the local Galaxy. To understand this counter-intuitive result, and avoid misleading interpretation due to beam dilution effects at such large distances, spatially resolved molecular emission maps are needed to link correctly measured abundances and local physical properties. We observed several organic molecules with the Atacama Large Millimeter Array towards WB89-671, the source with the largest Galactocentric distance (23.4~kpc) of the project "CHEMical complexity in star-forming regions of the OUTer Galaxy" (CHEMOUT), at a resolution of 15000~au. We compared the observed molecular abundances with chemical model predictions. We detected emission of c-C3H2, C4H, CH3OH, H2CO, HCO, H13CO+, HCS+, CS, HN13C, and SO. The emission morphology is complex, extended, and different in each tracer. The most intense emission in H13CO+, H2CO and c-C3H2 arises from two millimeter continuum, infrared-bright cores. The most intense CH3OH and SO emission arises predominantly from the part of the filament with no continuum sources. The narrow linewidths across the filament indicate quiescent gas, despite the two embedded protostars. Derived molecular column densities are comparable with those in local star-forming regions, and suggest anti-correlation between hydrocarbons, ions, HCO, and H2CO on one side, and CH3OH and SO on the other. Static chemical models that best match the observed column densities favour low energetic conditions, expected at large Galactocentric radii, but carbon elemental abundances 3 times higher than that derived extrapolating the [C/H] Galactocentric gradient at 23~kpc. This would indicate a flatter [C/H] trend at large Galactocentric radii, in line with a flat abundance of organics., Comment: accepted for publication in Astronomy & Astrophysics, 17 pages, 14 figures, 6 tables, 3 Appendices

Massive stars play an important role in the Universe. Unlike low-mass stars, the formation of these objects located at great distances is still unclear. It is expected to be governed by some combination of self-gravity, turbulence, and magnetic fields. Our aim is to study of the chemical and physical conditions of dense clumps in several high-mass star-forming regions. We performed observations towards 5 high-mass star-forming regions (L1287, S187, S231, DR21(OH), NGC7538) with the IRAM 30 m telescope. We covered the 2-3 and 4 mm wavelength band and analyzed the lines of HCN, HNC, HCO$^+$, HC$_3$N, HNCO, OCS, CS, SiO, SO$_2$ and SO. Using astrodendro algorithm on the 850 $\mu$m dust emission data from the SCUBA Legacy catalogue, we identified dense gas clumps and determined their masses, H$_2$ column densities and sizes. Furthermore, the kinetic temperatures, molecular abundances and dynamical state were obtained. The Red Midcourse Space Experiment Source survey (RMS) was used to determine the clump types. We identified $\sim$20 clumps. We found no significant correlation between line width and size, but the linewidth-mass and mass-size relationships are strongly correlated. Virial analysis indicated that the clumps with HII regions and young stellar objects (YSOs) are gravitationally bound. Furthermore, it was suggested that significant magnetic fields provide additional support for clump stability. The molecular abundances decrease from YSOs to submm and HII regions., Comment: 6 pages, 2 figures, Proceedings of the conference "Modern astronomy: from the early Universe to exoplanets and black holes" (25-31 August 2024, SAO RAS, Nizhny Arkhyz, Russia)

Hub-filament systems are considered as natural sites for high-mass star formation. Kinematic analysis of the surroundings of hub-filaments is essential to better understand high-mass star formation within such systems. In this work, we present a detailed study of the massive Galactic protocluster IRAS 15394$-$5358, using continuum and molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming Regions (ATOMS) survey. The 3~mm dust continuum map reveals the fragmentation of the massive ($\rm M=843~M_{\odot}$) clump into six cores. The core C-1A is the largest (radius = 0.04~pc), the most massive ($\rm M=157~M_{\odot}$), and lies within the dense central region, along with two smaller cores ($\rm M=7~and~3~M_{\odot}$). The fragmentation process is consistent with the thermal Jeans fragmentation mechanism and virial analysis shows that all the cores have small virial parameter values ($\rm \alpha_{vir}<<2$), suggesting that the cores are gravitationally bound. The mass vs. radius relation indicates that three cores can potentially form at least a single massive star. The integrated intensity map of $\rm H^{13}CO^{+}$ shows that the massive clump is associated with a hub-filament system, where the central hub is linked with four filaments. A sharp velocity gradient is observed towards the hub, suggesting a global collapse where the filaments are actively feeding the hub. We discuss the role of global collapse and the possible driving mechanisms for the massive star formation activity in the protocluster., Comment: 23 pages, 19 figures, accepted for publication in MNRAS

The present survey represents a continuation of our study of high mass star forming regions in the lines of deuterated molecules, the first results of which were published in Trofimova et al. (2020). This paper presents the results of observations of 50 objects in the line of ortho modification of singly deuterated ammonia NH$_2$D $1_{11}^s - 1_{01}^a$ at frequency 85.9 GHz, carried out with the 20-m radio telescope of the Onsala Space Observatory (Sweden). This line is detected in 29 sources. The analysis of obtained data, as well as the fact that gas density in the investigated sources, according to independent estimates, is significantly lower than the critical density for this NH$_2$D transition, indicate non-LTE excitation of NH$_2$D. Based on non-LTE modeling, estimates of the relative content of the NH$_2$D molecule and the degree of deuterium enrichment were obtained, and the dependencies of these parameters on temperature and velocity dispersion were analyzed with and without taking into account detection limits assuming the same gas density in all sources. An anti-correlation between the NH$_2$D relative abundances and the kinetic temperature is revealed in the temperature range 15-50K. At the same time, significant decrease in the ratio of the NH$_2$D/NH$_3$ abundances with increasing temperature, predicted by the available chemical models, is not observed under the adopted assumptions. An anti-correlation was also revealed between the relative content of the main isotopologue of ammonia NH$_3$ and the velocity dispersion, while no statistically significant correlation with the kinetic temperature of sources in the same temperature range was found., Comment: 36 pages, 9 figures, accepted for publication in Astronomy Reports

We argue that the North Polar Spur (NPS) and many less prominent structures are formed by gaseous metal-rich plumes associated with star-forming regions (SFRs). The SFRs located at the tangent to the 3-5~kpc rings might be particularly relevant to NPS. A multi-temperature mixture of gaseous components and cosmic rays rises above the Galactic disk under the action of their initial momentum and buoyancy. Eventually, the plume velocity becomes equal to that of the ambient gas, which rotates with different angular speed than the stars in the disk. As a result, the plumes acquire characteristic bent shapes. An ad hoc model of plumes' trajectories shows an interesting resemblance to the morphology of structures seen in the radio continuum and X-rays., Comment: A&A (accepted, unformatted version)

The contribution of the magnetic field to the formation of high-mass stars is poorly understood. We report the high-angular resolution ($\sim0.3^{\prime\prime}$, 870 au) map of the magnetic field projected on the plane of the sky (B$_\mathrm{POS}$) towards the high-mass star forming region G333.46$-$0.16 (G333), obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.2 mm as part of the Magnetic Fields in Massive Star-forming Regions (MagMaR) survey. The B$_\mathrm{POS}$ morphology found in this region is consistent with a canonical ``hourglass'' which suggest a dynamically important field. This region is fragmented into two protostars separated by $\sim1740$ au. Interestingly, by analysing H$^{13}$CO$^{+}$ ($J=3-2$) line emission, we find no velocity gradient over the extend of the continuum which is consistent with a strong field. We model the B$_\mathrm{POS}$, obtaining a marginally supercritical mass-to-flux ratio of 1.43, suggesting an initially strongly magnetized environment. Based on the Davis-Chandrasekhar-Fermi method, the magnetic field strength towards G333 is estimated to be 5.7 mG. The absence of strong rotation and outflows towards the central region of G333 suggests strong magnetic braking, consistent with a highly magnetized environment. Our study shows that despite being a strong regulator, the magnetic energy fails to prevent the process of fragmentation, as revealed by the formation of the two protostars in the central region.

The formation of the massive stars, and in particular, the role that the magnetic fields play in their early evolutionary phase is still far from being completely understood. Here, we present Atacama Large Millimeter/Submillimeter Array (ALMA) 1.2 mm full polarized continuum, and H$^{13}$CO$^+$(3$-$2), CS(5$-$4), and HN$^{13}$C(3$-$2) line observations with a high angular resolution ($\sim$0.4$''$ or 1100 au). In the 1.2 mm continuum emission, we reveal a dusty envelope surrounding the massive protostars, IRAS16547-E and IRAS16547-W, with dimensions of $\sim$10,000 au. This envelope has a bi-conical structure likely carved by the powerful thermal radio jet present in region. The magnetic fields vectors follow very-well the bi-conical envelope. The polarization fraction is $\sim$2.0\% in this region. Some of these vectors seem to converge to IRAS 16547-E, and IRAS 16547-W, the most massive protostars. Moreover, the velocity fields revealed from the spectral lines H$^{13}$CO$^+$(3$-$2), and HN$^{13}$C(3$-$2) show velocity gradients with a good correspondence with the magnetic fields, that maybe are tracing the cavities of molecular outflows or maybe in some parts infall. We derived a magnetic field strength in some filamentary regions that goes from 2 to 6.1\,mG. We also find that the CS(5$-$4) molecular line emission reveals multiple outflow cavities or bow-shocks with different orientations, some of which seem to follow the NW-SE radio thermal jet., Comment: Accepted by the Astrophysical Journal, 13 pages

In stars, metallicity is usually traced using Fe, while in nebulae, O serves as the preferred proxy. Both elements have different nucleosynthetic origins and are not directly comparable. Additionally, in ionized nebulae, Fe is heavily depleted onto dust grains. We investigate the distribution of Fe gas abundances in a sample of 452 star-forming nebulae with \feiii~$\lambda 4658$ detections and their relationship with O and N. Additionally, we analyze the depletion of Fe onto dust grains in photoionized environments. We homogeneously determine the chemical abundances with direct determinations of electron temperature ($T_e$), considering the effect of possible internal variations of this parameter. We adopt a sample of 300 Galactic stars to interpret the nebular findings. We find a moderate linear correlation ($r=-0.59$) between Fe/O and O/H. In turn, we report a stronger correlation ($r=-0.80$) between Fe/N and N/H. We interpret the tighter correlation as evidence of Fe and N being produced on similar timescales while Fe-dust depletion scales with the Fe availability. The apparently flat distribution between Fe/N and N/H in Milky Way stars supports this interpretation. We find that when 12+log(O/H)<7.6, the nebulae seem to reach a plateau value around $\text{log(Fe/O)} \approx -1.7$. If this trend is confirmed, it would be consistent with a very small amount of Fe-dust in these systems, similar to what is observed in high-z galaxies discovered by the James Webb Space Telescope (JWST). We derive a relationship that allows us to approximate the fraction of Fe trapped into dust in ionized nebulae. If the O-dust scales in the same way, its possible contribution in low metallicity nebulae would be negligible. After analyzing the Fe/O abundances in J0811+4730 and J1631+4426, we do not see evidence of the presence of very massive stars with $M_\text{init}>300M_{\odot}$ in these systems., Comment: Accepted for publication in A&A

We used the dendrogram algorithm to decompose the surface density distributions of stars into hierarchical structures. These structures were tied to the multiscale structures of star clusters. A similar power-law for the mass-size relation of star clusters measured at different scales suggests a self-similar structure of star clusters. We used the minimum spanning tree method to measure the separations between clusters and gas clumps in each massive star-forming region. The separations between clusters, between clumps, and between clusters and clumps were comparable, which indicates that the evolution from clump to embedded cluster proceeds in isolation and locally, and does not affect the surrounding objects significantly. By comparing the mass functions of the ATLASGAL clumps and the identified embedded clusters, we confirm that a constant star formation efficiency of $\approx$ 0.33 can be a typical value for the ATLASGAL clumps., Comment: 6 pages, 8 figures, Accepted for publication in A&A

We present the variations in far-ultraviolet (FUV) and H$\alpha$ star formation rates (SFR), SFR$_{UV}$ and SFR$_{H\alpha}$, respectively, at sub-kpc scales in 11 galaxies as part of the Deciphering the Interplay between the Interstellar Medium, Stars, and the Circumgalactic medium (DIISC) survey. Using archival GALEX FUV imagery and H$\alpha$+[NII] narrowband images obtained with the Vatican Advanced Technology Telescope, we detect a total of 1335 (FUV-selected) and 1474 (H$\alpha$-selected) regions of recent high-mass star formation, respectively. We find the H$\alpha$-to-FUV SFR ratios tend to be lower primarily for FUV-selected regions, where SFR$_{H\alpha}$ generally underestimates the SFR by an average factor of 2-3, for SFR $<$ 10$^{-4}$ M$_{\odot}$ yr$^{-1}$. In contrast, the SFRs are generally observed to be consistent for H$\alpha$-selected regions. This discrepancy arises from morphological differences between the two indicators. Extended FUV morphologies and larger areas covered by FUV-only regions, along with decreasing overlap between FUV clumps and compact HII regions with $R/R_{25}$ suggest that stochastic sampling of the initial mass function may be more pronounced in the outer regions of galaxies. Our observed H$\alpha$-to-FUV SFR ratios are also consistent with stochastic star formation model predictions. However, using larger apertures that include diffuse FUV emission results in an offset of 1 dex between SFR$_{H\alpha}$ and SFR$_{UV}$. This suggests that the observed low H$\alpha$-to-FUV SFR ratios in galaxies are likely caused by diffuse FUV emission, which can contribute $\sim$60-90\% to the total FUV flux., Comment: 12 figures, 18 pages

Dense outflowing gas, traced by transitions of molecules with large dipole moment, is important for understanding mass loss and feedback of massive star formation. HCN 3-2 and HCO$^+$ 3-2 are good tracers of dense outflowing molecular gas, which are closely related to active star formation. In this study, we present on-the-fly (OTF) mapping observations of HCN 3-2 and HCO$^+$ 3-2 toward a sample of 33 massive star-forming regions using the 10-m Submillimeter Telescope (SMT). With the spatial distribution of line wings of HCO$^+$ 3-2 and HCN 3-2, outflows are detected in 25 sources, resulting in a detection rate of 76$\%$. The optically thin H$^{13}$CN and H$^{13}$CO$^+$ 3-2 lines are used to identify line wings as outflows and estimate core mass. The mass $M_{out}$, momentum $P_{out}$, kinetic energy $E_{K}$, force $F_{out}$ and mass loss rate $\dot M_{out}$ of outflow and core mass, are obtained for each source. A sublinear tight correlation is found between the mass of dense molecular outflow and core mass, with an index of $\sim$ 0.8 and a correlation coefficient of 0.88., Comment: 22 pages, 5 figures, 4 tables, accepted in AJ

Stars form within dense cores composed of both gas and dust within molecular clouds. However, despite the crucial role that dust plays in the star formation process, its dynamics is frequently overlooked, with the common assumption being a constant, spatially uniform dust-to-gas ratio and grain size spectrum. In this study, we introduce a set of radiation-dust-magnetohydrodynamic simulations of star forming molecular clouds from the {\small STARFORGE} project. These simulations expand upon the earlier radiation MHD models, which included cooling, individual star formation, and feedback. Notably, they explicitly address the dynamics of dust grains, considering radiation, drag, and Lorentz forces acting on a diverse size spectrum of live dust grains. We find that once stars exceed a certain mass threshold ($\sim 2 M_{\odot}$), their emitted radiation can evacuate dust grains from their vicinity, giving rise to a dust-suppressed zone of size $\sim 100$ AU. This removal of dust, which interacts with gas through cooling, chemistry, drag, and radiative transfer, alters the gas properties in the region. Commencing during the early accretion stages and preceding the Main-sequence phase, this process results in a mass-dependent depletion in the accreted dust-to-gas (ADG) mass ratio within both the circumstellar disc and the star. We predict massive stars ($\gtrsim 10 M_{\odot}$) would exhibit ADG ratios that are approximately one order of magnitude lower than that of their parent clouds. Consequently, stars, their discs, and circumstellar environments would display notable deviations in the abundances of elements commonly associated with dust grains, such as carbon and oxygen., Comment: 14 pages, 9 figures, submitted to ApJ

Observations of star-forming regions provide snapshots in time of the star formation process, and can be compared with simulation data to constrain the initial conditions of star formation. In order to make robust inferences, different metrics must be used to quantify the spatial and kinematic distributions of stars. In this paper, we assess the suitability of the INDICATE (INdex to Define Inherent Clustering And TEndencies) method as a diagnostic to infer the initial conditions of star-forming regions that subsequently undergo dynamical evolution. We use INDICATE to measure the degree of clustering in N-body simulations of the evolution of star-forming regions with different initial conditions. We find that the clustering of individual stars, as measured by INDICATE, becomes significantly higher in simulations with higher initial stellar densities, and is higher in subvirial star-forming regions where significant amounts of dynamical mixing has occurred. We then combine INDICATE with other methods that measure the mass segregation, relative stellar surface density ratio and the morphology (Q-parameter) of star-forming regions, and show that the diagnostic capability of INDICATE increases when combined with these other metrics., Comment: 17 pages (including figures), 12 Figures. Accepted for publication in MNRAS

The Jeans criterion is one cornerstone in our understanding of gravitational fragmentation. A critical limitation of the Jeans criterion is that the background density is assumed to be a constant, which is often not true in dynamic conditions such as star-forming regions. For example, during the formation phase of the high-density gas filaments in a molecular cloud, a density increase rate $\dot \rho$ implies a mass accumulation time of $t_{\rm acc}= \rho / \dot \rho= - \rho (\nabla \cdot (\rho \vec{v}))^{-1}$. The system is non-stationary when the mass accumulation time becomes comparable to the free-fall time $t_{\rm ff} = 1 / \sqrt{G \rho}$. We study fragmentation in non-stationary settings, and find that accretion can significantly increase in the characteristic mass of gravitational fragmentation ( $\lambda_{\rm Jeans,\; aac}= \lambda_{\rm Jeans} (1 + t_{\rm ff} / t_{\rm acc})^{1/3}$, $m_{\rm Jeans,\, acc} = m_{\rm Jeans} (1 + t_{\rm ff} / t_{\rm acc})$). In massive star-forming regions, this mechanism of transport-driven super-jeans fragmentation can contribute to the formation of massive stars by causing order-of-magnitude increases in the mass of the fragments., Comment: accepted by MNRAS

For understanding the nature of gaseous flows in star-forming regions of nearby galaxies it is usually utilized the relation between surface brightness in H$\alpha$ line and velocity dispersion of ionized gas known as ''surface brightness -- velocity dispersion diagram''. Using the three-dimensional gasdynamic simulations we consider the evolution of the synthetic diagrams for supershells driven by multiple supernova explosions in stellar cluster located in the galactic disk. We investigate the shape and structure of the diagrams depending on gaseous density and metallicity, disk scale height values. We show that there are several loci of values typical for young bubbles evolved in dense or rarefied gas at large heights above the disk midplane. We find that the structure of the diagram is depended on age of a supershell and physical conditions of a gas in the disk. We argue that the diagrams obtained for the nearby dwarf galaxies can be interpreted by only using the dynamics of bubbles driven by multiple supernova explosions in low-mass stellar clusters of different age., Comment: 12 pages, 4 figures

The variations in metallicity and spatial patterns within star-forming regions of galaxies result from diverse physical processes unfolding throughout their evolutionary history, with a particular emphasis in recent events. Analysing MaNGA and \textsc{eagle} galaxies, we discovered an additional dependence of the mass-metallicity relation (MZR) on metallicity gradients ($\nabla_{{\rm (O/H)}}$). Two regimes emerged for low and high stellar mass galaxies, distinctly separated at approximately ${\rm M_{\star}} >10^{9.75}$. Low-mass galaxies with strong positive $\nabla_{{\rm (O/H)}}$ appear less enriched than the MZR median, while those with strong negative gradients are consistently more enriched in both simulated and observed samples. Interestingly, low-mass galaxies with strong negative $\nabla_{{\rm (O/H)}}$ exhibit high star-forming activity, regardless of stellar surface density or $\nabla_{{\rm (O/H)}}$. In contrast, a discrepancy arises for massive galaxies between MaNGA and \textsc{eagle} datasets. The latter exhibit a notable anticorrelation between specific star formation rate and stellar surface density, independent of $\nabla_{{\rm (O/H)}}$, while MaNGA galaxies show this trend mainly for strong positive $\nabla_{{\rm (O/H)}}$. Further investigation indicates that galaxies with strong negative gradients tend to host smaller central black holes in observed datasets, a trend not replicated in simulations. These findings suggest disparities in metallicity recycling and mixing history between observations and simulations, particularly in massive galaxies with varying metallicity gradients. These distinctions could contribute to a more comprehensive understanding of the underlying physics.

Results of the long-term monitoring observations of the 6.7 GHz Class II methanol masers associated with the four high-mass star-forming regions by Hitachi 32-m radio telescope are presented. We detected periodic flux variability in G06.795-0.257, G10.472+0.027, G12.209-0.102, and G13.657-0.599 with the periods of 968, 1624, 1272, and 1266 d, respectively, although the detected period is tentative due to the short monitoring term relative to the estimated period. The facts that the flux variation patterns show the symmetric sine curves and that the luminosities of the central protostar and periods of maser flux variation are consistent with the expected period-luminosity (PL) relation suggest that the mechanisms of maser flux variability of G10.472+0.027 and G12.209-0.102 can be explained by protostellar pulsation instability. From the PL relation, central stars of these two sources are expected to be very high-mass protostars with a mass of 40 M and to have the mass accretion rate of 2*10-2 M yr-1. On the other hand, G06.795-0.257 and G13.657-0.599 have the intermittent variation patterns and have luminosities that are an order of magnitude smaller than those expected from PL relation, suggesting variation mechanisms of these sources originated from binary system. Since almost all the maser features vary with the same period regardless of its geometry, periodic accretion model may be appropriate mechanisms for flux variability in G06.795-0.257 and G13.657-0.599., Comment: arXiv admin note: text overlap with arXiv:2301.01432

To investigate the effects of stellar feedback on the gravitational state of giant molecular clouds (GMCs), we study $^{12}$CO and $^{13}$CO ALMA maps of nine GMCs distributed throughout the Large Magellanic Cloud (LMC), the nearest star-forming galaxy to our own. We perform noise and resolution matching on the sample, working at a common resolution of 3.5 arcseconds (0.85 pc at the LMC distance of 50 kpc), and use the \textit{SCIMES} clustering algorithm to identify discrete substructure, or "clumps." We supplement these data with three tracers of recent star formation: $8\mu$m surface brightness, continuum-subtracted H$\alpha$ flux, and interstellar radiation field energy density inferred from dust emission. The $^{12}$CO clumps identified cover a range of 3.6 dex in luminosity-based mass and 2.4 dex in average $8\mu$m surface brightness, representative of the wide range of conditions of the interstellar medium in the LMC. Our observations suggest evidence for increased turbulence in these clouds. While the turbulent linewidths are correlated with clump surface density, in agreement with previous observations, we find even better correlation with the three star formation activity tracers considered, suggesting stellar energy injection plays a significant role in the dynamical state of the clumps. The excess linewidths we measure do not appear to result from opacity broadening. $^{12}$CO clumps are found to be typically less gravitationally bound than $^{13}$CO clumps, with some evidence of the kinetic-to-gravitational potential energy ratio increasing with star-formation tracers. Further multi-line analysis may better constrain the assumptions made in these calculations., Comment: 23 pages, 12 figures, accepted for publication in the Astrophysical Journal

The application of silicon monoxide (SiO) as a shock tracer arises from its propensity to occur in the gas phase as a result of shock-induced phenomena, including outflow activity and interactions between molecular clouds and expanding HII regions or supernova remnants. We searched for indications of shocks toward 366 massive star-forming regions by observing the ground rotational transition of SiO ($v=0$, $J=1-0$) at 43 GHz with the Korean VLBI Network (KVN) 21 m telescopes to extend our understanding on the origins of SiO in star-forming regions. We detected SiO emission toward 104 regions that consist of 57 IRDCs, 21 HMPOs, and 26 UCHIIs. The determined median SiO column density, $N$(SiO), and abundance, $X$(SiO), relative to $N$(H$_2$) are $8.12\times10^{12}$ cm$^{-2}$ and $1.28\times10^{-10}$, respectively. These values are similar to those obtained toward other star-forming regions and also consistent with predicted values from shock models with low-velocity shocks ($\lesssim$10 - 15 km s$^{-1}$). While the $X$(SiO) does not exhibit any strong correlation with the evolutionary stages of their host clumps, $L_{\rm SiO}$ is highly correlated with dust clump mass, and $L_{\rm SiO}/L_{\rm bol}$ also has a strong negative correlation with $T_{\rm dust}$. This shows that colder and younger clumps have high $L_{\rm SiO}/L_{\rm bol}$ suggestive of an evolutionary trend. This trend is not due to excess emission at higher velocities, such as SiO wing features, as the colder sources with high $L_{\rm SiO}/L_{\rm bol}$ ratios lack wing features. Comparing SiO emission with H$_2$O and Class I CH$_3$OH masers, we find a significant correlation between $L_{\rm SiO}$/$L_{\rm bol}$ and $L_{\rm CH_3OH}/L_{\rm bol}$ ratios, whereas no similar correlation is seen for the H$_2$O maser emission. This suggests a similar origin for the SiO and Class I CH$_3$OH emission in these sources., Comment: 23 pages, 16 figures, accepted for A&A for publication

To understand the relation between NH$_2$D and its physical environment, we mapped ortho-NH$_2$D $1_{11}^s-1_{01}^a$ at 85.9 GHz toward 24 Galactic late-stage massive star-forming regions with Institut de Radioastronomie Millim$ \'e$trique (IRAM) 30-m telescope. Ortho-NH$_2$D $1_{11}^s-1_{01}^a$ was detected in 18 of 24 sources. Comparing with the distribution of H$^{13}$CN 1-0 as a dense gas tracer and radio recombination line H42$\alpha$, ortho-NH$_2$D $1_{11}^s-1_{01}^a$ present complex and diverse spatial distribution in these targets. 11 of the 18 targets, present a different distribution between ortho-NH$_2$D $1_{11}^s-1_{01}^a$ and H$^{13}$CN 1-0, while no significant difference between these two lines can be found in the other 7 sources, mainly due to limited spatial resolution and sensitivity. Moreover, with H42$\alpha$ tracing massive young stellar objects, ortho-NH$_2$D $1_{11}^s-1_{01}^a$ seems to show a relatively weak emission near the massive young stellar objects., Comment: 30 pages, 20 figures, 4 tables. Accepted to MNRAS

Using images collected with the WFC3 camera on board of the Hubble Space Telescope, we detect stellar clumps in continuum-subtracted $H\alpha$ and ultraviolet (F275W filter), such clumps are often embedded in larger regions (star-forming complexes) detected in the optical (F606W filter). We model the photometry of these objects using BAGPIPES to obtain their stellar population parameters. The median mass-weighted stellar ages are 27 Myr for $H\alpha$ clumps and 39 Myr for F275W clumps and star-forming complexes, the oldest stars in the complexes can be older than $\sim$300 Myr which indicates that star-formation is sustained for long periods of time. Stellar masses vary from 10$^{3.5}$ to 10$^{7.1}$ $M_\odot$, with star-forming complexes being more massive objects in the sample. Clumps and complexes found further away from the host galaxy are younger, less massive and less obscured by dust. We interpret these trends as due to the effect of ram-pressure in different phases of the interstellar medium. $H\alpha$ clumps form a well-defined sequence in the stellar mass--SFR plane with slope 0.73. Some F275W clumps and star-forming complexes follow the same sequence while others stray away from it and passively age. The difference in stellar age between a complex and its youngest embedded clump scales with the distance between the clump and the center of the complex, with the most displaced clumps being hosted by the most elongated complexes. This is consistent with a fireball-like morphology, where star-formation proceeds in a small portion of the complex while older stars are left behind producing a linear stellar population gradient. The stellar masses of star-forming complexes are consistent with the ones of globular clusters, but stellar mass surface densities are lower by 2 dex, and their properties are more consistent with the population of dwarf galaxies in clusters., Comment: Accepted for publication in A&A

The chemical inventory of protoplanetary midplanes is the basis for forming planetesimals. Among them, solid-state reactions based on CO/CO$_2$ toward molecular complexity on interstellar dust grains have been studied in theoretical and laboratory work. In this work, the erosion of C dust grains induced by X-ray irradiation of H$_2$O ice was systematically investigated for the first time. The work aims to provide a better understanding of the reaction mechanism using selectively isotope-labeled oxygen/carbon species in kinetic analysis. Ultrahigh vacuum experiments were performed to study the interstellar ice analog on sub-$\mu$m thick C dust at $\sim$13~K. H$_2$O or O$_2$ ice was deposited on the pre-synthesized amorphous C dust and exposed to soft X-ray photons (250--1250~eV). Fourier-transform infrared spectroscopy was used to monitor in situ the newly formed species as a function of the incident photon fluence. Field emission scanning electron microscopy was used to monitor the morphological changes of (non-)eroded carbon samples. The X-ray processing of the ice/dust interface leads to the formation of CO$_2$, which further dissociates and forms CO. Carbonyl groups are formed by oxygen addition to grain surfaces and are confirmed as intermediate species in the formation process. The yields of CO and CO$_2$ were found to be dependent on the thickness of the carbon layer. The astronomical relevance of the experimental findings is discussed., Comment: 12 pages, 9 figures, 1 table

We present and publicly release a new star-forming regions emission library TODDLERS (Time evolution of Observables including Dust Diagnostics and Line Emission from Regions containing young Stars) for the publicly available radiative transfer code SKIRT. The library generation involves the spherical evolution of a homogeneous gas cloud around a young stellar cluster that accounts for stellar feedback processes including stellar winds, supernovae, and radiation pressure, as well as the gravitational forces on the gas. The semi-analytical evolution model is coupled with the photoionization code Cloudy to calculate time-dependent UV-mm spectral energy distributions (SEDs) from star-forming regions of varying metallicity, star-formation efficiency, birth-cloud density, and mass. The calculated SEDs include the stellar, nebular, and dust continuum emission along with a wide range of emission lines originating from H ii, photo-dissociation, and molecular gas regimes tabulated at high resolution. The SEDs incorporated in SKIRT are generated by calculating a stellar-mass normalized luminosity, which assumes that each emission source is composed of a power-law population of star-forming clouds. When compared to the previous treatment of star-forming regions in SKIRT, TODDLERS shows a better agreement with low-redshift observational data in the IR wavelength range while offering a more comprehensive line-emission support. This paves the way for a variety of applications using simulated galaxies at low and high redshift., Comment: Accepted for publication in MNRAS on 27/09/23

Optical depths of dense molecular gas are commonly used in Galactic and extragalactic studies to constrain the dense gas mass of the clouds or galaxies. The optical depths are often obtained based on spatially unresolved data, especially in galaxies, which may affect the reliability of such measurements. We examine such effects in spatially resolved Galactic massive star-forming regions. Using the 10-m SMT telescope, we mapped HCN and H13CN 3-2, HCO+, and H13CO+ 3-2 towards 51 Galactic massive star-forming regions, 30 of which resulted in robust determination of spatially resolved optical depths. Conspicuous spatial variations of optical depths have been detected within each source. We first obtained opacities for each position and calculated an optical-thick line intensity-weighted average, then averaged all the spectra and derived a single opacity for each region. The two were found to agree extremely well, with a linear least square correlation coefficient of 0.997 for the whole sample., Comment: 41 pages, 33 figures, 5 tables, publication in MNRAS

The effect of dust attenuation on a galaxy's light depends on a number of physical properties, such as geometry and dust composition, both of which can vary across the faces of galaxies. To investigate this variation, we continue analysis on star-forming regions in 29 galaxies studied previously. We analyse these regions using Swift/UVOT and WISE images, as well as SDSS/MaNGA emission line maps to constrain the relationship between the infrared excess (IRX) and the UV spectral index (beta) for each star forming region. This relationship can be used to constrain which dust attenuation law is appropriate for the region. We find that the value of Dn(4000) for a region is correlated with both IRX and beta, and that the gas-phase metallicity is strongly correlated with the IRX. This correlation between metallicity and IRX suggests that regardless of aperture, metal rich regions have steeper attenuation curves. We also find that integrated galactic light follows nearly the same IRX-beta relationship as that found for kiloparsec-sized star forming regions. This similarity may suggest that the attenuation law followed by the galaxy is essentially the same as that followed by the regions, although the relatively large size of our star forming regions complicates this interpretation because optical opacity and attenuation curves have been observed to vary within individual galaxies., Comment: Accepted to MNRAS September 5, 2023

We present a study aiming at detecting VMS in local star-forming region from the imprint they leave on the integrated UV and optical light. We analyzed a sample of 27 star-forming regions and galaxies in the local Universe. We selected sources with a metallicity close to that of the LMC. We defined empirical criteria to distinguish sources dominated by VMS and Wolf-Rayet stars (WR), using template spectra of VMS- and WR-dominated regions. We subsequently built population synthesis models with an updated treatment of VMS. We show that the UV range alone is not sufficient to distinguish between VMS- and WR-dominated sources. The region of the WR bumps in the optical breaks the degeneracy. In particular, the morphology of the blue bump at 4640-4686 A is a key diagnostic. Beyond the prototypical R136 region we identify two galaxies showing clear signatures of VMS. In two other galaxies or regions the presence of VMS can be suspected, as already discussed in the literature. The stellar population is clearly dominated by WR stars in seven other sources. The most recent BPASS population synthesis models can neither account for the strong HeII 1640 emission, nor for the shape of the blue bump in VMS- and WR-dominated sources. Our models that include VMS more realistically reproduce the UV-optical spectra of VMS-dominated sources. We conclude that VMS are present in some local star-forming regions, but that separating them from WR-dominated populations requires optical spectroscopy with a high signal-to-noise ratio. A high equivalent width of HeII 1640 is not a sufficient condition for identifying VMS. Populations synthesis models need to take VMS into account by incorporating not only evolutionary tracks, but also dedicated spectral libraries. Finally, we stress that the treatment of WR stars needs to be improved as well., Comment: 16 pages, 10 figures + appendix. Accepted in Astronomy and Astrophysics

The destiny of complex organic molecules (COMs) in star-forming regions is interlinked with various evolutionary phases. Therefore, identifying these species in diversified environments of identical star-forming regions would help to comprehend their physical and chemical heritage. We identified multiple COMs utilizing the Large Program `Astrochemical Surveys At IRAM' (ASAI) data, dedicated to chemical surveys in Sun-like star-forming regions with the IRAM 30 m telescope. It was an unbiased survey in the millimetre regime, covering the prestellar core, protostar, outflow region, and protoplanetary disk phase. Here, we have reported some transitions of seven COMs, namely, methanol (CH3OH), acetaldehyde (CH3CHO), methyl formate (CH3OCHO), ethanol (C2H5OH), propynal (HCCCHO), dimethyl ether (CH3OCH3), and methyl cyanide (CH3CN) in some sources L1544, B1-b, IRAS4A, and SVS13A. We found a trend among these species from the derived abundances using the rotational diagram method and MCMC fit. We have found that the abundances of all of the COMs, except for HCCCHO, increase from the L1544 (prestellar core) and peaks at IRAS16293-2422 (class 0 phase). It is noticed that the abundance of these molecules correlate with the luminosity of the sources. The obtained trend is also visible from the previous interferometric observations and considering the beam dilution effect., Comment: 44 pages, 25 figures, and 12 tables. Accepted for the publication in the Astrophysical Journal

Simulations of the effects of stellar fly-bys on planetary systems in star-forming regions show a strong dependence on subtle variations in the initial spatial and kinematic substructure of the regions. For similar stellar densities, the more substructured star-forming regions disrupt up to a factor of two more planetary systems. We extend this work to look at the effects of substructure on stellar binary populations. We present $N$-body simulations of substructured, and non-substructured (smooth) star-forming regions in which we place different populations of stellar binaries. We find that for binary populations that are dominated by close ($<$100au) systems, a higher proportion are destroyed in substructured regions. However, for wider systems ($>$100au), a higher proportion are destroyed in smooth regions. The difference is likely due to the hard-soft, or fast-slow boundary for binary destruction. Hard (fast/close) binaries are more likely to be destroyed in environments with a small velocity dispersion (kinematically substructured regions), whereas soft (slow/wide) binaries are more likely to be destroyed in environments with higher velocity dispersions (non-kinematically substructured regions). Due to the vast range of stellar binary semimajor axes in star-forming regions ($10^{-2} - 10^4$au) these differences are small and hence unlikely to be observable. However, planetary systems have a much smaller initial semimajor axis range (likely $\sim$1 -- 100au for gas giants) and here the difference in the fraction of companions due to substructure could be observed if the star-forming regions that disrupt planetary systems formed with similar stellar densities., Comment: 9 pages, 5 figures, accepted for publication in MNRAS

We use spectra and maps of the $J=1-0$ and $J=2-1$ DCO$^{+}$, DCN, DNC, $\rm N_2D^+$ lines and $1_{11}-1_{01}$ ortho- and para-NH$_{2}$D lines, obtained with the IRAM-30m telescope, as well as observations of their hydrogenated isotopologues to study deuteration processes in five high-mass star-forming regions. The temperature was estimated from CH$_3$CCH lines, also observed with the IRAM-30m telescope, and from NH$_3$ lines, observed with the 100-m radio telescope in Effelsberg, as well as using the integrated intensity ratios of the $J=1-0$ H$^{13}$CN and HN$^{13}$C lines and their main isotopologues. Applying a non-LTE radiative transfer model with RADEX, the gas density and the molecular column densities were estimated. D/H ratios are $0.001-0.05$ for DCO$^{+}$, $0.001-0.02$ for DCN, $0.001-0.05$ for DNC and $0.02-0.4$ for NH$_{2}$D. The D/H ratios decrease with increasing temperature in the range of $\rm 20-40 \,K$ and slightly vary at densities $n(\rm H_2) \sim 10^4-10^6\, cm^{-3}$. The deuterium fraction of $\rm N_2H^{+}$ is $0.008-0.1$ at temperatures in the range of $\rm 20-25\, K$ and at a density of $\sim 10^5\, \rm cm^{-3}$. We also estimate relative abundances and find $ \sim 10^{-11}-10^{-9}$ for DCO$^{+}$ and DNC, $ \sim 10^{-11}-10^{-10}$ for $\rm N_2D^+$ and $ \sim 10^{-10}-10^{-8}$ for NH$_{2}$D. The relative abundances of these species decrease with increasing temperature. However, the DCN/H$_2$ ratio is almost constant ($\sim 10^{-10}$). The observational results agree with the predictions of chemical models (although in some cases there are significant differences)., Comment: 13 pages + appendix (25 pages total), 11 figures. Accepted for publication in MNRAS

The study of molecules and their chemistry in star-forming regions is fundamental to understand the physical process occurring in such regions. The HCN and HNC J=1-0 emissions were used to derive their integrated line intensities (I), to probe a relation recently appeared in the literature between the kinetic temperatures (T$_{K}$) and the isomeric (I) ratio, and to obtain the isomers abundances (X) in 55 high-mass star-forming regions. These last ones are classified, according to the evolutive stage, as infrared dark clouds, high-mass protostellar objects, hot molecular cores, and ultracompact HII regions. It is inferred that the T$_{K}$ obtained from the isomeric integrated intensity ratio (I$^{HCN/HNC}$) are underestimated, and hence we suggest that this relation cannot be employed as an universal thermometer in the interstellar medium. The isomers abundances show a behavior that can be explained from the chemistry occurring as the temperature and the UV radiation increase according to the evolutive stage. We found that the abundance ratio (X$^{HCN/HNC}$) hardly could be used as a chemical clock, and we suggest that it can be approximated by I$^{HCN/HNC}$. This work is part of an on-going study of multiple molecules that stand in the sample of analyzed regions which intends to contribute in the chemical knowledge of high-mass star formation., Comment: to appear in Bolet\'in de la Asociaci\'on Argentina de Astronom\'ia

Astrochemical modelling of the interstellar medium typically makes use of complex computational codes with parameters whose values can be varied. It is not always clear what the exact nature of the relationship is between these input parameters and the output molecular abundances. In this work, a feature importance analysis is conducted using SHapley Additive exPlanations (SHAP), an interpretable machine learning technique, to identify the most important physical parameters as well as their relationship with each output. The outputs are the abundances of species and ratios of abundances. In order to reduce the time taken for this process, a neural network emulator is trained to model each species' output abundance and this emulator is used to perform the interpretable machine learning. SHAP is then used to further explore the relationship between the physical features and the abundances for the various species and ratios we considered. \ce{H2O} and CO's gas phase abundances are found to strongly depend on the metallicity. \ce{NH3} has a strong temperature dependence, with there being two temperature regimes (< 100 K and > 100K). By analysing the chemical network, we relate this to the chemical reactions in our network and find the increased temperature results in increased efficiency of destruction pathways. We investigate the HCN/HNC ratio and show that it can be used as a cosmic thermometer, agreeing with the literature. This ratio is also found to be correlated with the metallicity. The HCN/CS ratio serves as a density tracer, but also has three separate temperature-dependence regimes, which are linked to the chemistry of the two molecules., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 20 pages, 20 figures and 5 tables

Aims: We introduce a new deep-learning approach for the reconstruction of 3D dust density and temperature distributions from multi-wavelength dust emission observations on the scale of individual star-forming cloud cores (<0.2pc). Methods: We construct a training data set by processing cloud cores from the Cloud Factory simulations with the POLARIS radiative transfer code to produce synthetic dust emission observations at 23 wavelengths between 12 and 1300 $\mu$m. We simplify the task by reconstructing the cloud structure along individual lines of sight and train a conditional invertible neural network (cINN) for this purpose. The cINN belongs to the group of normalising flow methods and is able to predict full posterior distributions for the target dust properties. We test different cINN setups, ranging from a scenario that includes all 23 wavelengths down to a more realistically limited case with observations at only seven wavelengths. We evaluate the predictive performance of these models on synthetic test data. Results: We report an excellent reconstruction performance for the 23-wavelengths cINN model, achieving median absolute relative errors of about 1.8% in $\log(n/m^{-3})$ and 1% in $\log(T_{dust}/K)$, respectively. We identify trends towards overestimation at the low end of the density range and towards underestimation at the high end of both the density and temperature values, which may be related to a bias in the training data. Limiting our coverage to a combination of only seven wavelengths, we still find a satisfactory performance with average absolute relative errors of about 3.3% and 2.5% in $\log(n/m^{-3})$ and $\log(T_{dust}/K)$. Conclusions: This proof-of-concept study shows that the cINN-based approach for 3D reconstruction of dust density and temperature is very promising and even compatible with a more realistically constrained wavelength coverage., Comment: 38 pages, 25 figures, accepted for publication in Astronomy & Astrophysics on 26.01.2024

A pair of 6.0 and 6.9 $\mu$m absorption features are frequently observed in Milky-Way (MW) molecular-clouds and YSOs; they also occur in the $z=0.886$ rest-frame of a molecule-rich spiral galaxy obscuring blazar PKS 1830-211. I calibrate $\chi^2$-fitting methods which match observations with two or three laboratory spectra. The 6.0-$\mu$m component is dominated by H$_2$O ice, as expected. Included MW sources were selected using opacity criteria which limit the range of explored H$_2$O-ice column densities to 1.6--$2.4 \times 10^{18}$ molecules cm$^{-2}$, while the H$_2$O-ice density in the galaxy absorber is $(2.7\pm 0.5)\times 10^{18}$ molecules cm$^{-2}$. CH$_3$OH ice and / or small (< 0.1-$\mu$m-sized) Ca- and Mg-bearing carbonates contribute at 6.9 $\mu$m. The 41 % CH$_3$OH : H$_2$O molecular ratio in the PKS 1830-211 absorber is significantly higher than in the molecular cloud towards Taurus-Elias 16 (<7.5 %) and similar to the highest value in MW YSOs (35 % in AFGL 989). Fitted carbonate (-CO$_3$) : H$_2$O ratios in the galaxy absorber of 0.091 % are low in comparison to most of the ratios detected in the MW sample (0.2 - 0.4 %; $\sim 0$ % in AFGL 989). Inorganic carbonates could explain the increased oxygen depletion at the diffuse-medium-to-molecular-cloud transition which Jones \& Ysard associated with unobserved organic carbonates or materials with a C:O ratio of 1:3., Comment: 14 pages, 7 figures. Accepted June 22 2023. version before copy editing

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

The $^{18}$O$/^{17}$O abundance ratio is, in principle, a powerful tool to estimate the relative contributions of massive stars and low- to intermediate-mass stars to the chemical enrichment of galaxies. We present $^{18}$O$/^{17}$O ratios derived from simultaneous observations of C$^{18}$O and C$^{17}$O 1-0 toward fifty-one massive star forming regions with the Institut de Radioastronomie Millim\'etrique (IRAM) 30 meter telescope. Simultaneous observations of HC$^{18}$O$^{+}$ 1-0 and HC$^{17}$O$^{+}$ 1-0 with the Yebes 40m telescope toward five sources from this sample were also done to test the consistency of $^{18}$O$/^{17}$O ratios derived from different isotopic pairs. From our improved measurements, resulting in smaller errors than previous work in the literature, we obtain a clear trend of increasing $^{18}$O$/^{17}$O ratio with increasing galactocentric distance (D$_{GC}$), which provides a significant constraint on Galactic chemical evolution (GCE) models. Current GCE models have to be improved in order to explain the observed C$^{18}$O/C$^{17}$O 1-0 gradient., Comment: 15 pages, 5 figures, 1 table, Accepted for publication in MNRAS

Low-energy cosmic rays ($<1$ TeV) are a pivotal source of ionisation of the interstellar medium, where they play a central role in determining the gas chemical composition and drastically influence the formation of stars and planets. Over the last few decades, H$_3^+$ absorption lines observations in diffuse clouds have provided reliable estimates of the cosmic ray ionisation rate relative to H$_2$ ($\zeta^{\rm ion}_{{\rm H}_2}$). However, in denser clouds, where stars and planets form, this method is often inefficient due to the lack of H$_3^+$ rotational transitions. The $\zeta^{\rm ion}_{{\rm H}_2}$ estimates are, therefore, still provisional in this context and represent one of the least understood components when it comes to defining general models of star and planet formation. In this Letter, we present the first high-resolution maps of the $\zeta^{\rm ion}_{{\rm H}_2}$ in two high-mass clumps obtained with a new analytical approach recently proposed to estimate the $\zeta^{\rm ion}_{{\rm H}_2}$ in the densest regions of molecular clouds. We obtain $\langle\zeta^{\rm ion}_{{\rm H}_2}\rangle$ that span from $3\times10^{-17}$ to $10^{-16}\rm~s^{-1}$, depending on the different distribution of the main ion carriers, in excellent agreement with the most recent cosmic rays propagation models. The cores belonging to the same parental clump show comparable $\zeta^{\rm ion}_{{\rm H}_2}$, suggesting that the ionisation properties of prestellar regions are determined by global rather than local effects. These results provide important information for the chemical and physical modelling of star-forming regions., Comment: 10 pages, 4 figures, 2 tables. Accepted for publication in The Astrophysical Journal Letters (ApJL)

Protostellar outflows are considered a signpost of star formation. These outflows can cause shocks in the molecular gas and are typically traced by the line wings of certain molecules. HNCO (4--3) has been regarded as a shock tracer because of the high abundance in shocked regions. Here we present the first imaging results of HNCO (4--3) line wings toward nine sources in a sample of twenty three massive star-forming regions using the IRAM 30\,m telescope. We adopt the velocity range of the full width of HC$_{3}$N (10--9) and H$^{13}$CO$^+$ (1--0) emissions as the central emission values, beyond which the emission from HNCO (4--3) is considered to be from line wings. The spatial distributions of the red- and/or blue-lobes of HNCO (4--3) emission nicely associate with those lobes of HCO$^{+}$ (1--0) in most of the sources. High intensity ratios of HNCO (4--3) to HCO$^+$ (1--0) are obtained in the line wings. The derived column density ratios of HNCO to HCO$^+$ are consistent with those previously observed towards massive star-forming regions. These results provide direct evidence that HNCO could trace outflow in massive star-forming regions. This work also implies that the formation of some HNCO molecules is related to shock, either on the grain surface or within the shocked gas., Comment: 18 pages, 4 tables, 4 figures, and accepted for publication in ApJ