Search

We present a catalogue of filamentary structures identified in the SARAO (South African Radio Astronomy Observatory) MeerKAT 1.3 GHz Galactic Plane Survey (SMGPS). We extract 933 filaments across the survey area, 803 of which (⁠|${\sim }86~{{\ \rm per\ cent}}$|⁠) are associated with extended radio structures (e.g. supernova remnants and H  ii regions), while 130 (⁠|${\sim }14~{{\ \rm per\ cent}}$|⁠) are largely isolated. We classify filaments as thermal or non-thermal via their associated mid-infrared emission and find that 77/130 (⁠|${\sim }59~{{\ \rm per\ cent}}$|⁠) of the isolated sources are likely to be non-thermal, and are therefore excellent candidates for the first isolated, non-thermal radio filaments observed outside of the Galactic Centre (GC). Comparing the morphological properties of these non-thermal candidates to the non-thermal filaments observed towards the GC, we find that the GC filaments are on the whole angularly narrower and shorter than those across the SMGPS, potentially an effect of distance. The SMGPS filaments have flux densities similar to those of the GC; however, the distribution of the latter extends to higher flux densities. If the SMGPS filaments were closer than the GC population, it would imply a more energetic population of cosmic ray electrons in the GC. We find that the filament position angles in the SMGPS are uniformly distributed, implying that the local magnetic field traced by the filaments does not follow the large-scale Galactic field. Finally, although we have clearly shown that filaments are not unique to the GC, the GC nevertheless has the highest density of filaments in the Milky Way. [ABSTRACT FROM AUTHOR]

The relationship between B-field orientation and density structure in molecular clouds is often assessed using the histogram of relative orientations (HRO). We perform a plane-of-the-sky geometrical analysis of projected B-fields, by interpreting HROs in dense, spheroidal, pre-stellar, and protostellar cores. We use James Clerk Maxwell Telescope POL-2 850 |$\mu$| m polarization maps and Herschel column density maps to study dense cores in the Ophiuchus molecular cloud complex. We construct two-dimensional core models, assuming Plummer column density profiles and modelling both linear and hourglass B-fields. We find that high-aspect ratio ellipsoidal cores produce strong HRO signals, as measured using the shape parameter |$\xi$|⁠. Cores with linear fields oriented |$<\!\! 45 ^{\circ }$| from their minor axis produce constant HROs with |$-1 \lt \xi \lt 0$|⁠ , indicating that fields are preferentially parallel to column density gradients. Fields parallel to the core minor axis produce the most negative value of |$\xi$|⁠. For low-aspect ratio cores, |$\xi \approx 0$| for linear fields. Hourglass fields produce a minimum in |$\xi$| at intermediate densities in all cases, converging to the minor-axis-parallel linear field value at high and low column densities. We create HROs for six dense cores in Ophiuchus. |$\rho$| Oph A and IRAS 16293 have high aspect ratios and preferentially negative HROs, consistent with moderately strong field behaviour. |$\rho$| Oph C, L1689A, and L1689B have low aspect ratios, and |$\xi \approx 0$|⁠. |$\rho$| Oph B is too complex to be modelled using a simple spheroidal field geometry. We see no signature of hourglass fields, agreeing with previous findings that dense cores generally exhibit linear fields on these size scales. [ABSTRACT FROM AUTHOR]

We extend our time-dependent hydrogen ionization simulations of diffuse ionized gas to include metals important for collisional cooling and diagnostic emission lines. The combination of heating from supernovae and time-dependent collisional and photoionization from mid-plane OB stars produces emission line intensities (and emission line ratios) that follow the trends observed in the Milky Way and other edge-on galaxies. The long recombination times in low-density gas result in persistent large volumes of ions with high ionization potentials, such as O  iii and Ne  iii. In particular, the vertically extended layers of Ne  iii in our time-dependent simulations result in [Ne iii ] 15  |$\mu$| m/[Ne ii ] 12  |$\mu$| m emission line ratios in agreement with observations of the edge-on galaxy NGC 891. Simulations adopting ionization equilibrium do not allow for the persistence of ions with high ionization states and therefore cannot reproduce the observed emission lines from low-density gas at high altitudes. [ABSTRACT FROM AUTHOR]

Dense clumps distributed along filaments are the immediate medium for star formation. Kinematic properties of the clumps, such as velocity gradient and angular momentum, combined with filament orientation, provide important clues to the formation mechanism of filament-clump configurations and the role of filaments in star formation. By cross-matching the Milky Way atlas for linear filaments and the structure, excitation and dynamics of the inner Galactic interstellar medium (SEDIGISM) |$^{13}$| CO (2-1) data, we aim to derive the velocity gradient and its direction, the specific angular momentum (⁠|$J/M$|⁠), and the ratio (⁠|$\beta$|⁠) between the rotational energy and gravitational energy of clumps, as well as to investigate the alignment between clump rotation and filament orientation. We found a monotonic increase in |$J/M$| as a function of clump size (R), following a power-law relation |$J/M~\propto ~R^{1.5\pm 0.2}$|⁠. The ratio |$\beta$| ranges from 1.1 |$\times$| 10 |$^{-5}$| to 0.1, with a median value 1.0 |$\times$| 10 |$^{-3}$|⁠ , suggesting that clump rotation provides insignificant support against gravitational collapse. The distribution of the angle between clump rotation and natal filament orientation is random, indicating that the clumps' rotational axes have no discernible correlation with the orientation of their hosting filaments. Counting only the most massive clump in each filament also finds no alignment between clump rotation and filament orientation. [ABSTRACT FROM AUTHOR]

The extraplanar diffuse ionised gas (eDIG) is a key component for understanding the feedback processes that connect galactic discs and their halos. In this paper, we present the second study of the Bidimensional Exploration of the warm-Temperature Ionised Gas (BETIS) project, the aim of which is to explore the possible ionisation mechanisms and characteristics of the eDIG. We use a sample of eight edge-on galaxies observed with the Multi-Unit Spectroscopic Explorer (MUSE) integral field spectrograph (IFS) and apply the methodology developed in the first paper of the BETIS project for obtaining binned emission line maps. We find that the vertical and radial profiles of the [N II]/Hα, [S II]/Hα, [O III]/Hβ, and [O I]/Hα ratios depict a complex ionisation structure within galactic halos – which is influenced by the spatial distribution of H II regions across the galactic plane as observed from our line of sight–, with Lyman continuum photon leakage from OB associations constituting the main ionisation source. Moreover, the electron temperature and S+/S ionisation ratio also exhibit a dependency on the distribution of H II regions within the galactic discs. Our analysis excludes low-mass, hot, and evolved stars (HOLMES) as viable candidates for secondary ionisation sources to elucidate the unusual behaviour of the line ratios at greater distances from the galactic midplane. In contrast, we ascertain that shocks induced in the interstellar medium by star formation(SF)-related feedback mechanisms represent a promising secondary ionisation source of the eDIG. We present a suite of models integrating ionisation mechanisms arising from fast shocks and photoionisation associated with star formation. When applied to the classical Baldwin–Phillips–Terlevich (BPT) diagrams, these models reveal that the ionisation budget of the eDIG ranges from 20% to 50% across our sample, with local variations of up to 20% within individual galaxy halos. This contribution correlates with the presence of filaments and other structural components observed within galaxy halos. The presence of shocks is additionally supported by the observation of a high density of high [O I]/Hα ratios, which is characteristic of shock-compressed ionised gas, and is likely induced by feedback from regions of intense SF within the galactic disc. These results demonstrate consistency across all galaxies analysed in this sample. [ABSTRACT FROM AUTHOR]

