The faintest Local Group galaxies found lurking in and around the Milky Way halo provide a unique test bed for theories of structure formation and evolution on small scales. Deep Subaru and Hubble Space Telescope photometry demonstrates that their stellar populations are old, and that the star formation activity did not last longer than 2 Gyr in these systems. A few mechanisms that may lead to such a rapid quenching have been investigated by means of hydrodynamic simulations, without providing any final assessment so far. This is the first in a series of papers aimed at analysing the roles of stellar feedback, ram pressure stripping, host-satellite tidal interactions and reionization in cleaning the lowest-mass Milky Way companions of their cold gas, by using high-resolution, three-dimensional hydrodynamic simulations. We simulate an isolated ultrafaint dwarf galaxy loosely modeled after Bootes I, and examine whether or not stellar feedback alone could drive a substantial fraction of the ambient gas out from the shallow potential well. In contrast to simple analytical estimates, but in agreement with previous hydrodynamical studies, we find that most of the cold gas reservoir is retained. Conversely, a significant fraction of the metal-enriched stellar ejecta crosses the boundaries of the computational box with velocities exceeding the local escape velocity and is, thus, likely lost from the system. Although the total energy output from multiple supernova explosions exceeds the binding energy of the gas, no galactic-scale outflow develops in our simulations and as such, most of the ambient medium remains trapped within the weak potential well of the model galaxy. It seems thus unavoidable that, in order to explain the dearth of gas in ultrafaint dwarf galaxies, we will have to resort to environmental effects. This will be the subject of a forthcoming paper.

By means of grid-based, 3D hydrodynamical simulations we study the formation of second generation (SG) stars in a young globular cluster (GC) of mass 10^7 Msun, the possible progenitor of an old GC with a present mass ~(1-5) * 10^6 Msun. The cluster accretes external gas as its first generation (FG) asymptotic giant branch (AGB) stars release their ejecta and SG stars form. We consider two models characterised by different densities of the external gas. In both cases, we find that a very compact SG subsystem with central density >10^5M sun/pc^3 forms in the innermost regions of the cluster. The low-density model forms a population of extreme SG stars with high helium enhancement, followed by the formation of another SG group out of a mix of pristine gas and AGB ejecta and characterised by a modest helium enhancement. On the other hand, the high-density model forms in prevalence SG stars with modest helium enhancement. Our simulations illustrate the dynamical processes governing the formation of SG populations in GCs and shed light on the structural properties emerging at the end of this phase. The newly born SG groups have different concentrations, with more extreme SG stars more centrally concentrated than those with less extreme chemical abundances. The very high density of the SG subsystems implies that SG massive stars, if formed, might suffer frequent close encounters, collisions and gas stripping, thus possibly contributing further gas to the SG formation.

We present the results of a survey of N-body simulations aimed at studying the effects of the long-term dynamical evolution on the stellar mass function (MF) of multiple stellar populations in globular clusters. Our simulations show that if first-(1G) and second-generation (2G) stars have the same initial MF (IMF), the global MFs of the two populations are affected similarly by dynamical evolution and no significant differences between the 1G and the 2G MFs arise during the cluster's evolution. If the two populations have different IMFs, dynamical effects do not completely erase memory of the initial differences. Should observations find differences between the global 1G and 2G MF, these would reveal the fingerprints of differences in their IMFs. Irrespective of whether the 1G and 2G populations have the same global IMF or not, dynamical effects can produce differences between the local (measured at various distances from the cluster centre) 1G and 2G MFs; these differences are a manifestation of the process of mass segregation in populations with different initial structural properties. In dynamically old and spatially mixed clusters, however, differences between the local 1G and 2G MFs can reveal differences between the 1G and 2G global MFs. In general, for clusters with any dynamical age, large differences between the local 1G and 2G MFs are more likely to be associated with differences in the global MF. Our study also reveals a dependence of the spatial mixing rate on the stellar mass, another dynamical consequence of the multiscale nature of multiple-population clusters.

We addressed the so far unexplored issue of outflows induced by exponentially growing power sources, focusing on early supermassive black holes (BHs). We assumed that these objects grow to $10^9\;M_{\odot}$ by z=6 by Eddington-limited accretion and convert 5% of their bolometric output into a wind. We first considered the case of energy-driven and momentum-driven outflows expanding in a region where the gas and total mass densities are uniform and equal to the average values in the Universe at $z>6$. We derived analytic solutions for the evolution of the outflow, finding that, for an exponentially growing power with e-folding time $t_{Sal}$, the late time expansion of the outflow radius is also exponential, with e-folding time of $5t_{Sal}$ and $4t_{Sal}$ in the energy-driven and momentum-driven limit, respectively. We then considered energy-driven outflows produced by QSOs at the center of early dark matter halos of different masses and powered by BHs growing from different seeds. We followed the evolution of the source power and of the gas and dark matter density profiles in the halos from the beginning of the accretion until $z=6$. The final bubble radius and velocity do not depend on the seed BH mass but are instead smaller for larger halo masses. At z=6, bubble radii in the range 50-180 kpc and velocities in the range 400-1000 km s$^{-1}$ are expected for QSOs hosted by halos in the mass range $3\times10^{11}-10^{13}\;M_{\odot}$. By the time the QSO is observed, we found that the total thermal energy injected within the bubble in the case of an energy-driven outflow is $E_{th}\sim5 \times 10^{60}$ erg. This is in excellent agreement with the value of $E_{th}=(6.2\pm 1.7)\times 10^{60}$ erg measured through the detection of the thermal Sunyaev-Zeldovich effect around a large population of luminous QSOs at lower redshift. [abridged]

We report on five compact, extremely young (<10Myr) and blue (\beta_UV<-2.5, F_\lambda =\lambda^\beta) objects observed with VLT/MUSE at redshift 3.1169, 3.235, in addition to three objects at z=6.145. These sources are magnified by the Hubble Frontier Field galaxy clusters MACS~J0416 and AS1063. Their de-lensed half light radii (Re) are between 16 to 140pc, the stellar masses are ~1-20 X 10^6 Msun, the magnitudes are m_uv=28.8 - 31.4 (-17<Muv<-15) and specific star formation rates can be as large as ~800Gyr^-1. Multiple images of these systems are widely separated in the sky (up to 50'') and individually magnified by factors 3-40. Remarkably, the inferred physical properties of two objects are similar to those expected in some globular cluster formation scenarios, representing the best candidate proto-globular clusters (proto-GC) discovered so far. Rest-frame optical high dispersion spectroscopy of one of them at z=3.1169 yields a velocity dispersion \sigma_v~20km/s, implying a dynamical mass dominated by the stellar mass. Another object at z=6.145, with de-lensed Muv ~ -15.3 (m_uv ~ 31.4), shows a stellar mass and a star-formation rate surface density consistent with the values expected from popular GC formation scenarios. An additional star-forming region at z=6.145, with de-lensed m_uv ~ 32, a stellar mass of 0.5 X 10^6 Msun and a star formation rate of 0.06 Msun/yr is also identified. These objects currently represent the faintest spectroscopically confirmed star-forming systems at z>3, elusive even in the deepest blank fields. We discuss how proto-GCs might contribute to the ionization budget of the universe and augment Lya visibility during reionization. This work underlines the crucial role of JWST in characterizing the rest-frame optical and near-infrared properties of such low-luminosity high-z objects.

We study the interaction of the early spherical GC wind powered by Type II supernovae (SNe II) with the surrounding ambient medium consisting of the gaseous disk of a star forming galaxy at redshift z ~> 2. The bubble formed by the wind eventually breaks out of the disk, and most of the wind moves directly out of the galaxy and is definitively lost. The fraction of the wind moving nearly parallel to the galactic plane carves a hole in the disk which will contract after the end of the SN activity. During the interval of time between the end of the SN explosions and the "closure" of the hole, very O-poor stars (the Extreme population) can form out of the super-AGB (asymptotic giant branch) ejecta collected in the GC center. Once the hole contracts, the AGB ejecta mix with the pristine gas, allowing the formation of stars with an oxygen abundance intermediate between that of the very O-poor stars and that of the pristine gas. We show that this mechanism may explain why Extreme populations are present only in massive clusters, and can also produce a correlation between the spread in helium and the cluster mass. Finally, we also explore the possibility that our proposed mechanism can be extended to the case of multiple populations showing bimodality in the iron content, with the presence of two populations characterized by a small difference in [Fe/H]. Such a result can be obtained taking into account the contribution of delayed SN II.

