We extend our previous work on simulations with the code RAMSES on accretion driven turbulence by including self-gravity and study the effects of core formation and collapse. We show that radial accretion onto filaments drives turbulent motions which are not isotropic but radially dominated. In contrast to filaments without gravity, the velocity dispersion of self-gravitating filaments does not settle in an equilibrium. Despite showing similar amounts of driven turbulence, they continually dissipate their velocity dispersion until the onset of core formation. This difference is connected to the evolution of the radius as it determines the dissipation rate. In the non-gravitational case filament growth is not limited and its radius grows linearly with time. In contrast, there is a maximum extent in the self-gravitational case resulting in an increased dissipation rate. Furthermore, accretion driven turbulence shows a radial profile which is anti-correlated with density. This leads to a constant turbulent pressure throughout the filament. As the additional turbulent pressure does not have a radial gradient it does not contribute to the stability of filaments and does not increase the critical line-mass. However, this radial turbulence does affect the radius of a filament, adding to the extent and setting its maximum value. Moreover, the radius evolution also affects the growth timescale of cores which compared to the timescale of collapse of an accreting filament limits core formation to high line-masses.

Forming high-mass stars have a significant effect on their natal environment. Their feedback pathways, including winds, outflows, and ionising radiation, shape the evolution of their surroundings which impacts the formation of the next generation of stars. They create or reveal dense pillars of gas and dust towards the edges of the cavities they clear. They are modelled in feedback simulations, and the sizes and shapes of the pillars produced are consistent with those observed. However, these models predict measurably different kinematics which provides testable discriminants. Here we present the first ALMA Compact Array (ACA) survey of 13 pillars in Carina, observed in $^{12}$CO, $^{13}$CO and C$^{18}$O J=2-1, and the 230 GHz continuum. The pillars in this survey were chosen to cover a wide range in properties relating to the amount and direction of incident radiation, proximity to nearby irradiating clusters and cloud rims, and whether they are detached from the cloud. With these data, we are able to discriminate between models. We generally find pillar velocity dispersions of $<$ 1 km s$^{-1}$ and that the outer few layers of molecular emission in these pillars show no significant offsets from each other, suggesting little bulk internal motions within the pillars. There are instances where the pillars are offset in velocity from their parental cloud rim, and some with no offset, hinting at a stochastic development of these motions.

Observations of prestellar cores in star-forming filaments show two distinct morphologies. While molecular line measurements often show broad cores, submillimeter continuum observations predominantly display pinched cores compared to the bulk of the filament gas. In order to explain how different morphologies arise, we use the gravitational instability model where prestellar cores form by growing density perturbations. The radial extent at each position is set by the local line-mass. We show that the ratio of core radius to filament radius is determined by the initial line-mass of the filament. Additionally, the core morphology is independent of perturbation length scale and inclination, which makes it an ideal diagnostic for observations. Filaments with a line-mass of less than half its critical value should form broad cores, whereas filaments with more than half its critical line-mass value should form pinched cores. For filaments embedded in a constant background pressure, the dominant perturbation growth times significantly differ for low and high line-mass filaments. Therefore, we predict that only one population of cores is present if all filaments within a region begin with similar initial perturbations.

We study accretion driven turbulence for different inflow velocities in star forming filaments using the code ramses. Filaments are rarely isolated objects and their gravitational potential will lead to radially dominated accretion. In the non-gravitational case, accretion by itself can already provoke non-isotropic, radially dominated turbulent motions responsible for the complex structure and non-thermal line widths observed in filaments. We find that there is a direct linear relation between the absolute value of the total density weighted velocity dispersion and the infall velocity. The turbulent velocity dispersion in the filaments is independent of sound speed or any net flow along the filament. We show that the density weighted velocity dispersion acts as an additional pressure term supporting the filament in hydrostatic equilibrium. Comparing to observations, we find that the projected non-thermal line width variation is generally subsonic independent of inflow velocity.

