Smoothed Particle Hydrodynamics simulations of expanding HII regions. I. Numerical methods and tests

We describe a new algorithm for including the dynamical effects of ionizing radiation in SPH simulations, and we present several examples of how the algorithm can be applied to problems in star formation. We use the HEALPix software to tessellate the sky and to solve the equation of ionization equilibrium along a ray towards each of the resulting tesserae. We exploit the hierarchical nature of HEALPix to make the algorithm adaptive, so that fine angular resolution is invoked only where it is needed, and the computational cost is kept low. We present simulations of (i) the spherically symmetric expansion of an HII region inside a uniform-density, non--self-gravitating cloud; (ii) the spherically symmetric expansion of an HII region inside a uniform-density, self-gravitating cloud; (iii) the expansion of an off-centre HII region inside a uniform-density, non--self-gravitating cloud, resulting in rocket acceleration and dispersal of the cloud; and (iv) radiatively driven compression and ablation of a core overrun by an HII region. The new algorithm provides the means to explore and evaluate the role of ionizing radiation in regulating the efficiency and statistics of star formation.Comment: 12 pages, 16 figures, simulation movies available at http://galaxy.ig.cas.cz/~richard/HIIregion

ALMA observations of atomic carbon in z~4 dusty star-forming galaxies

We present ALMA [CI]($1-0$) (rest frequency 492 GHz) observations for a sample of 13 strongly-lensed dusty star-forming galaxies originally discovered at 1.4mm in a blank-field survey by the South Pole Telescope. We compare these new data with available [CI] observations from the literature, allowing a study of the ISM properties of $\sim 30$ extreme dusty star-forming galaxies spanning a redshift range $2 < z < 5$. Using the [CI] line as a tracer of the molecular ISM, we find a mean molecular gas mass for SPT-DSFGs of $6.6 \times 10^{10}$ M$_{\odot}$. This is in tension with gas masses derived via low-$J$ $^{12}$CO and dust masses; bringing the estimates into accordance requires either (a) an elevated CO-to-H$_2$ conversion factor for our sample of $\alpha_{\rm CO} \sim 2.5$ and a gas-to-dust ratio $\sim200$, or (b) an high carbon abundance $X_{\rm CI} \sim 7\times10^{-5}$. Using observations of a range of additional atomic and molecular lines (including [CI], [CII], and multiple transitions of CO), we use a modern Photodissociation Region code (3D-PDR) to assess the physical conditions (including the density, UV radiation field strength, and gas temperature) within the ISM of the DSFGs in our sample. We find that the ISM within our DSFGs is characterised by dense gas permeated by strong UV fields. We note that previous efforts to characterise PDR regions in DSFGs may have significantly underestimated the density of the ISM. Combined, our analysis suggests that the ISM of extreme dusty starbursts at high redshift consists of dense, carbon-rich gas not directly comparable to the ISM of starbursts in the local Universe.Comment: 21 pages, 12 figures. Accepted for publication in MNRA

Comparing simulations of ionisation triggered star formation and observations in RCW 120

Massive clumps within the swept-up shells of bubbles, like that surrounding the galactic HII region RCW 120, have been interpreted in terms of the Collect and Collapse (C&C) mechanism for triggered star formation. The cold, dusty clumps surrounding RCW 120 are arranged in an almost spherical shell and harbour many young stellar objects. By performing high-resolution, three-dimensional SPH simulations of HII regions expanding into fractal molecular clouds, we investigate whether the formation of massive clumps in dense, swept-up shells necessarily requires the C&C mechanism. In a second step, we use RADMC-3D to compute the synthetic dust continuum emission from our simulations, in order to compare them with observations of RCW 120 made with APEX-LABOCA at 870 micron. We show that a distribution of clumps similar to the one seen in RCW 120 can readily be explained by a non-uniform initial molecular cloud structure. Hence, a shell-like configuration of massive clumps does not imply that the C&C mechanism is at work. Rather, we find a hybrid form of triggering, which combines elements of C&C and Radiatively Driven Implosion (RDI). In addition, we investigate the reliability of deriving clump masses from their 870 micron emission. We find that for clumps with more than 100 M_sun the observational estimates are accurate to within a factor of two and that, even at these long wavelengths, it is important to account for the radiative heating from triggered, embedded protostars.Comment: 10 pages, 7 figures, submitted to MNRA

