Massive and very massive stars can play important roles in stellar populations by ejecting strong stellar winds and exploding in energetic phenomena. It is therefore of high importance that their behaviour is properly accounted for in synthetic model populations. Here we present nine grids of stellar evolutionary model sequences, together with finely resolved interpolated sequences and synthetic populations, of stars with 9-500 Msun and with metallicities ranging from Solar down to 1/250 Solar. The stellar models were computed with the 'Bonn' evolutionary code (covering core-hydrogen- and core-helium-burning phases, both complete). Post-processing for publication has been done using optimized methods, developed by our team, for massive and very massive stars. Interpolation and population synthesis were also performed on the models by our newly developed routine synStars. Eight of the grids represent slowly rotating massive stars with normal/classical evolution, while one grid represents fast rotating, chemically-homogeneously evolving models. Apart from the common stellar parameters such as mass, radius, surface temperature, luminosity and mass loss rate, we present stellar wind properties such as estimated wind velocity and kinetic energy of the wind. Additionally, we provide complete chemical yields of 34 isotopes, and estimates for the masses of the compact object remnants. The 'Bonn' Optimized Stellar Tracks (BoOST) project is published as simple tables - including stellar models, interpolated tracks and synthetic populations - thus ideal for further scientific applications. For example, star-formation studies could be done with BoOST to cover broad metallicity ranges, and so could be simulations of high-redshift galaxies. Additionally, gravitational-wave event rate predictions could be refined using BoOST by accounting for very massive stars at low-metallicity.

Although many galaxies in the Virgo cluster are known to have lost significant amounts of HI gas, only about a dozen features are known where the HI extends significantly outside its parent galaxy. Previous numerical simulations have predicted that HI removed by ram pressure stripping should have column densities far in excess of the sensitivity limits of observational surveys. We construct a simple model to try and quantify how many streams we might expect to detect. This accounts for the expected random orientation of the streams in position and velocity space as well as the expected stream length and mass of stripped HI. Using archival data from the Arecibo Galaxy Environment Survey, we search for any streams which might previously have been missed in earlier analyses. We report the confident detection of ten streams as well as sixteen other less sure detections. We show that these well-match our analytic predictions for which galaxies should be actively losing gas, however the mass of the streams is typically far below the amount of missing HI in their parent galaxies, implying that a phase change and/or dispersal renders the gas undetectable. By estimating the orbital timescales we estimate that dissolution rates of 1-10 Msolar/yr are able to explain both the presence of a few long, massive streams and the greater number of shorter, less massive features.

Following the observational and theoretical evidence that points at core collapse supernovae as major producers of dust, here we calculate the hydrodynamics of the matter reinserted within young and massive super stellar clusters under the assumption of gas and dust radiative cooling. The large supernova rate expected in massive clusters allows for a continuous replenishment of dust immersed in the high temperature thermalized reinserted matter and warrants a stationary presence of dust within the cluster volume during the type II supernova era. We first show that such a balance determines the range of dust to gas mass ratio and this the dust cooling law. We then search for the critical line that separates stationary cluster winds from the bimodal cases in the cluster mechanical luminosity (or cluster mass) vs cluster size parameter space. In the latter, strong radiative cooling reduces considerably the cluster wind mechanical energy output and affects particularly the cluster central regions, leading to frequent thermal instabilities that diminish the pressure and inhibit the exit of the reinserted matter. Instead matter accumulates there and is expected to eventually lead to gravitational instabilities and to further stellar formation with the matter reinserted by former massive stars. The main outcome of the calculations is that the critical line is almost two orders of magnitude or more, depending on the assumed value of the adiabatic wind terminal speed, lower than when only gas radiative cooling is applied. And thus, many massive clusters are predicted to enter the bimodal regime.

By means of 3-D hydrodynamical simulations, here we evaluate the impact that supernova explosions occurring within wind-driven bubbles have on the survival or destruction of dust grains. We consider both, the dust generated within the ejecta and the dust initially present in the ambient gas and later locked-up in the surrounding wind-driven shell. The collision of the supernova blast wave with the wind-driven shell leads to a transmitted shock that moves into the shell and a reflected shock into the ejecta. The transmitted shock is capable of destroying large amounts of the dust locked in the shell, but only if the mass of the wind-driven shell is small, less than a few tens the ejected mass. Conversely, massive wind-driven shells, with several times the ejected mass, lead upon the interaction to strong radiative cooling, which inhibits the Sedov-Taylor phase and weakens the transmitted shock, making it unable to traverse the wind-driven shell. In such a case, the destruction/disruption of the ambient dust is largely inhibited. On the other hand, the SNRs grow rapidly in the very tenuous region excavated by the stellar winds, and thus a large fraction of the dust generated within the ejecta is not efficiently destroyed by the supernova reverse shock, nor by the reflected shock. Our calculations favor a scenario in which core-collapse supernovae within sufficiently massive wind-driven shells supply more dust to the ISM than what they are able to destroy.

