By adopting empirical estimates of the Helium enhancement (Delta Y) between consecutive stellar generations for a sample of Galactic globular clusters (GGC), we uniquely constraint the star formation efficiency of each stellar generation in these stellar systems. In our approach, the star formation efficiency is the central factor that links stellar generations as it defines both their stellar mass and the remaining mass available for further star formation, fixing also the amount of matter required to contaminate the next stellar generation. In this way, the star formation efficiency is here shown to be fully defined by the He enhancement between successive stellar generations in a GC. Our approach has also an impact on the evolution of clusters and thus considers the possible loss of stars through evaporation, tidal interactions and stellar evolution. We focus on the present mass ratio between consecutive stellar generations and the present total mass of Galactic globular clusters. Such considerations suffice to determine the relative proportion of stars of consecutive generations that remain today in globular clusters. The latter is also shown to directly depend on the values of Delta Y and thus the He enhancement between consecutive stellar generations in GGC places major constraints on models of star formation and evolution of GC.

The existence of long (> 100 kpc) HI streams and small (< 20 kpc) free-floating HI clouds is well-known. While the formation of the streams has been investigated extensively, and the isolated clouds are often purported to be interaction debris, little research has been done on the formation of optically dark HI clouds that are not part of a larger stream. One possibility is that such features result from the fragmentation of more extended streams, while another idea is that they are primordial, optically dark galaxies. We test the validity of the fragmentation scenario (via harassment) using numerical simulations. In order to compare our numerical models with observations, we present catalogues of both the known long HI streams (42 objects) and free-floating HI clouds suggested as dark galaxy candidates (51 objects). In particular, we investigate whether it is possible to form compact features with high velocity widths (> 100 km/s), similar to observed clouds which are otherwise intriguing dark galaxy candidates. We find that producing such features is possible but extremely unlikely, occurring no more than 0.2% of the time in our simulations. In contrast, we find that genuine dark galaxies could be extremely stable to harassment and remain detectable even after 5 Gyr in the cluster environment (with the important caveat that our simulations only explore harassment and do not yet include the intracluster medium, heating and cooling, or star formation). We also discuss the possibility that such objects could be the progenitors of recently discovered ultra diffuse galaxies.

Context. Bubble N107 was discovered in the infrared emission of dust in the Galactic Plane observed by the Spitzer Space Telescope (GLIMPSE survey: l ~ 51.0 deg, b ~ 0.1 deg). The bubble represents an example of shell-like structures found all over the Milky Way Galaxy. Aims. We aim to analyse the atomic and molecular components of N107, as well as its radio continuum emission. With the help of numerical simulations, we aim to estimate the bubble age and other parameters which cannot be derived directly from observations. Methods. From the observations of the HI (I-GALFA) and 13CO (GRS) lines we derive the bubble's kinematical distance and masses of the atomic and molecular components. With the algorithm DENDROFIND, we decompose molecular material into individual clumps. From the continuum observations at 1420 MHz (VGPS) and 327 MHz (WSRT), we derive the radio flux density and the spectral index. With the numerical code ring, we simulate the evolution of stellar-blown bubbles similar to N107. Results. The total HI mass associated with N107 is 5.4E3 Msun. The total mass of the molecular component (a mixture of cold gasses of H2, CO, He and heavier elements) is 1.3E5 Msun, from which 4.0E4 Msun is found along the bubble border. We identified 49 molecular clumps distributed along the bubble border, with the slope of the clump mass function of -1.1. The spectral index of -0.30 of a strong radio source located apparently within the bubble indicates nonthermal emission, hence part of the flux likely originates in a supernova remnant, not yet catalogued. The numerical simulations suggest N107 is likely less than 2.25 Myr old. Since first supernovae explode only after 3 Myr or later, no supernova remnant should be present within the bubble. It may be explained if there is a supernova remnant in the direction towards the bubble, however not associated with it.