The Milky Way Project: A statistical study of massive star formation associated with infrared bubbles

The Milky Way Project citizen science initiative recently increased the number of known infrared bubbles in the inner Galactic plane by an order of magnitude compared to previous studies. We present a detailed statistical analysis of this dataset with the Red MSX Source catalog of massive young stellar sources to investigate the association of these bubbles with massive star formation. We particularly address the question of massive triggered star formation near infrared bubbles. We find a strong positional correlation of massive young stellar objects (MYSOs) and H II regions with Milky Way Project bubbles at separations of < 2 bubble radii. As bubble sizes increase, a statistically significant overdensity of massive young sources emerges in the region of the bubble rims, possibly indicating the occurrence of triggered star formation. Based on numbers of bubble-associated RMS sources we find that 67+/-3% of MYSOs and (ultra)compact H II regions appear associated with a bubble. We estimate that approximately 22+/-2% of massive young stars may have formed as a result of feedback from expanding H II regions. Using MYSO-bubble correlations, we serendipitously recovered the location of the recently discovered massive cluster Mercer 81, suggesting the potential of such analyses for discovery of heavily extincted distant clusters.Comment: 16 pages, 17 figures. Accepted for publication in ApJ, comments welcome. Milky Way Project public data release available at http://www.milkywayproject.org/dat

STARBENCH: the D-type expansion of an HII region

STARBENCH is a project focused on benchmarking and validating different star formation and stellar feedback codes. In this first STARBENCH paper we perform a comparison study of the D-type expansion of an H II region. The aim of this work is to understand the differences observed between the 12 participating numerical codes against the various analytical expressions examining the D-type phase of H II region expansion. To do this, we propose two well-defined tests which are tackled by 1D and 3D grid- and smoothed particle hydrodynamics-based codes. The first test examines the ‘early phase’ D-type scenario during which the mechanical pressure driving the expansion is significantly larger than the thermal pressure of the neutral medium. The second test examines the ‘late phase’ D-type scenario during which the system relaxes to pressure equilibrium with the external medium. Although they are mutually in excellent agreement, all 12 participating codes follow a modified expansion law that deviates significantly from the classical Spitzer solution in both scenarios. We present a semi-empirical formula combining the two different solutions appropriate to both early and late phases that agrees with high-resolution simulations to ≲ 2 per cent. This formula provides a much better benchmark solution for code validation than the Spitzer solution. The present comparison has validated the participating codes and through this project we provide a data set for calibrating the treatment of ionizing radiation hydrodynamics codes

Evolution of prolate molecular clouds at H II boundaries - II. Formation of BRCs of asymmetrical morphology

A systematic investigation on the evolution of a prolate cloud at an H II boundary is conducted using smoothed particle hydrodynamics in order to understand the mechanism for a variety of irregular morphological structures found at the boundaries of various H II regions. The prolate molecular clouds in this investigation are set with their semimajor axes at inclinations between 0° and 90° to a plane-parallel ionizing radiation flux. A set of four parameters, the number density n, the ratio of major to minor axis γ, the inclination angle Φ and the incident flux FEUV, are used to define the initial state of the simulated clouds. The dependence of the evolution of a prolate cloud under radiation-driven implosion (RDI) on each of the four parameters is investigated. It is found that (i) in addition to the well-studied standard type A, B or C bright-rimmed clouds (BRCs), many other types such as asymmetrical BRCs, filamentary structures and irregular horse-head structures could also be developed at H II boundaries with only simple initial conditions; (ii) the final morphological structures are very sensitive to the four initial parameters, especially to the initial density and the inclination; (iii) the previously defined ionizing radiation penetration depth can still be used as a good indicator of the final morphology. Based on the simulation results, the formation time-scales and masses of the early RDI-triggered star formation from clouds of different initial conditions are also estimated. Finally a unified mechanism for the various morphological structures found in many different H II boundaries is suggested