Sources of X-rays such as active galactic nuclei and X-ray binaries are often variable by orders of magnitude in luminosity over timescales of years. During and after these flares the surrounding gas is out of chemical and thermal equilibrium. We introduce a new implementation of X-ray radiative transfer coupled to a time-dependent chemical network for use in 3D magnetohydrodynamical simulations. A static fractal molecular cloud is irradiated with X-rays of different intensity, and the chemical and thermal evolution of the cloud are studied. For a simulated $10^5$ M$_\odot$ fractal cloud an X-ray flux $<0.01$ erg cm$^{-2}$ s$^{-1}$ allows the cloud to remain molecular, whereas most of the CO and H$_2$ are destroyed for a flux of $>1$ erg cm$^{-2}$ s$^{-1}$. The effects of an X-ray flare, which suddenly increases the X-ray flux by $10^5 \times$ are then studied. A cloud exposed to a bright flare has 99% of its CO destroyed in 10-20 years, whereas it takes $>10^3$ years for 99% of the H$_2$ to be destroyed. CO is primarily destroyed by locally generated far-UV emission from collisions between non-thermal electrons and H$_2$; He$^+$ only becomes an important destruction agent when the CO abundance is already very small. After the flare is over, CO re-forms and approaches its equilibrium abundance after $10^3-10^5$ years. This implies that molecular clouds close to Sgr A$^*$ in the Galactic Centre may still be out of chemical equilibrium, and we predict the existence of clouds near flaring X-ray sources in which CO has been mostly destroyed but H is fully molecular.

Multiple stellar populations are observed in almost all globular-clusters, but the origin of this phenomenon is still debated. We investigate the role cool supergiants may have played. To do this, we combine two investigative methods: state-of-the-art massive stellar evolution and calculations of the hydrodynamic structure of the cluster-gas. This approach allows us to study how star-formation in young massive clusters depends on the energy- and mass-input of the first-generation of stars, while predicting the chemical composition of the second-generation. We find that the presence of massive (9-500 M$_{\odot}$) metal-poor supergiants in the young cluster leads to a star-formation episode within the first 4 Myr of the cluster's lifetime, that is, before the first core-collapse supernovae explode or the gas is expelled. The stellar winds accumulate in the cluster center, forming the second-generation there. Its composition is predicted to show variations in O & Na abundances, consistently with observations. The abundance of helium is, similarly to other scenarios involving massive stars, higher than what is referred from observations. Supposing dynamical removal of stars from the outskirts of the cluster, or applying a top-heavy initial-mass-function, we can predict a number ratio of the second-generation as high as 20-80%. The effect of metallicity is shown to be important, as the most luminous supergiants are only predicted at low-metallicity, thus limiting $-$ but not excluding $-$ the extent of a polluted second-generation at high-metallicity. These massive stars becoming black-holes suggests globular-clusters hosting gravitational-wave progenitors. Our scenario predicts a correlation between the mass of the cluster and the extent of the multiple population phenomenon.

The C$^{+}$ ion is an important coolant of interstellar gas, and so the [CII] fine structure line is frequently observed in the interstellar medium. However, the physical and chemical properties of the [CII]-emitting gas are still unclear. We carry out non-LTE radiative transfer simulations with RADMC-3D to study the [CII] line emission from a young, turbulent molecular cloud before the onset of star formation, using data from the SILCC-Zoom project. The [CII] emission is optically thick over 40% of the observable area with $I_{[\textrm{CII}]} > 0.5$ K km s$^{-1}$. To determine the physical properties of the [CII] emitting gas, we treat the [CII] emission as optically thin. We find that the [CII] emission originates primarily from cold, moderate density gas ($40 \lesssim T \lesssim 65$ K and $50 \lesssim n \lesssim 440$ cm$^{-3}$), composed mainly of atomic hydrogen and with an effective visual extinction between $\sim 0.50$ and $\sim 0.91$. Gas dominated by molecular hydrogen contributes only $\lesssim$20% of the total [CII] line emission. Thus, [CII] is not a good tracer for CO-dark H$_2$ at this early phase in the cloud's lifetime. We also find that the total gas, H and C$^+$ column densities are all correlated with the integrated [CII] line emission, with power law slopes ranging from 0.5 to 0.7. Further, the median ratio between the total column density and the [CII] line emission is $Y_{{\rm CII}}\approx 1.1 \times 10^{21}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$, and $Y_{{\rm CII}}$ scales with $I_{[\textrm{CII}]}^{-0.3}$. We expect $Y_{{\rm CII}}$ to change in environments with a lower or higher radiation field than simulated here.

