Ionization Feedback Research Articles

The Impact of Stellar Radiative Feedback on Formation of Young Massive Clusters via Fast H i Gas Collisions

Young massive clusters (YMCs) are dense aggregates of young stars and are often speculated as potential precursors to globular clusters. However, the formation mechanism of massive and compact gas clumps that precede YMCs remains unknown. In this paper, we study the formation of such massive clumps via fast H i gas collisions (∼100 km s−1) as suggested by recent observations and their subsequent evolution into YMCs by using three-dimensional magnetohydrodynamics simulations involving self-gravity and detailed thermal/chemical processes. In particular, the impact of ionization feedback from stellar radiation is included in an approximate fashion where the temperature within the H ii regions is elevated to 10,000 K, while supernova feedback is not included. We examine whether the resulting massive clumps can survive this ionization feedback and evolve into YMCs. Our simulations reveal the emergence of gas clumps that not only possess substantial mass (∼105 M ⊙) but also sufficient compactness (∼5 pc). Notably, these clumps exhibit significantly higher escape velocities compared to the sound speed of the H ii region, indicating effective gravitational retention of gas against feedback-induced evaporation. Consequently, these conditions foster efficient star formation within the massive gas clumps, ultimately leading to their evolution into YMCs. We also perform simulations involving lower-velocity gas collisions, approximately 15 km s−1, typical shock velocities induced by galactic superbubbles. In contrast to the high-velocity collisions, we find that molecular cloud formation does not occur in the case of 1 cm−3 gas collision, while YMC formation is observed in the presence of denser gas of 10 cm−3. However, the formation of YMCs requires compression periods exceeding 10 Myr in these cases, indicating a potential preference for gas collisions driven by intergalactic interactions rather than galactic superbubbles for YMC formation.

Simulating ionization feedback from young massive stars: impact of numerical resolution

ABSTRACT Modelling galaxy formation in hydrodynamic simulations has increasingly adopted various radiative transfer methods to account for photoionization feedback from young massive stars. However, the evolution of H ii regions around stars begins in dense star-forming clouds and spans large dynamical ranges in both space and time, posing severe challenges for numerical simulations in terms of both spatial and temporal resolution that depends strongly on gas density (∝n−1). In this work, we perform a series of idealized H ii region simulations using the moving-mesh radiation-hydrodynamic code arepo-rt to study the effects of numerical resolution. The simulated results match the analytical solutions and the ionization feedback converges only if the Strömgren sphere is resolved by at least 10–100 resolution elements and the size of each time integration step is smaller than 0.1 times the recombination time-scale. Insufficient spatial resolution leads to reduced ionization fraction but enhanced ionized gas mass and momentum feedback from the H ii regions, as well as degrading the multiphase interstellar medium into a diffuse, partially ionized, warm (∼8000 K) gas. On the other hand, insufficient temporal resolution strongly suppresses the effects of ionizing feedback. This is because longer time-steps are not able to resolve the rapid variation of the thermochemistry properties of the gas cells around massive stars, especially when the photon injection and thermochemistry are performed with different cadences. Finally, we provide novel numerical implementations to overcome the above issues when strict resolution requirements are not achievable in practice.

ATOMS: ALMA three-millimeter observations of massive star-forming regions – XIII. Ongoing triggered star formation within clump-fed scenario found in the massive (∼1500 M⨀) clump

ABSTRACT Whether ionization feedback triggers the formation of massive stars is highly debated. Using ALMA 3-mm observations with a spatial resolution of ∼0.05 pc and a mass sensitivity of 1.1 $\rm M_\odot$ per beam at 20 K, we investigate the star formation and gas flow structures within the ionizing feedback-driven structure, a clump-scale massive (≳ 1500 $\rm M_\odot$) bright-rimmed cloud (BRC) associated with IRAS 18290–0924. This BRC is bound only if external compression from ionized gas is considered. A small-scale (≲ 1 pc) age sequence along the direction of ionizing radiation is revealed for the embedded cores and protostars, which suggests triggered star formation via radiation-driven implosion (RDI). Furthermore, filamentary gas structures converge towards the cores located in the BRC’s centre, indicating that these filaments are fueling mass towards cores. The local core-scale mass infall rate derived from H13CO+ J = 1 − 0 blue profile is of the same order of magnitude as the filamentary mass inflow rate, approximately 1 $\rm M_\odot$ kyr−1. A photodissociation region (PDR) covering the irradiated clump surface is detected in several molecules, such as CCH, HCO+, and CS whereas the spatial distribution stratification of these molecules is indistinct. CCH spectra of the PDR possibly indicate a photoevaporation flow leaving the clump surface with a projected velocity of ∼2 km s−1. Our new observations show that RDI accompanied by a clump-fed process is operating in this massive BRC. Whether this combined process works in other massive BRCs is worth exploring with dedicated surveys.

