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.

We present the first high-resolution map of the cold molecular gas distribution, as traced by CO(2-1) emission with ALMA, in a prominent ram pressure stripped tail. The Norma cluster galaxy ESO 137-001 is undergoing a strong interaction with the surrounding intra-cluster medium and is one of the nearest jellyfish galaxies with a long multi-phase tail. We have mapped the full extent of the tail at 1" (350 pc) angular resolution and found a rich distribution of mostly compact CO regions extending to nearly 60 kpc in length and 25 kpc in width. In total about 10^9 M_sun of molecular gas was detected. The CO features are found predominantly at the heads of numerous small-scale (~ 1.5 kpc) fireballs (i.e., star-forming clouds with linear streams of young stars extending toward the galaxy) but also of large-scale (~ 8 kpc) super-fireballs, and double-sided fireballs that have additional diffuse ionized gas tails extending in the direction opposite to the stellar tails. The new data help to shed light on the origin of the molecular tail - CO filaments oriented in the direction of the tail with only diffuse associated Halpha emission are likely young molecular features formed in situ, whereas other large CO features tilted with respect to the tail may have originated from the densest gas complexes that were pushed gradually away from the disk. The ALMA observations of ESO 137-001, together with observations from HST, Chandra and VLT/MUSE, offer the most complete view of a spectacular ram pressure stripped tail to date.

HI shells and supershells, found in discs of many galaxies including our own, are formed by the activity of young and massive stars (supernova explosions and stellar winds), but the formation of these structures may be linked to other energetic events, such as interactions of high-velocity clouds with the galactic disc. The larger structures in particular significantly influence their surroundings; their walls are often places where molecular clouds reside and where star formation happens. We explore the HI supershell GS242-03+37, a large structure in the outer Milky Way. Its size and position make it a good case for studying the effects of large shells on their surrounding. We perform numerical simulations of the structure with the simplified hydrodynamical code RING, which uses the thin-shell approximation. The best fit is found by a comparison with the HI data and then we compare our model with the distribution of star clusters near this supershell. The best model of GS242-03+37 requires, contrary to previous estimates, a relatively low amount of energy, and it has an old age of $\sim$ 100 Myr. We also find that the distribution of young star clusters (with ages $<$ 120 Myr) is correlated with walls of the supershell, while the distribution of older clusters is not. Clusters that have the highest probability of being born in the wall of the supershell show an age sequence along the wall. GS242-03+37 is a relatively old structure, shaped by the differential rotation, and its wall is a birthplace of several star clusters. The star formation started at a time when the supershell was not already supersonically expanding; it was a result of the density increase due to the galactic shear and oscillations perpendicular to the disc of the Milky Way.

The removal of gas by ram pressure stripping of galaxies is treated by a purely kinematic description. The solution has two asymptotic limits: if the duration of the ram pressure pulse exceeds the period of vertical oscillations perpendicular to the galactic plane, the commonly used quasi-static criterion of Gunn & Gott is obtained which uses the maximum ram pressure that the galaxy has experienced along its orbit. For shorter pulses the outcome depends on the time-integrated ram pressure. This parameter pair fully describes the gas mass fraction that is stripped from a given galaxy. This approach closely reproduces results from SPH simulations. We show that typical galaxies follow a very tight relation in this parameter space corresponding to a pressure pulse length of about 300 Myr. Thus, the Gunn & Gott criterion provides a good description for galaxies in larger clusters. Applying the analytic description to a sample of 232 Virgo galaxies from the GoldMine database, we show that the ICM provides indeed the ram pressures needed to explain the deficiencies. We also can distinguish current and past strippers, including objects whose stripping state was unknown.

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 have discovered large amounts of molecular gas, as traced by CO emission, in the ram pressure stripped gas tail of the Coma cluster galaxy D100 (GMP 2910), out to large distances of about 50 kpc. D100 has a 60 kpc long, strikingly narrow tail which is bright in X-rays and H{\alpha}. Our observations with the IRAM 30m telescope reveal in total ~ 10^9 M_sun of H_2 (assuming the standard CO-to-H_2 conversion) in several regions along the tail, thus indicating that molecular gas may dominate its mass. Along the tail we measure a smooth gradient in the radial velocity of the CO emission that is offset to lower values from the more diffuse H{\alpha} gas velocities. Such a dynamic separation of phases may be due to their differential acceleration by ram pressure. D100 is likely being stripped at a high orbital velocity >2200 km/s by (nearly) peak ram pressure. Combined effects of ICM viscosity and magnetic fields may be important for the evolution of the stripped ISM. We propose D100 has reached a continuous mode of stripping of dense gas remaining in its nuclear region. D100 is the second known case of an abundant molecular stripped-gas tail, suggesting that conditions in the ICM at the centers of galaxy clusters may be favorable for molecularization. From comparison with other galaxies, we find there is a good correlation between the CO flux and the H{\alpha} surface brightness in ram pressure stripped gas tails, over about 2 dex.

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.