Following the current debate on the fate of SN-condensed dust grains, here we explore by means of three-dimensional hydrodynamical simulations the interaction of dusty supernova remnants (SNRs) with the shocked winds of neighboring massive stars within young massive stellar clusters (SSCs). As a comparison, we have also explored the evolution of supernova remnants in the diffuse ISM with constant density. Since the hydrodynamics of SNRs is intimately related to the properties of their immediate environment, the lifecycle of dust grains in SNRs within SSCs is radically different from that in the diffuse ISM. Moreover, off-centered SNRs evolving in the steep density gradient established due to a star cluster wind experience a blowout phase: shell fragmentation due to protruding Rayleigh-Taylor instabilities and the venting of SN ejecta. Our main finding is that clustered SN explosions will cause a net increase in the amount of dust in the surroundings of young massive stellar clusters. Our analysis considers the multiple dust processing resulting from the passage of the SN reverse shock, including its reflection at the SNR's center, the injection of shocked stellar winds within the respective remnant's volume and the effect of secondary forward shocks produced in sequential SN explosions. In the simulations, we have on-the-fly calculated the rates of thermal sputtering and dust-induced radiative cooling provided an initial distribution of grain sizes and dust content. Fast-moving elongated dusty SN ejecta resemble mushroom clouds violently ascending in a stratified atmosphere after volcanic super-eruptions, where the pyroclasts carried by the clouds are wind-driven and eventually accumulate into the vast surroundings.

Magnetic fields are ubiquitously observed in the interstellar medium (ISM) of present-day star-forming galaxies with dynamically relevant energy densities. Using three-dimensional magneto-hydrodynamic (MHD) simulations of the supernova (SN) driven ISM in the flux-freezing approximation (ideal MHD) we investigate the impact of the magnetic field on the chemical and dynamical evolution of the gas, fragmentation and the formation of molecular clouds. We follow the chemistry with a network of six species (H$^{+}$, H, H$_2$, C$^+$, CO, free electrons) including local shielding effects. We find that magnetic fields thicken the disc by a factor of a few to a scale height of $\sim100\,\mathrm{pc}$, delay the formation of dense (and molecular) gas by $\sim25\,\mathrm{Myr}$ and result in differently shaped gas structures. The magnetised gas fragments into fewer clumps, which are initially at subcritical mass-to-flux ratios, $M/\Phi\approx0.3(M/\Phi)_\mathrm{crit}$, and accrete gas preferentially parallel to the magnetic field lines until supercritial mass-to-flux ratios of up to order 10 are reached. The accretion rates onto molecular clouds scale with $\dot{M}\propto M^{1.5}$. The median of the inter-cloud velocity dispersion is $\sim2-5\,\mathrm{km\,s}^{-1}$ and lower than the internal velocity dispersion in the clouds ($\sim3-7\,\mathrm{km\,s}^{-1}$). However, individual cloud-cloud collisions occur at speeds of a few $10\,\mathrm{km\,s}^{-1}$.

Most detected neutral atomic hydrogen (HI) at low redshift is associated with optically bright galaxies. However, a handful of HI clouds are known which appear to be optically dark and have no nearby potential progenitor galaxies, making tidal debris an unlikely explanation. In particular, 6 clouds identified by the Arecibo Galaxy Environment Survey are interesting due to the combination of their small size, isolation, and especially their broad line widths atypical of other such clouds. A recent suggestion is that these clouds exist in pressure equilibrium with the intracluster medium, with the line width arising from turbulent internal motions. Here we explore that possibility by using the FLASH code to perform a series of 3D hydro simulations. Our clouds are modelled using spherical Gaussian density profiles, embedded in a hot, low-density gas representing the intracluster medium. The simulations account for heating and cooling of the gas, and we vary the structure and strength of their internal motions. We create synthetic HI spectra, and find that none of our simulations reproduce the observed cloud parameters for longer than about 100 Myr : the clouds either collapse, disperse, or experience rapid heating which would cause ionisation and render them undetectable to HI surveys. While the turbulent motions required to explain the high line widths generate structures which appear to be inherently unstable, making this an unlikely explanation for the observed clouds, these simulations demonstrate the importance of including the intracluster medium in any model seeking to explain the existence of these objects.

Photoionizing radiation and stellar winds from massive stars deposit energy and momentum into the interstellar medium (ISM). They might disperse the local ISM, change its turbulent multi-phase structure, and even regulate star formation. Ionizing radiation dominates the massive stars' energy output, but the relative effect of winds might change with stellar mass and the properties of the ambient ISM. We present simulations of the interaction of stellar winds and ionizing radiation of 12, 23, and 60 M$_{\odot}$ stars within a cold neutral (CNM, $n_{0}$ = 100 cm$^{-3}$), warm neutral (WNM, $n_{0}$ = 1, 10 cm$^{-3}$) or warm ionized (WIM, $n_{0}$ = 0.1 cm$^{-3}$) medium. The FLASH simulations adopt the novel tree-based radiation transfer algorithm TreeRay. With the On-the-Spot approximation and a temperature-dependent recombination coefficient, it is coupled to a chemical network with radiative heating and cooling. In the homogeneous CNM, the total momentum injection ranges from 1.6$\times$10$^{4}$ to 4$\times$10$^{5}$ M$_{\odot}$ km s$^{-1}$ and is always dominated by the expansion of the ionized H$_{\text{II}}$ region. In the WIM, stellar winds dominate (2$\times$10$^{2}$ to 5$\times$10$^{3}$ M$_{\odot}$ km s$^{-1}$), while the input from radiation is small ($\sim$ 10$^{2}$ M$_{\odot}$ km s$^{-1}$). The WNM ($n_{0}$ = 1 cm$^{-3}$) is a transition regime. Energetically, stellar winds couple more efficiently to the ISM ($\sim$ 0.1 percent of wind luminosity) than radiation ($<$ 0.001 percent of ionizing luminosity). For estimating the impact of massive stars, the strongly mass-dependent ratios of wind to ionizing luminosity and the properties of the ambient medium have to be considered.