Non-ideal magnetohydrodynamic simulations of the first star formation: the effect of ambipolar diffusion

ABSTRACT In the present-day universe, magnetic fields play such essential roles in star formation as angular momentum transport and outflow driving, which control circumstellar disc formation/fragmentation and also the star formation efficiency. While only a much weaker field has been believed to exist in the early universe, recent theoretical studies find that strong fields can be generated by turbulent dynamo during the gravitational collapse. Here, we investigate the gravitational collapse of a cloud core ($\sim 10^{3}\ \rm cm^{-3}$) up to protostar formation ($\sim 10^{20}\ \rm cm^{-3}$) by non-ideal magnetohydrodynamic simulations considering ambipolar diffusion (AD), the dominant non-ideal effects in the primordial-gas. We systematically study rotating cloud cores either with or without turbulence and permeated with uniform fields of different strengths. We find that AD can slightly suppress the field growth by dynamo especially on scales smaller than the Jeans-scale at the density range $10^{10}\!-\!10^{14}\ \rm cm^{-3}$, while we could not see the AD effect on the temperature evolution, since the AD heating rate is always smaller than compression heating. The inefficiency of AD makes the field as strong as $10^{3}\!-\!10^{5} \rm \ G$ near the formed protostar, much stronger than in the present-day cases, even in cases with initially weak fields. The magnetic field affects the inflow motion when amplified to the equipartition level with turbulence on the Jeans-scale, although disturbed fields do not launch winds. This might suggest that dynamo amplified fields have smaller impact on the dynamics in the later accretion phase than other processes such as ionization feedback.

The Role of Magnetic Fields in Triggered Star Formation of RCW 120

We report on the near-infrared polarimetric observations of RCW 120 with the 1.4 m IRSF telescope. The starlight polarization of the background stars reveals for the first time the magnetic field of RCW 120. The global magnetic field of RCW 120 is along the direction of 20°, parallel to the Galactic plane. The field strength on the plane of the sky is 100 ± 26 μG. The magnetic field around the eastern shell shows evidence of compression by the H ii region. The external pressure (turbulent pressure + magnetic pressure) and the gas density of the ambient cloud are minimum along the direction where RCW 120 breaks out, which explains the observed elongation of RCW 120. The dynamical age of RCW 120, depending on the magnetic field strength, is ∼1.6 Myr for field strength of 100 μG, older than the hydrodynamic estimates. In direction perpendicular to the magnetic field, the density contrast of the western shell is greatly reduced by the strong magnetic field. The strong magnetic field in general reduces the efficiency of triggered star formation, in comparison with the hydrodynamic estimates. Triggered star formation via the “collect and collapse” mechanism could occur in the direction along the magnetic field. Core formation efficiency (CFE) is found to be higher in the southern and eastern shells of RCW 120 than in the infrared dark cloud receiving little influence from the H ii region, suggesting increase in the CFE related to triggering from ionization feedback.

Breaking Orion's Veil with fossil outflows

