Recent results from the James Webb Space Telescope show that nearby spiral galaxies are dominated by the presence of H I and H II bubbles that strongly shape their surrounding medium. These bubbles result from the feedback of high-mass stars at different stages of their life cycle. However, early (pre-supernova) feedback from high-mass stars is still poorly quantified. Recent results from numerical simulations suggest that the impact of high-mass star early feedback (photoionization, wind) on star formation properties is complex, time-dependent, and strongly depends on physical conditions, including the magnetic field properties. In our Galaxy, ionized (H II) regions observed in different evolution stages show a high diversity of star formation in their associated photo-dissociation regions (PDRs). However, the way in which the low- to high-density interstellar medium evolves to this situation remains elusive. Quantifying the impact of early feedback from high-mass stars on star formation properties and star formation laws (star formation rate, star formation efficiency versus gas surface density, Sigma(gas)) will allow for a better understanding of the evolution of star formation laws in external galaxies, the laws that are key ingredients of galaxy evolution models. PRIMA, with its high sensitivity, large mapping efficiency, and polarimetric capabilities, offers a unique opportunity to address the way radiative feedback and magnetic field control star formation in the Milky Way. (c) 2025 Society of Photo-Optical Instrumentation Engineers (SPIE)

The stellar initial mass function (sIMF) describes the distribution of stellar masses formed in a single star formation event in a molecular cloud clump. It is fundamental to astrophysics and cosmology, shaping our understanding of unresolved stellar populations, galactic chemical enrichment and habitable zones, and black hole growth. This White Paper reviews studies on the core mass function, stellar multiplicity, and dynamical processes affecting sIMF determinations, as well as the link between star-forming clumps and the galaxy-wide IMF (gIMF). The evidence gleaned from observed systems for the dependency of the sIMF on the metallicity and density of the clump is portrayed. We examine evidence from gravitational lensing, stellar and gas kinematics, and spectral diagnostics to assess environmental dependencies of the gIMF. Theoretical perspectives provide further insights into the sIMF's variability. Beyond summarizing current knowledge, this work aims to establish a shared framework and define strategies for studying a variable IMF in the era of near-infrared integral-field spectroscopy, 30m-class telescopes and major space-based observatories.

Dense filaments are believed to be representative of the initial conditions of star formation in molecular clouds. We have used the MIRI instrument on JWST to image the massive filament NGC 6334M at d ∼ 1.3 kpc with unprecedented resolution and dynamic range at 7.7 and 25.5 μ m. Our observations reveal the fine structure of the filament in absorption against mid-infrared background emission. From the absorption data, we derive high-resolution column density maps and perform a detailed analysis of the structure of NGC 6334M. We find a median filament width of 0.12 ± 0.02 pc at both wavelengths, resolved by almost two orders of magnitude with MIRI, and consistent with the typical half-power width of Herschel filaments in nearby ( d < 500 pc) clouds. The JWST data also reveal the presence of a quasi-periodic series of side filaments with a similar projected spacing of 0.125 ± 0.015 pc. Combining our JWST results with Spitzer and APEX/Herschel data, we perform a study of cloud structure over four orders of magnitude in linear scale. A convergence test shows that our width estimates for NGC 6334M are robust and reflect the presence of a true characteristic scale. While there is evidence of a Kolmogorov-like spectrum of small-scale fluctuations down to the 1.6 × 10 ^−3 pc resolution of the JWST observations, we identify a break in the power spectrum of column density fluctuations at a scale ∼0.1–0.4 pc comparable to the width of NGC 6334M and its side filaments. This characteristic scale ∼0.1 pc has important implications for the origin of the star formation efficiency in dense gas and the initial mass function.

We report 850 $\mu$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from north to east. Field strengths estimates and a virial analysis for the major clumps indicate that NGC 2264C is globally dominated by gravity while in 2264D magnetic, gravitational, and kinetic energies are roughly balanced. We present an analysis scheme that utilizes the locally resolved magnetic field structures, together with the locally measured gravitational vector field and the extracted filamentary network. From this, we infer statistical trends showing that this network consists of two main groups of filaments oriented approximately perpendicular to one another. Additionally, gravity shows one dominating converging direction that is roughly perpendicular to one of the filament orientations, which is suggestive of mass accretion along this direction. Beyond these statistical trends, we identify two types of filaments. The type-I filament is perpendicular to the magnetic field with local gravity transitioning from parallel to perpendicular to the magnetic field from the outside to the filament ridge. The type-II filament is parallel to the magnetic field and local gravity. We interpret these two types of filaments as originating from the competition between radial collapsing, driven by filament self-gravity, and the longitudinal collapsing, driven by the region's global gravity.