We describe an OctTree algorithm for the MPI-parallel, adaptive mesh-refinement code {\sc FLASH}, which can be used to calculate the gas self-gravity, and also the angle-averaged local optical depth, for treating ambient diffuse radiation. The algorithm communicates to the different processors only those parts of the tree that are needed to perform the tree walk locally. The advantage of this approach is a relatively low memory requirement, important in particular for the optical depth calculation, which needs to process information from many different directions. This feature also enables a general tree-based radiation transport algorithm that will be described in a subsequent paper, and delivers excellent scaling up to at least 1500 cores. Boundary conditions for gravity can be either isolated or periodic, and they can be specified in each direction independently, using a newly developed generalisation of the Ewald method. The gravity calculation can be accelerated with the {\em adaptive block update} technique by partially re-using the solution from the previous time-step. Comparison with the {\sc Flash} internal multi-grid gravity solver shows that tree based methods provide a competitive alternative, particularly for problems with isolated or mixed boundary conditions. We evaluate several multipole acceptance criteria (MACs) and identify a relatively simple APE MAC which provides high accuracy at low computational cost. The optical depth estimates are found to agree very well with those of the {\sc RADMC-3D} radiation transport code, with the tree solver being much faster. Our algorithm is available in the standard release of the {\sc FLASH} code in version 4.0 and later.

Primordial clouds are supposed to host the so-called population III stars. These stars are very massive and completely metal-free. The final stage of the life of population III stars with masses between 130 and 260 solar masses is a very energetic hypernova explosion. A hypernova drives a shock, behind which a spherically symmetric very dense supershell forms, which might become gravitationally unstable, fragment, and form stars. In this paper we study under what conditions can an expanding supershell become gravitationally unstable and how the feedback of these supershell stars (SSSs) affects its surroundings. We simulate, by means of a 1-D Eulerian hydrocode, the early evolution of the primordial cloud after the hypernova explosion, the formation of SSSs, and the following evolution, once the SSSs start to release energy and heavy elements into the interstellar medium. Our results indicate that a shell, enriched with nucleosynthetic products from SSSs, propagates inwards, towards the center of the primordial cloud. In a time span of a few Myr, this inward-propagating shell reaches a distance of only a few parsec away from the center of the primordial cloud. Its density is extremely high and its temperature very low, thus the conditions for a new episode of star formation are achieved. We study what fraction of these two distinct populations of stars can remain bound and survive until the present day. We study also under what conditions can this process repeat and form multiple stellar populations. We extensively discuss whether the proposed scenario can help to explain some open questions of the formation mechanism of globular clusters.

We present 3D "zoom-in" simulations of the formation of two molecular clouds out of the galactic interstellar medium. We model the clouds - identified from the SILCC simulations - with a resolution of up to 0.06 pc using adaptive mesh refinement in combination with a chemical network to follow heating, cooling, and the formation of H$_2$ and CO including (self-) shielding. The two clouds are assembled within a few million years with mass growth rates of up to $\sim$ 10$^{-2}$ M$_\mathrm{sun}$ yr$^{-1}$ and final masses of $\sim$ 50 000 M$_\mathrm{sun}$. A spatial resolution of $\lesssim$ 0.1 pc is required for convergence with respect to the mass, velocity dispersion, and chemical abundances of the clouds, although these properties also depend on the cloud definition such as based on density thresholds, H$_2$ or CO mass fraction. To avoid grid artefacts, the progressive increase of resolution has to occur within the free-fall time of the densest structures (1 - 1.5 Myr) and $\gtrsim$ 200 time steps should be spent on each refinement level before the resolution is progressively increased further. This avoids the formation of spurious, large-scale, rotating clumps from unresolved turbulent flows. While CO is a good tracer for the evolution of dense gas with number densities $n \geq$ 300 cm$^{-3}$, H$_2$ is also found for $n \lesssim 30$ cm$^{-3}$ due to turbulent mixing and becomes dominant at column densities around 30 - 50 M$_\mathrm{sun}$ pc$^{-2}$. The CO-to-H$_2$ ratio steadily increases within the first 2 Myr whereas $X_\mathrm{CO}\simeq$ 1 - 4 $\times$ 10$^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$ is approximately constant since the CO(1-0) line quickly becomes optically thick.