Context. The role of feedback in the self-regulation of star formation is a fundamental question in astrophysics. The Orion Nebula is the nearest site of ongoing and recent massive star formation. It is a unique laboratory for the study of stellar feedback. Recent SOFIA [C II] 158 μm observations have revealed an expanding bubble, the Veil shell, that is powered by stellar winds and ionization feedback. Aims. We have identified a protrusion-like substructure in the northwestern portion of the Orion Veil shell that may indicate additional, highly directional feedback mechanisms. Our goal is to investigate the origin of the protrusion by quantifying its possible driving mechanisms. Methods. We use the [C II] 158 μm map of the Orion Nebula obtained with the upGREAT instrument on board SOFIA. The spectral and spatial resolution of the observations are 0.3 km s−1 and 16", respectively. The velocity-resolved [C II] observations allow us to construct position-velocity (pv) diagrams to measure the morphology and the expansion velocity of the protrusion. For the morphology, we also use new observations of 12CO and 13CO J = 2-1 (to trace molecular gas), Spitzer 8 μm observations (to trace the far-ultraviolet-illuminated surfaces of photodissociation regions), and Hα observations (to trace ionized gas). For the kinematics, we perform a line-profile analysis of [C II], 13CO, and 12CO at 12 positions, covering the entire protrusion. To quantify the stellar feedback, we estimate the mass of the protrusion by fitting the dust thermal emission. We compare the kinetic energy with the stellar wind of θ1 Ori C and the momentum of the outflows of massive protostars to investigate the driving mechanism of the protrusion. Results. The pv diagrams reveal two half-shells expanding at velocities of +6 km s−1 and +12 km s−1. We find that the protrusion has a diameter of ~1.3 pc with a ~45 M⊙ shell expanding at +12 km s−1 at the northwestern rim of the Veil. The thickness of the expanding shell is ~0.1 pc. Using the mass in the limb-brightened shell and the maximum expansion velocity, we calculate the kinetic energy and the momentum of the protrusion to be ~7 × 1046 erg and 540 M⊙ km s−1, respectively. We consider three possible origins for this protrusion: fossil outflow cavities created by jets and outflows during the protostellar accretion phase, preexisting "clumpiness" in the OMC-1 core, and the stellar wind during the main-sequence phase. Based on the energetics and the morphology, we conclude that the northwestern part of the preexisting cloud was locally perturbed by outflows ejected from massive protostars in the Trapezium cluster. This suggests that the protrusion of the Veil is the result of mechanical rather than radiative feedback. Furthermore, we argue that the location of the protrusion is a suitable place to break Orion's Veil owing to the photo-ablation from the walls of the protrusion. Conclusions. We conclude that the outflows of massive protostars can influence the morphology of the future H II region and even cause breakages in the ionization front. Specifically, the interaction of stellar winds of main-sequence stars with the molecular core preprocessed by the protostellar jet is important.

H II regions and high-mass starless clump candidates

Context. The ionization feedback from H II regions modifies the properties of high-mass starless clumps (HMSCs, of several hundred to a few thousand solar masses with a typical size of 0.1–1 pc), such as dust temperature and turbulence, on the clump scale. The question of whether the presence of H II regions modifies the core-scale (~0.025 pc) fragmentation and star formation in HMSCs remains to be explored. Aims. We aim to investigate the difference of 0.025 pc-scale fragmentation between candidate HMSCs that are strongly impacted by H II regions and less disturbed ones. We also search for evidence of mass shaping and induced star formation in the impacted candidate HMSCs. Methods. Using the ALMA 1.3 mm continuum, with a typical angular resolution of 1.3′′, we imaged eight candidate HMSCs, including four impacted by H II regions and another four situated in the quiet environment. The less-impacted candidate HMSCs are selected on the basis of their similar mass and distance compared to the impacted ones to avoid any possible bias linked to these parameters. We carried out a comparison between the two types of candidate HMSCs. We used multi-wavelength data to analyze the interaction between H II regions and the impacted candidate HMSCs. Results. A total of 51 cores were detected in eight clumps, with three to nine cores for each clump. Within our limited sample, we did not find a clear difference in the ~0.025 pc-scale fragmentation between impacted and non-impacted candidate HMSCs, even though H II regions seem to affect the spatial distribution of the fragmented cores. Both types of candidate HMSCs present a thermal fragmentation with two-level hierarchical features at the clump thermal Jeans length λJ,clumpth and 0.3λJ,clumpth. The ALMA emission morphology of the impacted candidate HMSCs AGAL010.214-00.306 and AGAL018.931-00.029 sheds light on the capacities of H II regions to shape gas and dust in their surroundings and possibly to trigger star formation at ~0.025 pc-scale in candidate HMSCs. Conclusions. The fragmentation at ~0.025 pc scale for both types of candidate HMSCs is likely to be thermal-dominant, meanwhile H II regions probably have the capacity to assist in the formation of dense structures in the impacted candidate HMSCs. Future ALMA imaging surveys covering a large number of impacted candidate HMSCs with high turbulence levels are needed to confirm the trend of fragmentation indicated in this study.

