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.

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.

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.

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

Gould's Belt is a flat local system composed of young OB stars, molecular clouds and neutral hydrogen within 500 pc from the Sun. It is inclined about 20 degrees to the galactic plane and its velocity field significantly deviates from rotation around the distant center of the Milky Way. We discuss possible models of its origin: free expansion from a point or from a ring, expansion of a shell, or a collision of a high velocity cloud with the plane of the Milky Way. Currently, no convincing model exists. Similar structures are identified in HI and CO distribution in our and other nearby galaxies.

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.

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.

Here we consider the strong evolution experienced by the matter reinserted by massive stars, both in giant star forming regions driven by a constant star formation rate, and in massive and coeval superstar clusters. In both cases we take into consideration the changes induced by stellar evolution on the number of massive stars, the number of ionizing photons and the integrated mechanical luminosity of the star forming regions. The latter is at all times compared with the critical luminosity that defines, for a given size, the lower mechanical luminosity limit above which the matter reinserted via strong winds and supernova explosions suffers frequent and recurrent thermal instabilities that reduce its temperature and pressure and inhibit its exit as part of a global wind. Instead, the unstable reinserted matter is compressed by the pervasive hot gas, and photoionization maintains its temperature at T $\sim$ 10$^4$ K. As the evolution proceeds, more unstable matter accumulates and the unstable clumps grow in size. Here we evaluate the possible self-shielding of thermally unstable clumps against the UV radiation field. Self shielding allows for a further compression of the reinserted matter which rapidly develops a high density neutral core able to absorb in its outer skin the incoming UV radiation. Under such conditions the cold (T $\sim$ 10 K) neutral cores soon surpass the Jeans limit and become gravitationally unstable, causing a new stellar generation with the matter reinserted by former massive stars. We present the results of several calculations of this positive star formation feedback scenario promoted by strong radiative cooling and mass loading.

The local velocity patterns of star forming regions, young OB stars, nearby OB associations, atomic and molecular gas are confronted with models of an expanding region. We test free expansion from a point or from a ring, expanding 2D shell, and expanding 3D belt with abrupt or gradual energy injection snow-plowing the ambient medium with or without the drag forces including fragmentation and porosity of the medium. There is no agreement on the expansion time, which varies from 30 - 100 Myr. The inclination of the Gould belt is not explained by the above models of expansion. An oblique impact of a high velocity cloud may explain it, but the observed velocity pattern is difficult to reproduce. The Gould's belt may be one of the many structures resulting from shell-shell collisions in the galactic plane. The origin of the Gould's belt may be connected to instabilities in the curling gas flows downstream from the Galaxy spiral arms, forming ISM clouds and star formation complexes.

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.

The evolution of galaxies results from a combination of internal and external processes. The star formation is an internal process transforming cold and dense cores of molecular clouds to stars. It may be triggered internally by expanding shells, or externally, e. g. by galaxy collisions. The gas accretion and galaxy merging events are external contributors to galaxy evolution. They compete with another internal process of galaxy evolution, which is the secular evolution redistributing the mass and angular momentum inside of galaxies as a consequence of bar and spiral arms formation. As a nearby example of gas accretion we mention the interacting system Milky Way - Large and Small Magellanic Clouds. Stripping of ISM in galaxy clusters is reviewed as another example of interaction of galaxies with their environment.