We present the results of a survey of N-body simulations aimed at exploring the evolution of compact binaries in multiple-population globular clusters.We show that as a consequence of the initial differences in the structural properties of the first-generation (FG) and the second-generation (SG) populations and the effects of dynamical processes on binary stars, the SG binary fraction decreases more rapidly than that of the FG population. The difference between the FG and SG binary fraction is qualitatively similar to but quantitatively smaller than that found for wider binaries in our previous investigations.The evolution of the radial variation of the binary fraction is driven by the interplay between binary segregation, ionization and ejection. Ionization and ejection counteract in part the effects of mass segregation but for compact binaries the effects of segregation dominate and the inner binary fraction increases during the cluster evolution. We explore the variation of the difference between the FG and the SG binary fraction with the distance from the cluster centre and its dependence on the binary binding energy and cluster structural parameters. The difference between the binary fraction in the FG and the SG populations found in our simulations is consistent with the results of observational studies finding a smaller binary fraction in the SG population.

We explain the multiple populations recently found in the 'prototype' Globular Cluster (GC) NGC 2808 in the framework of the asymptotic giant branch (AGB) scenario. The chemistry of the five -or more- populations is approximately consistent with a sequence of star formation events, starting after the supernovae type II epoch, lasting approximately until the time when the third dredge up affects the AGB evolution (age ~90-120Myr), and ending when the type Ia supernovae begin exploding in the cluster, eventually clearing it from the gas. The formation of the different populations requires episodes of star formation in AGB gas diluted with different amounts of pristine gas. In the nitrogen-rich, helium-normal population identified in NGC 2808 by the UV Legacy Survey of GCs, the nitrogen increase is due to the third dredge up in the smallest mass AGB ejecta involved in the star formation of this population. The possibly-iron-rich small population in NGC 2808 may be a result of contamination by a single type Ia supernova. The NGC 2808 case is used to build a general framework to understand the variety of 'second generation' stars observed in GCs. Cluster-to-cluster variations are ascribed to differences in the effects of the many processes and gas sources which may be involved in the formation of the second generation. We discuss an evolutionary scheme, based on pollution by delayed type II supernovae, which accounts for the properties of s-Fe-anomalous clusters.

Globular clusters are considerably more complex structures than previously thought, harbouring at least two stellar generations which present clearly distinct chemical abundances. Scenarios explaining the abundance patterns in globular clusters mostly assume that originally the clusters had to be much more massive than today, and that the second generation of stars originates from the gas shed by stars of the first generation (FG). The lack of metallicity spread in most globular clusters further requires that the supernova-enriched gas ejected by the FG is completely lost within ~30 Myr, a hypothesis never tested by means of three-dimensional hydrodynamic simulations. In this paper, we use 3D hydrodynamic simulations including stellar feedback from winds and supernovae, radiative cooling and self-gravity to study whether a realistic distribution of OB associations in a massive proto-GC of initial mass M_tot ~ 10^7 M_sun is sufficient to expel its entire gas content. Our numerical experiment shows that the coherence of different associations plays a fundamental role: as the bubbles interact, distort and merge, they carve narrow tunnels which reach deeper and deeper towards the innermost cluster regions, and through which the gas is able to escape. Our results indicate that after 3 Myr, the feedback from stellar winds is responsible for the removal of ~40% of the pristine gas, and that after 14 Myr, ~ 99% of the initial gas mass has been removed.

Horizontal Branch stars belong to an advanced stage in the evolution of the oldest stellar galactic population, occurring either as field halo stars or grouped in globular clusters. The discovery of multiple populations in these clusters, that were previously believed to have single populations gave rise to the currently accepted theory that the hottest horizontal branch members (the blue hook stars, which had late helium-core flash ignition, followed by deep mixing) are the progeny of a helium-rich "second generation" of stars. It is not known why such a supposedly rare event (a late flash followed by mixing) is so common that the blue hook of {\omega} Cen contains \sim 30% of horizontal branch stars 10 , or why the blue hook luminosity range in this massive cluster cannot be reproduced by models. Here we report that the presence of helium core masses up to \sim 0.04 solar masses larger than the core mass resulting from evolution is required to solve the luminosity range problem. We model this by taking into account the dispersion in rotation rates achieved by the progenitors, whose premain sequence accretion disc suffered an early disruption in the dense environment of the cluster's central regions where second-generation stars form. Rotation may also account for frequent late-flash-mixing events in massive globular clusters.

The discovery of multiple stellar populations in globular clusters has implications for all the aspects of the study of these stellar systems. In this paper, by means of N-body simulations, we study the evolution of binary stars in multiple-population clusters and explore the implications of the initial differences in the spatial distribution of different stellar populations for the evolution and survival of their binary stars. Our simulations show that initial differences between the spatial distribution of first-generation (FG) and second-generation (SG) stars can leave a fingerprint in the current properties of the binary population. SG binaries are disrupted more efficiently than those of the FG population resulting in a global SG binary fraction smaller than that of the FG. As for surviving binaries, dynamical evolution produces a difference between the SG and the FG binary binding energy distribution with the SG population characterized by a larger fraction of high binding energy (more bound) binaries. We have also studied the dependence of the binary properties on the distance from the cluster centre. Although the global binary fraction decreases more rapidly for the SG population, the local binary fraction measured in the cluster inner regions may still be dominated by SG binaries. The extent of the differences between the surviving FG and SG binary binding energy distribution also varies radially within the cluster and is larger in the cluster inner regions.

We investigated the impact of supernova feedback in gas-rich dwarf galaxies experiencing a low-to-moderate star formation rate, typical of relatively quiescent phases between starbursts. We calculated the long term evolution of the ISM and the metal-rich SN ejecta using 3D hydrodynamic simulations, in which the feedback energy is deposited by SNeII exploding in distinct OB associations. We found that a circulation flow similar to galactic fountains is generally established, with some ISM lifted at heights of one to few kpc above the galactic plane. This gas forms an extra-planar layer, which falls back to the plane in about $10^8$ yr, once the star formation stops. Very little or no ISM is expelled outside the galaxy system for the considered SFRs, even though in the most powerful model the SN energy is comparable to the gas binding energy. The metal-rich SN ejecta is instead more vulnerable to the feedback and we found that a significant fraction (25-80\%) is vented in the intergalactic medium, even for low SN rate ($7\times 10^{-5}$ - $7\times 10^{-4}$ yr$^{-1}$). About half of the metals retained by the galaxy are located far ($z >$ 500 pc) from the galactic plane. Moreover, our models indicate that the circulation of the metal-rich gas out from and back to the galactic disk is not able to erase the chemical gradients imprinted by the (centrally concentrated) SN explosions.

Context: The study of chemical abundance patterns in globular clusters is of key importance to constrain the different candidates for intra-cluster pollution of light elements. Aims: We aim at deriving accurate abundances for a large range of elements in the globular cluster 47 Tucanae (NGC 104) to add new constraints to the pollution scenarios for this particular cluster, expanding the range of previously derived element abundances. Methods: Using tailored 1D LTE atmospheric models together with a combination of equivalent width measurements, LTE, and NLTE synthesis we derive stellar parameters and element abundances from high-resolution, high signal-to-noise spectra of 13 red giant stars near the tip of the RGB. Results: We derive abundances of a total 27 elements (O, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Y, Zr, Mo, Ru, Ba, La, Ce, Pr, Nd, Eu, Dy). Departures from LTE were taken into account for Na, Al and Ba. We find a mean [Fe/H] = $-0.78\pm0.07$ and $[\alpha/{\rm Fe}]=0.34\pm0.03$ in good agreement with previous studies. The remaining elements show good agreement with the literature, but the inclusion of NLTE for Al has a significant impact on the behaviour of this key element. Conclusions: We confirm the presence of an Na-O anti-correlation in 47 Tucanae found by several other works. Our NLTE analysis of Al shifts the [Al/Fe] to lower values, indicating that this may be overestimated in earlier works. No evidence for an intrinsic variation is found in any of the remaining elements.

We investigate the viability of a model in which the chemical anomalies among Globular Cluster stars are due to accretion of gas onto the protostellar discs of low mass stars. This model has been suggested as a way to reduce the large initial cluster masses required by other models for the formation of multiple stellar generations. We numerically follow the evolution of the accreting stars, and we show that the structure of the seed star does not remain fully convective for the whole duration of the accretion phase. Stellar populations showing discrete abundances of helium in the core, that seem to be present in some clusters, might be formed with this mechanism only if accretion occurs before the core of the stars become radiative (within 2-3Myr) or if a thermohaline instability is triggered, to achieve full mixing after the accretion phase ends. We also show that the lithium abundances in accreted structures may vary by orders of magnitude in equal masses obtained by accreting different masses. In addition, the same thermohaline mixing which could provide a homogeneous helium distribution down to the stellar center, would destroy any lithium surviving in the envelope, so that both helium homogeneity and lithium survival require that the accretion phase be limited to the first couple of million years of the cluster evolution. Such a short accretion phase strongly reduces the amount of processed matter available, and reintroduces the requirement of an extremely large initial mass for the protocluster.