Infrared shells and bubbles are ubiquitous in the Galaxy and can generally be associated with HII regions formed around young, massive stars. In this paper, we use high-resolution 3D SPH simulations to explore the effect of a single O7 star emitting photons at 10^49 1/s and located at the centre of a molecular cloud with mass 10^4 M_sun and radius 6.4 pc; the internal structure of the cloud is characterised by its fractal dimension, D (with 2.0 <= D <= 2.8), and the variance of its (log-normal) density distribution, sigma_0^2 (with 0.36 <= sigma_0^2 <= 1.42). Our study focuses on the morphology of the swept-up cold gas and the distribution and statistics of the resulting star formation. If the fractal dimension is low, the border of the HII region is dominated by extended shell-like structures, and these break up into a small number of massive high-density clumps which then spawn star clusters; star formation occurs relatively quickly, and delivers somewhat higher stellar masses. Conversely, if the fractal dimension is high, the border of the HII region is dominated by a large number of pillars and cometary globules, which contain compact dense clumps and tend to spawn single stars or individual multiple systems; star formation occurs later, the stellar masses are somewhat lower, and the stars are more widely distributed.

Hydrodynamic models for spherically-symmetric winds driven by young stellar clusters with a generalized Schuster stellar density profile are explored. For this we use both semi-analytic models and 1D numerical simulations. We determine the properties of quasi-adiabatic and radiative stationary winds and define the radius at which the flow turns from subsonic into supersonic for all stellar density distributions. Strongly radiative winds diminish significantly their terminal speed and thus their mechanical luminosity is strongly reduced. This also reduces their potential negative feedback into their host galaxy ISM. The critical luminosity above which radiative cooling becomes dominant within the clusters, leading to thermal instabilities which make the winds non-stationary, is determined, and its dependence on the star cluster density profile, core radius and half mass radius is discussed.

This paper presents a hydrodynamical model describing the evolution of the gas reinserted by stars within a rotating young nuclear star cluster (NSC). We explicitly consider the impact of the stellar component to the flow by means of a uniform insertion of mass and energy within the stellar cluster. The model includes the gravity force of the stellar component and a central supermassive black hole (SMBH), and accounts for the heating from the central source of radiation and the radiative cooling of the thermalized gas. By using a set of parameters typical for NSCs and SMBHs in Seyfert galaxies our simulations show that a filamentary/clumpy structure is formed in the inner part of the cluster. This "torus" is Compton thick and covers a large fraction of the sky (as seen from the SMBH). In the outer parts of the cluster a powerful wind is produced, that inhibits the infall of matter from larger scales and thus the NSC-SMBH interplay occurs in isolation.

The role of feedback from massive stars is believed to be a key element in the evolution of molecular clouds. We use high-resolution 3D SPH simulations to explore the dynamical effects of a single O7 star located at the centre of a molecular cloud with mass 10^4M_sun and radius 6.4pc. The initial internal structure of the cloud is characterised by its fractal dimension, D=2.0 - 2.8, and its log-normal density PDF. (i) As regards star formation, in the short term ionising feedback is positive, in the sense that star formation occurs much more quickly in gas that is compressed by the high pressure of the ionised gas. However, in the long term ionising feedback is negative, in the sense that most of the cloud is dispersed with an outflow rate of up to ~0.01M_sun/yr, on a timescale comparable with the sound-crossing time for the ionised gas (~1-2Myr), and triggered star formation is therefore limited to a few percent of the cloud's mass. (ii) As regards the morphology of the ionisation fronts (IFs) bounding the HII region and the systematics of outflowing gas, we distinguish two regimes. For low D<=2.2, the initial cloud is dominated by large-scale structures, so the neutral gas tends to be swept up into a few extended coherent shells, and the ionised gas blows out through a few large holes between these shells; we term these HII regions "shell-dominated". Conversely, for high D>=2.6, the initial cloud is dominated by small-scale structures, and these are quickly overrun by the advancing IF, thereby producing neutral pillars whilst the ionised gas blows out through a large number of small holes between the pillars; we term these HII regions "pillar-dominated". (iii) As regards the injection of bulk kinetic energy, by ~1Myr, the expansion of the HII region has delivered a rms velocity of ~6km/s; this represents less than 0.1% of the total energy radiated by the O7 star.