We present a physically-motivated model involving the different processes affecting supernova dust grains as they are incorporated into the thermalized medium within young massive star clusters. The model is used to explain the near- to mid-infrared (NIR-MIR) excess found in such clusters and usually modeled as a blackbody with temperature $\sim (400-1000)$ K. In our approach, dust grains are efficiently produced in the clumpy ejecta of core-collapse supernovae, fragmented into small pieces ($\lesssim 0.05$ $\mu$m) as they are incorporated into the hot and dense ISM, heated via frequent collisions with electrons and the absorption of energetic photons. Grains with small sizes can more easily acquire the high temperatures ($\sim 1000$ K) required to produce a NIR-MIR excess with respect to the emission of foreground PAHs and starlight. However, the extreme conditions inside young massive clusters make difficult for these small grains to have a persistent manifestation at NIR-MIR wavelengths as they are destroyed by efficient thermal sputtering. Nevertheless, the chances for a persistent manifestation are increased by taking into account that small grains become increasingly transparent to their impinging ions as their size decreases. For an individual SN event, we find that the NIR-MIR excess last longer if the time required to incorporate all the grains into the thermalized medium is also longer, and in some cases, comparable to the characteristic interval between supernova explosions. Our models, can successfully explain the near-infrared excesses found in the star clusters observed in M33 (Rela\~no et al. 2016) assuming a low heating efficiency and mass-loading. In this scenario, the presence of the NIR-MIR excess is an indication of efficient dust production in SNe and its subsequent destruction

Isolated HI clouds with no optical counterparts are often taken as evidence for galaxy-galaxy interactions, though an alternative hypothesis is that these are primordial 'dark galaxies' which have not formed stars. Similarly, certain kinematic features in HI streams are also controversial, sometimes taken as evidence of dark galaxies but also perhaps explicable as the result of harassment. We numerically model the passage of a galaxy through the gravitational field of cluster. The galaxy consists of SPH particles for the gas and n-bodies for the stars and dark matter, while the cluster includes the gravitational effects of substructure using 400 subhalos (the effects of the intracluster medium are ignored). We find that harassment can indeed produce long HI streams and these streams can include kinematic features resembling dark galaxy candidates such as VIRGOHI21. We also show that apparent clouds with diameter < 20 kpc and velocity widths < 50 km/s are almost invariably produced in these simulations, making tidal debris a highly probable explanation. In contrast, we show that the frequency of isolated clouds of the same size but velocity width > 100 km/s is negligible - making this a very unlikely explanation for the observed clouds in the Virgo cluster with these properties.