The role of Galactic H II regions in the formation of filaments

Context. Massive stars and their associated ionized (H II) regions could play a key role in the formation and evolution of filaments that host star formation. However, the properties of filaments that interact with H II regions are still poorly known. Aims. To investigate the impact of H II regions on the formation of filaments, we imaged the Galactic H II region RCW 120 and its surroundings where active star formation takes place and where the role of ionization feedback on the star formation process has already been studied. Methods. We used the large-format bolometer camera ArTéMiS on the APEX telescope and combined the high-resolution ArTéMiS data at 350 and 450 μm with Herschel-SPIRE/HOBYS data at 350 and 500 μm to ensure good sensitivity to a broad range of spatial scales. This allowed us to study the dense gas distribution around RCW 120 with a resolution of 8′′ or 0.05 pc at a distance of 1.34 kpc. Results. Our study allows us to trace the median radial intensity profile of the dense shell of RCW 120. This profile is asymmetric, indicating a clear compression from the H II region on the inner part of the shell. The profile is observed to be similarly asymmetric on both lateral sides of the shell, indicating a homogeneous compression over the surface. On the contrary, the profile analysis of a radial filament associated with the shell, but located outside of it, reveals a symmetric profile, suggesting that the compression from the ionized region is limited to the dense shell. The mean intensity profile of the internal part of the shell is well fitted by a Plummer-like profile with a deconvolved Gaussian full width at half maximum of 0.09 pc, as observed for filaments in low-mass star-forming regions. Conclusions. Using ArTéMiS data combined with Herschel-SPIRE data, we found evidence for compression from the inner part of the RCW 120 ionized region on the surrounding dense shell. This compression is accompanied with a significant (factor 5) increase of the local column density. This study suggests that compression exerted by H II regions may play a key role in the formation of filaments and may further act on their hosted star formation. ArTéMiS data also suggest that RCW 120 might be a 3D ring, rather than a spherical structure.

H II regions and high-mass starless clump candidates

Context. The role of ionization feedback on high-mass (>8 M⊙) star formation is still highly debated. Questions remain concerning the presence of nearby H II regions changes the properties of early high-mass star formation and whether H II regions promote or inhibit the formation of high-mass stars. Aims. To characterize the role of H II regions on the formation of high-mass stars, we study the properties of a sample of candidates high-mass starless clumps (HMSCs), of which about 90% have masses larger than 100 M⊙. These high-mass objects probably represent the earliest stages of high-mass star formation; we search if (and how) their properties are modified by the presence of an H II region. Methods. We took advantage of the recently published catalog of HMSC candidates. By cross matching the HMSCs and H II regions, we classified HMSCs into three categories: (1) the HMSCs associated with H II regions both in the position in the projected plane of the sky and in velocity; (2) HMSCs associated in the plane of the sky, but not in velocity; and (3) HMSCs far away from any H II regions in the projected sky plane. We carried out comparisons between associated and nonassociated HMSCs based on statistical analyses of multiwavelength data from infrared to radio. Results. We show that there are systematic differences of the properties of HMSCs in different environments. Statistical analyses suggest that HMSCs associated with H II regions are warmer, more luminous, more centrally-peaked and turbulent. We also clearly show, for the first time, that the ratio of bolometric luminosity to envelope mass of HMSCs (L∕M) could not be a reliable evolutionary probe for early massive star formation due to the external heating effects of the H II regions. Conclusions. We show HMSCs associated with H II regions present statistically significant differences from HMSCs far away from H II regions, especially for dust temperature and L∕M. More centrally peaked and turbulent properties of HMSCs associated with H II regions may promote the formation of high-mass stars by limiting fragmentation. High-resolution interferometric surveys toward HMSCs are crucial to reveal how H II regions impact the star formation process inside HMSCs.

Probing star formation and feedback using CCOSMA and archival data in the CFG028.68–0.28 quasi-sinusoidal filament