$^{13}$CO(J=2--1) and C$^{18}$O(J=2--1) observations of the molecular cloud G285.90+4.53 (Cloud~16) in the Carina Flare supershell (GSH287+04-17) with the APEX telescope are presented. With an algorithm DENDROFIND we identify 51 fragments and compute their sizes and masses. We discuss their mass spectrum and interpret it as being the result of the shell fragmentation process described by the pressure assisted gravitational instability - PAGI. We conclude that the explanation of the clump mass function needs a combination of gravity with pressure external to the shell.

We study the evolution of Super Star Cluster (SSC) winds driven by stellar winds and supernova (SN) explosions. Time-dependent rates at which mass and energy are deposited into the cluster volume, as well as the time-dependent chemical composition of the re-inserted gas, are obtained from the population synthesis code Starburst99. These results are used as input for a semi-analytic code which determines the hydrodynamic properties of the cluster wind as a function of cluster age. Two types of winds are detected in the calculations. For the quasi-adiabatic solution, all of the inserted gas leaves the cluster in the form of a stationary wind. For the bimodal solution, some of the inserted gas becomes thermally unstable and forms dense warm clumps which accumulate inside the cluster. We calculate the evolution of the wind velocity and energy flux and integrate the amount of accumulated mass for clusters of different mass, radius and initial metallicity. We consider also conditions with low heating efficiency of the re-inserted gas or mass loading of the hot thermalized plasma with the gas left over from star formation. We find that the bimodal regime and the related mass accumulation occur if at least one of the two conditions above is fulfilled.

We present simulations of initially stable isothermal clouds exposed to ionising 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 and an HEALPix-based photo-ionisation algorithm to simulate the propagation of the ionising radiation and the resulting dynamical evolution of the cloud. We find that the incident ionising flux, $\Phi_{_{\rm LyC}}$, is the critical parameter determining the cloud evolution. At moderate fluxes, a large fraction of the cloud mass is converted into stars. As the flux is increased, the fraction of the cloud mass that is converted into stars and the mean masses of the individual stars both decrease. Very high fluxes simply disperse the cloud. Newly-formed stars tend to be concentrated along the central axis of the cloud (i.e. the axis pointing in the direction of the incident flux). For given cloud parameters, the time, $t_{_\star}$, at which star formation starts is proportional to $\Phi_{_{\rm LyC}}^{-1/3}$. The pattern of star formation found in the simulations is similar to that observed in bright-rimmed clouds.

We use SPH simulations to investigate the gravitational fragmentation of expanding shells through the linear and non--linear regimes. The results are analysed using spherical harmonic decomposition to capture the initiation of structure during the linear regime; the potential-based method of Smith et al. (2009) to follow the development of clumps in the mildly non-linear regime; and sink particles to capture the properties of the final bound objects during the highly non-linear regime. In the early, mildly non--linear phase of fragmentation, we find that the clump mass function still agrees quite well with the mass function predicted by the analytic model. However, the sink mass function is quite different, in the sense of being skewed towards high-mass objects. This is because, once the growth of a condensation becomes non-linear, it tends to be growing non-competitively from its own essentially separate reservoir; we call this Oligarchic Accretion.

We study analytically the development of gravitational instability in an expanding shell having finite thickness. We consider three models for the radial density profile of the shell: (i) an analytic uniform-density model, (ii) a semi-analytic model obtained by numerical solution of the hydrostatic equilibrium equation, and (iii) a 3D hydrodynamic simulation. We show that all three profiles are in close agreement, and this allows us to use the first model to describe fragments in the radial direction of the shell. We then use non-linear equations describing the time-evolution of a uniform oblate spheroid to derive the growth rates of shell fragments having different sizes. This yields a dispersion relation which depends on the shell thickness, and hence on the pressure confining the shell. We compare this dispersion relation with the dispersion relation obtained using the standard thin-shell analysis, and show that, if the confining pressure is low, only large fragments are unstable. On the other hand, if the confining pressure is high, fragments smaller than predicted by the thin-shell analysis become unstable. Finally, we compare the new dispersion relation with the results of 3D hydrodynamic simulations, and show that the two are in good agreement.