HI shells, which may be formed by the activity of young and massive stars, or connected to energy released by interactions of high-velocity clouds with the galactic disk, may be partly responsible both for the destruction of CO clouds and for the creation of others. It is not known which effect prevails. We study the relation between HI shells and CO in the outer parts of the Milky Way, using HI and CO surveys and a catalogue of previously identified HI shells. For each individual location, the distance to the nearest HI shell is calculated and it is specified whether it lies in the interior of an HI shell, in its walls, or outside an HI shell. The method takes into account irregular shapes of HI shells. We find a lack of CO clouds in the interiors of HI shells and their increased occurrence in walls. Properties of clouds differ for different environments: interiors of HI shells, their walls, and unperturbed medium. CO clouds found in the interiors of HI shells are those that survived and were robbed of their more diffuse gas. Walls of HI shells have a high molecular content, indicative of an increased rate of CO formation. Comparing the CO fractions within HI shells and outside in the unperturbed medium, we conclude that HI shells are responsible for approx. 20 % increase in the total amount of CO in the outer Milky Way.

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.

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.

We present IRAM 30m sensitive upper limits on CO emission in the ram pressure stripped dwarf Virgo galaxy IC3418 and in a few positions covering HII regions in its prominent 17 kpc UV/Ha gas-stripped tail. In the central few arcseconds of the galaxy, we report a possible marginal detection of about 1x10^6 M_sun of molecular gas (assuming a Galactic CO-to-H_2 conversion factor) that could correspond to a surviving nuclear gas reservoir. We estimate that there is less molecular gas in the main body of IC3418, by at least a factor of 20, than would be expected from the pre-quenching UV-based star formation rate assuming the typical gas depletion timescale of 2 Gyr. Given the lack of star formation in the main body, we think the H_2-deficiency is real, although some of it may also arise from a higher CO-to-H_2 factor typical in low-metallicity, low-mass galaxies. The presence of HII regions in the tail of IC3418 suggests that there must be some dense gas; however, only upper limits of < 1x10^6 M_sun were found in the three observed points in the outer tail. This yields an upper limit on the molecular gas content of the whole tail < 1x10^7 M_sun, which is an amount similar to the estimates from the observed star formation rate over the tail. We also present strong upper limits on the X-ray emission of the stripped gas in IC3418 from a new Chandra observation. The measured X-ray luminosity of the IC3418 tail is about 280 times lower than that of ESO 137-001, a spiral galaxy in a more distant cluster with a prominent ram pressure stripped tail. Non-detection of any diffuse X-ray emission in the IC3418 tail may be due to a low gas content in the tail associated with its advanced evolutionary state and/or due to a rather low thermal pressure of the surrounding intra-cluster medium.

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.

We analyse the all-sky Leiden/Argentina/Bonn HI survey, where we identify shells belonging to the Milky Way. We used an identification method based on the search of continuous regions of a low brightness temperature that are compatible with given properties of HI shells. We found 333 shells in the whole Galaxy. The size distribution of shells in the outer Galaxy is fitted by a power law with the coefficient of 2.6 corresponding to the index 1.8 in the distribution of energy sources. Their surface density decreases exponentially with a scale length of 2.8 kpc. The surface density of shells with radii >= 100 pc in the solar neighbourhood is around 4 per kpc^2 and the 2D porosity is approximately 0.7.