We use the combination of photometric and spectroscopic data of 47 Tuc stars to reconstruct the possible formation of a second generation of stars in the central regions of the cluster, from matter ejected from massive Asymptotic Giant Branch stars, diluted with pristine gas. The yields from massive AGB stars with the appropriate metallicity (Z=0.004, i.e. [Fe/H]=-0.75) are compatible with the observations, in terms of extension and slope of the patterns observed, involving oxygen, nitrogen, sodium and aluminium. Based on the constraints on the maximum helium of 47 Tuc stars provided by photometric investigations, and on the helium content of the ejecta, we estimate that the gas out of which second generation stars formed was composed of about one-third of gas from intermediate mass stars, with M>= 5Mo and about two-thirds of pristine gas. We tentatively identify the few stars whose Na, Al and O abundances resemble the undiluted AGB yields with the small fraction of 47 Tuc stars populating the faint subgiant branch. From the relative fraction of first and second generation stars currently observed, we estimate that the initial FG population in 47 Tuc was about 7.5 times more massive than the cluster current total mass.

We analyze the evolution of colour gradients predicted by the hydrodynamical models of early type galaxies (ETGs) in Pipino et al. (2008), which reproduce fairly well the chemical abundance pattern and the metallicity gradients of local ETGs. We convert the star formation (SF) and metal content into colours by means of stellar population synthetic model and investigate the role of different physical ingredients, as the initial gas distribution and content, and eps_SF, i.e. the normalization of SF rate. From the comparison with high redshift data, a full agreement with optical rest-frame observations at z < 1 is found, for models with low eps_SF, whereas some discrepancies emerge at 1 < z < 2, despite our models reproduce quite well the data scatter at these redshifts. To reconcile the prediction of these high eps_SF systems with the shallower colour gradients observed at lower z we suggest intervention of 1-2 dry mergers. We suggest that future studies should explore the impact of wet galaxy mergings, interactions with environment, dust content and a variation of the Initial Mass Function from the galactic centers to the peripheries.

We examine the photometric data for Fornax clusters, focussing our attention on their horizontal branch color distribution and, when available, on the RR Lyr variables fraction and period distribution. Based on our understanding of the HB morphology in terms of varying helium content in the context of multiple stellar generations, we show that clusters F2, F3 and F5 must contain substantial fractions of second generation stars (~54-65%). On the basis of a simple chemical evolution model we show that the helium distribution in these clusters can be reproduced by models with cluster initial masses ranging from values equal to ~4 to ~10 times larger than the current masses. Models with a very short second generation star formation episode can also reproduce the observed helium distribution but require larger initial masses up to about twenty times the current mass. While the lower limit of this range of possible initial GC masses is consistent with those suggested by the observations of the low metallicity field stars, we also discuss the possibility that the metallicity scale of field stars (based on CaII triplet spectroscopy) and the metallicities derived for the clusters in Fornax may not be consistent with each other. The reproduction of the HB morphology in F2,F3,F5 requires two interesting hypotheses: 1) the first generation HB stars lie all at "red" colours. According to this interpretation, the low metallicity stars in the field of Fornax, populating the HB at colours bluer than the blue side ((V-I)o<=0.3 or (B-V)o<=0.2) of the RR Lyrs, should be second generation stars born in the clusters;a preliminary analysis of available colour surveys of Fornax field provides a fraction ~20% of blue HB stars, in the low metallicity range; 2) the mass loss from individual second generation red giants is a few percent of a solar mass larger than the mass loss from first generation stars.

In this paper we study the long-term dynamical evolution of multiple-population clusters, focusing on the evolution of the spatial distributions of the first- (FG) and second-generation (SG) stars.In previous studies we have suggested that SG stars formed from the ejecta of FG AGB stars are expected initially to be concentrated in the cluster inner regions. Here, by means of N-body simulations, we explore the time scales and the dynamics of the spatial mixing of the FG and the SG populations and their dependence on the SG initial concentration.Our simulations show that, as the evolution proceeds, the radial profile of the SG/FG number ratio, NSG/NFG, is characterized by three regions: 1) a flat inner part; 2) a declining part in which FG stars are increasingly dominant; and 3) an outer region where the NSG/NFG profile flattens again (the NSG/NFG profile may rise slightly again in the outermost cluster regions). The radial variation of NSG/NFG implies that the fraction of SG stars determined by observations covering a limited range of radial distances is not, in general, equal to the SG global fraction, (NSG/NFG)glob. The distance at which NSG/NFG equals (NSG/NFG)glob is approximately between 1 and 2 cluster half-mass radii. The results of our simulations suggest that in many Galactic globular clusters the SG should still be more spatially concentrated than the FG.[abridged]

We follow the scenario of formation of second generation stars in globular clusters by matter processed by hot bottom burning (HBB) in massive asymptotic giant branch (AGB) stars and Super AGB stars (SAGB). In the cluster NGC 2419 we assume the presence of an extreme population directly formed from the AGB and SAGB ejecta, so we can directly compare the yields for a metallicity Z=0.0003 with the chemical inventory of the cluster NGC 2419. At such a low metallicity, the HBB temperatures (well above 108K) allow a very advanced nucleosynthesis. Masses of about 6Mo deplete Mg and synthesize Si, going beyond Al, so this latter element results only moderately enhanced; sodium can not be enhanced. The models are consistent with the observations, although the predicted Mg depletion is not as strong as in the observed stars. We predict that the oxygen abundance must be depleted by a huge factor in the Mg poor stars. The HBB temperatures are close to the region where other p capture reactions on heavier nuclei become possible. We show that high potassium abundance found in Mg poor stars can be achieved during HBB, by p captures on the argon nuclei, if the relevant cross section(s) are larger than listed in the literature or if the HBB temperature is higher. Finally, we speculate that some calcium production is occurring owing to proton capture on potassium. We emphasize the importance of a strong effort to measure a larger sample of abundances in this cluster.

Globular cluster stars show chemical abundance patterns typical of hot-CNO processing. Lithium is easily destroyed by proton capture in stellar environments, so its abundance may be crucial to discriminate among different models proposed to account for multiple populations. In order to reproduce the observed O-Na anticorrelation and other patterns typical of multiple populations, the formation of second generation stars must occur from the nuclearly processed stellar ejecta, responsible of the chemical anomalies, diluted with pristine gas having the composition of first generation stars. The lithium abundance in the unprocessed gas -which is very likely to be equal to the lithium abundance emerging from the Big Bang- affects the lithium chemical patterns among the cluster stars. This paper focuses on a scenario in which processed gas is provided by asymptotic giant branch (AGB) stars. We examine the predictions of this scenario for the lithium abundances of multiple populations. We study the role of the non-negligible lithium abundance in the ejecta of massive AGB (A(Li)~2), and, at the same time, we explore how our models can constrain the extremely large ---and very model dependent--- lithium yields predicted by recent super--AGB models. We show that the super--AGB yields may be tested by examining the lithium abundances in a large set of blue main sequence stars in wCen and/or NGC2808. In addition, we examine the different model results obtained by assuming for the pristine gas either the Big Bang abundance predicted by the standard models (A(Li)=2.6-2.7), or the abundance detected at the surface of population II stars (A(Li)=2.2-2.3). Once a chemical model is well constrained, the O--Li distribution could perhaps be used to shed light on the primordial lithium abundance.

In order to account for the chemical composition of a stellar second generation (SG), Globular Clusters (GCs) evolution models based on the asymptotic giant branch (AGB) scenario so far included only the yields available for the massive AGB stars, while the possible role of super-AGB ejecta was either extrapolated or not considered. In this work, we explore the role of super-AGB ejecta using yields recently calculated by Ventura and D'Antona. Models of clusters showing extended Na-O anticorrelations, like NGC 2808, indicate that a SG formation history similar to that outlined in our previous work is required: formation of an Extreme population with very large helium content from the pure ejecta of super-AGB stars, followed by formation of an Intermediate population by dilution of stellar ejecta with pristine gas. The very O-poor Na-rich Extreme stars can be accounted for once deep-mixing is assumed in SG giants forming in a gas with helium abundance Y> 0.34, which significantly reduces the atmospheric oxygen content, while preserving the sodium abundance. On the other hand, for clusters showing a mild O-Na anticorrelation, like M 4, the use of the new yields broadens the range of SG formation routes leading to abundance patterns consistent with observations. It is shown that models in which SG stars form only from super-AGB ejecta promptly diluted with pristine gas can reproduce the observations. We discuss the variety of small helium variations occurring in this model and its relevance for the horizontal branch morphology. In some of these models the duration of the SG formation episode can be as short as \sim10 Myr; the formation time of the SG is thus compatible with the survival of a cooling flow in the GC, previous to the explosion of the SG core collapse supernovae. We also explore models with formation of multiple populations in individual bursts, each lasting no longer than \sim10 Myr.