Context. Magnetic turbulence is classified as weak or strong based on the relative amplitude of the magnetic field fluctuations compared to the mean field. These two classes have different energy transport properties. Aims. The purpose of this study is to analyze turbulence in the interstellar medium (ISM) based on this classification. Specifically, we examined the ISM of simulated galaxies to detect evidence of strong magnetic turbulence and provide statistics on the associated magnetic coherent structures (MCoSs), such as current sheets, that arise in this context. Methods. We analyzed magnetohydrodynamic galaxy simulations with different initial magnetic field structures (either completely ordered or completely random) and recorded statistics on the magnetic field fluctuations (δB/B0) and the MCoSs, which are defined here as regions where the current density surpasses a certain threshold. We also studied the MCoS sizes and kinematics. Results. The magnetic field disturbances in both models follow a log-normal distribution, peaking at values close to unity; this distribution turns into a power law at large values (δB/B0 > 1), which is consistent with strong magnetic turbulence The current densities are widely distributed, with non-power-law deviations from a log-normal at the largest values. These deviating values of the current density define MCoSs. We find that, in both models, MCoSs are fractally distributed in space, with a typical volume-filling factor of about 10%, and tend to coincide with peaks of star formation density. Their fractal dimension is close to unity on sub-kiloparsec scales, and between 2 and 3 on larger scales. These values are consistent with MCoSs having a sheet-like or filament-like morphology. Conclusions. Our work challenges the prevailing paradigm of weak magnetic turbulence in the ISM by demonstrating that strong magnetic disturbances can occur even when the initial magnetic field is completely ordered. This strong magnetic turbulence arises self-consistently from differential rotation and supernova feedback. Our findings provide a foundation for a magnetic turbulence description of the galactic ISM that includes strong fluctuations of the magnetic field. [ABSTRACT FROM AUTHOR]

Context. Magnetohydrodynamic (MHD) turbulence plays a critical role in many key astrophysical processes, such as star formation, acceleration of cosmic rays, and heat conduction. However, its properties are still poorly understood. Aims. In this work, we explore how to extract the intermittency of compressible MHD turbulence from synthetic and real observations. Methods. We used three statistical methods, namely the probability distribution function, kurtosis, and scaling exponent of the multi-order structure function, to reveal the intermittency of MHD turbulence. Results. Our numerical results demonstrate that: (1) the synchrotron polarization intensity statistics can be used to probe the intermittency of magnetic turbulence, by which we can distinguish different turbulence regimes; (2) the intermittency of MHD turbulence is dominated by the slow mode in the sub-Alfvénic turbulence regime; and (3) the Galactic interstellar medium (ISM) in the low latitude region corresponds to the sub-Alfvénic and supersonic turbulence regime. Conclusions.We have successfully measured the intermittency of the Galactic ISM from synthetic and realistic observations. [ABSTRACT FROM AUTHOR]

Context. IRAS 16293E is a rare case of a prestellar core being subjected to the effects of at least one outflow. Aims. We want to disentangle the actual structure of the core from the outflow impact and evaluate the evolutionary stage of the core. Methods. Prestellar cores being cold and depleted, the best tracers of their central regions are the two isotopologues of the trihydrogen cation that are observable from the ground: ortho-H2D+ and para-D2H+. We used the Atacama Pathfinder EXperiment (APEX) telescope to map the para-D2H+ emission in IRAS 16293E and collected James Clerk Maxwell Telescope (JCMT) archival data of ortho-H2D+. We compared their emission to that of other tracers, including dust emission, and analysed their abundance with the help of a 1D radiative transfer tool. The ratio of the abundances of ortho-H2D+ to para-D2H+ can be used to estimate the stage of the chemical evolution of the core. Results. We have obtained the first complete map of para-D2H+ emission in a prestellar core. We compare it to a map of ortho-H2D+ and show their partial anti-correlation. This reveals a strongly evolved core with a para-D2H+/ortho-H2D+ abundance ratio towards the centre for which we obtain a conservative lower limit from 3.9 (at 12 K) to 8.3 (at 8 K), while the high extinction of the core is indicative of a central temperature below 10 K. This ratio is higher than predicted by the known chemical models found in the literature. Para-D2H+ (and indirectly ortho-H2D+) is the only species that reveals the true centre of this core, while the emission of other molecular tracers and dust are biased by the temperature structure that results from the impact of the outflow. Conclusions. This study is an invitation to reconsider the analysis of previous observations of this source and possibly questions the validity of the deuteration chemical models or of the reaction and inelastic collisional rate coefficients of the H+3 H 3 + isotopologue family. This could impact the deuteration clock predictions for all sources. [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]