As part of the B-fields In Star-forming Region Observations survey, we present James Clerk Maxwell Telescope (JCMT) 850 mu m polarimetric observations toward the Orion integral-shaped filament (ISF) that covers three portions known as OMC-1, OMC-2, and OMC-3. The magnetic field threading the ISF seen in the JCMT POL-2 map appears as a tale of three: pinched for OMC-1, twisted for OMC-2, and nearly uniform for OMC-3. A multiscale analysis shows that the magnetic field structure in OMC-3 is very consistent at all the scales, whereas the field structure in OMC-2 shows no correlation across different scales. In OMC-1, the field retains its mean orientation from large to small scales but shows some deviations at small scales. Histograms of relative orientations between the magnetic field and filaments reveal a bimodal distribution for OMC-1, a relatively random distribution for OMC-2, and a distribution with a predominant peak at 90 degrees for OMC-3. Furthermore, the magnetic fields in OMC-1 and OMC-3 both appear to be aligned perpendicular to the fibers, which are denser structures within the filament, but the field in OMC-2 is aligned along with the fibers. All these suggest that gravity, turbulence, and magnetic field are each playing a leading role in OMC-1, 2, and 3, respectively. While OMC-2 and 3 have almost the same gas mass, density, and nonthermal velocity dispersion, there are on average younger and fewer young stellar objects in OMC-3, providing evidence that a stronger magnetic field will induce slower and less efficient star formation in molecular clouds.

Context. Recent studies suggest that filamentary structures are representative of the initial conditions of star formation in molecular clouds and support a filament paradigm for star formation, potentially accounting for the origin of the stellar initial mass function (IMF). The detailed, local physical properties of molecular filaments remain poorly characterized, however. Aims. Using Herschel imaging observations of the California giant molecular cloud, we aim to further investigate the filament paradigm for low- to intermediate-mass star formation and to better understand the exact role of filaments in the origin of stellar masses. Methods. Using the multiscale, multiwavelength extraction method getsf, we identify starless cores, protostars, and filaments in the Herschel data set and separate these components from the background cloud contribution to determine accurate core and filament properties. Results. We find that filamentary structures contribute approximately 20% of the overall mass of the California cloud, while compact sources such as dense cores contribute a mere 2% of the total mass. Considering only dense gas (defined as gas with A(V,bg) > 4.5-7), filaments and cores contribute similar to 66-73% and 10-14% of the dense gas mass, respectively. The transverse half-power diameter measured for California molecular filaments has a median undeconvolved value of 0.18 pc, consistent within a factor of 2 with the typical similar to 0.1 pc width of nearby filaments from the Herschel Gould Belt survey. A vast majority of identified prestellar cores (similar to 82-90%) are located within similar to 0.1 pc of the spines of supercritical filamentary structures. Both the prestellar core mass function (CMF) and the distribution of filament masses per unit length or filament line mass function (FLMF) are consistent with power-law distributions at the high-mass end, Delta N/Delta logM proportional to M-1.4 +/- 0.2 at M > 1 M-circle dot for the CMF and Delta N/Delta log M-line proportional to M-line(-1.5 +/- 0.2) for the FLMF at M-line > 10 M-circle dot pc(-1), which are both consistent with the Salpeter power-law IMF. Based on these results, we propose a revised model for the origin of the CMF in filaments, whereby the global prestellar CMF in a molecular cloud arises from the integration of the CMFs generated by individual thermally supercritical filaments within the cloud. Conclusions. Our findings support the existence a tight connection between the FLMF and the CMF/IMF and suggests that filamentary structures represent a critical evolutionary step in establishing a Salpeter-like mass function.