The discovery of multiple stellar populations in globular clusters raises fundamental questions concerning the formation and dynamical history of these systems. In a previous study aimed at exploring the formation of second-generation (SG) stars from the ejecta of first-generation (FG) AGB stars, and the subsequent dynamical evolution of the cluster, we showed that SG stars are expected to form in a dense subsystem concentrated in the inner regions of the FG cluster. In this paper we explore the implications of the structural properties of multiple-population clusters, and in particular the presence of the inner SG subsystem, for the disruption of binary stars. We quantify the enhancement of the binary disruption rate due to the presence of the central SG subsystem for a number of different initial conditions. Our calculations show that SG binaries, which are assumed to be more concentrated in the cluster inner regions, are disrupted at a substantially larger rate than FG binaries. Assuming a similar initial fraction of FG and SG binaries, our dynamical study indicates that the SG population is now expected to contain a significantly smaller binary fraction than the FG population.

Recent exam of large samples of omega Cen giants shows that it shares with mono-metallic globular clusters the presence of the sodium versus oxygen anticorrelation, within each subset of stars with iron content in the range -1.9<~[Fe/H]<~-1.3. These findings suggest that, while the second generation formation history in omega Cen is more complex than that of mono-metallic clusters, it shares some key steps with those simpler cluster. In addition, the giants in the range -1.3<[Fe/H]<~-0.7 show a direct O--Na correlation, at moderately low O, but Na up to 20 times solar. These peculiar Na abundances are not shared by stars in other environments often assumed to undergo a similar chemical evolution, such as in the field of the Sagittarius dwarf galaxy. These O and Na abundances match well the yields of the massive asymptotic giant branch stars (AGB) in the same range of metallicity, suggesting that the stars at [Fe/H]>-1.3 in omega Cen are likely to have formed directly from the pure ejecta of massive AGBs of the same metallicities. This is possible if the massive AGBs of [Fe/H]>-1.3 in the progenitor system evolve when all the pristine gas surrounding the cluster has been exhausted by the previous star formation events, or the proto--cluster interaction with the Galaxy caused the loss of a significant fraction of its mass, or of its dark matter halo, and the supernova ejecta have been able to clear the gas out of the system. The absence of dilution in the metal richer populations lends further support to a scenario of the formation of second generation stars in cooling flows from massive AGB progenitors. We suggest that the entire formation of omega Cen took place in a few 10^8yr, and discuss the problem of a prompt formation of s--process elements.

The star-to-star differences in the abundance of light elements observed in the globular clusters (GCs) can be explained assuming that a second generation (SG) of stars form in the gas ejected by the asymptotic giant branch (AGB) stars belonging to a first stellar generation. However, while Na and O appear to be anticorrelated in the cluster stars, from the stellar models they turn out to be correlated into the AGB ejecta. In order to reconcile the stellar theory with the observational findings, all the GC models invoke an early dilution of AGB ejecta with pristine gas occurring during the SG formation. Despite a vast consensus about the occurrence of such dilution, the physical process behind it is still unknown. In the present paper we set some general constraints on the pristine gas dynamics and on the possible amount of pristine gas involved in the SG formation, making use of a one zone chemical model. We find that such a dilution is a necessary ingredient in the SG star formation to explain the observed abundance patterns. We confirm the conclusion of our previous works showing that clusters must have been initially much more massive. We also show that models assuming that clusters had an initial mass similar to their current one, and adopting a large fraction of pristine gas to form SG stars, fail to reproduce the observed Na-O anticorrelation and are not viable. We finally show that the dilution event should be restricted in time, rather than extended for the all duration of the SG formation.

AGN heating, through massive subrelativistic outflows, might be the key to solve the long-lasting `cooling flow problem' in cosmological systems. In a previous paper, we showed that cold accretion feedback and, to a lesser degree, Bondi self-regulated models are in fact able to quench cooling rates for several Gyr, at the same time preserving the mainc ool core features, like observed density and temperature profiles. Is it true also for lighter systems, such as galaxy groups? The answer is globally yes, although with remarkable differences. Adopting a modified version of the AMR code FLASH 3.2, we found that successful 3D simulations with cold and Bondi models are almost convergent in the galaxy group environment, with mechanical efficiencies in the range 5.e-4 - 1.e-3 and 5.e-2 - 1.e-1, respectively. The evolutionary storyline of galaxy groups is dominated by a quasi-continuous gentle injection with sub-Eddington outflows (with mechanical power and velocity around 1.e44 erg/s and 1.e4 km/s). The cold and hybrid accretion models present, in addition, very short quiescence periods, followed by moderate outbursts (10 times the previous phase), which generate a series of 10-20 kpc size cavities with high density contrast, temperatures similar to the ambient medium and cold rims. After shock heating, a phase of turbulence promotes gas mixing and diffusion of metals, which peak along jet-axis (up to 40 kpc) during active phases. At this stage the tunnel, produced by the enduring outflow (hard to detect in the mock SBx maps), is easily fragmented, producing tiny buoyant bubbles, typically a few kpc in size. In contrast to galaxy clusters, the AGN self-regulated feedback has to be persistent, with a `delicate touch', rather than rare and explosive strokes. This evolutionary difference dictates in the end that galaxy groups are not scaled-down versions of clusters.