Evolution of prolate molecular clouds at H II boundaries - I. Formation of fragment-core structures

The evolution of a prolate cloud at an H II boundary is investigated using smoothed particle hydrodynamics. The prolate molecular clouds in our investigation are set with their semi-major axis perpendicular to the radiative direction of a plane-parallel ionizing extreme ultraviolet (EUV) flux. Simulations on three high-mass prolate clouds reveal that EUV radiation can trigger distinctive high-density core formation embedded in a final linear structure. This contrasts with results of the previous work in which only an isotropic far-ultraviolet interstellar background flux was applied. A systematic investigation on a group of prolate clouds of equal mass but different initial densities and geometric shapes finds that the distribution of the cores over the final linear structure changes with the initial conditions of the prolate cloud and the strength of the EUV radiation flux. These highly condensed cores may either scatter over the full length of the final linear structure or form two groups of high-density cores at two foci, depending on the value of the ionizing radiation penetration depth dEUV, the ratio of the physical ionizing radiation penetration depth to the minor axis of the cloud. Data analysis on the total mass of the high-density cores and the core formation time finds that the potential for EUV radiation triggered star formation efficiency is higher in prolate clouds with shallow ionization penetration depth and intermediate major-to-minor axial ratio, for the physical environments investigated. Finally, it is suggested that the various fragment-core structures observed at H II boundaries may result from the interaction between ionizing radiation and pre-existing prolate clouds of different initial geometrical and physical conditions

Ionization compression impact on dense gas distribution and star formation: probability density functions around H II regions as seen by <i>Herschel</i>

Aims. Ionization feedback should impact the probability distribution function (PDF) of the column density of cold dust around the ionized gas. We aim to quantify this effect and discuss its potential link to the core and initial mass function (CMF/IMF). Methods. We used Herschel column density maps of several regions observed within the HOBYS key program in a systematic way: M 16, the Rosette and Vela C molecular clouds, and the RCW 120 H II region. We computed the PDFs in concentric disks around the main ionizing sources, determined their properties, and discuss the effect of ionization pressure on the distribution of the column density. Results. We fitted the column density PDFs of all clouds with two lognormal distributions, since they present a "double-peak" or an 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. Such a double peak is not visible for all clouds associated with ionization fronts, but it depends on the relative importance of ionization pressure and turbulent ram pressure. A power-law tail is present for higher column densities, which are generally ascribed to the effect of gravity. 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 that is able to disentangle triggered star formation from pre-existing star formation. Conclusions. In the context of the gravo-turbulent scenario for the origin of the CMF/IMF, the double-peaked or enlarged shape of the PDF may affect the formation of objects at both the low-mass and the high-mass ends of the CMF/IMF. In particular, a broader PDF is required by the gravo-turbulent scenario to fit the IMF properly 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 said to broaden the PDF, the relative importance of the different effects remains an open question

Physical Processes in Star Formation

© 2020 Springer-Verlag. The final publication is available at Springer via https://doi.org/10.1007/s11214-020-00693-8.Star formation is a complex multi-scale phenomenon that is of significant importance for astrophysics in general. Stars and star formation are key pillars in observational astronomy from local star forming regions in the Milky Way up to high-redshift galaxies. From a theoretical perspective, star formation and feedback processes (radiation, winds, and supernovae) play a pivotal role in advancing our understanding of the physical processes at work, both individually and of their interactions. In this review we will give an overview of the main processes that are important for the understanding of star formation. We start with an observationally motivated view on star formation from a global perspective and outline the general paradigm of the life-cycle of molecular clouds, in which star formation is the key process to close the cycle. After that we focus on the thermal and chemical aspects in star forming regions, discuss turbulence and magnetic fields as well as gravitational forces. Finally, we review the most important stellar feedback mechanisms.Peer reviewedFinal Accepted Versio

Towards sustainable credit union development A research report

SIGLEAvailable from British Library Document Supply Centre-DSC:f99/3819 / BLDSC - British Library Document Supply CentreGBUnited Kingdo