Stars form in dense cores within molecular clouds, and newly formed stars influence their natal environments. How stellar feedback impacts core properties and evolution has been the subject of extensive investigation. We performed a hierarchical clustering (dendrogram) analysis of a STARFORGE (STAR FORmation in Gaseous Environments) simulation, modelling a giant molecular cloud to identify gas overdensities (cores) and study changes in their radius, mass, velocity dispersion, and virial parameter with respect to stellar feedback. We binned these cores on the basis of the fraction of gas affected by protostellar outflows, stellar winds, and supernovae and analysed the property distributions for each feedback bin. We find that cores that experience more feedback influence are smaller. Feedback notably enhances the velocity dispersion and virial parameter of the cores, more so than it reduces their radius. This is also evident in the linewidth–size relation, according to which cores in higher-feedback bins exhibit higher velocities than their similarly sized pristine counterparts. We conclude that stellar feedback mechanisms, which impart momentum to the molecular cloud, simultaneously compress and disperse the dense molecular gas. [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. Recent studies of the optical depth comparing [12C II] and [13C II] line profiles in Galactic star-forming regions have revealed strong self-absorption in [12C II] by low excitation foreground material. This implies a high column density for C+, corresponding to equivalent AV values of a few (up to about 10) mag. Aims. As the nature and origin of such a great column of cold C+ foreground gas are difficult to determine, it is essential to constrain the physical conditions of this material. Methods. We conducted high-resolution observations of [O I] 63 μm and [O I] 145 μm lines in M17 SW and Mon R2. The [O I] 145 μm transition traces warm PDR-material, while the [O I] 63 μm line traces the foreground material, as manifested by the absorption dips. Results. A comparison of both [O I] line profiles with [C II] isotopic lines confirm warm PDR-origin background emission and a significant column of cold foreground material, causing the self-absorption to be visible in the [12C II] and [O I] 63 μm profiles. In M17 SW, the C+ and O0 column densities are comparable for both layers. Mon R2 exhibits larger O0 columns compared to C+, indicating additional material where the carbon is neutral or in molecular form. Small-scale spatial variations in the foreground absorption profiles and the large column density (~1018 cm−2) of the foreground material suggest the emission is coming from high-density regions associated with the cloud complex – and not a uniform diffuse foreground cloud. Conclusions. The analysis confirms that the previously detected intense [C II] foreground absorption is attributable to a large column of low-excitation dense atomic material, where carbon is ionized and oxygen is in a neutral atomic form. [ABSTRACT FROM AUTHOR]

Context. Predicting the star formation rate (SFR) in galaxies is crucial to understand their evolution and morphology. To do so requires a fine understanding of how dense structures of gas are created and collapse. In that, turbulence and gravity play a major role. Aims. Within the gravo-turbulent framework, we assume that turbulence shapes the interstellar medium (ISM), creating density fluctuations that, if gravitationally unstable, will collapse and form stars. The goal of this work is to quantify how different regimes of turbulence, characterized by the strength and compressibility of the driving, shape the density field. We are interested in the outcome in terms of SFR and how it compares with existing analytical models for the SFR. Methods. We ran a series of hydrodynamical simulations of turbulent gas. The simulations were first conducted without gravity, so that the density and velocity were shaped by the turbulence driving. Gravity was then switched on, and the SFR was measured and compared with analytical models. The physics included in these simulations was very close to the one assumed in the classical gravo-turbulent SFR analytical models, which makes the comparison straightforward. Results. We found that the existing analytical models convincingly agree with simulations at low Mach number, but we measure a much lower SFR in the simulation with a high Mach number. We develop, in a companion paper, an updated physically motivated SFR model that reproduces well the inefficient high Mach regime of the simulations. Conclusions. Our work demonstrates that accurate estimations of the turbulent-driven replenishment time of dense structures and the dense gas spatial distribution are necessary to correctly predict the SFR in the high Mach regime. The inefficient high-Mach regime is a possible explanation for the low SFR found in dense and turbulent environments such as the centre of our Milky Way and other galaxies. [ABSTRACT FROM AUTHOR]

Context. Magnetic fields permeate the diffuse interstellar medium (ISM) of the Milky Way, and are essential to explain the dynamical evolution and current shape of the Galaxy. Magnetic fields reveal themselves via their influence on the surrounding matter, and as such are notoriously hard to measure independently of other tracers. Aims. In this work, we attempt to disentangle an all-sky map of the line-of-sight (LoS)-parallel component of the Galactic magnetic field from the Faraday effect, utilizing several tracers of the Galactic electron density, ne. Additionally, we aim to produce a Galactic electron dispersion measure map and quantify several tracers of the structure of the ionized medium of the Milky Way. Methods. The method developed to reach these aims is based on information field theory, a Bayesian inference framework for fields, which performs well when handling noisy and incomplete data and constraining high-dimensional-parameter spaces. We rely on compiled catalogs of extragalactic Faraday rotation measures and Galactic pulsar dispersion measures, a well as data on bremsstrahlung and the hydrogen α spectral line to trace the ionized medium of the Milky Way. Results. We present the first full sky map of the LoS-averaged Galactic magnetic field. Within this map, we find LoS-parallel and LoS-averaged magnetic field strengths of up to 4 µG, with an all-sky root mean square of 1.1 µG, which is consistent with previous local measurements and global magnetic field models. Additionally, we produce a detailed electron dispersion measure map that agrees with existing parametric models at high latitudes but suffers from systematic effects in the disk. Further analysis of our results with regard to the 3D structure of ne reveals that it follows a Kolmogorov-type turbulence for most of the sky. From the reconstructed dispersion measure and emission measure maps, we construct several tracers of variability in ne along the LoS. Conclusions. This work demonstrates the power of consistent joint statistical analysis including multiple datasets and physical quantities and defines a road map toward a full three-dimensional joint reconstruction of the Galactic magnetic field and the ionized ISM. [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]

Magnetic fields play an important role in shaping and regulating star formation in molecular clouds. Here, we present one of the first studies examining the relative orientations between magnetic (B) fields and the dust emission, gas column density, and velocity centroid gradients on the 0.02 pc (core) scales, using the BISTRO and VLA+GBT observations of the NGC 1333 star-forming clump. We quantified these relative orientations using the Project Rayleigh Statistic (PRS) and found preferential global parallel alignment between the B field and dust emission gradients, consistent with large-scale studies with Planck. No preferential global alignments, however, are found between the B field and velocity gradients. Local PRS calculated for subregions defined by either dust emission or velocity coherence further revealed that the B field does not preferentially align with dust emission gradients in most emission-defined subregions, except in the warmest ones. The velocity-coherent structures, on the other hand, also showed no preferred B field alignments with velocity gradients, except for one potentially bubble-compressed region. Interestingly, the velocity gradient magnitude in NGC 1333 ubiquitously features prominent ripple-like structures that are indicative of magnetohydrodynamic (MHD) waves. Finally, we found B field alignments with the emission gradients to correlate with dust temperature and anticorrelate with column density, velocity dispersion, and velocity gradient magnitude. The latter two anticorrelations suggest that alignments between gas structures and B fields can be perturbed by physical processes that elevate velocity dispersion and velocity gradients, such as infall, accretions, and MHD waves. [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]

Three-dimensional dust density maps are crucial for understanding the structure of the interstellar medium of the Milky Way and the processes that shape it. However, constructing these maps requires large data sets and the methods used to analyse them are computationally expensive and difficult to scale up. As a result, it has only recently become possible to map kiloparsec-scale regions of our Galaxy at parsec-scale grid sampling. We present all-sky three-dimensional dust density and extinction maps of the Milky Way out to 2.8 kpc in distance from the Sun using the fast and scalable Gaussian Process algorithm Dustribution. The sampling of the three-dimensional map is l, b, d  = 1° × 1° × 1.7 pc. The input extinction and distance catalogue contains 120 million stars with photometry and astrometry from Gaia DR2, 2MASS and AllWISE. This combines the strengths of optical and infrared data to probe deeper into the dusty regions of the Milky Way. We compare our maps with other published 3D dust maps. All maps quantitatively agree at the 0.001 mag pc−1 scale with many qualitatively similar features, although each map also has its own features. We recover Galactic features previously identified in the literature. Moreover, we also see a large under-density that may correspond to an inter-arm or -spur gap towards the Galactic Centre. [ABSTRACT FROM AUTHOR]

We use smoothed particle hydrodynamics simulations of isolated Milky Way-mass disc galaxies that include cold, interstellar gas to test subgrid prescriptions for star formation (SF). Our fiducial model combines a Schmidt law with a gravitational instability criterion, but we also test density thresholds and temperature ceilings. While SF histories are insensitive to the prescription for SF, the Kennicutt–Schmidt (KS) relations between SF rate and gas surface density can discriminate between models. We show that our fiducial model, with an SF efficiency per free-fall time of 1 per cent, agrees with spatially resolved and azimuthally averaged observed KS relations for neutral, atomic, and molecular gas. Density thresholds do not perform as well. While temperature ceilings selecting cold, molecular gas can match the data for galaxies with solar metallicity, they are unsuitable for very low-metallicity gas and hence for cosmological simulations. We argue that SF criteria should be applied at the resolution limit rather than at a fixed physical scale, which means that we should aim for numerical convergence of observables rather than of the properties of gas labelled as star-forming. Our fiducial model yields good convergence when the mass resolution is varied by nearly 4 orders of magnitude, with the exception of the spatially resolved molecular KS relation at low surface densities. For the gravitational instability criterion, we quantify the impact on the KS relations of gravitational softening, the SF efficiency, and the strength of supernova feedback, as well as of observable parameters such as the inclusion of ionized gas, the averaging scale, and the metallicity. [ABSTRACT FROM AUTHOR]

We present an analysis of the diffuse ionized gas (DIG) in a high-resolution simulation of an isolated Milky Way-like galaxy, incorporating on-the-fly radiative transfer and non-equilibrium thermochemistry. We utilize the Monte-Carlo radiative transfer code colt to self-consistently obtain ionization states and line emission in post-processing. We find a clear bimodal distribution in the electron densities of ionized gas (⁠|$n_{\rm e}$|⁠), allowing us to define a threshold of |$n_{\rm e}=10\, \mathrm{cm}^{-3}$| to differentiate DIG from |${\rm H\, {\small II}}$| regions. The DIG is primarily ionized by stars aged 5 – 25 Myr, which become exposed directly to low-density gas after |${\rm H\, {\small II}}$| regions have been cleared. Leakage from recently formed stars (⁠|$\lt 5$|  Myr) is only moderately important for DIG ionization. We forward model local observations and validate our simulated DIG against observed line ratios in [ |${\rm S\, {\small II}}$| ]/H |$\alpha$|⁠ , [ |${\rm N\, {\small II}}$| ]/H |$\alpha$|⁠ , [ |${\rm O\, {\small I}}$| ]/H |$\alpha$|⁠ , and [ |${\rm O\, {\small III}}$| ]/H |$\beta$| against |$\Sigma _{\rm H\alpha }$|⁠. The mock observations not only reproduce observed correlations, but also demonstrate that such trends are related to an increasing temperature and hardening ionizing radiation field with decreasing |$n_{\rm e}$|⁠. The hardening of radiation within the DIG is caused by the gradual transition of the dominant ionizing source with decreasing |$n_{\rm e}$| from 0 to 25 Myr stars, which have progressively harder intrinsic ionizing spectra primarily due to the extended Wolf–Rayet phase caused by binary interactions. Consequently, the DIG line ratio trends can be attributed to ongoing star formation, rather than secondary ionization sources, and therefore present a potent test for stellar feedback and stellar population models. [ABSTRACT FROM AUTHOR]

To advance our understanding of the evolution of the interstellar medium (ISM) of our Galaxy, numerical models of Milky Way (MW) type galaxies are widely used. However, most models only vaguely resemble the MW (e.g. in total mass), and often use imposed analytic potentials (which cannot evolve dynamically). This poses a problem in asserting their applicability for the interpretation of observations of our own Galaxy. The goal of this work is to identify a numerical model that is not only an MW-type galaxy, but one that can mimic some of the main observed structures of our Galaxy, using dynamically evolving potentials, so that it can be used as a base model to study the ISM cycle in a galaxy like our own. This paper introduces a suite of 15 MW-type galaxy models developed using the arepo numerical code, that are compared to Galactic observations of |$^{12}$| CO and H  i emission via longitude–velocity plots, from where we extract and compare the skeletons of major galactic features and the terminal gas velocities. We found that our best-fitting model to the overall structure, also reproduces some of the more specific observed features of the MW, including a bar with a pattern speed of |$30.0 \pm 0.2$| km s |$^{-1}$|  kpc |$^{-1}$|⁠ , and a bar half-length of |$3.2 \pm 0.8$|  kpc. Our model shows large streaming motions around spiral arms, and strong radial motions well beyond the inner bar. This model highlights the complex motions of a dynamic MW-type galaxy and has the potential to offer valuable insight into how our Galaxy regulates the ISM and star formation. [ABSTRACT FROM AUTHOR]

Recent advancements in 3D dust mapping have transformed our understanding of the Milky Way's local interstellar medium, enabling us to explore its structure in three spatial dimensions for the first time. In this Letter, we use the most recent 3D dust map by Edenhofer et al. to study the well-known Chameleon, Musca, and Coalsack cloud complexes, located about 200 pc from the Sun. We find that these three complexes are not isolated but rather connect to form a surprisingly well-defined half-ring, constituting a single C-shaped cloud with a radius of about 50 pc, a thickness of about 45 pc, and a total mass of about 5 × 104M⊙, or 9 × 104M⊙ if including everything in the vicinity of the C-shaped cloud. Despite the absence of an evident feedback source at its center, the dynamics of young stellar clusters associated with the C structure suggest that a single supernova explosion about 4 Myr–10 Myr ago likely shaped this structure. Our findings support a single origin story for these cloud complexes, suggesting that they were formed by feedback-driven gas compression, and offer new insights into the processes that govern the birth of star-forming clouds in feedback-dominated regions, such as the Scorpius-Centaurus association. [ABSTRACT FROM AUTHOR]

Local gravitational instability (LGI) is considered crucial for regulating star formation and gas turbulence in galaxy discs, especially at high redshift. Instability criteria usually assume infinitesimally thin discs or rely on approximations to include the stabilising effect of the gas disc thickness. We test a new 3D instability criterion for rotating gas discs that are vertically stratified in an external potential. This criterion reads Q3D < 1, where Q3D is the 3D analogue of the Toomre parameter Q. The advantage of Q3D is that it allows us to study LGI in and above the galaxy midplane in a rigorous and self-consistent way. We apply the criterion to a sample of 44 star-forming galaxies at 0 ≲ z ≲ 5 hosting rotating discs of cold gas. The sample is representative of galaxies on the main sequence at z ≈ 0 and includes massive star-forming and starburst galaxies at 1 ≲ z ≲ 5. For each galaxy, we first apply the Toomre criterion for infinitesimally thin discs, finding ten unstable systems. We then obtain maps of Q3D from a 3D model of the gas disc derived in the combined potential of dark matter, stars and the gas itself. According to the 3D criterion, two galaxies with Q < 1 show no evidence of instability and the unstable regions that are 20% smaller than those where Q < 1. No unstable disc is found at 0 ≲ z ≲ 1, while ≈60% of the systems at 2 ≲ z ≲ 5 are locally unstable. In these latter, a relatively small fraction of the total gas (≈30%) is potentially affected by the instability. Our results disfavour LGI as the main regulator of star formation and turbulence in moderately star-forming galaxies in the present-day Universe. LGI likely becomes important at high redshift, but the input by other mechanisms seems required in a significant portion of the disc. We also estimate the expected mass of clumps in the unstable regions, offering testable predictions for observations. [ABSTRACT FROM AUTHOR]

We perform an infrared (IR) spectral and visible morphological study of two young planetary nebulae (YPNe) Hen 2-73 and Hen 2-96 using archival Spitzer Space Telescope and Hubble Space Telescope (HST) observations to understand their dust properties and nebular structures. High-resolution HST images of these nebulae show several bipolar lobes and ionised tori in the central regions of both objects. The presence of these multi-lobe structures suggests that the formation process of these nebulae is complex. To search for a possible link between the central sources and multipolar appearances of these objects, the Transiting Exoplanet Survey Satellite (TESS) observations are used to examine whether their central stars (CSs) exhibit periodic photometric variations. In the TESS observations, the CS light curve of Hen 2-96 shows a photometric variation with a period of 2.23 h. The IR spectra of these two YPNe suggest that the nebulae have mixed dust environments, which are associated with the presence of dense tori created by central binary interactions in these objects. Two three-dimensional models are constructed to study the complex nebular structures of the YPNe. These simulations suggest that the number of multipolar YPNe may be larger than observed. In addition, we analyse the spectral energy distributions of these nebulae to study their gas, dust, and photospheric components. [ABSTRACT FROM AUTHOR]

Aims. Recent simulations have shown that asymmetries in the ejecta distribution of supernova remnants (SNRs) may be a reflection of asymmetries left over from the initial supernova explosion. Thus, SNR studies provide a vital means for testing and constraining model predictions in relation to the distribution of heavy elements, which are key to improving our understanding of the explosion mechanisms in Type Ia supernovae. Methods. The use of a novel blind source separation method applied to the megasecond X-ray observations of the historic Kepler and Tycho supernova remnants has revealed maps of the ejecta distribution. These maps are endowed with an unprecedented level of detail and clear separations from the continuum emission. Our method also provides a three-dimensional (3D) view of the ejecta by individually disentangling red- and blueshifted spectral components associated with images of the Si, S, Ar, Ca, and Fe emission. This approach provides insights into the morphology of the ejecta distribution in those two remnants. Results. Those mappings have allowed us to thoroughly investigate the asymmetries in the intermediate-mass elements and Fe distribution in two Type Ia supernova remnants. We also compared the results with the core-collapse Cassiopeia A remnant, which we had studied previously. The images obtained confirm, as expected for Type Ia SNRs, that the Fe distribution is mostly closer to the core than that of intermediate-mass elements. They also highlight peculiar features in the ejecta distribution, such as the Fe-rich southeastern knot in Tycho. [ABSTRACT FROM AUTHOR]

Faraday tomography observations with the Low Frequency Array (LOFAR) have unveiled a remarkable network of structures in polarized synchrotron emission at high Galactic latitudes. The observed correlation between LOFAR structures, dust polarization, and HI emission suggests a connection to the neutral interstellar medium (ISM). We investigated this relationship by estimating the rotation measure (RM) of the warm neutral (partially ionized) medium (WNM) in the local ISM. Our work combines UV spectroscopy from FUSE and dust polarization observations from Planck with LOFAR data. We derived electron column densities from UV absorption spectra toward nine background stars, within the field of published data from the LOFAR two-meter sky survey. The associated RMs were estimated using a local magnetic field model fitted to the dust polarization data of Planck. A comparison with Faraday spectra at the position of the stars suggests that LOFAR structures delineate a slab of magnetized WNM and synchrotron emission, located ahead of the bulk of the warm ionized medium. This conclusion establishes an astrophysical framework for exploring the link between Faraday structures and the dynamics of the magnetized multiphase ISM. It will be possible to test it on a larger sample of stars when maps from the full northern sky survey of LOFAR become available. [ABSTRACT FROM AUTHOR]

We present the Bidimensional Exploration of the warm-Temperature Ionised gaS (BETIS) project, designed for the spatial and spectral study of the diffuse ionised gas (DIG) in a selection of nearby spiral galaxies observed with the MUSE integral-field spectrograph. Our primary objective is to investigate the various ionisation mechanisms at play within the DIG. We analysed the distribution of high- and low-ionisation species in the optical spectra of the sample on a spatially resolved basis. We introduced a new methodology for spectroscopically defining the DIG, optimised for galaxies of different resolutions. Firstly, we employed an innovative adaptive binning technique on the observed datacube based on the spectroscopic signal-to-noise ratio (S/N) of the collisional [S II] line to increase the S/N of the rest of the lines including [O III], [O I], and He I. Subsequently, we created a DIG mask by eliminating the emissions associated with both bright and faint H II regions. We also examined the suitability of using Hα equivalent width (EWHα) as a proxy for defining the DIG and its associated ionisation regime. Notably, for EWHα < 3 Å – the expected emission from hot low-mass evolved stars (HOLMES) – the measured value is contingent on the chosen population synthesis technique performed. Our analysis of the showcase sample reveals a consistent cumulative DIG fraction across all galaxies in the sample, averaging around 40%–70%. The average radial distribution of the [N II]/Hα, [S II]/Hα, [O I]/Hα, and [O III]/Hβ ratios are enhanced in the DIG regimes (up to 0.2 dex). It follows similar trends between the DIG regime and the H II regions, as well as the Hα surface brightness (ΣHα), indicating a correlation between the ionisation of these species in both the DIG and the H II regions. The DIG loci in typical diagnostic diagrams are found, in general, within the line ratios that correspond to photoionisation due to the star formation. There is a noticeable offset correspondent to ionisation due to fast shocks. However, an individual diagnosis performed for each galaxy reveals that all the DIG in these galaxies can be attributed to photoionisation from star formation. The offset is primarily due to the contribution of Seyfert galaxies in our sample, which is closely aligned with models of ionisation from fast shocks and galactic outflows, thus mimicking the DIG emission. Our results indicate that galaxies exhibiting active galactic nucleus (AGN) activity should be considered separately when conducting a general analysis of the DIG ionisation mechanisms, since this emission is indistinguishable from high-excitation DIG. [ABSTRACT FROM AUTHOR]

Context. We conducted a follow-up study on the analysis of the LOw Frequency ARray (LOFAR) Two-metre Sky Survey (LoTSS) mosaic in the high-latitude outer Galaxy presented in the first paper of this series. Here, we focus on the search for alignment between the magnetic field traced by dust, HI filaments, starlight optical linear polarisation, and linear depolarised structures (depolarisation canals) observed in low-frequency synchrotron polarisation. This alignment was previously found in several smaller fields observed with LOFAR, offering valuable insights into the nature of the interstellar medium and the 3D spatial distribution of the diffuse ionised medium. Aims. We aim to determine whether the alignment of the interstellar medium (ISM) phases observed through multiple tracers is a common occurrence or an exception. Additionally, in areas where depolarisation canals align with the magnetic field, we use starlight polarisation to constrain the distance to the structures associated with the observed canals. Methods. We employed the Rolling Hough Transform (RHT) and projected Rayleigh statistics (PRS) to identify and quantify the alignment between the different tracers. We used these tools to detect linear features in the data and quantify the significance of the orientation trends between pairs of tracers. Results. On the scale of the whole mosaic, we did not find any evidence of a universal alignment among the three tracers. However, in one particular area, the western region (Dec between 29° and 70° and RA between 7h44m and 9h20m), we do find a significant alignment between the magnetic field, depolarisation canals, and HI filaments. Based on this alignment, we used the starlight polarisation of stars with known parallax distances to estimate that the minimum distance to the structures observed by LOFAR in this region lies within the range of 200 to 240 pc. We associate these structures with the edge of the Local Bubble. [ABSTRACT FROM AUTHOR]

The dynamic spectra of pulsars frequently exhibit diverse interference patterns, often associated with parabolic arcs in the Fourier-transformed (secondary) spectra. Our approach differs from previous ones in two ways: first, we extend beyond the traditional Fresnel–Kirchhoff method by using the Green's function of the Helmholtz equation, i.e. we consider spherical waves originating from three-dimensional space, not from a two-dimensional screen. Secondly, the discrete structures observed in the secondary spectrum result from discrete scatterer configurations, namely plasma concentrations in the interstellar medium, and not from the selection of points by the stationary phase approximation. Through advanced numerical techniques, we model both the dynamic and secondary spectra, providing a comprehensive framework that describes all components of the latter spectra in terms of physical quantities. Additionally, we provide a thorough analytical explanation of the secondary spectrum. [ABSTRACT FROM AUTHOR]

We suggest that filaments in star-forming regions undergo frequent mergers. As stellar nurseries, filaments play a vital role in understanding star formation and mergers could pave the way for understanding the formation of more complex filamentary systems, such as networks and hubs. We compare the physical properties derived from hydrodynamic ramses simulations of merging filaments to those obtained from ALMA observations towards the LDN 1641-North region in Orion. We find similarities in the distributions of line-mass, column density, and velocity dispersion. Such common features support the hypothesis of filament mergers shaping the structure of the interstellar medium. [ABSTRACT FROM AUTHOR]

Widespread, high altitude, diffuse ionized gas with scale heights of around a kiloparsec is observed in the Milky Way and other star-forming galaxies. Numerical radiation-magnetohydrodynamic simulations of a supernova-driven turbulent interstellar medium show that gas can be driven to high altitudes above the galactic mid-plane, but the degree of ionization is often less than inferred from observations. For computational expediency, ionizing radiation from massive stars is often included as a post-processing step assuming ionization equilibrium. We extend our simulations of an Milky Way-like interstellar medium to include the combined effect of supernovae and photoionization feedback from mid-plane OB stars and a population of hot evolved low mass stars. The diffuse ionized gas has densities below 0.1 |${\rm \,cm^{-3}}$|⁠ , so recombination time-scales can exceed millions of years. Our simulations now follow the time-dependent ionization and recombination of low density gas. The long recombination time-scales result in diffuse ionized gas that persists at large altitudes long after the deaths of massive stars that produce the vast majority of the ionized gas. The diffuse ionized gas does not exhibit the large variability inherent in simulations that adopt ionization equilibrium. The vertical distribution of neutral and ionized gas is close to what is observed in the Milky Way. The volume filling factor of ionized gas increases with altitude resulting in the scale height of free electrons being larger than that inferred from H  |$\alpha$| emission, thus reconciling the observations of ionized gas made in H  |$\alpha$| and from pulsar dispersion measurements. [ABSTRACT FROM AUTHOR]

Recent advances in identifying giant molecular filaments in Galactic surveys allow us to study the interstellar material and its dense, potentially star forming phase on scales comparable to resolved extragalactic clouds. Two large filaments detected in the 13CO/C18O(J  = 3–2) Heterodyne Inner Milky Way Plane Survey (CHIMPS) survey, one in the Sagittarius-arm and one in an interarm region, were mapped with dense gas tracers inside a 0.06 square degrees area and with a spatial resolution of around 0.4 and 0.65 pc at the distance of the targets using the 30 m telescope of the Institut de Radioastronomie Millimétrique (IRAM) to investigate the environmental dependence of the dense gas fraction. The N2H+(1 − 0) transition, an excellent tracer of the dense gas, was detected in parsec-scale, elliptical clumps and with a filling factor of around 8.5  per cent in our maps. The N2H+-emitting areas appear to have higher dense gas fraction (e.g. the ratio of N2H+ and 13CO emission) in the interarm than in the arm which is opposite to the behaviour found by previous studies, using dust emission rather than N2H+ as a tracer of dense gas. However, the arm filament is brighter in 13CO and the infrared emission of dust, and the dense gas fraction determined as above is governed by the 13CO brightness. We caution that measurements regarding the distribution and fraction of dense gas on these scales may be influenced by many scale- and environment-dependent factors, as well as the chemistry and excitation of the particular tracers, then consider several scenarios that can reproduce the observed effect. [ABSTRACT FROM AUTHOR]

Aims. W43-main is a massive molecular complex undergoing starburst activities, located at the interaction of the Scutum arm and the Galactic bar. We aim to investigate the gas dynamics, in particular, the prevailing shock signatures from cloud to clump scales. We also look to assess the impact of shocks on the formation of dense gas and early-stage cores in OB cluster formation processes. Methods. We carried out NOEMA and IRAM-30 m observations at 3 mm towards five molecular gas clumps in W43 main located within large-scale interacting gas components. We used CH3CCH and H2CS lines to trace the extended gas temperature and CH3OH lines to probe the volume density of the dense gas components (≳105 cm−3). We adopted multiple tracers that are sensitive to different gas density regimes to reflect the global gas motions. The density enhancements constrained by CH3OH and a population of NH2D cores are correlated (in the spatial and velocity domains) with SiO emission, which is a prominent indicator of shock processing in molecular clouds. Results. The emission of SiO (2–1) is extensive across the region (~4 pc) and it is contained within a low-velocity regime, hinting at a large-scale origin for the shocks. Position-velocity maps of multiple tracers show systematic spatio-kinematic offsets supporting the cloud-cloud collision-merging scenario. We identified an additional extended velocity component in the CCH emission, which coincides with one of the velocity components of the larger scale 13CO (2−1) emission, likely representing an outer, less-dense gas layer in the cloud merging process. We find that the 'V-shaped', asymmetric SiO wings are tightly correlated with localised gas density enhancements, which is direct evidence of dense gas formation and accumulation in shocks. The dense gas that is formed in this way may facilitate the accretion of the embedded, massive pre-stellar and protostellar cores. We resolved two categories of NH2D cores: those exhibiting only subsonic to transonic velocity dispersions and those with an additional supersonic velocity dispersion. The centroid velocities of the latter cores are correlated with the shock front seen via SiO. The kinematics of the ~0.1 pc NH2D cores are heavily imprinted by shock activities and may represent a population of early-stage cores forming around the shock interface. [ABSTRACT FROM AUTHOR]

Context. High-resolution 3D maps of interstellar dust are critical for probing the underlying physics shaping the structure of the interstellar medium, and for foreground correction of astrophysical observations affected by dust. Aims. We aim to construct a new 3D map of the spatial distribution of interstellar dust extinction out to a distance of 1.25 kpc from the Sun. Methods. We leveraged distance and extinction estimates to 54 million nearby stars derived from the Gaia BP/RP spectra. Using the stellar distance and extinction information, we inferred the spatial distribution of dust extinction. We modeled the logarithmic dust extinction with a Gaussian process in a spherical coordinate system via iterative charted refinement and a correlation kernel inferred in previous work. In total, our posterior has over 661 million degrees of freedom. We probed the posterior distribution using the variational inference method MGVI. Results. Our 3D dust map has an angular resolution of up to 14′ (Nside = 256), and we achieve parsec-scale distance resolution, sampling the dust in 516 logarithmically spaced distance bins spanning 69 pc to 1250 pc. We generated 12 samples from the variational posterior of the 3D dust distribution and release the samples alongside the mean 3D dust map and its corresponding uncertainty. Conclusions. Our map resolves the internal structure of hundreds of molecular clouds in the solar neighborhood and will be broadly useful for studies of star formation, Galactic structure, and young stellar populations. It is available for download in a variety of coordinate systems online and can also be queried via the publicly available dustmaps Python package. [ABSTRACT FROM AUTHOR]

Magnetic fields play a crucial role in star formation, yet tracing them becomes particularly challenging, especially in the presence of outflow feedback in protostellar systems. We targeted the star-forming region L1551, notable for its apparent outflows, to investigate the magnetic fields. These fields were probed using polarimetry observations from the Planck satellite at 353 GHz/849 μm, the Stratospheric Observatory for Infrared Astronomy's (SOFIA) High-resolution Airborne Wide-band Camera (HAWC+) measurement at 214 μm, and the James Clerk Maxwell Telescope's (JCMT) Submillimetre Common-User POLarimeter (SCUPOL) 850 μm survey. Consistently, all three measurements show that the magnetic fields twist towards the protostar IRS 5. Additionally, we utilized the velocity gradient technique on the 12CO (J = 1–0) emission data to distinguish the magnetic fields directly associated with the protostellar outflows. These were then compared with the polarization results. Notably, in the outskirts of the region, these measurements generally align. However, as one approaches the centre of IRS 5, the measurements tend to yield mostly perpendicular relative orientations. This suggests that the outflows might be dynamically significant from a scale of ∼0.2 pc, causing the velocity gradient to change direction by 90°. Furthermore, we discovered that the polarization fraction p and the total intensity I correlate as p  ∝  I −α. Specifically, α is approximately 1.044 ± 0.06 for SCUPOL and around 0.858 ± 0.15 for HAWC+. This indicates that the outflows could significantly impact the alignment of dust grains and magnetic fields in the L1551 region. [ABSTRACT FROM AUTHOR]

The fragmentation of gas to form stars in molecular clouds is intrinsically linked to the turbulence within them. These internal motions are set at the birth of the cloud and may vary with galactic environment and as the cloud evolves. In this paper, we introduce a new suite of 15 high-resolution 3D molecular cloud simulations using the moving mesh code arepo to investigate the role of different decaying turbulent modes (mixed, compressive, and solenoidal) and virial ratios on the evolution of a |$10^4\, \mathrm{M}_{\odot }$| molecular cloud. We find that diffuse regions maintain a strong relic of the initial turbulent mode, whereas the initial gravitational potential dominates dense regions. Solenoidal seeded models thus give rise to a diffuse cloud with filament-like morphology, and an excess of brown dwarf mass fragments. Compressive seeded models have an early onset of star-formation, centrally condensed morphologies and a higher accretion rate, along with overbound clouds. 3D filaments identified using disperse and analysed through a new python toolkit we develop and make publicly available with this work called fiesta , show no clear trend in lengths, masses and densities between initial turbulent modes. Overbound clouds, however, produce more filaments and thus have more mass in filaments. The hubs formed by converging filaments are found to favour star-formation, with surprisingly similar mass distributions independent of the number of filaments connecting the hub. [ABSTRACT FROM AUTHOR]

Young Stellar Objects (YSO) are newly formed stars from molecular clouds. They stay close to where they were born and serve as good tracers to study gas and star formation. During cloud evolution, young massive stars can disrupt the surrounding gas through stellar feedback, changing the gas distribution. We study the distribution of the gas around a sample of YSO associations located at |$d \lesssim 3.5 \,\,\rm kpc$| from the Sun by comparing the location and morphology between 12CO (J  = 1–0) emission, Planck 870  |$\mu$| m maps and YSO associations. Based on the spatial distribution of the gas compared to that of the YSOs, we classify the YSO associations into six types: direct, close, bubble, complex, diffuse, and clean associations. The complex associations are large structures consisting of both gas-rich and gas-poor segments. We study the velocity dispersion-size relation towards different association types. From the ratio between different types, we estimate a feedback time of ≈ 1.7 Myr in the solar neighbourhood. The sample sets a solid foundation to explore the relationship between interstellar medium evolution, star formation, and Galaxy structure. [ABSTRACT FROM AUTHOR]

Geostatistical methods are powerful tools for understanding the spatial structure of the metallicity distribution of galaxies, and enable construction of accurate predictive models of the 2D metallicity distribution. However, so far these methods have only been applied to very high spatial resolution metallicity maps, leaving it uncertain if they will work on lower quality data. In this study, we apply geostatistical techniques to high-resolution spectroscopic maps of three local galaxies convolved to eight different spatial resolutions ranging from ∼40 pc to ∼1 kpc per pixel. We fit a geostatistical model to the data at all resolutions, and find that for metallicity maps where small-scale structure is visible by eye (with ≳10 resolution elements per R e), all parameters, including the metallicity correlation scale, can be recovered accurately. At all resolutions tested, we find that point metallicity predictions from such a geostatistical model outperform a circularly symmetric metallicity gradient model. We also explore dependence on the number of data points, and find that N ≳ 100 spatially resolved metallicity values are sufficient to train a geostatistical model that yields more accurate metallicity predictions than a radial gradient model. Finally, we investigate the potential detrimental effects of having spaxels smaller than an individual H  ii region by repeating our analysis with metallicities integrated over H  ii regions. We see that spaxel-based measurements have more noise, as expected, but the underlying spatial metallicity distribution can be recovered regardless of whether spaxels or integrated regions are used. [ABSTRACT FROM AUTHOR]

We present the preliminary results of 4 yr of observations of the scattering of pulsar radiation in the interstellar medium using the Low Frequency Array (LOFAR) PL611 station located in Lazy near Krakow. In this work we show the initial results for six pulsars from our observing campaign. We used the HBA antennas of the station, with a central frequency of 154 MHz and a 72 MHz bandwidth and we were able to detect the frequency-dependent change in the pulse profiles. Splitting the bandwidth into a number of separate sub-bands we obtained independent scatter time estimates, which allowed for the estimation of the scattering frequency scaling slope for each individual observation. Our project has been accumulating data for more than 4 yr and as a result we are in a position to study the time variability of the scattering parameters over this period. We detected significant changes in the scatter time and its frequency scaling for at least two of pulsars. The average values of the frequency scaling index for the six pulsars selected for this work are below the range of values predicted by the simple thin screen model of interstellar scattering. This is in accordance with previous results shown for both higher dispersion measure pulsars observed in the past, as well as the more recent LOFAR observations and other projects similar in the observing frequency. We also discuss the advantages of using individual LOFAR stations (or similar instruments) for this kind of research. [ABSTRACT FROM AUTHOR]

The Sun is located close to the Galactic mid-plane, meaning that we observe the Galaxy through significant quantities of dust. Moreover, the vast majority of the Galaxy's stars also lie in the disc, meaning that dust has an enormous impact on the massive astrometric, photometric and spectroscopic surveys of the Galaxy that are currently underway. To exploit the data from these surveys we require good three-dimensional maps of the Galaxy's dust. We present a new method for making such maps in which we form the best linear unbiased predictor of the extinction at an arbitrary point based on the extinctions for a set of observed stars. This method allows us to avoid the artificial inhomogeneities (so-called 'fingers of God') and resolution limits that are characteristic of many published dust maps. Moreover, it requires minimal assumptions about the statistical properties of the interstellar medium. In fact, we require only a model of the first and second moments of the dust density field. The method is suitable for use with directly measured extinctions, such as those provided by the Rayleigh–Jeans colour excess method, and inferred extinctions, such as those provided by hierarchical Bayesian models like StarHorse. We test our method by mapping dust in the region of the giant molecular cloud Orion A. Our results indicate a foreground dust cloud at a distance of 350 pc, which has been identified in work by another author. [ABSTRACT FROM AUTHOR]

Cosmic rays (CRs) with energies ≪ TeV comprise a significant component of the interstellar medium (ISM). Major uncertainties in CR behaviour on observable scales (much larger than CR gyroradii) stem from how magnetic fluctuations scatter CRs in pitch angle. Traditional first-principles models, which assume these magnetic fluctuations are weak and uniformly scatter CRs in a homogeneous ISM, struggle to reproduce basic observables such as the dependence of CR residence times and scattering rates on rigidity. We therefore explore a new category of 'patchy' CR scattering models, wherein CRs are pre-dominantly scattered by intermittent strong scattering structures with small volume-filling factors. These models produce the observed rigidity dependence with a simple size distribution constraint, such that larger scattering structures are rarer but can scatter a wider range of CR energies. To reproduce the empirically inferred CR scattering rates, the mean free path between scattering structures must be |$\ell _{\rm mfp}\sim 10\, {\rm pc}$| at GeV energies. We derive constraints on the sizes, internal properties, mass/volume-filling factors, and the number density any such structures would need to be both physically and observationally consistent. We consider a range of candidate structures, both large scale (e.g. H ii  regions) and small scale (e.g. intermittent turbulent structures, perhaps even associated with radio plasma scattering) and show that while many macroscopic candidates can be immediately ruled out as the primary CR scattering sites, many smaller structures remain viable and merit further theoretical study. We discuss future observational constraints that could test these models. [ABSTRACT FROM AUTHOR]

Magnetic fields are a defining yet enigmatic aspect of the interstellar medium, with their three-dimensional (3D) mapping posing a substantial challenge. In this study, we harness the innovative velocity gradient technique (VGT), underpinned by magnetohydrodynamic turbulence theories, to map the magnetic field structure by applying it to the atomic neutral hydrogen (H  i) emission line and the molecular tracer 12CO. We construct the tomography of the magnetic field in the low-mass star-forming region L1688, utilizing two approaches: (1) VGT-H  i combined with the Galactic rotational curve, and (2) stellar polarization paired with precise star parallax measurements. Our analysis reveals that the magnetic field orientations deduced from stellar polarization undergo a distinct directional change in the vicinity of L1688, providing evidence that the misalignment between VGT-H  i and stellar polarization stems from the influence of the molecular cloud's magnetic field on the polarization of starlight. When comparing VGT-12CO to stellar polarization and Planck polarization data, we observe that VGT-12CO effectively reconciles the misalignment noted with VGT-H  i , showing statistical alignment with Planck polarization measurements. This indicates that VGT-12CO could be integrated with VGT-H  i , offering vital insights into the magnetic fields of molecular clouds, thereby enhancing the accuracy of our 3D magnetic field reconstructions. [ABSTRACT FROM AUTHOR]

To understand physical properties of the interstellar medium (ISM) on various scales, we investigate it at parsec resolution on the kiloparsec scale. Here, we report on the sub-kpc scale gas density histogram (GDH) of the Milky Way. The GDH is a density probability distribution function (PDF) of the gas volume density. Using this method, we are free from an identification of individual molecular clouds and their spatial structures. We use survey data of 12CO and 13CO (J = 1–0) emission in the Galactic plane (l = 10○–50○) obtained as a part of the FOREST Unbiased Galactic plane Imaging survey with the Nobeyama 45m telescope (FUGIN). We make a GDH for every channel map of 2○ × 2○ area including the blank sky component, and without setting cloud boundaries. This is a different approach from previous works for molecular clouds. The GDH fits well to a single or double lognormal distribution, which we name the low-density lognormal (L-LN) and high-density lognormal (H-LN) components, respectively. The multi-lognormal components suggest that the L-LN and H-LN components originate from two different stages of structure formation in the ISM. Moreover, we find that both the volume ratios of H-LN components to total (f H) and the width of the L-LN along the gas density axis (σL) show coherent structure in the Galactic-plane longitude-velocity diagram. It is possible that these GDH parameters are related to strong galactic shocks and other weak shocks in the Milky Way. [ABSTRACT FROM AUTHOR]

We present an analysis of the outer Galaxy giant molecular filament (GMF) G214.5-1.8 (G214.5) using IRAM 30m data of 12CO, 13CO, and C18O. We find that the 12CO (1-0) and (2-1) derived excitation temperatures are near identical and are very low, with a median of 8.2 K, showing that the gas is extremely cold across the whole cloud. Investigating the abundance of 13CO across G214.5, we find that there is a significantly lower abundance along the entire 13 pc spine of the filament, extending out to a radius of ∼0.8 pc, corresponding to Av ≳ 2 mag and Tdust ≲ 13.5 K. Due to this, we attribute the decrease in abundance to CO freeze-out, making G214.5 the largest scale example of freeze-out yet. We construct an axisymmetric model of G214.5's 13CO volume density considering freeze-out and find that to reproduce the observed profile significant depletion is required beginning at low volume densities, n ≳ 2000 cm−3. Freeze-out at this low number density is possible only if the cosmic-ray ionization rate is ∼1.9 × 10−18 s−1, an order of magnitude below the typical value. Using time scale arguments, we posit that such a low ionization rate may lead to ambipolar diffusion being an important physical process along G214.5's entire spine. We suggest that if low cosmic-ray ionization rates are more common in the outer Galaxy, and other quiescent regions, cloud-scale CO freeze-out occurring at low column and number densities may also be more prevalent, having consequences for CO observations and their interpretation. [ABSTRACT FROM AUTHOR]

Star formation takes place in filamentary molecular clouds which arise by physical processes that take place in the cold neutral medium (CNM). We address the necessary conditions for this diffuse (n ≈ 30 cm−3), cold (T ≈ 60 K), magnetized gas undergoing shock waves, and supersonic turbulence, to produce filamentary structures capable of fragmenting into cluster forming regions. Using ramses and a magnetized CNM environment as our initial conditions, we simulate a 0.5 kpc turbulent box to model a uniform gas with magnetic field strength of 7 μG, varying the 3D velocity dispersion via decaying turbulence. We use a surface density of 320 M⊙ pc−2, representative of the inner 4.0 kpc central molecular zone of the Milky Way and typical luminous galaxies. Filamentary molecular clouds are formed dynamically via shocks within a narrow range of velocity dispersions in the CNM of 5–10 km s−1 with a preferred value at 8 km s−1. Cluster sink particles appear in filaments which exceed their critical line mass, occurring optimally for velocity dispersions of 8 km s−1. Tracking the evolution of magnetic fields, we find that they lead to double the dense star-forming gas than in purely hydro runs. Perpendicular orientations between magnetic field and filaments can increase the accretion rates onto filaments and hence their line masses. Because magnetic fields help support gas, magnetohydrodynamic runs result in average temperatures an order of magnitude higher than unmagnetized counterparts. Finally, we find magnetic fields delay the onset of cluster formation by ∝ 0.4 Myr. [ABSTRACT FROM AUTHOR]

Observing 3D magnetic fields, including orientation and strength, within the interstellar medium is vital but notoriously difficult. However, recent advances in our understanding of anisotropic magnetohydrodynamic (MHD) turbulence demonstrate that MHD turbulence and 3D magnetic fields leave their imprints on the intensity features of spectroscopic observations. Leveraging these theoretical frameworks, we propose a novel Convolutional Neural Network (CNN) model to extract this embedded information, enabling the probe of 3D magnetic fields. This model examines the plane-of-the-sky magnetic field orientation (ϕ), the magnetic field's inclination angle (γ) relative to the line-of-sight, and the total magnetization level (M |$_{\rm A}^{-1}$|⁠) of the cloud. We train the model using synthetic emission lines of 13CO (J  = 1–0) and C18O (J  = 1–0), generated from 3D MHD simulations that span conditions from sub-Alfvénic to super-Alfvénic molecular clouds. Our tests confirm that the CNN model effectively reconstructs the 3D magnetic field topology and magnetization. The median uncertainties are under 5° for both ϕ and γ, and less than 0.2 for M A in sub-Alfvénic conditions (M A ≈ 0.5). In super-Alfvénic scenarios (M A ≈ 2.0), they are under 15° for ϕ and γ, and 1.5 for M A. We applied this trained CNN model to the L1478 molecular cloud. Results show a strong agreement between the CNN-predicted magnetic field orientation and that derived from Planck 353 GHz polarization. The CNN approach enabled us to construct the 3D magnetic field map for L1478, revealing a global inclination angle of ≈76° and a global M A of ≈1.07. [ABSTRACT FROM AUTHOR]

We propose to employ emission line luminosities obtained via optical spectroscopy to estimate the content of neutral hydrogen (H  i) in galaxies. We use the optical spectroscopy data from the Mapping of Nearby Galaxies at APO (MaNGA) survey released in the frames of public DR17 of the Sloan Digital Sky Survey (SDSS). We compare the H  i mass estimated by us for a large sample of SDSS/MaNGA galaxies with direct H  i measurements from the Arecibo Legacy Fast ALFA survey and find a tight correlation between the masses with the correlation coefficient (CC) of 0.91 and the rms scatter of 0.15 dex for the logarithmic mass. The obtained relationship is verified via another sample of MaNGA galaxies with H  i masses measured with the Green Bank Telescope. Despite the coarser angular resolution of the radio data, the relation between the estimated and measured directly H  i mass is tight as well – in this case CC = 0.74 and the rms is 0.29 dex. The established relations allow us to estimate the total mass of neutral hydrogen as well as the spatial distribution of H  i surface density in galaxies from optical spectroscopic observations only in a simple and efficient way. [ABSTRACT FROM AUTHOR]

ABSTRACT We present the stability analysis of two regions, OMC-3 and OMC-4, in the massive and long molecular cloud complex of Orion A. We obtained 214 $\mu$m HAWC + /SOFIA polarization data, and we make use of archival data for the column density and C18O (1–0) emission line. We find clear depolarization in both observed regions and that the polarization fraction is anticorrelated with the column density and the polarization-angle dispersion function. We find that the filamentary cloud and dense clumps in OMC-3 are magnetically supercritical and strongly subvirial. This region should be in the gravitational collapse phase and is consistent with many young stellar objects (YSOs) forming in the region. Our histogram of relative orientation (HRO) analysis shows that the magnetic field is dynamically sub-dominant in the dense gas structures of OMC-3. We present the first polarization map of OMC-4. We find that the observed region is generally magnetically subcritical except for an elongated dense core, which could be a result of projection effect of a filamentary structure aligned close to the line of sight. The relative large velocity dispersion and the unusual positive shape parameters at high column densities in the HROs analysis suggest that our viewing angle may be close to axes of filamentary substructures in OMC-4. The dominating strong magnetic field in OMC-4 is unfavourable for star formation and is consistent with much fewer YSOs than in OMC-3.

Optical and infrared polarization mapping and recent Planck observations of the filametary cloud L1495 in Taurus show that the large-scale magnetic field is approximately perpendicular to the long axis of the cloud. We use the HAWC + polarimeter on SOFIA to probe the complex magnetic field in the B211 part of the cloud. Our results reveal a dispersion of polarization angles of 36 °, about five times that measured on a larger scale by Planck. Applying the Davis-Chandrasekhar-Fermi (DCF) method with velocity information obtained from Institut de Radioastronomie Millimétrique 30 m C 18 O(1-0) observations, we find two distinct sub-regions with magnetic field strengths differing by more than a factor 3. The quieter sub-region is magnetically critical and sub-Alfvé nic; the field is comparable to the average field measured in molecular clumps based on Zeeman observations. The more chaotic, super-Alfvé nic sub-region shows at least three velocity components, indicating interaction among multiple substructures. Its field is much less than the average Zeeman field in molecular clumps, suggesting that the DCF value of the field there may be an underestimate. Numerical simulation of filamentary cloud formation shows that filamentary substructures can strongly perturb the magnetic field. DCF and true field values in the simulation are compared. Pre-stellar cores are observed in B211 and are seen in our simulation. The appendices give a deri v ation of the standard DCF method that allows for a dispersion in polarization angles that is not small, present an alternate deri v ation of the structure function version of the DCF method, and treat fragmentation of filaments.

We present a new physically motivated model for estimating the molecular line emission in active galaxies. The model takes into account (i) the internal density structure of giant molecular clouds (GMCs), (ii) the heating associated with both stars and the active galactic nuclei (AGNs), respectively, producing photodissociation regions (PDRs) and X-ray-dominated regions (XDRs) within the GMCs, and (iii) the mass distribution of GMCs within the galaxy volume. The model needs, as input parameters, the radial profiles of molecular mass, far-UV flux and X-ray flux for a given galaxy, and it has two free parameters: the CO-to-H2 conversion factor αCO, and the X-ray attenuation column density N H. We test this model on a sample of 24 local (z ≤ 0.06) AGN-host galaxies, simulating their carbon monoxide spectral-line energy distribution (CO SLED). We compare the results with the available observations and calculate, for each galaxy, the best (αCO, N H) with a Markov chain Monte Carlo algorithm, finding values consistent with those present in the literature. We find a median αCO = 4.8 M⊙ (K km s−1 pc2)−1 for our sample. In all the modelled galaxies, we find the XDR component of the CO SLED to dominate the CO luminosity from J upp ≥ 4. We conclude that, once a detailed distribution of molecular gas density is taken into account, PDR emission at mid-/high- J becomes negligible with respect to XDR. [ABSTRACT FROM AUTHOR]