We have performed a multi-wavelength study toward a quasi-sinusoidal filament (CFG028.68–0.28). A new large-scale 12CO J = 3 – 2 map was obtained from the China-Cologne Observation for SubMillimeter Astronomy (CCOSMA) 3m radio telescope. Based on the ATLASGAL catalog, we have identified 27 dust clumps in the filament. Through the relationship between the mass and radius of these clumps, 67% of these clumps are dense and massive enough to potentially form massive stars. The obtained CFE is ∼11% in the filament. The filament has a linear mass density of ∼305 M⊙ pc−1, which is smaller than its critical mass to length ratio. This suggests that the external pressure from the neighboring H ii regions may help prevent the filament from dispersing under the effects of turbulence. Comparing the energy injection from outflows and H ii regions in the filament, the ionization feedback from the H ii regions can help maintain the observed turbulence.

The effects of ionization feedback on star formation: a case study of the M 16 H II region

Aims. We aim to investigate the impact of the ionized radiation from the M 16 H II region on the surrounding molecular cloud and on its hosted star formation. Methods. To present comprehensive multi-wavelength observations towards the M 16 H II region, we used new CO data and existing infrared, optical, and submillimeter data. The 12CO J = 1−0, 13CO J = 1−0, and C18O J = 1−0 data were obtained with the Purple Mountain Observatory (PMO) 13.7 m radio telescope. To trace massive clumps and extract young stellar objects (YSOs) associated with the M 16 H II region, we used the ATLASGAL and GLIMPSE I catalogs, respectively. Results. From CO data, we discern a large-scale filament with three velocity components. Because these three components overlap with each other in both velocity and space, the filament may be made of three layers. The M 16 ionized gas interacts with the large-scale filament and has reshaped its structure. In the large-scale filament, we find 51 compact cores from the ATLASGAL catalog, 20 of them being quiescent. The mean excitation temperature of these cores is 22.5 K, while this is 22.2 K for the quiescent cores. This high temperature observed for the quiescent cores suggests that the cores may be heated by M 16 and do not experience internal heating from sources in the cores. Through the relationship between the mass and radius of these cores, we obtain that 45% of all the cores are massive enough to potentially form massive stars. Compared with the thermal motion, the turbulence created by the nonthermal motion is responsible for the core formation. For the pillars observed towards M 16, the H II region may give rise to the strong turbulence.

The SOMA Radio Survey. I. Comprehensive SEDs of High-mass Protostars from Infrared to Radio and the Emergence of Ionization Feedback

Abstract We study centimeter continuum emission of eight high- and intermediate-mass protostars that are part of the SOFIA Massive Star Formation Survey, thus building extended spectral energy distributions (SEDs) from the radio to the infrared. We discuss the morphology seen in the centimeter continuum images, which are mostly derived from archival Very Large Array data, and the relation to infrared morphology. We use the SEDs to test new models of high-mass star formation including radiative and disk-wind feedback and associated free–free and dust continuum emission. We show that interferometric data of the centimeter continuum flux densities provide additional, stringent tests of the models by constraining the ionizing luminosity of the source; they also help to break degeneracies encountered when modeling the infrared-only SEDs, especially for the protostellar mass. Our derived parameters are consistent with physical parameters estimated by other methods, such as dynamical protostellar masses. We find a few examples of additional stellar sources in the vicinity of the high-mass protostars, which may be low-mass young stellar objects. However, the stellar multiplicity of the regions, at least as traced by radio continuum emission, appears to be relatively low.

Self-consistent semi-analytic models of the first stars

We have developed a semi-analytic framework to model the large-scale evolution of the first Population III (Pop III) stars and the transition to metal-enriched star formation. Our model follows dark matter halos from cosmological N-body simulations, utilizing their individual merger histories and three-dimensional positions, and applies physically motivated prescriptions for star formation and feedback from Lyman-Werner (LW) radiation, hydrogen ionizing radiation, and external metal enrichment due to supernovae winds. This method is intended to complement analytic studies, which do not include clustering or individual merger histories, and hydrodynamical cosmological simulations, which include detailed physics, but are computationally expensive and have limited dynamic range. Utilizing this technique, we compute the cumulative Pop III and metal-enriched star formation rate density (SFRD) as a function of redshift at $z \geq 20$. We find that varying the model parameters leads to significant qualitative changes in the global star formation history. The Pop III star formation efficiency and the delay time between Pop III and subsequent metal-enriched star formation are found to have the largest impact. The effect of clustering (i.e. including the three-dimensional positions of individual halos) on various feedback mechanisms is also investigated. The impact of clustering on LW and ionization feedback is found to be relatively mild in our fiducial model, but can be larger if external metal enrichment can promote metal-enriched star formation over large distances.