There is growing observational evidence that giant planet formation happens early, within a million years of the coalescence of the protoplanetary disk. Ionization rate is one of the most important parameters controlling both the chemical and dynamical processes in these disks. What few observational constrains on ionization currently exists suggest overall low ionization, limiting the processes able to take place. This is seemingly in conflict with chemical models which demonstrate the importance of ionization for the chemical processing of volatile carbon and observations which suggest such processing is ubiquitous and happens quickly. I will present new NOEMA observations which, when combined with chemical modeling, are indicative of enhanced ionization rates in the envelopes of three Class I protostars. I will then discuss the potential impact of this early enhancement on the chemical composition of the material available to forming planets.

We present 850 mu m polarization observations of the IC 348 star-forming region in the Perseus molecular cloud as part of the B-fields In STar-forming Region Observation survey. We study the magnetic properties of two cores (HH 211 MMS and IC 348 MMS) and a filamentary structure of IC 348. We find that the overall field tends to be more perpendicular than parallel to the filamentary structure of the region. The polarization fraction decreases with intensity, and we estimate the trend by power law and the mean of the Rice distribution fittings. The power indices for the cores are much smaller than 1, indicative of possible grain growth to micron size in the cores. We also measure the magnetic field strengths of the two cores and the filamentary area separately by applying the Davis-Chandrasekhar-Fermi method and its alternative version for compressed medium. The estimated mass-to-flux ratios are 0.45-2.20 and 0.63-2.76 for HH 211 MMS and IC 348 MMS, respectively, while the ratios for the filament are 0.33-1.50. This result may suggest that the transition from subcritical to supercritical conditions occurs at the core scale (similar to 0.05 pc) in the region. In addition, we study the energy balance of the cores and find that the relative strength of turbulence to the magnetic field tends to be stronger for IC 348 MMS than for HH 211 MMS. The result could potentially explain the different configurations inside the two cores: a single protostellar system in HH 211 MMS and multiple protostars in IC 348 MMS.

We study the H ii regions associated with the NGC 6334 molecular cloud observed in the submillimeter and taken as part of the B -fields In STar-forming Region Observations Survey. In particular, we investigate the polarization patterns and magnetic field morphologies associated with these H ii regions. Through polarization pattern and pressure calculation analyses, several of these bubbles indicate that the gas and magnetic field lines have been pushed away from the bubble, toward an almost tangential (to the bubble) magnetic field morphology. In the densest part of NGC 6334, where the magnetic field morphology is similar to an hourglass, the polarization observations do not exhibit observable impact from H ii regions. We detect two nested radial polarization patterns in a bubble to the south of NGC 6334 that correspond to the previously observed bipolar structure in this bubble. Finally, using the results of this study, we present steps (incorporating computer vision; circular Hough transform) that can be used in future studies to identify bubbles that have physically impacted magnetic field lines.

The total mass estimate of molecular clouds suffers from the uncertainty in the H$_2$-CO conversion factor, the so-called $X_{\rm CO}$ factor, which is used to convert the $^{12}$CO (1--0) integrated intensity to the H$_2$ column density. We demonstrate the machine learning's ability to predict the H$_2$ column density from the $^{12}$CO, $^{13}$CO, and C$^{18}$O (1--0) data set of four star-forming molecular clouds; Orion A, Orion B, Aquila, and M17. When the training is performed on a subset of each cloud, the overall distribution of the predicted column density is consistent with that of the Herschel column density. The total column density predicted and observed is consistent within 10\%, suggesting that the machine learning prediction provides a reasonable total mass estimate of each cloud. However, the distribution of the column density for values $> \sim 2 \times 10^{22}$ cm$^{-2}$, which corresponds to the dense gas, could not be predicted well. This indicates that molecular line observations tracing the dense gas are required for the training. We also found a significant difference between the predicted and observed column density when we created the model after training the data on different clouds. This highlights the presence of different $X_{\rm CO}$ factors between the clouds, and further training in various clouds is required to correct for these variations. We also demonstrated that this method could predict the column density toward the area not observed by Herschel if the molecular line and column density maps are available for the small portion, and the molecular line data are available for the larger areas.