We investigate, using 3D hydrodynamic simulations, the fragmentation of pressure-confined, vertically stratified, self-gravitating gaseous layers. The confining pressure is either thermal pressure acting on both surfaces, or thermal pressure acting on one surface and ram-pressure on the other. In the linear regime of fragmentation, the dispersion relation we obtain agrees well with that derived by Elmegreen & Elmegreen (1978), and consequently deviates from the dispersion relations based on the thin shell approximation (Vishniac 1983) or pressure assisted gravitational instability (W\"unsch et al. 2010). In the non-linear regime, the relative importance of the confining pressure to the self-gravity is a crucial parameter controlling the qualitative course of fragmentation. When confinement of the layer is dominated by external pressure, self- gravitating condensations are delivered by a two-stage process: first the layer fragments into gravitationally bound but stable clumps, and then these clumps coalesce until they assemble enough mass to collapse. In contrast, when external pressure makes a small contribution to confinement of the layer, the layer fragments monolithically into gravitationally unstable clumps and there is no coalescence. This dichotomy persists whether the external pressure is thermal or ram. We apply these results to fragments forming in a shell swept up by an expanding H II region, and find that, unless the swept up gas is quite hot or the surrounding medium has low density, the fragments have low-mass ( ~< 3 M_Sun ), and therefore they are unlikely to spawn stars that are sufficiently massive to promote sequential self-propagating star formation.

We study a model of rapidly cooling shocked stellar winds in young massive clusters and estimate the circumstances under which secondary star formation, out of the reinserted winds from a first stellar generation (1G), is possible. We have used two implementations of the model: a highly idealized computationally inexpensive spherically symmetric semi-analytic model, and a complex three-dimensional radiation-hydrodynamic simulations, and they are in a good mutual agreement. The results confirm our previous findings that in a cluster with 1G mass $10^7$ M$_\odot$ and half-mass radius $2.38$ pc, the shocked stellar winds become thermally unstable, collapse into dense gaseous structures that partially accumulate inside the cluster, self-shield against ionizing stellar radiation and form the second generation (2G) of stars. We have used the semi-analytic model to explore a subset of the parameter space covering a wide range of the observationally poorly constrained parameters: the heating efficiency, $\eta_\mathrm{he}$, and the mass loading, $\eta_\mathrm{ml}$. The results show that the fraction of the 1G stellar winds accumulating inside the cluster can be larger than $50$ % if $\eta_\mathrm{he} \lesssim 10$ % which is suggested by the observations. Furthermore, for low $\eta_\mathrm{he}$, the model provides a self-consistent mechanism predicting 2G stars forming only in the central zones of the cluster. Finally, we have calculated the accumulated warm gas emission in the H30$\alpha$ recombination line, analyzed its velocity profile and estimated its intensity for super star clusters in interacting galaxies NGC4038/9 (Antennae) showing that the warm gas should be detectable with ALMA.

We present three-dimensional magneto-hydrodynamical simulations of the self-gravitating interstellar medium (ISM) in a periodic (256 pc)$^3$ box with a mean number density of 0.5 cm$^{-3}$. At a fixed supernova rate we investigate the multi-phase ISM structure, H$_{2}$ molecule formation and density-magnetic field scaling for varying initial magnetic field strengths (0, $6\times 10^{-3}$, 0.3, 3 $\mu$G). All magnetic runs saturate at mass weighted field strengths of $\sim$ 1 $-$ 3 $\mu$G but the ISM structure is notably different. With increasing initial field strengths (from $6\times 10^{-3}$ to 3 $\mu$G) the simulations develop an ISM with a more homogeneous density and temperature structure, with increasing mass (from 5% to 85%) and volume filling fractions (from 4% to 85%) of warm (300 K $<$ T $<$ 8000 K) gas, with decreasing volume filling fractions (VFF) from $\sim$ 35% to $\sim$ 12% of hot gas (T $> 10^5$ K) and with a decreasing H$_{2}$ mass fraction (from 70% to $<$ 1%). Meanwhile the mass fraction of gas in which the magnetic pressure dominates over the thermal pressure increases by a factor of 10, from 0.07 for an initial field of $6\times 10^{-3}$ $\mu$G to 0.7 for a 3 $\mu$G initial field. In all but the simulations with the highest initial field strength self-gravity promotes the formation of dense gas and H$_{2}$, but does not change any other trends. We conclude that magnetic fields have a significant impact on the multi-phase, chemical and thermal structure of the ISM and discuss potential implications and limitations of the model.

We present three-dimensional radiation-hydrodynamical simulations of the impact of stellar winds, photoelectric heating, photodissociating and photoionising radiation, and supernovae on the chemical composition and star formation in a stratified disc model. This is followed with a sink-based model for star clusters with populations of individual massive stars. Stellar winds and ionising radiation regulate the star formation rate at a factor of ~10 below the simulation with only supernova feedback due to their immediate impact on the ambient interstellar medium after star formation. Ionising radiation (with winds and supernovae) significantly reduces the ambient densities for most supernova explosions to rho < 10^-25 g cm^-3, compared to 10^-23 g cm^-3 for the model with only winds and supernovae. Radiation from massive stars reduces the amount of molecular hydrogen and increases the neutral hydrogen mass and volume filling fraction. Only this model results in a molecular gas depletion time scale of 2 Gyr and shows the best agreement with observations. In the radiative models, the Halpha emission is dominated by radiative recombination as opposed to collisional excitation (the dominant emission in non-radiative models), which only contributes ~1-10 % to the total Halpha emission. Individual massive stars (M >= 30 M_sun) with short lifetimes are responsible for significant fluctuations in the Halpha luminosities. The corresponding inferred star formation rates can underestimate the true instantaneous star formation rate by factors of ~10.

We study the impact of stellar winds and supernovae on the multi-phase interstellar medium using three-dimensional hydrodynamical simulations carried out with FLASH. The selected galactic disc region has a size of (500 pc)$^2$ x $\pm$ 5 kpc and a gas surface density of 10 M$_{\odot}$/pc$^2$. The simulations include an external stellar potential and gas self-gravity, radiative cooling and diffuse heating, sink particles representing star clusters, stellar winds from these clusters which combine the winds from indi- vidual massive stars by following their evolution tracks, and subsequent supernova explosions. Dust and gas (self-)shielding is followed to compute the chemical state of the gas with a chemical network. We find that stellar winds can regulate star (cluster) formation. Since the winds suppress the accretion of fresh gas soon after the cluster has formed, they lead to clusters which have lower average masses (10$^2$ - 10$^{4.3}$ M$_{\odot}$) and form on shorter timescales (10$^{-3}$ - 10 Myr). In particular we find an anti-correlation of cluster mass and accretion time scale. Without winds the star clusters easily grow to larger masses for ~5 Myr until the first supernova explodes. Overall the most massive stars provide the most wind energy input, while objects beginning their evolution as B-type stars contribute most of the supernova energy input. A significant outflow from the disk (mass loading $\gtrsim$ 1 at 1 kpc) can be launched by thermal gas pressure if more than 50% of the volume near the disc mid-plane can be heated to T > 3x10$^5$ K. Stellar winds alone cannot create a hot volume-filling phase. The models which are in best agreement with observed star formation rates drive either no outflows or weak outflows.

Massive stars in young massive clusters insert tremendous amounts of mass and energy into their surroundings in the form of stellar winds and supernova ejecta. Mutual shock-shock collisions lead to formation of hot gas, filling the volume of the cluster. The pressure of this gas then drives a powerful cluster wind. However, it has been shown that if the cluster is massive and dense enough, it can evolve in the so--called bimodal regime, in which the hot gas inside the cluster becomes thermally unstable and forms dense clumps which are trapped inside the cluster by its gravity. We will review works on the bimodal regime and discuss the implications for the formation of subsequent stellar generations. The mass accumulates inside the cluster and as soon as a high enough column density is reached, the interior of the clumps becomes self-shielded against the ionising radiation of stars and the clumps collapse and form new stars. The second stellar generation will be enriched by products of stellar evolution from the first generation, and will be concentrated near the cluster center.

Young and massive star clusters above a critical mass form thermally unstable clumps reducing locally the temperature and pressure of the hot 10$^{7}$~K cluster wind. The matter reinserted by stars, and mass loaded in interactions with pristine gas and from evaporating circumstellar disks, accumulate on clumps that are ionized with photons produced by massive stars. We discuss if they may become self-shielded when they reach the central part of the cluster, or even before it, during their free fall to the cluster center. Here we explore the importance of heating efficiency of stellar winds.

We present a hydrodynamical simulation of the turbulent, magnetized, supernova (SN)-driven interstellar medium (ISM) in a stratified box that dynamically couples the injection and evolution of cosmic rays (CRs) and a self-consistent evolution of the chemical composition. CRs are treated as a relativistic fluid in the advection-diffusion approximation. The thermodynamic evolution of the gas is computed using a chemical network that follows the abundances of H+, H, H2, CO, C+, and free electrons and includes (self-)shielding of the gas and dust. We find that CRs perceptibly thicken the disk with the heights of 90% (70%) enclosed mass reaching ~1.5 kpc (~0.2 kpc). The simulations indicate that CRs alone can launch and sustain strong outflows of atomic and ionized gas with mass loading factors of order unity, even in solar neighborhood conditions and with a CR energy injection per SN of 10^50 erg, 10% of the fiducial thermal energy of an SN. The CR-driven outflows have moderate launching velocities close to the midplane (~100 km/s) and are denser (\rho~1e-24 - 1e-26 g/cm^3), smoother, and colder than the (thermal) SN-driven winds. The simulations support the importance of CRs for setting the vertical structure of the disk as well as the driving of winds.

The SILCC project (SImulating the Life-Cycle of molecular Clouds) aims at a more self-consistent understanding of the interstellar medium (ISM) on small scales and its link to galaxy evolution. We present three-dimensional (magneto)hydrodynamic simulations of the ISM in a vertically stratified box including self-gravity, an external potential due to the stellar component of the galactic disc, and stellar feedback in the form of an interstellar radiation field and supernovae (SNe). The cooling of the gas is based on a chemical network that follows the abundances of H+, H, H2, C+, and CO and takes shielding into account consistently. We vary the SN feedback by comparing different SN rates, clustering and different positioning, in particular SNe in density peaks and at random positions, which has a major impact on the dynamics. Only for random SN positions the energy is injected in sufficiently low-density environments to reduce energy losses and enhance the effective kinetic coupling of the SNe with the gas. This leads to more realistic velocity dispersions (\sigma_HI ~ 0.8\sigma_(300-8000K) ~ 10-20km/s, \sigma_H\alpha ~ 0.6\sigma_(8000-3e5K) ~ 20-30km/s), and strong outflows with mass loading factors of up to 10 even for solar neighbourhood conditions. Clustered SNe abet the onset of outflows compared to individual SNe but do not influence the net outflow rate. The outflows do not contain any molecular gas and are mainly composed of atomic hydrogen. The bulk of the outflowing mass is dense (\rho ~ 1e-25-1e-24g/cc) and slow (v ~ 20-40km/s) but there is a high-velocity tail of up to v ~ 500km/s with \rho ~ 1e-28-1e-27g/cc.

The halo of the Milky Way contains a hot plasma with a surface brightness in soft X-rays of the order $10^{-12}$erg cm$^{-2}$ s$^{-1}$ deg$^{-2}$. The origin of this gas is unclear, but so far numerical models of galactic star formation have failed to reproduce such a large surface brightness by several orders of magnitude. In this paper, we analyze simulations of the turbulent, magnetized, multi-phase interstellar medium including thermal feedback by supernova explosions as well as cosmic-ray feedback. We include a time-dependent chemical network, self-shielding by gas and dust, and self-gravity. Pure thermal feedback alone is sufficient to produce the observed surface brightness, although it is very sensitive to the supernova rate. Cosmic rays suppress this sensitivity and reduce the surface brightness because they drive cooler outflows. Self-gravity has by far the largest effect because it accumulates the diffuse gas in the disk in dense clumps and filaments, so that supernovae exploding in voids can eject a large amount of hot gas into the halo. This can boost the surface brightness by several orders of magnitude. Although our simulations do not reach a steady state, all simulations produce surface brightness values of the same order of magnitude as the observations, with the exact value depending sensitively on the simulation parameters. We conclude that star formation feedback alone is sufficient to explain the origin of the hot halo gas, but measurements of the surface brightness alone do not provide useful diagnostics for the study of galactic star formation.

We use hydrodynamical simulations in a $(256\;{\rm pc})^3$ periodic box to model the impact of supernova (SN) explosions on the multi-phase interstellar medium (ISM) for initial densities $n = 0.5-30$ cm$^{-3}$ and SN rates $1-720$ Myr$^{-1}$. We include radiative cooling, diffuse heating, and the formation of molecular gas using a chemical network. The SNe explode either at random positions, at density peaks, or both. We further present a model combining thermal energy for resolved and momentum input for unresolved SNe. Random driving at high SN rates results in hot gas ($T\gtrsim 10^6$ K) filling $> 90$% of the volume. This gas reaches high pressures ($10^4 < P/k_\mathrm{B} < 10^7$ K cm$^{-3}$) due to the combination of SN explosions in the hot, low density medium and confinement in the periodic box. These pressures move the gas from a two-phase equilibrium to the single-phase, cold branch of the cooling curve. The molecular hydrogen dominates the mass ($>50$%), residing in small, dense clumps. Such a model might resemble the dense ISM in high-redshift galaxies. Peak driving results in huge radiative losses, producing a filamentary ISM with virtually no hot gas, and a small molecular hydrogen mass fraction ($\ll 1$%). Varying the ratio of peak to random SNe yields ISM properties in between the two extremes, with a sharp transition for equal contributions. The velocity dispersion in HI remains $\lesssim 10$ km s$^{-1}$ in all cases. For peak driving the velocity dispersion in H$_\alpha$ can be as high as $70$ km s$^{-1}$ due to the contribution from young, embedded SN remnants.

The SILCC project (SImulating the Life-Cycle of molecular Clouds) aims at a more self-consistent understanding of the interstellar medium (ISM) on small scales and its link to galaxy evolution. We simulate the evolution of the multi-phase ISM in a 500 pc x 500 pc x 10 kpc region of a galactic disc, with a gas surface density of $\Sigma_{_{\rm GAS}} = 10 \;{\rm M}_\odot/{\rm pc}^2$. The Flash 4.1 simulations include an external potential, self-gravity, magnetic fields, heating and radiative cooling, time-dependent chemistry of H$_2$ and CO considering (self-) shielding, and supernova (SN) feedback. We explore SN explosions at different (fixed) rates in high-density regions (peak), in random locations (random), in a combination of both (mixed), or clustered in space and time (clustered). Only random or clustered models with self-gravity (which evolve similarly) are in agreement with observations. Molecular hydrogen forms in dense filaments and clumps and contributes 20% - 40% to the total mass, whereas most of the mass (55% - 75%) is in atomic hydrogen. The ionised gas contributes <10%. For high SN rates (0.5 dex above Kennicutt-Schmidt) as well as for peak and mixed driving the formation of H$_2$ is strongly suppressed. Also without self-gravity the H$_2$ fraction is significantly lower ($\sim$ 5%). Most of the volume is filled with hot gas ($\sim$90% within $\pm$2 kpc). Only for random or clustered driving, a vertically expanding warm component of atomic hydrogen indicates a fountain flow. Magnetic fields have little impact on the final disc structure. However, they affect dense gas ($n\gtrsim 10\;{\rm cm}^{-3}$) and delay H$_2$ formation. We highlight that individual chemical species, in particular atomic hydrogen, populate different ISM phases and cannot be accurately accounted for by simple temperature-/density-based phase cut-offs.

We present simulations of initially stable isothermal clouds exposed to ionizing radiation from a discrete external source, and identify the conditions that lead to radiatively driven implosion and star formation. We use the Smoothed Particle Hydrodynamics code SEREN (Hubber et al. 2010) and the HEALPix-based photoionization algorithm described in Bisbas et al. (2009). We find that the incident ionizing flux is the critical parameter determining the evolution: high fluxes simply disperse the cloud, whereas low fluxes trigger star formation. We find a clear connection between the intensity of the incident flux and the parameters of star formation.

We present 1D numerical simulations, which consider the effects of radiative cooling and gravity on the hydrodynamics of the matter reinserted by stellar winds and supernovae within young nuclear starbursts with a central supermassive black hole (SMBH). The simulations confirm our previous semi-analytic results for low energetic starbursts, evolving in a quasi-adiabatic regime, and extend them to more powerful starbursts evolving in the catastrophic cooling regime. The simulations show a bimodal hydrodynamic solution in all cases. They present a quasi-stationary accretion flow onto the black hole, defined by the matter reinserted by massive stars within the stagnation volume and a stationary starburst wind, driven by the high thermal pressure acquired in the region between the stagnation and the starburst radii. In the catastrophic cooling regime, the stagnation radius rapidly approaches the surface of the starburst region, as one considers more massive starbursts. This leads to larger accretion rates onto the SMBH and concurrently to powerful winds able to inhibit interstellar matter from approaching the nuclear starburst. Our self-consistent model thus establishes a direct physical link between the SMBH accretion rate and the nuclear star formation activity of the host galaxy and provides a good upper limit to the accretion rate onto the central black hole.

The mass reinserted by young stars of an emerging massive compact cluster shows a bimodal hydrodynamic behaviour. In the inner part of the cluster, it is thermally unstable, while in its outer parts it forms an out-blowing wind. The chemical homogeneity/inhomogeneity of low/high mass clusters demonstrates the relevance of this solution to the presence of single/multiple stellar populations. We show the consequences that the thermal instability of the reinserted mass has to the galactic super-winds.