We present a detailed study of the hydrodynamics of the matter reinserted by massive stars via stellar winds and supernovae explosions in young assembling galaxies. We show that the interplay between the thermalization of the kinetic energy provided by massive stars, radiative cooling of the thermalized plasma and the gravitational pull of the host galaxy, lead to three different hydrodynamic regimes. These are: a) The quasi-adiabatic supergalactic winds. b) The bimodal flows, with mass accumulation in the central zones and gas expulsion from the outer zones of the assembling galaxy. c) The gravitationally bound regime, for which all of the gas returned by massive stars remains bound to the host galaxy and is likely to be reprocessed into futher generations of stars. Which of the three possible solutions takes place, depends on the mass of the star forming region its mechanical luminosity (or star formation rate) and its size. The model predicts that massive assembling galaxies with large star formation rates similar to those detected in SCUBA sources ($\sim 1000$ M$_\odot$ yr$^{-1}$) are likely to evolve in a positive star-formation feedback condition, either in the bimodal, or in the gravitationally bound regime. This implies that star formation in these sources may have little impact on the intergalactic medium and result instead into a fast interstellar matter enrichment, as observed in high redshift quasars.

The origin of supersonic infrared and radio recombination nebular lines often detected in young and massive superstar clusters are discussed. We suggest that these arise from a collection of repressurizing shocks (RSs), acting effectively to re-establish pressure balance within the cluster volume and from the cluster wind which leads to an even broader although much weaker component. The supersonic lines are here shown to occur in clusters that undergo a bimodal hydrodynamic solution (Tenorio-Tagle et al. 2007), that is within clusters that are above the threshold line in the mechanical luminosity or cluster mass vs the size of the cluster (Silich et al. 2004). The plethora of repressurizing shocks is due to frequent and recurrent thermal instabilities that take place within the matter reinserted by stellar winds and supernovae. We show that the maximum speed of the RSs and of the cluster wind, are both functions of the temperature reached at the stagnation radius. This temperature depends only on the cluster heating efficiency ($\eta$). Based on our two dimensional simulations (Wunsch et al. 2008) we calculate the line profiles that result from several models and confirm our analytical predictions. From a comparison between the predicted and observed values of the half-width zero intensity of the two line components we conclude that the thermalization efficiency in SSC's above the threshold line must be lower than 20%.

We investigate the gravitational fragmentation of expanding shells in the context of the linear thin--shell analysis. We make use of two very different numerical schemes; the FLASH Adaptive Mesh Refinement code and a version of the Benz Smoothed Particle Hydrodynamics code. We find that the agreement between the two codes is excellent. We use our numerical results to test the thin--shell approximation and we find that the external pressure applied to the shell has a strong effect on the fragmentation process. In cases where shells are not pressure--confined, the shells thicken as they expand and hydrodynamic flows perpendicular to the plane of the shell suppress fragmentation at short wavelengths. If the shells are pressure--confined internally and externally, so that their thickness remains approximately constant during their expansion, the agreement with the analytical solution is better.

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.

Using the hydrodynamic code ZEUS, we perform 2D simulations to determine the fate of the gas ejected by massive stars within super star clusters. It turns out that the outcome depends mainly on the mass and radius of the cluster. In the case of less massive clusters, a hot high velocity ($\sim 1000$ km s$^{-1}$) stationary wind develops and the metals injected by supernovae are dispersed to large distances from the cluster. On the other hand, the density of the thermalized ejecta within massive and compact clusters is sufficiently large as to immediately provoke the onset of thermal instabilities. These deplete, particularly in the central densest regions, the pressure and the pressure gradient required to establish a stationary wind, and instead the thermally unstable parcels of gas are rapidly compressed, by a plethora of re-pressurizing shocks, into compact high density condensations. Most of these are unable to leave the cluster volume and thus accumulate to eventually feed further generations of star formation. The simulations cover an important fraction of the parameter-space, which allows us to estimate the fraction of the reinserted gas which accumulates within the cluster and the fraction that leaves the cluster as a function of the cluster mechanical luminosity, the cluster size and heating efficiency.