$^{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 studied the impact of the revisited values for the LSR circular velocity of the Milky Way (Reid et al. 2004) on the formation of the Magellanic Stream. The LSR circular velocity was varied within its observational uncertainties as a free parameter of the interaction between the Large (LMC) and the Small (SMC) Magellanic Clouds and the Galaxy. We have shown that the large-scale morphology and kinematics of the Magellanic Stream may be reproduced as tidal features, assuming the recent values of the proper motions of the Magellanic Clouds (Kallivayalil et al. 2006). Automated exploration of the entire parameter space for the interaction was performed to identify all parameter combinations that allow for modeling the Magellanic Stream. Satisfactory models exist for the dynamical mass of the Milky Way within a wide range of 0.6*10^12Msun to 3.0*10^12Msun and over the entire 1-sigma errors of the proper motions of the Clouds. However, the successful models share a common interaction scenario. The Magellanic Clouds are satellites of the Milky Way, and in all cases two close LMC-SMC encounters occurred within the last 4Gyr at t<-2.5Gyr and t approx. -150Myr, triggering the formation of the Stream and of the Magellanic Bridge, respectively. The latter encounter is encoded in the observed proper motions and inevitable in any model of the interaction. We conclude that the tidal origin of the Magellanic Stream implies the previously introduced LMC/SMC orbital history, unless the parameters of the interaction are revised substantially.

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.

Ram pressure stripping of galaxies in clusters can yield gas deficient disks. Previous numerical simulations based on various approaches suggested that, except for near edge-on disk orientations, the amount of stripping depends very little on the inclination angle. Following our previous study of face-on stripping, we extend the set of parameters with the disk tilt angle and explore in detail the effects of the ram pressure on the interstellar content (ISM) of tilted galaxies that orbit in various environments of clusters, with compact or extended distributions of the intra-cluster medium (ICM). We further study how results of numerical simulations could be estimated analytically. A grid of numerical simulations with varying parameters is produced using the tree/SPH code GADGET with a modified method for calculating the ISM-ICM interaction. These SPH calculations extend the set of existing results obtained from different codes using various numerical techniques. The simulations confirm the general trend of less stripping at orientations close to edge-on. The dependence on the disk tilt angle is more pronounced for compact ICM distributions, however it almost vanishes for strong ram pressure pulses. Although various hydrodynamical effects are present in the ISM-ICM interaction, the main quantitative stripping results appear to be roughly consistent with a simple scenario of momentum transfer from the encountered ICM. This behavior can also be found in previous simulations. To reproduce the numerical results we propose a fitting formula depending on the disk tilt angle and on the column density of the encountered ICM. Such a dependence is superior to that on the peak ram pressure used in previous simple estimates.

Recently, Kallivayalil et al. derived new values of the proper motion for the Large and Small Magellanic Clouds (LMC and SMC, respectively). The spatial velocities of both Clouds are unexpectedly higher than their previous values resulting from agreement between the available theoretical models of the Magellanic System and the observations of neutral hydrogen (HI) associated with the LMC and the SMC. Such proper motion estimates are likely to be at odds with the scenarios for creation of the large-scale structures in the Magellanic System suggested so far. We investigated this hypothesis for the pure tidal models, as they were the first ones devised to explain the evolution of the Magellanic System, and the tidal stripping is intrinsically involved in every model assuming the gravitational interaction. The parameter space for the Milky Way (MW)-LMC-SMC interaction was analyzed by a robust search algorithm (genetic algorithm) combined with a fast restricted N-body model of the interaction. Our method extended the known variety of evolutionary scenarios satisfying the observed kinematics and morphology of the Magellanic large-scale structures. Nevertheless, assuming the tidal interaction, no satisfactory reproduction of the HI data available for the Magellanic Clouds was achieved with the new proper motions. We conclude that for the proper motion data by Kallivayalil et al., within their 1-sigma errors, the dynamical evolution of the Magellanic System with the currently accepted total mass of the MW cannot be explained in the framework of pure tidal models. The optimal value for the western component of the LMC proper motion was found to be pm_w(LMC) > -1.3 mas/yr in case of tidal models. It corresponds to the reduction of the Kallivayalil et al. value for pm_w(LMC) by approx. 40% in its magnitude.

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.

The influence of a time-varying ram pressure on spiral galaxies in clusters is explored with a new simulation method based on the N-body SPH/tree code GADGET. We have adapted the code to describe the interaction of two different gas phases, the diffuse hot intracluster medium (ICM) and the denser and colder interstellar medium (ISM). Both the ICM and ISM components are introduced as SPH particles. As a galaxy arrives on a highly radial orbit from outskirts to cluster center, it crosses the ICM density peak and experiences a time-varying wind. Depending on the duration and intensity of the ISM-ICM interaction, early and late type galaxies in galaxy clusters with either a large or small ICM distribution are found to show different stripping efficiencies, amounts of reaccretion of the extra-planar ISM, and final masses. We compare the numerical results with analytical approximations of different complexity and indicate the limits of the Gunn & Gott simple stripping formula. Our investigations emphasize the role of the galactic orbital history to the stripping amount. We discuss the contribution of ram pressure stripping to the origin of the ICM and its metallicity. We propose gas accumulations like tails, filaments, or ripples to be responsible for stripping in regions with low overall ICM occurrence.

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.

Different hydrodynamic regimes for the gaseous outflows generated by multiple supernovae explosions and stellar winds occurring within compact and massive star clusters are discussed. It is shown that there exists the threshold energy that separates clusters whose outflows evolve in the quasi-adiabatic or radiative regime from those within which catastrophic cooling and a positive feedback star-forming mode sets in. The role of the surrounding ISM and the observational appearance of the star cluster winds evolving in different hydrodynamic regimes are also discussed.

We performed an extended analysis of the parameter space for the interaction of the Magellanic System with the Milky Way (MW). The varied parameters cover the phase space parameters, the masses, the structure, and the orientation of both Magellanic Clouds, as well as the flattening of the dark matter halo of the MW. The analysis was done by a specially adopted optimization code searching for the best match between numerical models and the detailed HI map of the Magellanic System by Bruens et al. (2005). The applied search algorithm is a genetic algorithm combined with a code based on the fast, but approximative restricted N-body method. By this, we were able to analyze more than 10^6 models, which makes this study one of the most extended ones for the Magellanic System. Here we focus on the flattening q of the axially symmetric MW dark matter halo potential, that is studied within the range 0.74<=q<=1.20. We show that creation of a trailing tail (Magellanic Stream) and a leading stream (Leading Arm) is quite a common feature of the Magellanic System-MW interaction, and such structures were modeled across the entire range of halo flattening values. However, important differences exist between the models, concerning density distribution and kinematics of HI, and also the dynamical evolution of the Magellanic System. Detailed analysis of the overall agreement between modeled and observed distribution of neutral hydrogen shows that the models assuming an oblate (q<1.0) dark matter halo of the Galaxy allow for better satisfaction of HI observations than models with other halo configurations.

We present results of a method for an automatic search for HI shells in 3D data cubes and apply it to the Leiden-Dwingeloo HI survey of the northern Milky Way. In the 2nd Galactic quadrant, where identifications of structures are not substantially influenced by overlapping, we find nearly 300 structures. The Galactic distribution of shells has an exponential profile in the radial direction with a scale length of 3 kpc. In the z direction, one half of the shells are found at distances smaller than 500 pc. We also calculate the energies necessary to create the shells: there are several structures with energies greater than 10 E_SN but only one with an energy exceeding 100 E_SN. Their size distribution, corrected for distance effects, is approximated by a power-law with an index 2.1. Our identifications provide a lower limit to the filling factor of shells in the outer Milky Way: f_2D = 0.4 and f_3D = 0.05.

A review of gravitational and hydrodynamical processes during formation of clusters and evolution of galaxies is given. Early, at the advent of N-body computer simulations, the importance of tidal fields in galaxy encounters has been recognized. Orbits are crowded due to tides along spiral arms, where the star formation is enhanced. Low relative velocity encounters lead to galaxy mergers. The central dominating galaxies in future clusters form before the clusters in a merging process in galaxy groups. Galaxy clusters are composed in a hierarchical scenario due to relaxation processes between galaxies and galaxy groups. As soon as the overall cluster gravitational potential is built, high speed galaxy versus galaxy encounters start to play a role. These harassment events gradually thicken and shorten spiral galaxy disks leading to the formation of S0 galaxies and ellipticals. Another aspect of the high speed motion in the hot and diluted Intracluster Medium (ICM) is the ram pressure exerted on the Interstellar Matter (ISM) leading to stripping of the ISM from parent spirals. The combinations of tides and ram pressure stripping efficiently removes the gas from spirals, quenching the star formation in galactic disks, while triggering it in the tidal arms and at the leading edge of gaseous disk. Gas stripping from disks transports the metals to the ICM. In some cases, the gas extracted from the galactic disks becomes self-gravitating forming tidal dwarf galaxies.

Here we model a star forming factory in which the continuous creation of stars results in a highly concentrated, massive (globular cluster-like) stellar system. We show that under very general conditions a large-scale gravitational instability in the ISM, which triggers the collapse of a massive cloud, leads with the aid of a spontaneous first generation of massive stars, to a standing, small-radius, cold and dense shell. Eventually, as more of the collapsing matter is processed and incorporated, the shell becomes gravitationally unstable and begins to fragment, allowing the formation of new stars, while keeping its location. This is due to a detailed balance established between the ram pressure from the collapsing cloud which, together with the gravitational force exerted on the shell by the forming cluster, acts against the mechanical energy deposited by the collection of new stars. We analyze the mass spectrum of fragments that result from the continuous fragmentation of the standing shell and show that its shape is well approximated at the high mass end by a power law with slope -2.25, very close to the value that fits the universal IMF. The self-contamination resultant from the continuous generation of stars is shown to lead to a large metal spread in massive ($\sim$ 10$^6$ M$_\odot$) clusters, while clusters with a mass similar to 10$^5$ M$_\odot$ or smaller, simply reflect the initial metalicity of the collapsing cloud. This is in good agreement with the data available for globular clusters in the Galaxy.

Criteria for gravitational collapse of expanding shells in rotating, shearing galaxy disks were determined using three-dimensional numerical simulations in the thin shell approximation. The simulations were run over a grid of 7 independent variables, and the resultant probabilities for triggering and unstable masses were determined as functions of 8 dimensionless parameters. When the ratio of the midplane gas density to the midplane total density is small, an expanding shell reaches the disk scale height and vents to the halo before it collapses. When the Toomre instability parameter Q, or a similar shear parameter, Q_A, are large, Coriolis forces and shear stall or reverse the collapse before the shell accumulates enough mass to be unstable. With large values of C=c_sh/(GL)^0.2, for rms velocity dispersion c_sh in the swept-up matter and shell-driving luminosity L, the pressure in the accumulated gas is too large to allow collapse during the expansion time. Considering ~5000 models covering a wide range of parameter space, the common properties of shell collapse as a mechanism for triggered star formation are: (1) the time scale is 4*sqrt(C/2 pi G rho) for ambient midplane density rho, (2) the total fragment mass is ~2x10^7 Msun, of which only a small fraction is likely to be molecular, (3) the triggering radius is ~2 times the scale height, and the triggering probability is ~0.5 for large OB associations. Star formation triggered by shell collapse should be most common in gas-rich galaxies, such as young galaxies or those with late Hubble types.

The gravitational instability of expanding shells triggering the formation of clouds and stars is analyzed. Disks with different scale-heights, ambient and shell velocity dispersions, mid-plane densities, rotation rates and shear rates are explored with three dimensional numerical simulations in the thin shell approximation. Three conditions for the shell collapse are specified: the first is that it happens before a significant blow-out, the second requires that the shell collapses before it is distorted by Coriolis forces and shear, and the third requires that the internal pressure in the accumulated gas is small and the fragmentation is achieved within the expansion time. The gas-rich and slowly rotating galaxies are the best sites of the triggered star formation, concluding that its importance has been much larger at the times of galaxy formation compared to the present epoch.

We discuss fragmentation processes which induce star formation in dense walls of expanding shells. The influence of the energy input, the ISM scale-height and speed of sound in the ambient medium is tested. We formulate the condition for the gravitational fragmentation of expanding shells: if the total surface density of the disc is higher than a certain critical value, shells are unstable. The value of the critical density depends on the energy of the shell and the sound speed in the ISM.

We present an N-body computer code - aimed at studies of galactic dynamics - with a CPU-efficient algorithm for a continuous (i.e. time-dependent) stellar mass-loss. First, we summarize available data on stellar mass-loss and derive the long-term (20 Gyr) dependence of mass-loss rate of a coeval stellar population. We then implement it, through a simple parametric form, into a particle-mesh code with stellar and gaseous particles. We perform several tests of the algorithm reliability and show an illustrative application: a 2D simulation of a disk galaxy, starting as purely stellar but evolving as two-component due to gradual mass-loss from initial stars and due to star formation. In a subsequent paper we will use the code to study what changes are induced in galactic disks by the continuous gas recycling compared to the instantaneous recycling approximation, especially the changes in star formation rate and radial inflow of matter.

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.

We describe HI observations of a 4 x 4 degrees field in the Milky Way centered on l=61 deg, b=0 deg made by the Effelsberg radiotelescope. The field contains one previously identified HI supershell, GS061+00+51 (Heiles, 1979); apart from it we find several new structures. We also study the HI distribution in the vicinity of four HII regions, S86, S87, S88 and S89. We confirm the existence of the shell GS061+00+51, and we find that it has two smaller neighbours, spherical shells with a radius around 30 pc. We identify at least one more regular shell at v_LSR = -18 km/s; and one blown-out shell at v_LSR = -54 km/s. In two cases we are able to connect HII regions with features in the HI distribution (S86 and S87), in two other cases no connection is found. Apart from quite regular HI shells we see a number of non-coherent objects, which are probably a result of the turbulence in the interstellar medium.

The star formation triggered in dense walls of expanding shells will be discussed. The fragmentation process is studied using the linear and non-linear perturbation theory. The influence of the energy input, the ISM distribution and the speed of sound is examined analytically and by numerical simulations. We formulate the condition for the gravitational fragmentation of expanding shells: if the total surface density of the disc is higher than a certain critical value, shells are unstable. This value depends on the energy of the shell and the sound speed in the ISM. As an example the formation of OB associations near the Sun will be discussed. We trace their orbits in the Milky Way to see where they have been born: 10 - 12 Myr ago progenitors of Scorpius-Centaurus OB associations and the Orion OB association resided together within a sheet-like region elongated in the $l = 20-200degrees direction, showing that the local OB associations may be formed as fragments of an expanding supershell.

We discuss the induced star formation in dense walls of expanding shells. The fragmentation process is studied using the linear perturbation theory. The influence of the energy input, the ISM distribution and the ISM speed of sound is examined analytically and by numerical simulations. We formulate the universal condition for the gravitational fragmentation of expanding shells: if the total surface density of the disk is higher than a certain critical value, shells are unstable. The value of the critical density depends on the energy of the shell and the sound speed in the ISM.

We discuss differences in shapes, expansion velocities and fragmentation times of structures created by an energy deposition from a single Gamma Ray Burst (GRB) or an OB association to the ISM. After the initial inflation, supershells produced by GRBs are almost static or slowly expanding, contrary to more rapidly expanding supershells created by OB associations. We discuss the position of the energy source relative to the symmetry plane of the galaxy: observed arc-like structures can be the most dense parts of structures formed by an expansion from a source above or below the galactic plane. Arcs may also form, if the expansion takes place inside a giant HI cloud. We try to reproduce the size, the age and the average distance between OB associations in the Sextant region at the edge of LMC 4.

Conditions for the fragmentation of expanding shells due to gravitational instability are discussed. The self-similar analytical solution is compared with the results of 3-dimensional computer simulations for the expansion into homogeneous media. For realistic galactic disks we show that the formation of fragments is influenced by the amount of energy supply from the final number of young stars in an OB association, the value of the sound speed, the stratification and density of the ambient medium, the galactic differential rotation and the gravitational force perpendicular to the galactic plane. The typical size of gravitationally unstable shells is 1 kpc for an ambient gas density n=1 cm^-3. In thick disk galaxies the fragmentation occurs in nearly the whole shell while in thin disks it is restricted to the galactic equator. Unstable fragments may become molecular and trigger the formation of molecular clouds, and finally new star formation. We conclude that in dwarf galaxies the star formation may propagate in all directions turning the system into a star- burst. Contrary to that, the star formation in spiral galaxies propagates only in some directions in a thin strip near the symmetry plane, basically at the tips of the expanding shell.