We study the stability of filaments in equilibrium between gravity and internal as well as external pressure using the grid based AMR-code RAMSES. A homogeneous, straight cylinder below a critical line mass is marginally stable. However, if the cylinder is bent, e.g. with a slight sinusoidal perturbation, an otherwise stable configuration starts to oscillate, is triggered into fragmentation and collapses. This previously unstudied behavior allows a filament to fragment at any given scale, as long as it has slight bends. We call this process `geometrical fragmentation'. In our realization the spacing between the cores matches the wavelength of the sinusoidal perturbation, whereas up to now, filaments were thought to be only fragmenting on the characteristical scale set by the mass-to-line ratio. Using first principles, we derive the oscillation period as well as the collapse timescale analytically. To enable a direct comparison with observations, we study the line-of-sight velocity for different inclinations. We show that the overall oscillation pattern can hide the infall signature of cores.

We present an analysis of the effect of feedback from O- and B-type stars with data from the integral field spectrograph MUSE mounted on the Very Large Telescope of pillar-like structures in the Carina Nebular Complex, one of the most massive star-forming regions in the Galaxy. For the observed pillars, we compute gas electron densities and temperatures maps, produce integrated line and velocity maps of the ionised gas, study the ionisation fronts at the pillar tips, analyse the properties of the single regions, and detect two ionised jets originating from two distinct pillar tips. For each pillar tip we determine the incident ionising photon flux $Q_\mathrm{0,pil}$ originating from the nearby massive O- and B-type stars and compute the mass-loss rate $\dot{M}$ of the pillar tips due to photo-evaporation caused by the incident ionising radiation. We combine the results of the Carina data set with archival MUSE data of a pillar in NGC 3603 and with previously published MUSE data of the Pillars of Creation in M16, and with a total of 10 analysed pillars, find tight correlations between the ionising photon flux and the electron density, the electron density and the distance from the ionising sources, and the ionising photon flux and the mass-loss rate. The combined MUSE data sets of pillars in regions with different physical conditions and stellar content therefore yield an empirical quantification of the feedback effects of ionising radiation. In agreement with models, we find that $\dot{M}\propto Q_\mathrm{0,pil}^{1/2}$.

The origin of the stellar initial mass function (IMF) is a fundamental issue in the theory of star formation. It is generally fit with a composite power law. Some clues on the progenitors can be found in dense starless cores that have a core mass function (CMF) with a similar shape. In the low-mass end, these mass functions increase with mass, albeit the sample may be somewhat incomplete; in the high-mass end, the mass functions decrease with mass. There is an offset in the turn-over mass between the two mass distributions. The stellar mass for the IMF peak is lower than the corresponding core mass for the CMF peak in the Pipe Nebula by about a factor of three. Smaller offsets are found between the IMF and the CMFs in other nebulae. We suggest that the offset is likely induced during a starburst episode of global star formation which is triggered by the formation of a few O/B stars in the multi-phase media, which naturally emerged through the onset of thermal instability in the cloud-core formation process. We consider the scenario that the ignition of a few massive stars photoionizes the warm medium between the cores, increases the external pressure, reduces their Bonnor?Ebert mass, and triggers the collapse of some previously stable cores. We quantitatively reproduce the IMF in the low-mass end with the assumption of additional rotational fragmentation.

Integral field unit (IFU) data of the iconic Pillars of Creation in M16 are presented. The ionisation structure of the pillars was studied in great detail over almost the entire visible wavelength range, and maps of the relevant physical parameters, e.g. extinction, electron density, electron temperature, line-of-sight velocity of the ionised and neutral gas are shown. In agreement with previous authors, we find that the pillar tips are being ionised and photo-evaporated by the massive members of the nearby cluster NGC 6611. They display a stratified ionisation structure where the emission lines peak in a descending order according to their ionisation energies. The IFU data allowed us to analyse the kinematics of the photo-evaporative flow in terms of the stratified ionisation structure, and we find that, in agreement with simulations, the photo-evaporative flow is traced by a blueshift in the position-velocity profile. The gas kinematics and ionisation structure have allowed us to produce a sketch of the 3D geometry of the Pillars, positioning the pillars with respect to the ionising cluster stars. We use a novel method to detect a previously unknown bipolar outflow at the tip of the middle pillar and suggest that it has an embedded protostar as its driving source. Furthermore we identify a candidate outflow in the leftmost pillar. With the derived physical parameters and ionic abundances, we estimate a mass loss rate due to the photo-evaporative flow of 70 M$_{\odot}$ Myr$^{-1}$ which yields an expected lifetime of approximately 3 Myr.

Proplyds and stars inside HII-regions are a well studied phenomenon. It is possible that they were triggered by the expansion of the HII-region itself. Here, we present calculations on the dynamics of HII-regions. We show that the triggered stars that form in the expanding shell of swept up material around the HII region rarely return into the HII regions on timescales that are inferred for the proplyds and observed young stars. However, in very dense environments like Orion, the triggered stars return in time. Thus, our model can explain why proplyds are barely observed in other HII regions. We propose that the properties of young stellar objects in HII regions in general depend critically on the distance from the massive, ionizing central star cluster. Closest in, there are proplyds, where the disk of a young star interacts directly with the feedback of the massive star. Further out are Class II protostars, where the ionization already removed the envelope. Even further away, one should find Class I stars, which either have been triggered by the ionizing radiation or pre-existed and have not lost their envelope yet. This radial sequence is not necessarily an age sequence but rather a result of the dwindling importance of stellar winds and ionizing radiation with distance. We investigate the observational signature of triggered star formation and find that the stellar distribution for ionization triggered star formation shows a distinct feature, a peak at the current position of the ionization front. Therefore, it is generally possible to tell triggered and in situ distributions of stars apart.

In the widely adopted LambdaCDM scenario for galaxy formation, dwarf galaxies are the building blocks of larger galaxies. Since they formed at relatively early epochs when the background density was relatively high, they are expected to retain their integrity as satellite galaxies when they merge to form larger entities. Although many dwarf spheroidal galaxies (dSphs) are found in the galactic halo around the Milky Way, their phase space density (or velocity dispersion) appears to be significantly smaller than that expected for satellite dwarf galaxies in the LambdaCDM scenario. In order to account for this discrepancy, we consider the possibility that they may have lost a significant fraction of their baryonic matter content during the first infall at the Hubble expansion turnaround. Such mass loss arises naturally due to the feedback by relatively massive stars which formed in their centers briefly before the maximum contraction. Through a series of N-body simulations, we show that the timely loss of a significant fraction of the dSphs initial baryonic matter content can have profound effects on their asymptotic half-mass radius, velocity dispersion, phase-space density, and the mass fraction between residual baryonic and dark matter.

We propose that the Pipe Nebula is an HII region shell swept up by the B2 IV beta Cephei star theta Ophiuchi. After reviewing the morphological evidence by recent observations, we perform a series of analytical calculations. We use realistic HII region parameters derived with the radiative transfer code Cloudy from observed stellar parameters. We are able to show that the current size, mass and pressure of the region can be explained in this scenario. We investigate the configuration today and come to the conclusion that the Pipe Nebula can be best described by a three phase medium in pressure equilibrium. The pressure support is provided by the ionized gas and mediated by an atomic component to confine the cores at the observed current pressure. In the future, star formation in these cores is likely to be either triggered by feedback of the most massive, gravitationally bound cores as soon as they collapse or by the supernova explosion of theta Ophiuchi itself.

We investigate the enrichment of the pre-solar cloud core with short lived radionuclides (SLRs), especially 26Al. The homogeneity and the surprisingly small spread in the ratio 26Al/27Al observed in the overwhelming majority of calcium-aluminium-rich inclusions (CAIs) in a vast variety of primitive chondritic meteorites places strong constraints on the formation of the the solar system. Freshly synthesized radioactive 26Al has to be included and well mixed within 20kyr. After discussing various scenarios including X-winds, AGB stars and Wolf-Rayet stars, we come to the conclusion that triggering the collapse of a cold cloud core by a nearby supernova is the most promising scenario. We then narrow down the vast parameter space by considering the pre-explosion survivability of such a clump as well as the cross-section necessary for sufficient enrichment. We employ numerical simulations to address the mixing of the radioactively enriched SN gas with the pre-existing gas and the forced collapse within 20kyr. We show that a cold clump of 10Msun at a distance of 5pc can be sufficiently enriched in 26Al and triggered into collapse fast enough - within 18kyr after encountering the supernova shock - for a range of different metallicities and progenitor masses, even if the enriched material is assumed to be distributed homogeneously in the entire supernova bubble. In summary, we envision an environment for the birth place of the Solar System 4.567Gyr ago similar to the situation of the pillars in M16 nowadays, where molecular cloud cores adjacent to an HII region will be hit by a supernova explosion in the future. We show that the triggered collapse and formation of the Solar System as well as the required enrichment with radioactive 26Al are possible in this scenario.

We present synthetic observations of star-forming interstellar medium structures obtained by hydrodynamic calculations of a turbulent box under the influence of an ionising radiation field. The morphological appearance of the pillar-like structures in optical emission lines is found to be very similar to observations of nearby star forming regions. We calculate line profiles as a function of position along the pillars for collisionally excited [OIII]5007, [NII]6584 and [SII]6717, which show typical FWHM of 2--4 km/s. Spatially resolved emission line diagnostic diagrams are also presented which show values in general agreement with observations of similar regions. The diagrams, however, also highlight significant spatial variations in the line ratios, including values that would be classically interpreted as shocked regions based on one-dimensional photoionisation calculations. These values tend to be instead the result of lines of sight intersecting which intersect for large portions of their lenghts the ionised-to-neutral transition regions in the gas. We caution therefore against a straightforward application of classical diagnostic diagrams and one--dimensional photoionisation calculations to spatially resolved observations of complex three-dimensional star forming regions.

We compare the three-dimensional gas temperature distributions obtained by a dedicated radiative transfer and photoionisation code, MOCASSIN, against those obtained by the recently-developed Smooth Particle Hydrodynamics (SPH) plus ionisation code iVINE for snapshots of an hydrodynamical simulation of a turbulent interstellar medium (ISM) irradiated by a nearby O star. Our tests demonstrate that the global ionisation properties of the region are correctly reproduced by iVINE, hence validating further application of this code to the study of feedback in star forming regions. However we highlight potentially important discrepancies in the detailed temperature distribution. In particular we show that in the case of highly inhomogenous density distributions the commonly employed on-the-spot (OTS) approximation yields unrealistically sharp shadow regions which can affect the dynamical evolution of the system. We implement a simple strategy to include the effects of the diffuse field in future calculations, which makes use of physically motivated temperature calibrations of the diffuse-field dominated regions and can be readily applied to similar codes. We find that while the global qualitative behaviour of the system is captured by simulations with the OTS approximation, the inclusion of the diffuse field in iVINE (called DiVINE) results in a stronger confinement of the cold gas, leading to denser and less coherent structures. This in turn leads to earlier triggering of star formation. We confirm that turbulence is being driven in simulations that include the diffuse field, but the efficiency is slightly lower than in simulations that use the OTS approximation.

Feedback from photoionisation may dominate on parsec scales in massive star-forming regions. Such feedback may inhibit or enhance the star formation efficiency and sustain or even drive turbulence in the parent molecular cloud. Photoionisation feedback may also provide a mechanism for the rapid expulsion of gas from young clusters' potentials, often invoked as the main cause of 'infant mortality'. There is currently no agreement, however, with regards to the efficiency of this process and how environment may affect the direction (positive or negative) in which it proceeds. The study of the photoionisation process as part of hydrodynamical simulations is key to understanding these issues, however, due to the computational demand of the problem, crude approximations for the radiation transfer are often employed. We will briefly review some of the most commonly used approximations and discuss their major drawbacks. We will then present the results of detailed tests carried out using the detailed photoionisation code MOCASSIN and the SPH+ionisation code iVINE code, aimed at understanding the error introduced by the simplified photoionisation algorithms. This is particularly relevant as a number of new codes have recently been developed along those lines. We will finally propose a new approach that should allow to efficiently and self-consistently treat the photoionisation problem for complex radiation and density fields.

We study the structural evolution of turbulent molecular clouds under the influence of ionizing radiation emitted from a nearby massive star by performing a high resolution parameter study with the iVINE code. The temperature is taken to be 10K or 100K, the mean number density is either 100cm^3 or 300cm^3. Besides, the turbulence is varied between Mach 1.5 and Mach 12.5 and the main driving scale between 1pc and 8pc. We vary the ionizing flux by an order of magnitude. In our simulations the ionizing radiation enhances the initial turbulent density distribution and thus leads to the formation of pillar-like structures observed adjacent to HII regions in a natural way. Gravitational collapse occurs regularly at the tips of the structures. We find a clear correlation between the initial state of the turbulent cold cloud and the final morphology and physical properties of the structures formed. The most favorable regime for the formation of pillars is Mach 4-10. Structures and therefore stars only form if the initial density contrast between the high density unionized gas and the gas that is going to be ionized is lower than the temperature contrast between the hot and the cold gas. The density of the resulting pillars is determined by a pressure equilibrium between the hot and the cold gas. A thorough analysis of the simulations shows that the complex kinematical and geometrical structure of the formed elongated filaments reflects that of observed pillars to an impressive level of detail. In addition, we find that the observed line-of sight velocities allow for a distinct determination of different formation mechanisms. Comparing the current simulations to previous results and recent observations we conclude that e.g. the pillars of creation in M16 formed by the mechanism proposed here and not by the radiation driven implosion of pre-existing clumps.

Contradicting results have been reported in the literature with respect to the performance of the numerical techniques employed for the study of supersonic turbulence. We aim at characterising the performance of different particle-based and grid-based techniques on the modelling of decaying supersonic turbulence. Four different grid codes (ENZO, FLASH, TVD, ZEUS) and three different SPH codes (GADGET, PHANTOM, VINE) are compared. We additionally analysed two calculations denoted as PHANTOM A and PHANTOM B using two different implementations of artificial viscosity. Our analysis indicates that grid codes tend to be less dissipative than SPH codes, though details of the techniques used can make large differences in both cases. For example, the Morris & Monaghan viscosity implementation for SPH results in less dissipation (PHANTOM B and VINE versus GADGET and PHANTOM A). For grid codes, using a smaller diffusion parameter leads to less dissipation, but results in a larger bottleneck effect (our ENZO versus FLASH runs). As a general result, we find that by using a similar number of resolution elements N for each spatial direction means that all codes (both grid-based and particle-based) show encouraging similarity of all statistical quantities for isotropic supersonic turbulence on spatial scales k<N/32 (all scales resolved by more than 32 grid cells), while scales smaller than that are significantly affected by the specific implementation of the algorithm for solving the equations of hydrodynamics. At comparable numerical resolution, the SPH runs were on average about ten times more computationally intensive than the grid runs, although with variations of up to a factor of ten between the different SPH runs and between the different grid runs. (abridged)

We investigate the collapse and fragmentation of low-mass, trans-sonically turbulent prestellar cores, using SPH simulations. The initial conditions are slightly supercritical Bonnor-Ebert spheres, all with the same density profile, the same mass (M_O=6.1 Msun) and the same radius (R_O=17,000 AU), but having different initial turbulent velocity fields. Four hundred turbulent velocity fields have been generated, all scaled so that the mean Mach number is M=1. Then a subset of these, having a range of net angular momenta, j, has been evolved. The evolution of these turbulent cores is not strongly correlated with j. Instead it is moderated by the formation of filamentary structures due to converging turbulent flows. A high fraction (~ 82%) of the protostars forming from turbulent cores are attended by protostellar accretion discs, but only a very small fraction (~16%) of these discs is sufficiently cool and extended to develop non-linear gravitational instabilities and fragment.

Abridged: We use three-dimensional SPH simulations to investigate the collapse of low-mass prestellar cores and the formation and early evolution of protostellar discs. The initial conditions are slightly supercritical Bonnor-Ebert spheres in rigid rotation. The core mass and initial radius are held fixed at M_O=6.1 M_sun and R_O=17,000 AU, and the only parameter that we vary is the initial angular speed \Omega_O. Protostellar discs forming from cores with \Omega_O<1.35 10d-13 1/s have radii between 100 and 300 AU and are quite centrally concentrated; due to heating by gas infall onto the disc and accretion onto the central object, they are also quite warm, T>100 K, and therefore stable against gravitational fragmentation. In contrast, more rapidly rotating cores form discs which are less concentrated and cooler, and have radii between 400 and 1000 AU; as a consequence they are prone to gravitational fragmentation and the formation of multiple systems. We derive a criterion that predicts whether a rigidly rotating core having given M_O, R_O and \Omega_O will produce a protostellar disc which fragments whilst material is still infalling from the core envelope. We then apply this criterion to core samples for which M_O, R_O and \Omega_O have been estimated observationally. We conclude that the observed cores are stable against fragmentation at this stage, due to their low angular speeds and the heat delivered at the accretion shock where the infalling material hits the disc.

We present high resolution simulations on the impact of ionizing radiation of massive O-stars on the surrounding turbulent interstellar medium (ISM). The simulations are performed with the newly developed software iVINE which combines ionization with smoothed particle hydrodynamics (SPH) and gravitational forces. We show that radiation from hot stars penetrates the ISM, efficiently heats cold low density gas and amplifies over-densities seeded by the initial turbulence. The formation of observed pillar-like structures in star forming regions (e.g. in M16) can be explained by this scenario. At the tip of the pillars gravitational collapse can be induced, eventually leading to the formation of low mass stars. Detailed analysis of the evolution of the turbulent spectra shows that UV-radiation of O-stars indeed provides an excellent mechanism to sustain and even drive turbulence in the parental molecular cloud.

We present a three-dimensional, fully parallelized, efficient implementation of ionizing UV radiation for smoothed particle hydrodynamics (SPH) including self-gravity. Our method is based on the SPH/tree code VINE. We therefore call it iVINE (for Ionization + VINE). This approach allows detailed high-resolution studies of the effects of ionizing radiation from e.g. young massive stars on their turbulent parental molecular clouds. In this paper we describe the concept and the numerical implementation of the radiative transfer for a plain-parallel geometry and we discuss several test cases demonstrating the efficiency and accuracy of the new method. As a first application, we study the radiatively driven implosion of marginally stable molecular clouds at various distances of a strong UV source and show that they are driven into gravitational collapse. The resulting cores are very compact and dense exactly as it is observed in clustered environments. Our simulations indicate that the time of triggered collapse depends on the distance of the core from the UV source. Clouds closer to the source collapse several $10^5$ years earlier than more distant clouds. This effect can explain the observed age spread in OB associations where stars closer to the source are found to be younger. We discuss possible uncertainties in the observational derivation of shock front velocities due to early stripping of proto-stellar envelopes by ionizing radiation.

Recent observations with the Spitzer Space Telescope show clear evidence that star formation takes place in the surrounding of young massive O-type stars, which are shaping their environment due to their powerful radiation and stellar winds. In this work we investigate the effect of ionising radiation of massive stars on the ambient interstellar medium (ISM): In particular we want to examine whether the UV-radiation of O-type stars can lead to the observed pillar-like structures and can trigger star formation. We developed a new implementation, based on a parallel Smooth Particle Hydrodynamics code (called IVINE), that allows an efficient treatment of the effect of ionising radiation from massive stars on their turbulent gaseous environment. Here we present first results at very high resolution. We show that ionising radiation can trigger the collapse of an otherwise stable molecular cloud. The arising structures resemble observed structures (e.g. the pillars of creation in the Eagle Nebula (M16) or the Horsehead Nebula B33). Including the effect of gravitation we find small regions that can be identified as formation places of individual stars. We conclude that ionising radiation from massive stars alone can trigger substantial star formation in molecular clouds.

Plerionic supernova remnants exhibit radio emission with remarkably flat spectral indices ranging from $\alpha=0.0$ to $\alpha=-0.3$. The origin of very hard particle energy distributions still awaits an explanation, since shock waves generate particle distributions with synchrotron spectra characterized by $\alpha\le-0.5$. Acceleration of high energy leptons in magnetohydrodynamic turbulence instead may be responsible for the observed hard spectra. This process is studied by means of relativistic test particle calculations using electromagnetic fields produced by three-dimensional simulations of resistive magnetohydrodynamical turbulence. The particles receive power-law energy spectra $N(\gamma)\propto \gamma^{-s}$ with $s$ ranging from 1.2 to 1.6, i.e. particle spectra that are required to explain the radio emission of plerions.