We investigate numerically the orbital evolution of massive extrasolar planets within central cavities of their parent protoplanetary discs. Assuming that they arrive at the inner edge of the disc due to type II migration, we show that they spiral further in. We find that in magnetospheric cavities more massive planets stop migrating at a larger distance from the edge of the disc. This effect may qualitatively explain the correlation between masses and orbital periods found for massive planets with P shorter than 5 days.

We look for a simple analytic model to distinguish between stellar clusters undergoing a bimodal hydrodynamic solution from those able to drive only a stationary wind. Clusters in the bimodal regime undergo strong radiative cooling within their densest inner regions, which results in the accumulation of the matter injected by supernovae and stellar winds and eventually in the formation of further stellar generations, while their outer regions sustain a stationary wind. The analytic formulae are derived from the basic hydrodynamic equations. Our main assumption, that the density at the star cluster surface scales almost linearly with that at the stagnation radius, is based on results from semi-analytic and full numerical calculations. The analytic formulation allows for the determination of the threshold mechanical luminosity that separates clusters evolving in either of the two solutions. It is possible to fix the stagnation radius by simple analytic expressions and thus to determine the fractions of the deposited matter that clusters evolving in the bimodal regime blow out as a wind or recycle into further stellar generations.

We present semi-analytical and numerical models, accounting for the impact of radiative cooling on the hydrodynamics of the matter reinserted as strong stellar winds and supernovae within the volume occupied by young, massive and compact superstellar clusters. First of all we corroborate the location of the threshold line in the mechanical energy input rate vs the cluster size plane, found by Silich et al. (2004). Such a line separates clusters able to drive a quasi-adiabatic or a strongly radiative wind from clusters in which catastrophic cooling occurs within the star cluster volume. Then we show that the latter, clusters above the threshold line, undergo a bimodal behavior in which the central densest zones cool rapidly and accumulate the injected matter to eventually feed further generations of star formation, while the outer zones are still able to drive a stationary wind. The results are presented into a series of universal dimensionless diagrams from which one can infer: the size of the two zones, the fraction of the deposited mass that goes into each of them and the luminosity of the resultant winds, for clusters of all sizes and energy input rates, regardless the assumed adiabatic terminal speed V_A.

We study axisymmetric models of layered protoplanetary discs taking radiative transfer effects into account, and allowing for a residual viscosity in the dead zone. We also explore the effect of different viscosity prescriptions. In addition to the ring instability reported in the first paper of the series we find an oscillatory instability of the dead zone, accompanied by variations of the accretion rate onto the central star. We provide a simplified analytical description explaining the mechanism of the oscillations. Finally, we find that the residual viscosity enables stationary accretion in large regions of layered discs. Based on results obtained with the help of a simple 1-D hydrocode we identify these regions, and discuss conditions in which layered discs can give rise to FU~Orionis phenomena.

In this work we use the radiation hydrodynamic code TRAMP to perform a two-dimensional axially symmetric model of the layered disc. Using this model we follow the accumulation of mass in the dead zone due to the radially varying accretion rate. We found a new type of instability which causes the dead zone to split into rings. This "ring instability" works due to the positive feedback between the thickness of the dead zone and the mass accumulation rate. We give an analytical description of this instability, taking into account non-zero thickness of the dead zone and deviations from the Keplerian rotational velocity. The analytical model agrees reasonably well with results of numerical simulations. Finally, we speculate about the possible role of the ring instability in protoplanetary discs and in the formation of planets.

The gravitational instability of expanding shells is discussed. Linear and nonlinear terms are included in an analytical solution in the static and homogeneous medium. We discuss the interaction of modes and give the time needed for fragmentation. Masses of individual fragments are also estimated and their formation rates and the initial mass function are derived. Results of simulations are compared to observation.

A model of the thin shell expanding into a uniform ambient medium is developed. Density perturbations are described using equations with linear and quadratic terms, and the linear and the nonlinear solutions are compared. We follow the time evolution of the fragmentation process and separate the well defined fragments. Their mass spectrum is compared to observations and we also estimate their formation time.