Infrared dust bubble CS51 and its interaction with the surrounding interstellar medium

A multiwavelength investigation of the southern infrared dust bubble CS51 is presented in this paper. We probe the associated ionized, cold dust, molecular and stellar components. Radio continuum emission mapped at 610 and 1300 MHz, using the Giant Metrewave Radio Telescope, India, reveal the presence of three compact emission components (A, B, and C) apart from large-scale diffuse emission within the bubble interior. Radio spectral index map show the coexistence of thermal and non-thermal emission components. Modified blackbody fits to the thermal dust emission using Herschel PACS and SPIRE data is performed to generate dust temperature and column density maps. We identify five dust clumps associated with CS51 with masses and radius in the range 810 - 4600 M{\sun} and 1.0 - 1.9 pc, respectively. We further construct the column density probability distribution functions of the surrounding cold dust which display the impact of ionization feedback from high-mass stars. The estimated dynamical and fragmentation timescales indicate the possibility of collect and collapse mechanism in play at the bubble border. Molecular line emission from the MALT90 survey is used to understand the nature of two clumps which show signatures of expansion of CS51.

The super star cluster driven feedback in ESO338-IG04 and Haro 11

AbstractThe stellar content of young massive star clusters emit large amounts of Lyman continuum photons and inject momentum into the inter stellar medium (ISM) by the strong stellar winds of the most massive stars in the cluster. When the most massive stars explode as supernovae, large amounts of mechanical energy are injected in the ISM. A detailed study of the ISM around these massive cluster provides insights on the effect of cluster feedback.We present high quality integral field spectroscopy taken with VLT/MUSE of two starburst galaxies: ESO 338-IG04 and Haro 11. Both galaxies contain a significant number of super star clusters. The MUSE data provide us with an unprecedented view of the state and kinematics of the ionized gas in the galaxy allowing us to study the effect of stellar feedback on small and large spatial scales. We present our recent results on studying the ISM state of these two galaxies. The data of both galaxies show that the mechanical and ionization feedback of the super star clusters in the galaxy modify the state and kinematics of the ISM substancially by creating highly ionized bubbles around the cluster, making the central part of the galaxy highly ionized. This shows that the HII regions around the individual clusters are density bounded, allowing the ionizing photons to escape and ionize the ISM further out.

VLT/MUSE view of the highly ionized outflow cones in the nearby starburst ESO338-IG04

The Ly$\alpha$ line is an important diagnostic for star formation at high redshift, but interpreting its flux and line profile is difficult because of the resonance nature of Ly$\alpha$. Trends between the escape of Ly$\alpha$ photons and dust and properties of the interstellar medium (ISM) have been found, but detailed comparisons between Ly$\alpha$ emission and the properties of the gas in local high-redshift analogs are vital for understanding the relation between Ly$\alpha$ emission and galaxy properties. For the first time, we can directly infer the properties of the ionized gas at the same location and similar spatial scales of the extended Ly$\alpha$ halo around ESO 338-IG04. We obtained VLT/MUSE integral field spectra. We used ionization parameter mapping of the [SII]/[OIII] line ratio and the kinematics of H$\alpha$ to study the ionization state and kinematics of the ISM of ESO338-IG04. The velocity map reveals two outflows. The entire central area of the galaxy is highly ionized by photons leaking from the HII regions around the youngest star clusters. Three highly ionized cones have been identified, of which one is associated with an outflow. We propose a scenario where the outflows are created by mechanical feedback of the older clusters, while the highly ionized gas is caused by the hard ionizing photons emitted by the youngest clusters. A comparison with the Ly$\alpha$ map shows that the (approximately bipolar) asymmetries observed in the Ly$\alpha$ emission are consistent with the base of the outflows detected in H$\alpha$. No clear correlation with the ionization cones is found. The mechanical and ionization feedback of star clusters significantly changes the state of the ISM by creating ionized cones and outflows. The comparison with Ly$\alpha$ suggests that especially the outflows could facilitate the escape of Ly$\alpha$ photons [Abridged].