It is now widely accepted that heating processes play a fundamental role in galaxy clusters, struggling in an intricate but fascinating `dance' with its antagonist, radiative cooling. Last generation observations, especially X-ray, are giving us tiny hints about the notes of this endless ballet. Cavities, shocks, turbulence and wide absorption-lines indicate the central active nucleus is injecting huge amount of energy in the intracluster medium. However, which is the real dominant engine of self-regulated heating? One of the model we propose are massive subrelativistic outflows, probably generated by a wind disc or just the result of the entrainment on kpc scale by the fast radio jet. Using a modified version of AMR code FLASH 3.2, we explored several feedback mechanisms which self-regulate the mechanical power. Two are the best schemes that answer our primary question, id est quenching cooling flow and at the same time preserving a cool core appearance for a long term evolution (7 Gyr): one more explosive (with efficiencies 0.005 - 0.01), triggered by central cooled gas, and the other gentler, ignited by hot gas Bondi accretion (with efficiency 0.1). These three-dimensional simulations show that the total energy injected is not the key aspect, but the results strongly depend on how energy is given to the ICM. We follow the dynamics of best model (temperature, density, SB maps and profiles) and produce many observable predictions: buoyant bubbles, ripples, turbulence, iron abundance maps and hydrostatic equilibrium deviation. We present a deep discussion of merits and flaws of all our models, with a critical eye towards observational concordance.

Many observational studies have revealed the presence of multiple stellar generations in Galactic globular clusters. These studies suggest that second-generation stars make up a significant fraction of the current mass of globular clusters, with the second-generation mass fraction ranging from ~50 to 80 per cent in individual clusters. In this Letter we carry out hydrodynamical simulations to explore the dependence of the mass of second-generation stars on the initial mass and structural parameters and stellar initial mass function (IMF) of the parent cluster. We then use the results of these simulations to estimate the fraction, f_{SG,H}, of the mass of the Galactic stellar halo composed of second-generation stars that originated in globular clusters. We study the dependence of f_{SG,H} on the parameters of the initial mass function of the Galactic globular cluster system. For a broad range of initial conditions, we find that the fraction of mass of the Galactic stellar halo in second-generation stars is always small, f_{SG,H}<4-6 per cent for a Kroupa-1993 IMF and f_{SG,H}<7-9 per cent for a Kroupa-2001 IMF.

We investigate whether it is possible to explain the wide range of observed gradients in early type galaxies in the framework of monolithic models. To do so, we extend the set of hydrodynamical simulations by Pipino et al. (2008a) by including low-mass ellipticals and spiral (true) bulges. These models satisfy the mass-metallicity and the mass-[alpha/Fe] relations. The typical metallicity gradients predicted by our models have a slope of -0.3 dex per decade variation in radius, consistent with the mean values of several observational samples. However, we also find a few quite massive galaxies in which this slope is -0.5 dex per decade, in agreement with some recent data. In particular, we find a mild dependence from the mass tracers when we transform the stellar abundance gradients into radial variations of the Mg_2 line-strength index, but not in the Mg_b. We conclude that, rather than a mass- slope relation, is more appropriate to speak of an increase in the scatter of the gradient slope with the galactic mass. We can explain such a behaviour with different efficiencies of star formation in the framework of the revised monolithic formation scenario, hence the scatter in the observed gradients should not be used as an evidence of the need of mergers. Indeed, model galaxies that exhibit the steepest gradient slopes are preferentially those with the highest star formation efficiency at that given mass.

A large number of spectroscopic studies have provided evidence of the presence of multiple populations in globular clusters by revealing patterns in the stellar chemical abundances. This paper is aimed at studying the origin of these abundance patterns. We explore a model in which second generation (SG) stars form out of a mix of pristine gas and ejecta of the first generation of asymptotic giant branch stars. We first study the constraints imposed by the spectroscopic data of SG stars in globular clusters on the chemical properties of the asymptotic and super asymptotic giant branch ejecta. With a simple one-zone chemical model, we then explore the formation of the SG population abundance patterns focussing our attention on the Na-O, Al-Mg anticorrelations and on the helium distribution function. We carry out a survey of models and explore the dependence of the final SG chemical properties on the key parameters affecting the gas dynamics and the SG formation process. Finally, we use our chemical evolution framework to build specific models for NGC 2808 and M4, two Galactic globular clusters which show different patterns in the Na-O and Mg-Al anticorrelation and have different helium distributions. We find that the amount of pristine gas involved in the formation of SG stars is a key parameter to fit the observed O-Na and Mg-Al patterns. The helium distribution function for these models is in general good agreement with the observed one. Our models, by shedding light on the role of different parameters and their interplay in determining the final SG chemical properties, illustrate the basic ingredients, constraints and problems encountered in this self-enrichment scenario which must be addressed by more sophisticated chemical and hydrodynamic simulations.

We consider the horizontal branch (HB) of the Globular Cluster Terzan 5, recently shown to be split into two parts, the fainter one (delta M_K ~ 0.3mag) having a lower metallicity than the more luminous. Both features show that it contains at least two stellar populations. The separation in magnitude has been ascribed to an age difference of ~6 Gyr and interpreted as the result of an atypical evolutionary history for this cluster. We show that the observed HB morphology is also consistent with a model in which the bright HB is composed of second generation stars that are metal enriched and with a helium mass fraction larger (by delta Y ~ 0.07) than that of first generation stars populating the fainter part of the HB. Terzan 5 would therefore be anomalous, compared to most "normal" clusters hosting multiple populations, only because its second generation is strongly contaminated by supernova ejecta; the previously proposed prolonged period of star formation, however, is not required. The iron enrichment of the bright HB can be ascribed either to contamination from Type Ia supernova ejecta of the low-iron, helium rich, ejecta of the massive asympotic giant branch stars of the cluster, or to its mixing with gas, accreting on the cluster from the environment, that has been subject to fast metal enrichment due to its proximity with the galactic bulge. The model here proposed requires only a small age difference, of ~100Myr.

Spectroscopic and photometric observations show that many globular clusters host multiple stellar populations, challenging the common paradigm that globular clusters are "simple stellar populations" composed of stars of uniform age and chemical composition. The chemical abundances of second-generation (SG) stars constrain the sources of gas out of which these stars must have formed, indicating that the gas must contain matter processed through the high-temperature CNO cycle. First-generation massive Asymptotic Giant Branch (AGB) stars have been proposed as the source of this gas. In a previous study, by means of hydrodynamical and N-body simulations, we have shown that the AGB ejecta collect in a cooling flow in the cluster core, where the gas reaches high densities, ultimately forming a centrally concentrated subsystem of SG stars. In this Letter we show that the high gas density can also lead to significant accretion onto a pre-existing seed black hole. We show that gas accretion can increase the black hole mass by up to a factor of 100. The details of the gas dynamics are important in determining the actual black hole growth. Assuming a near-universal seed black hole mass and small cluster-to-cluster variations in the duration of the SG formation phase, the outcome of our scenario is one in which the present intermediate-mass black hole (IMBH) mass may have only a weak dependence on the current cluster properties. The scenario presented provides a natural mechanism for the formation of an IMBH at the cluster center during the SG star-formation phase.

We compute 3D gasdynamical models of jet outflows from the central AGN, that carry mass as well as energy to the hot gas in galaxy clusters and groups. These flows have many attractive attributes for solving the cooling flow problem: why the hot gas temperature and density profiles resemble cooling flows but show no spectral evidence of cooling to low temperatures. Subrelativistic jets, described by a few parameters, are assumed to be activated when gas flows toward or cools near a central SMBH. Using approximate models for a rich cluster (A1795), a poor cluster (2A 0336+096) and a group (NGC 5044), we show that mass-carrying jets with intermediate mechanical efficiencies ($\sim10^{-3}$) can reduce for many Gyr the global cooling rate to or below the low values implied by X-spectra, while maintaining $T$ and $\rho$ profiles similar to those observed, at least in clusters. Groups are much more sensitive to AGN heating and present extreme time variability in both profiles. Finally, the intermittency of the feedback generates multiple generations of X-ray cavities similar to those observed in Perseus cluster and elsewhere. Thus we also study the formation of buoyant bubbles and weak shocks in the ICM, along with the injection of metals by SNIa and stellar winds.

Using 3D hydrodynamical simulations, we studied in detail the fountain flow and its dependence with several factors, such as the Galactic rotation, the distance to the Galactic center, and the presence of a hot gaseous halo. We have considered the observed size-frequency distribution of young stellar clusters within the Galaxy in order to appropriately fuel the multiple fountains in our simulations. The present work confirms the localized nature of the fountain flows: the freshly ejected metals tend to fall back close to the same Galactocentric region where they are delivered. Therefore, the fountains do not change significantly the radial profile of the disk chemical abundance. The multiple fountains simulations also allowed to consistently calculate the feedback of the star formation on the halo gas. Finally, we have also considered the possibility of mass infall from the intergalactic medium and its interaction with the clouds that are formed by the fountains. Though our simulations are not suitable to reproduce the slow rotational pattern that is typically observed in the halos around the disk galaxies, they indicate that the presence of an external gas infall may help to slow down the rotation of the gas in the clouds and thus the amount of angular momentum that they transfer to the coronal gas, as previously suggested in the literature.

We review here the effects of supernovae (SNe) explosions on the environment of star-forming galaxies. Randomly distributed, clustered SNe explosions cause the formation of hot superbubbles that drive either galactic fountains or supersonic winds out of the galactic disk. In a galactic fountain, the ejected gas is re-captured by the gravitational potential and falls back onto the disk. From 3D non-equilibrium radiative cooling hydrodynamical simulations of these fountains, we find that they may reach altitudes smaller than 5 kpc in the halo and hence explain the formation of the so-called intermediate-velocity-clouds (IVCs). On the other hand, the high-velocity-clouds (HVCs) that are observed at higher altitudes (of up to 12 kpc) require another mechanism to explain their production. We argue that they could be formed either by the capture of gas from the intergalactic medium and/or by the action of magnetic fields that are carried out to the halo with the gas in the fountains. Due to angular momentum losses (of 10-15%) to the halo, we find that the fountain material falls back to smaller radii and is not largely spread over the galactic disk, as previously expected. This result is consistent with the metal distribution derived from recent chemical models of the galaxy. We also find that after about 150 Myr, the gas circulation between the halo and the disk in the fountains reaches a steady state regime (abridged).

Some globular clusters are observed to host a population of second generation (SG) stars which show chemical anomalies and must have formed from gas containing matter processed in the envelopes of first generation (FG) cluster stars. We study the SG formation process by means of 1D hydrodynamical simulations, assuming that the SG is formed by the gas ejected by AGB stars. This gas collects in a cooling flow into the cluster core, where it forms a SG star subsystem strongly concentrated in the cluster innermost regions with structural properties largely independent of the FG initial properties. We also present the results of a model in which pristine gas contributes to the SG formation. In this model a very helium-rich SG population and one with a moderate helium enrichment form; the resulting SG bimodal helium distribution resembles that observed for SG stars in NGC 2808. By means of N-body simulations, we study the two-population cluster dynamical evolution. In our simulations, a large fraction of FG stars are lost early in the cluster evolution due to the expansion and stripping of the cluster outer layers resulting from early mass loss associated with FG SN ejecta. The SG population is largely unscathed by this early mass loss, and this early evolution leads to values of the SG to FG number ratio consistent with observations. We also demonstrate possible evolutionary routes leading to the loss of most of the FG population, leaving an SG-dominated cluster. Until mixing of the two populations is complete, the radial profile of the SG to FG number ratio is characterized by a flat inner part and a declining portion in the outer cluster regions. (abridged)

We present three dimensional hydrodynamical simulations aimed at studying the dynamical and chemical evolution of the interstellar medium (ISM) in isolated dwarf spheroidal galaxies (dSphs). This evolution is driven by the explosion of Type II and Type Ia supernovae, whose different contribution on both the dynamics and chemical enrichment is taken into account. Radiative losses are effective in radiating away the huge amount of energy released by SNe explosions, and the dSph is able to retain most of the gas allowing a long period (>2-3 Gyr) of star formation, as usually observed in this kind of galaxies. We are able to reproduce the stellar metallicity distribution function (MDF) as well as the peculiar chemical properties of strongly O-depleted stars observed in several dSphs. The model also naturally predicts two different stellar populations, with an anti-correlation between [Fe/H] and velocity dispersion, similarly to what observed in the Sculptor and Fornax dSphs. These results derive from the inhomogeneous pollution of the SNe Ia, a distinctive characteristic of our model. We also applied the model to the peculiar globular cluster (GC) Omega Cen in the hypothesis that it is the remnant of a formerly larger stellar system, possibly a dSph.

The ejection of the gas out of the disk in late-type galaxies is related to star formation and is due mainly to Type II supernovae. In this paper we studied in detail the development of the Galactic fountains in order to understand their dynamical evolution and their influence in the redistribution of the freshly delivered metals over the disk. To this aim, we performed a number of 3D hydrodynamical radiative cooling simulations of the gas in the Milky Way where the whole Galaxy structure, the Galactic differential rotation and the supernovae explosions generated by a single OB association are considered. A typical fountain powered by 100 Type II supernovae may eject material up to $\sim 2$ kpc which than collapses back mostly in form of dense, cold clouds and filaments. The majority of the gas lifted up by the fountains falls back on the disk remaining within a radial distance $\Delta R=0.5$ kpc from the place where the fountain originated. This localized circulation of disk gas does not influence the radial chemical gradients on large scale, as required by the chemical models of the Milky Way which reproduce the metallicity distribution without invoking large fluxes of metals. Simulations of multiple fountains fuelled by Type II supernovae of different OB associations will be presented in a companion paper.

The scope of this paper is two-fold: i) to test and improve our previous models of an outside-in formation for the majority of ellipticals in the context of the SN-driven wind scenario, by means of a careful study of gas inflows/outflows; ii) to explain the observed slopes, either positive or negative, in the radial gradient of the mean stellar [alpha/Fe], and their apparent lack of any correlation with all the other observables. In order to pursue these goals we present a new class of hydrodynamical simulations for the formation of single elliptical galaxies in which we implement detailed prescriptions for the chemical evolution of H, He, O and Fe. We find that all the models which predict chemical properties (such as the central mass-weighted abundance ratios, the colours as well as the [<Fe/H>] gradient) within the observed ranges for a typical elliptical, also exhibit a variety of gradients in the [<alpha/Fe>] ratio, in agreement with the observations (namely positive, null or negative). All these models undergo an outside-in formation, in the sense that star formation stops earlier in the outermost than in the innermost regions, owing to the onset of a galactic wind. The typical [<Z/H>] gradients predicted by our models have a slope of -0.3 dex per decade variation in radius, consistent with the mean values of several observational samples. We can safely conclude that the history of star formation is fundamental for the creation of abundance gradients in ellipticals but that radial flows with different velocity in conjunction with the duration and efficiency of star formation in different galactic regions are responsible for the gradients in the [<alpha/Fe>] ratios.

Using 3-dimensional hydrodynamical simulations of isolated dwarf spheroidal galaxies (dSphs), we undertake an analysis of the chemical properties of their inner regions, identifying the respective roles played by Type Ia (SNe Ia) and Type II (SNe II) supernovae. The effect of inhomogeneous pollution from SNe Ia is shown to be prominent within two core radii, with the stars forming therein amounting to ~20% of the total. These stars are relatively iron-rich and alpha-element-depleted compared to the stars forming in the rest of the galaxy. At odds with the projected stellar velocity dispersion radial profile, the actual 3-dimensional one shows a depression in the central region, where the most metal-rich (ie. [Fe/H]-rich) stars are partly segregated. This naturally results in two different stellar populations, with an anti-correlation between [Fe/H] and velocity dispersion, in the same sense as that observed in the Sculptor and Fornax dSphs. Because the most iron-rich stars in our model are also the most alpha-depleted, a natural prediction and test of our model is that the same radial segregation effects should exist between [alpha/Fe] and velocity dispersion.

We present a new class of hydrodynamical models for the formation of bulges (either massive elliptical galaxies or classical bulges in spirals) in which we implement detailed prescriptions for the chemical evolution of H, He, O and Fe. Our results hint toward an outside-in formation in the context of the supernovae-driven wind scenario. The build-up of the chemical properties of the stellar populations inhabiting the galactic core is very fast. Therefore we predict a non significant evolution of both the mass-metallicity and the mass-[alpha/Fe] relations after the first 0.5 - 1 Gyr. In this framework we explain how the observed slopes, either positive or negative, in the radial gradient of the mean stellar [alpha/Fe], and their apparent lack of any correlation with all the other observables, can arise as a consequence of the interplay between star formation and metal-enhanced internal gas flows.

We present a hydrodynamical and chemical model for the globular cluster Omega Cen, under the assumption that it is the remnant of an ancient dwarf spheroidal galaxy (dSph), the bulk of which was disrupted and accreted by our Galaxy ~10 Gyr ago. We highlight the very different roles played by Type II and Type Ia supernovae (SNe) in the chemical enrichment of the inner regions of the putative parent dSph. While the SNe II pollute the interstellar medium rather uniformly, the SNe Ia ejecta may remain confined inside dense pockets of gas as long as succesive SNe II explosions spread them out. Stars forming in such pockets have lower alpha-to-iron ratios than the stars forming elsewhere. Owing to the inhomogeneous pollution by SNe Ia, the metal distribution of the stars in the central region differs substantially from that of the main population of the dwarf galaxy, and resembles that observed in Omega Cen. This inhomogeneous mixing is also responsible for a radial segregation of iron-rich stars with depleted [alpha/Fe] ratios, as observed in some dSphs. Assuming a star formation history of ~1.5 Gyr, our model succeeds in reproducing both the iron and calcium distributions observed in Omega Cen and the main features observed in the empirical alpha/Fe versus Fe/H plane. Finally, our model reproduces the overall spread of the color-magnitude diagram, but fails in reproducing the morphology of the SGB-a and the double morphology of the main sequence. However, the inhomogeneous pollution reduces (but does not eliminate) the need for a significantly enhanced helium abundance to explain the anomalous position of the blue main sequence. Further models taking into account the dynamical interaction of the parent dwarf galaxy with the Milky Way and the effect of AGB pollution will be required.

Elliptical galaxies probably host the most metal rich stellar populations in the Universe. The processes leading to both the formation and the evolution of such stars are discussed by means of a new gas dynamical model which implements detailed chemical evolution prescriptions. Moreover, the radial variations in the metallicity distribution of these stars are investigated by means of G-dwarf-like diagrams. By comparing model predictions with observations, we derive a picture of galaxy formation in which the higher is the mass of the galaxy, the shorter are the infall and the star formation timescales. The galaxies seem to have formed outside-in, namely the most external regions accrete gas, form stars and develop a galactic wind very quickly (a few Myr) compared to the central core, where the star formation can last up to 1 Gyr. We show for the first time a model able in reproducing the mass-metallicity and the color-magnitude relations as well as the radial metallicity gradient, and, at the same time, the observed either positive or negative slopes in the [alpha/Fe] abundace ratio gradient in stars.

We present 3D hydrodynamic simulations aimed at studying the dynamical and chemical evolution of the interstellar medium in dwarf spheroidal galaxies. This evolution is driven by the explosions of Type II and Type Ia supernovae, whose different contribution is explicitly taken into account in our models. We compare our results with avaiable properties of the Draco galaxy. Despite the huge amount of energy released by SNe explosions, in our model the galaxy is able to retain most of the gas allowing a long period ($> 3$ Gyr) of star formation, consistent with the star formation history derived by observations. The stellar [Fe/H] distribution found in our model matches very well the observed one. The chemical properties of the stars derive from the different temporal evolution between Type Ia and Type II supernova rate, and from the different mixing of the metals produced by the two types of supernovae. We reproduce successfully the observed [O/Fe]-[Fe/H] diagram.

We present 3D hydrodynamical simulations aimed to study the dynamical and chemical evolution of the interstellar medium in dwarf spheroidal galaxies. This evolution is driven by the explosions of Type II and Type Ia supernovae, whose different contribution is explicity taken into account in our models. We compare our results with detailed observations of the Draco galaxy. We assume star formation histories consisting of a number of instantaneous burst separated by quiescent periods. Because of the large effectiveness of the radiative losses and the extended dark matter halo, no galactic wind develops, despite the total energy released by the supernovae is much larger than the binding energy of the gas. This explains why the galaxy is able to form stars for a long period (> 3 Gyr), consistently with observations. In this picture, the end of the star formation and gas removal must result from external mechanisms, such as ram pressure and/or tidal interaction with the Galaxy. The metallicity distributions of the stars found in our models agree very well with the observed one. We find a mean value <[Fe/H]>=-1.65 with a spread of ~1.5 dex. The chemical properties of the stars derive by the different temporal evolution between Type Ia and Type II supernova rate, and by the different mixing of the metals produced by the two types of SNe. We reproduce successfully the observed [O/Fe]-[Fe/H] diagram. However, our interpretation of this diagram differs from that generally adopted by previous chemical models. In fact, we find that the chemical properties of the stars derive, besides the different temporal evolution of the SNe II and SNe Ia rates, from the spatial inhomogeneous chemical enrichment due to the different dynamical behaviour between the remnants of the two types of supernovae.

In this paper we present high-resolution hydrodynamical models of warm ionized clouds embedded in a superwind, and compare the OVI and soft X-ray properties to the existing observational data. These models include thermal conduction, which we show plays an important role in shaping both the dynamics and radiative properties of the resulting wind/cloud interaction. Heat conduction stabilizes the cloud by inhibiting the growth of K-H and R-T instabilities, and also generates a shock wave at the cloud's surface that compresses the cloud. This dynamical behaviour influences the observable properties. We find that while OVI emission and absorption always arises in cloud material at the periphery of the cloud, most of the soft X-ray arises in the region between the wind bow shock and the cloud surface, and probes either wind or cloud material depending on the strength of conduction and the relative abundances of the wind with respect to the cloud. In general only a small fraction (<1%) of the wind mechanical energy intersecting a cloud is radiated away at UV and X-ray wavelengths, with more wind energy going into accelerating the cloud. Models with heat conduction at Spitzer-levels are found to produce observational properties closer to those observed in superwinds than models with no thermal conduction, in particular in terms of the OVI-to-X-ray luminosity ratio, but cloud life times are uncomfortably short (<1Myr) compared to the dynamical ages of real winds. We experimented with reducing the thermal conductivity and found that even when we reduced conduction by a factor of 25 that the simulations retained the beneficial hydrodynamical stability and low O{\sc vi}-to-X-ray luminosity ratio found in the Spitzer-level conductive models, while also having reduced evaporation rates.

We present preliminary results of 3-D hydro simulations of the interstellar medium evolution in dwarf spheroidal galaxies undergoing star formation for the first time. The star formation is assumed to occur in a sequence of instantaneus bursts separated by quiescent periods. Different models are made changing the number and the intensity of the bursts in such a way that the final mass of the formed stars remains the same. We followed the enrichment of the ISM taking into account the contribution of both type Ia and II supernovae. The aim of our models is to find a star formation history compatible with the observed spread of stellar age and metallicity in such galaxies and to reproduce the observed mass-metallicity relation.

We performed high resolution hydrodynamical simulations of dense cool clouds embedded in supernova-driven galactic superwinds. Here we present preliminary results of our reference model in which the effect of heat conduction are taken into account. Significant dynamical differences occur between simulations with and without heat conduction. In absence of heat conduction the cloud fragments in few dynamical timescale. The inclusion of heat conduction has the effect to stabilize the cloud and inhibit the growth of Kelvin-Helmholtz and Reyleigh-Taylor instabilities. Furthermore in the conditions met in our simulations the strong heat flux at the cloud edge generates a converging shock which compresses the cloud. We also calculate the high energy emission (O{\sc vi} and soft X-ray) of the cloud and O{\sc vi} absorption line properties and compare the results with obvervations. Models in which heat conduction is taken into account seem to fit the observations much better. In general only a small fraction (0.1-0.4%) of the wind mechanical energy intersecting the cloud is radiated away. Finally some of our models are able to explain the low metallicity abundance seen in X-ray osservation of superwinds.

We study the dynamical and chemical evolution of a galaxy similar to IZw18 under the assumption of a continuous star formation during bursts. We adopt a 2-D hydrocode coupled with detailed chemical yields originating from SNeII, SNeIa and from single intermediate-mass stars. Different nucleosynthetic yields and different IMF slopes are tested. In most of the explored cases, a galactic wind develops, mostly carrying out of the galaxy the metal-enriched gas produced by the burst itself. The chemical species with the largest escape probabilities are Fe and N. Consequently, we predict that the [$\alpha$/Fe] and [$\alpha$/N] ratios outside the galaxy are lower than inside. In order to reproduce the chemical composition of IZw18, the best choice seems to be the adoption of the yields of Meynet & Maeder (2002) which take into account stellar rotation, although these authors do not follow the whole evolution of all the stars. Models with a flat IMF (x=0.5) seem to be able to better reproduce the chemical properties of IZw18, but they inject in the gas a much larger amount of energy and the resulting galactic wind is very strong, at variance with observations. We also predict the evolution of the abundances in the \hi medium and compare them with recent {\sl FUSE} observations.

We study the hydrodynamical evolution of galactic winds in disky dwarf galaxies moving through an intergalactic medium. In agreement with previous investigations,we find that when the ram pressure stripping does not disrupt the ISM, it usually has a negligible effect on the galactic wind dynamics. Only when the IGM ram pressure is comparable to the central ISM thermal pressure the stripping and the superwind influence each other increasing the gas removal rate. In this case several parameters regulate the ISM ejection process, as the original distribution of the ISM and the geometry of the IGM-galaxy interaction. When the ISM is not removed by the ram pressure or the wind, it loses memory of the starburst episode and recovers almost its pre-burst distribution in a timescale of 50-200 Myr. After this time another star formation episode becomes, in principle, possible. Evidently, galactic winds are consistent with a recurrent bursts star formation history. Contrary to the ISM content, the amount of the metal-rich ejecta retained by the galaxy is more sensitive to the ram pressure action. Part of the ejecta is first trapped in a low density, extraplanar gas produced by the IGM-ISM interaction, and then pushed back onto the galactic disc. The amount of trapped metals in a moving galaxy may be up to three times larger than in a galaxy at rest. This prediction may be tested comparing metallicity of dwarf galaxies in nearby poor clusters or groups, such as Virgo or Fornax, with the field counterpart. The sensitivity of the metal entrapment efficiency on the geometry of the interaction may explain part of the observed scatter in the metallicity-luminosity relation for dwarf galaxies.

We present 3D hydrodynamic simulations of ram pressure stripping in dwarf galaxies. Analogous studies on this subject usually deal with much higher ram pressures, typical of galaxy clusters, or mild ram pressure due to the gas halo of the massive galactic companions. We extend over previous investigations by considering flattened, rotating dwarf galaxies subject to ram pressures typical of poor galaxy groups. We study the ram pressure effects as a function of several parameters such as galactic mass and velocity, ambient gas density, and angle between the galactic plane and the direction of motion. It turns out that this latter parameter plays a role only when the gas pressure in the galactic centre is comparable to the ram pressure. Despite the low values of the ram pressure, some dwarf galaxies can be completely stripped after 1-2 hundred of million years. This pose an interesting question on the aspect of the descents and, more in general, on the morphological evolution of dwarf galaxies. In cases in which the gas is not completely stripped, the propagation of possible galactic wind may be influenced by the disturbed distribution of the interstellar matter. We also consider the modification of the ISM surface density induced by the ram pessure and find that the resulting compression may trigger star formation over long time spans.

We present 3D hydrodynamical simulations of galactic winds in dwarf, gas-rich galaxies. The galaxy is moving through the ICM of a small galaxy group at v=200 km/s. The ram pressure removes the galactic gas at large radii, but does not strongly influence the ISM near the center. A starburst generates a galactic wind. The newly produced metals are expelled in the ICM and carried to large distance from the galaxy by the ram pressure. 500 Myr after the starburst only a few percent of the heavy elements produced are present in the central region of the dwarf galaxy. A large collection of ram pressure + wind models will be presented in a forthcoming paper.

In this paper we present some results concerning the effects of two instantaneous starbursts, separated by a quiescent period, on the dynamical and chemical evolution of blue compact dwarf galaxies. In particular, we compare the model results to the galaxy IZw18, which is a very metal-poor, gas-rich dwarf galaxy, possibly experiencing its first or second burst of star formation. We follow the evolution of a first weak burst of star formation followed by a second more intense one occurring after several hundreds million years. We find that a galactic wind develops only during the second burst and that metals produced in the burst are preferentially lost relative to the hydrogen gas. We predict the evolution of several chemical abundances (H, He, C, N, O, \alpha-elements, Fe) in the gas inside and outside the galaxy, by taking into account in detail the chemical and energetical contributions from type II and Ia supernovae. We find that the abundances predicted for the star forming region are in good agreement with the HII region abundances derived for IZw18. We also predict the abundances of C, N and O expected for the HI gas to be compared with future FUSE abundance determinations. We conclude that IZw18 must have experienced two bursts of star formation, one occurred \sim 300 Myr ago and a present one with an age between 4-7 Myr. However, by taking into account also other independent estimates, such as the color-magnitude diagram and the spectral energy distribution of stars in IZw18, and the fact that real starbursts are not instantaneous, we suggest that it is more likely that the burst age is between 4 and 15 Myr.

In this paper we present some results of numerical simulations concerning the development of galactic winds in starburst galaxies. In particular, we focus on a galaxy similar to IZw18, the most metal-poor galaxy locally known. We compute the chemio-dynamical evolution of this galaxy, considering the energetic input and the chemical yields originating from Supernovae (SNe) of Type II and Ia and from intermediate-mass stars. We consider both single, instantaneous starburst and two starbursts separated by a quiescent period. In all considered cases a metal enriched winds develops and in particular the metals produced by Type Ia SNe are ejected more efficiently than the other metals. We suggest that two burst of star formation, the first being weaker and the last having an age of some tenth of Myr, can satisfactorily reproduce the abundances and abundance ratios found in literature for IZw18

We study the effect of a single, instantaneous starburst on the dynamical and chemical evolution of a gas-rich dwarf galaxy, whose potential well is dominated by a dark matter halo. We follow the dynamical and chemical evolution of the ISM by means of an improved 2-D hydrodynamical code coupled with detailed chemical yields originating from type II SNe, type Ia SNe and single low and intermediate mass stars (IMS). In particular we follow the evolution of the abundances of H, He, C, N, O, Mg, Si and Fe. We find that for a galaxy resembling IZw18, a galactic wind develops as a consequence of the starburst and it carries out of the galaxy mostly the metal-enriched gas. In addition, we find that different metals are lost differentially in the sense that the elements produced by type Ia SNe are more efficiently lost than others. As a consequence of that we predict larger [$\alpha$/Fe] ratios for the gas inside the galaxy than for the gas leaving the galaxy. A comparison of our predicted abundances of C, N, O and Si in the case of a burst occurring in a primordial gas shows a very good agreement with the observed abundances in IZw18 as long as the burst has an age of $\sim 31$ Myr and IMS produce some primary nitrogen. However, we cannot exclude that a previous burst of star formation had occurred in IZw18 especially if the preenrichment produced by the older burst was lower than $Z=0.01$ Z$_{\odot}$. Finally, at variance with previous studies, we find that most of the metals reside in the cold gas phase already after few Myr. This result is mainly due to the assumed low SNII heating efficiency, and justifies the generally adopted homogeneous and instantaneous mixing of gas in chemical evolution models.

We study the effect of a single, instantaneous starburst in a gas-rich dwarf galaxy on the dynamical and chemical evolution of its interstellar medium. We consider the energetic input and the chemical yields originating from SNeII, SNeIa and intermediate-mass stars. We find that a galaxy resembling IZw18 develops a galactic wind carrying out mostly the metal-rich gas. The various metals are lost differentially and the metals produced by the SNeIa are lost more efficiently than the others. As a consequence, we find larger [$\alpha$/Fe] ratios for the gas inside the galaxy than for the gas leaving the galaxy. Finally we find that a single burst occurring in primordial gas (without pre-enrichment), gives chemical abundances and dynamical structures in good agreement with what observed in IZw18 after $\sim$ 29 Myr from the beginning of star formation.

During a Hubble time, cluster galaxies may undergo several mutual encounters close enough to gravitationally perturb their hot, X-ray emitting gas flows. We ran several 2D, time dependent hydrodynamical models to investigate the effects of such perturbations on the gas flow inside elliptical galaxies. In particular, we studied in detail the modifications occurring in the scenario proposed by D'Ercole et al. (1989), in which the galactic interstellar medium produced by the aging galactic stellar population, is heated by SNIa at a decreasing rate. We find that, although the tidal interaction in our models lasts less than 1 Gyr, its effect extends over several Gyrs. The tidally induced turbulent flows create dense filaments which cool quickly and accrete onto the galactic center, producing large spikes in the global Lx. Once this mechanism starts, it is fed by gravity and amplified by SNIa. In cooling flow models without supernovae the amplitude of the Lx fluctuations due to the tidal interaction is substantially reduced. We conclude that, if SNIa significantly contribute to the energetics of the gas flows in ellipticals, then the observed spread in the Lx-Lb diagram may be caused, at least in part, by this mechanism. On the contrary, tidal interactions cannot be responsible for the observed spread if the pure cooling flow scenario applies (abridged).

We study, through 2D hydrodynamical simulations, the feedback of a starburst on the ISM of typical gas rich dwarf galaxies. The main goal is to address the circulation of the ISM and metals following the starburst. We assume a single-phase rotating ISM in equilibrium in the galactic potential generated by a stellar disk and a spherical dark halo. The starburst is assumed to occur in a small volume in the center of the galaxy, and it generates a mechanical power of 3.8e39 erg/s or 3.8e40 erg/s for 30 Myr. We found, consistently with previous investigations, that the galactic wind is not very effective in removing the ISM. The metal rich stellar ejecta, instead, may be efficiently expelled from the galaxy and dispersed in the intergalactic medium. Moreover, we found that the central region of the galaxy is always replenished with cold and dense gas after a few 100 Myr from the starbust, achieving the requisite for a new star formation event in 0.5-1 Gyr. The hydrodynamical evolution of galactic winds is thus consistent with the episodic star formation regime suggested by many chemical evolution studies. We also discuss the X-ray emission of these galaxies and find that the observable (emission averaged) abundance of the hot gas underestimates the real one if thermal conduction is effective. This could explain the very low hot gas metallicities estimated in starburst galaxies.

A chemical evolution model following the evolution of the abundances of H, He, C, N, O and Fe for dwarf irregular and blue compact galaxies is presented. This model takes into account detailed nucleosynthesis and computes in detail the rates of supernovae of type II and I. The star formation is assumed to have proceeded in short but intense bursts. The novelty relative to previous models is that the development of a galactic wind is studied in detail by taking into account the energy injected into the interstellar medium (ISM) from both supernovae and stellar winds from massive stars as well as the presence of dark matter halos. Both metal enriched and normal winds have been considered. Our main conclusions are: i) a substantial amount of dark matter (from 1 to 50 times larger than the luminous matter) is required in order to avoid the complete destruction of such galaxies during strong starbursts, and ii) the energy injected by stellar winds and type Ia supernovae into the ISM is negligible relative to the total thermal energy, and in particular to the type II supernovae, which in fact, dominate the energetics during starbursts.

A recent analysis of the "Einstein" sample of early-type galaxies has revealed that at any fixed optical luminosity Lb S0 galaxies have lower mean X-ray luminosity Lx per unit Lb than ellipticals. Following a previous analytical investigation of this problem (Ciotti & Pellegrini 1996), we have performed 2D numerical simulations of the gas flows inside S0 galaxies in order to ascertain the effectiveness of rotation and/or galaxy flattening in reducing the Lx/Lb ratio. The flow in models without SNIa heating is considerably ordered, and essentially all the gas lost by the stars is cooled and accumulated in the galaxy center. If rotation is present, the cold material settles in a disk on the galactic equatorial plane. Models with a time decreasing SNIa heating host gas flows that can be much more complex. After an initial wind phase, gas flows in energetically strongly bound galaxies tend to reverse to inflows. This occurs in the polar regions, while the disk is still in the outflow phase. In this phase of strong decoupling, cold filaments are created at the interface between inflowing and outflowing gas. Models with more realistic values of the dynamical quantities are preferentially found in the wind phase with respect to their spherical counterparts of equal Lb. The resulting Lx of this class of models is lower than in spherical models with the same Lb and SNIa heating. At variance with cooling flow models, rotation is shown to have only a marginal effect in this reduction, while the flattening is one of the driving parameters for such underluminosity, in accordance with the analytical investigation.