We present observations of polarized dust emission at 850 mu m from the L43 molecular cloud, which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense (N-H2 similar to 10(22) - 10(23) cm(-2)) complex molecular cloud with a submillimeter-bright starless core and two protostellar sources. There appears to be an evolutionary gradient along the isolated filament that L43 is embedded within, with the most evolved source closest to the Sco OB2 association. One of the protostars drives a CO outflow that has created a cavity to the southeast. We see a magnetic field that appears to be aligned with the cavity walls of the outflow, suggesting interaction with the outflow. We also find a magnetic field strength of up to similar to 160 +/- 30 mu G in the main starless core and up to similar to 90 +/- 40 mu G in the more diffuse, extended region. These field strengths give magnetically super- and subcritical values, respectively, and both are found to be roughly trans-Alfvenic. We also present a new method of data reduction for these denser but fainter objects like starless cores.

The study of the development of structures on multiple scales in the cold interstellar medium has experienced rapid expansion in the past decade, on both the observational and the theoretical front. Spectral line studies at (sub-)millimeter wavelengths over a wide range of physical scales have provided unique probes of the kinematics of dense gas in star-forming regions, and have been complemented by extensive, high dynamic range dust continuum surveys of the column density structure of molecular cloud complexes, while dust polarization maps have highlighted the role of magnetic fields. This has been accompanied by increasingly sophisticated numerical simulations including new physics (e.g., supernova driving, cosmic rays, non-ideal magneto-hydrodynamics, radiation pressure) and new techniques such as zoom-in simulations allowing multi-scale studies. Taken together, these new data have emphasized the anisotropic growth of dense structures on all scales, from giant ISM bubbles driven by stellar feedback on $\sim$50-100 pc scales through parsec-scale molecular filaments down to $<$0.1 pc dense cores and $<$1000 au protostellar disks. Combining observations and theory, we present a coherent picture for the formation and evolution of these structures and synthesize a comprehensive physical scenario for the initial conditions and early stages of star and disk formation.

We present BISTRO Survey 850 μm dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to the filaments in almost all cases. We also find that the large-scale magnetic field orientation measured by Planck is not correlated with any of the core or filament structures, except in the case of the lowest-density core. We propose a scenario for early prestellar evolution that is both an extension to, and consistent with, previous models, introducing an additional evolutionary transitional stage between field-dominated and matter-dominated evolution, observed here for the first time. In this scenario, the cloud collapses first to a sheet-like structure. Uniquely, we appear to be seeing this sheet almost face-on. The sheet fragments into filaments, which in turn form cores. However, the material must reach a certain critical density before the evolution changes from being field-dominated to being matter-dominated. We measure the sheet surface density and the magnetic field strength at that transition for the first time and show consistency with an analytical prediction that had previously gone untested for over 50 years (Mestel 1965).

[Abridged] The interstellar medium is observed to be organised in filamentary structures, as well as neutral (HI) and ionized (HII) bubbles. The expanding nature of these bubbles makes them shape their surroundings and possibly play a role in the formation and evolution of interstellar filaments. We present APEX ^13CO and C^18O(2-1) observations of the NGC 6334 molecular cloud. We investigate the gas velocity structure along and across the 50 pc-long cloud and towards 75 identified velocity-coherent-filaments (VCFs). We measure a wealth of velocity gradients along the VCFs. We derive the column density and velocity power spectra of the VCFs. These power spectra are well represented with power laws showing similar slopes for both quantities (with a mean of about -2), albeit some differ by up to a factor of two. The position velocity diagrams perpendicular to three VCFs show the V-shaped velocity pattern, corresponding to a bent structure in velocity space with the filament at the tip of the V surrounded by an extended structure connected to it with a velocity gradient. This velocity structure is qualitatively similar to that resulting from numerical simulations of filament formation from large-scale compression from propagating shock fronts. In addition, the radial profiles perpendicular to these VCFs hint to small-scale internal impacts from neighbouring HII bubbles. The observed opposite curvature in velocity space towards the VCFs points to various origins of large-scale external compressions from propagating HI bubbles. This suggests the plausible importance of multiple HI compressions, separated in space and time, in the formation and evolution of molecular clouds and their star formation history. These latter atomic compressions due to past and distant star formation events are complemented by the impact of HII bubbles from present time and local star formation activity.

We present 850 μm polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope as part of the B-fields In STar-forming Region Observations survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less-dense filamentary structures where NH2<0.93×1022 cm−2 (magnetic fields perpendicular to density gradients), while they are perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at NH2≈4.6×1022 cm−2. This can be interpreted as a signature of core formation. At NH2≈16×1022 cm−2, magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments (B POS = 60–300 μG)) using the Davis–Chandrasekhar–Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities.