Ionization compression impact on dense gas distribution and star formation

Ionization feedback should impact the probability distribution function (PDF) of the column density around the ionized gas. We aim to quantify this effect and discuss its potential link to the Core and Initial Mass Function (CMF/IMF). We used in a systematic way Herschel column density maps of several regions observed within the HOBYS key program: M16, the Rosette and Vela C molecular cloud, and the RCW 120 H ii region. We fitted the column density PDFs of all clouds with two lognormal distributions, since they present a double-peak or enlarged shape in the PDF. Our interpretation is that the lowest part of the column density distribution describes the turbulent molecular gas while the second peak corresponds to a compression zone induced by the expansion of the ionized gas into the turbulent molecular cloud. The condensations at the edge of the ionized gas have a steep compressed radial profile, sometimes recognizable in the flattening of the power-law tail. This could lead to an unambiguous criterion able to disentangle triggered from pre-existing star formation. In the context of the gravo-turbulent scenario for the origin of the CMF/IMF, the double peaked/enlarged shape of the PDF may impact the formation of objects at both the low-mass and the high-mass end of the CMF/IMF. In particular a broader PDF is required by the gravo-turbulent scenario to fit properly the IMF with a reasonable initial Mach number for the molecular cloud. Since other physical processes (e.g. the equation of state and the variations among the core properties) have already been suggested to broaden the PDF, the relative importance of the different effects remains an open question.

Ionization driven molecular outflow in K3-50A

Whether high mass stars continue to accrete material beyond the formation of an HII region is still an open question. Ionized infall and outflow have been seen in some sources, but their ties to the surrounding molecular gas are not well constrained. We aim to quantify the ionized and molecular gas dynamics in a high mass star forming region (K3-50A) and their interaction. We present CARMA observations of the 3mm continuum, H41\alpha, and HCO+ emission, and VLA continuum observations at 23 GHz and 14.7 GHz to quantify the gas and its dynamics in K3-50A. We find large scale dynamics consistent with previous observations. On small scales, we find evidence for interaction between the ionized and molecular gas which suggests the ionized outflow is entraining the molecular one. This is the first time such an outflow entrained by photo ionized gas has been observed. Accretion may be ongoing in K3-50A because an ionized bipolar outflow is still being powered, which is in turn entraining part of the surrounding molecular gas. This outflow scenario is similar to that predicted by ionization feedback models.

The birth and death of star clusters in the Milky Way

Simple though admittedly speculative considerations explain why most of our galaxy’s stellar nurseries are highly fragile but a few survive for a remarkably long time.

Formation of structures around HII regions: ionization feedback from massive stars

AbstractWe present a new model for the formation of dense clumps and pillars around HII regions based on shocks curvature at the interface between a HII region and a molecular cloud. UV radiation leads to the formation of an ionization front and of a shock ahead. The gas is compressed between them forming a dense shell at the interface. This shell may be curved due to initial interface or density modulation caused by the turbulence of the molecular cloud. Low curvature leads to instabilities in the shell that form dense clumps while sufficiently curved shells collapse on itself to form pillars. When turbulence is high compared to the ionized-gas pressure, bubbles of cold gas have sufficient kinetic energy to penetrate into the HII region and detach themselves from the parent cloud, forming cometary globules.Using computational simulations, we show that these new models are extremely efficient to form dense clumps and stable and growing elongated structures, pillars, in which star formation might occur (see Tremblin et al.2012a). The inclusion of turbulence in the model shows its importance in the formation of cometary globules (see Tremblin et al.2012b). Globally, the density enhancement in the simulations is of one or two orders of magnitude higher than the density enhancement of the classical “collect and collapse“ scenario. The code used for the simulation is the HERACLES code, that comprises hydrodynamics with various equation of state, radiative transfer, gravity, cooling and heating.Our recent observations with Herschel (see Schneider et al.2012a) and SOFIA (see Schneider et al.2012b) and additional Spitzer data archives revealed many more of these structures in regions where OB stars have already formed such as the Rosette Nebula, Cygnus X, M16 and Vela, suggesting that the UV radiation from massive stars plays an important role in their formation. We present a first comparison between the simulations described above and recent observations of these regions.