ABSTRACT 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 C18O(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$\acute{\rm e}$nic; the field is comparable to the average field measured in molecular clumps based on Zeeman observations. The more chaotic, super-Alfv$\acute{\rm e}$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 derivation of the standard DCF method that allows for a dispersion in polarization angles that is not small, present an alternate derivation of the structure function version of the DCF method, and treat fragmentation of filaments.

We present and analyze observations of polarized dust emission at 850 μ m toward the central 1 × 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields in Star-forming Region Observations survey. The orientations of the magnetic field follow the spiral structure of Mon R2, which are well described by an axisymmetric magnetic field model. We estimate the turbulent component of the magnetic field using the angle difference between our observations and the best-fit model of the underlying large-scale mean magnetic field. This estimate is used to calculate the magnetic field strength using the Davis–Chandrasekhar–Fermi method, for which we also obtain the distribution of volume density and velocity dispersion using a column density map derived from Herschel data and the C 18 O ( J = 3 − 2) data taken with HARP on the JCMT, respectively. We make maps of magnetic field strengths and mass-to-flux ratios, finding that magnetic field strengths vary from 0.02 to 3.64 mG with a mean value of 1.0 ± 0.06 mG, and the mean critical mass-to-flux ratio is 0.47 ± 0.02. Additionally, the mean Alfvén Mach number is 0.35 ± 0.01. This suggests that, in Mon R2, the magnetic fields provide resistance against large-scale gravitational collapse, and the magnetic pressure exceeds the turbulent pressure. We also investigate the properties of each filament in Mon R2. Most of the filaments are aligned along the magnetic field direction and are magnetically subcritical.

We present 850 μ m dust polarization observations of the massive DR21 filament from the B -fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the subfilaments, the magnetic fields are mainly parallel to the filamentary structures and smoothly connect to the magnetic fields of the main filament. We compare the POL-2 and Planck dust polarization observations to study the magnetic field structures of the DR21 filament on 0.1–10 pc scales. The magnetic fields revealed in the Planck data are well-aligned with those of the POL-2 data, indicating a smooth variation of magnetic fields from large to small scales. The plane-of-sky magnetic field strengths derived from angular dispersion functions of dust polarization are 0.6–1.0 mG in the DR21 filament and ∼0.1 mG in the surrounding ambient gas. The mass-to-flux ratios are found to be magnetically supercritical in the filament and slightly subcritical to nearly critical in the ambient gas. The alignment between column density structures and magnetic fields changes from random alignment in the low-density ambient gas probed by Planck to mostly perpendicular in the high-density main filament probed by James Clerk Maxwell Telescope. The magnetic field structures of the DR21 filament are in agreement with MHD simulations of a strongly magnetized medium, suggesting that magnetic fields play an important role in shaping the DR21 main filament and subfilaments.

We report the first high spatial resolution measurement of magnetic fields surrounding LkH α 101, part of the Auriga–California molecular cloud. The observations were taken with the POL-2 polarimeter on the James Clerk Maxwell Telescope within the framework of the B-fields In Star-forming Region Observations (BISTRO) survey. Observed polarization of thermal dust emission at 850 μ m is found to be mostly associated with the redshifted gas component of the cloud. The magnetic field displays a relatively complex morphology. Two variants of the Davis–Chandrasekhar–Fermi method, unsharp masking and structure function, are used to calculate the strength of magnetic fields in the plane of the sky, yielding a similar result of B POS ∼ 115 μ G. The mass-to-magnetic-flux ratio in critical value units, λ ∼ 0.3, is the smallest among the values obtained for other regions surveyed by POL-2. This implies that the LkH α 101 region is subcritical, and the magnetic field is strong enough to prevent gravitational collapse. The inferred δ B / B 0 ∼ 0.3 implies that the large-scale component of the magnetic field dominates the turbulent one. The variation of the polarization fraction with total emission intensity can be fitted by a power law with an index of α = 0.82 ± 0.03, which lies in the range previously reported for molecular clouds. We find that the polarization fraction decreases rapidly with proximity to the only early B star (LkH α 101) in the region. Magnetic field tangling and the joint effect of grain alignment and rotational disruption by radiative torques can potentially explain such a decreasing trend.