HARPS spectra with signal-to-noise ratios S/N > 600 at 6000 A were analysed with MARCS model atmospheres to obtain 1D LTE abundances of C, O, Na, Mg, Al, Si, Ca, Ti, Cr, Fe, Ni, Sr, and Y for 72 nearby solar-type stars with metallicities in the range -0.3 < [Fe/H] < +0.3 and ASTEC stellar models were used to determine stellar ages from effective temperatures, luminosities obtained via Gaia DR2 parallaxes, and heavy element abundances. The resulting age-metallicity distribution appears to consist of two distinct populations: a sequence of old stars with a steep rise of [Fe/H] to about +0.3 dex at an age of ~7 Gyr and a younger sequence with [Fe/H] increasing from -0.3 dex to about +0.2 dex over the last six Gyr. Furthermore, the trends of several abundance ratios, [O/Fe], [Na/Fe], [Ca/Fe], and [Ni/Fe], as a function of stellar age split into two corresponding sequences. The [Y/Mg]-age relation, on the other hand, shows no offset between the two age sequences and has no significant dependence on [Fe/H], but the components of a visual binary star, zeta Reticuli, have a large and puzzling deviation. The split of the age-metallicity distribution into two sequences may be interpreted as evidence of two episodes of accretion of gas onto the Galactic disk with a quenching of star formation in between. Some of the [X/Fe]-age relations support this scenario but other relations are not so easy to explain, which calls for a deeper study of systematic errors in the derived abundances as a function of [Fe/H], in particular 3D non-LTE effects.

The trends of chemical abundances and abundance ratios observed in stars of different ages, kinematics, and metallicities bear the imprints of several physical processes that concur to shape the host galaxy properties. By inspecting these trends, we get precious information on stellar nucleosynthesis, the stellar mass spectrum, the timescale of structure formation, the efficiency of star formation, as well as any inward or outward flows of gas. In this paper, we analyse recent determinations of carbon-to-iron and carbon-to-oxygen abundance ratios in different environments (the Milky Way and elliptical galaxies) using our latest chemical evolution models that implement up-to-date stellar yields and rely on the tight constraints provided by asteroseismic stellar ages (whenever available). A scenario where most carbon is produced by rotating massive stars, with yields largely dependent on the metallicity of the parent proto-star clouds, allows us to fit simultaneously the high-quality data available for the local Galactic components (thick and thin discs) and for microlensed dwarf stars in the Galactic bulge, as well as the abundance ratios inferred for massive elliptical galaxies. Yet, more efforts are needed from both observers and theoreticians in order to base these conclusions on firmer grounds.

We present a detailed study of the composition of 20 M giants in the Galactic center with 15 of them confirmed to be in the Nuclear Star Cluster. As a control sample we have also observed 7 M giants in the Milky Way Disk with similar stellar parameters. All 27 stars are observed using the NIRSPEC spectograph on the KECK II telescope in the K-band at a resolving power of R=23,000. We report the first silicon abundance trends versus [Fe/H] for stars in the Galactic center. While finding a disk/bulge like trend at subsolar metallicities, we find that [Si/Fe] is enhanced at supersolar metallicities. We speculate on possible enrichment scenarios to explain such a trend. However, the sample size is modest and the result needs to be confirmed by additional measurements of silicon and other \textalpha-elements. We also derive a new distribution of [Fe/H] and find the most metal rich stars at [Fe/H]=+0.5 dex, confirming our earlier conclusions that the Galactic center hosts no stars with extreme chemical composition.

The cosmic origin of fluorine is still not well constrained. Several nucleosynthetic channels at different phases of stellar evolution have been suggested, but these must be constrained by observations. For this, the fluorine abundance trend with metallicity spanning a wide range is required. Our aim is to determine stellar abundances of fluorine for -1.1<[Fe/H]<+0.4. We determine the abundances from HF lines in infrared K-band spectra (at 2.3 microns) of cool giants, observed with the IGRINS and Phoenix high-resolution spectrographs. We derive accurate stellar parameters for all our observed K giants, which is important since the HF lines are very temperature sensitive. We find that [F/Fe] is flat as a function of metallicity at [F/Fe]~0, but increases as the metallicity increases. The fluorine slope shows a clear secondary behavior in this metallicity range. We also find that the [F/Ce] ratio is relatively flat for -0.6<[Fe/H]<0, and that for two metal-poor ([Fe/H]<-0.8), s-process element enhanced giants, we do not detect an elevated fluorine abundance. We interpret all these observational constraints to indicate that several major processes are at play for the cosmic budget of fluorine over time; from those in massive stars at low metallicities, through the asymptotic giant branch-star contribution at -0.6<[Fe/H]<0, to processes with increasing yields with metallicity at super-solar metallicities. The origins of the latter, and whether or not Wolf-Rayet stars and/or novae could contribute at super-solar metallicities, is currently not known. To quantify these observational results, theoretical modelling is required. More observations in the metal-poor region are required to clarify the processes there.

In their evolution, star-forming galaxies are known to follow scaling relations between some fundamental physical quantities, such as the mass-metallicity and the main sequence relations. We aim at studying the evolution of galaxies that, at a given redshift, lie simultaneously on the mass-metallicity (MZ) and main sequence (MS) relations. To this aim, we use the analytical, 'leaky-box' chemical evolution model of Spitoni et al. (2017), in which galaxy evolution is described by an infall timescale $\tau$ and a wind efficiency $\lambda$. We provide a detailed analysis of the temporal evolution of galactic metallicity, stellar mass, mass-weighted age and gas fraction. The evolution of the galaxies lying on the MZ and MS at $z\sim0.1$ suggests that the average infall time-scale in two different bins of stellar masses ($M_{\star}<10^{10} M_{\odot}$ and $M_{\star}>10^{10} M_{\odot}$) decreases with decreasing redshift. This means that at each redshift, only the youngest galaxies can be assembled on the shortest timescales and still belong to the MS. In the lowest mass bin, a decrease of the median $\tau$ is accompanied by an increase of the median $\lambda$ value. This implies that systems which have formed at more recent times will need to eject a larger amount of mass to keep their metallicity at low values. Another important result is that galactic downsizing, as traced by the age-mass relation, is naturally recovered by imposing that local galaxies lie on both the MZ and MS. Finally, we study the evolution of the hosts of C IV-selected AGN, which at $z\sim 2$ follow a flat MZ, as found by Mignoli et al. (2019). If we impose that these systems lie on the MS, at lower redshifts we find an 'inverted' MZ, meaning that some additional processes must be at play in their evolution.

Since the onset of the `space revolution' of high-precision high-cadence photometry, asteroseismology has been demonstrated as a powerful tool for informing Galactic archaeology investigations. The launch of the NASA TESS mission has enabled seismic-based inferences to go full sky -- providing a clear advantage for large ensemble studies of the different Milky Way components. Here we demonstrate its potential for investigating the Galaxy by carrying out the first asteroseismic ensemble study of red giant stars observed by TESS. We use a sample of 25 stars for which we measure their global asteroseimic observables and estimate their fundamental stellar properties, such as radius, mass, and age. Significant improvements are seen in the uncertainties of our estimates when combining seismic observables from TESS with astrometric measurements from the Gaia mission compared to when the seismology and astrometry are applied separately. Specifically, when combined we show that stellar radii can be determined to a precision of a few percent, masses to 5-10% and ages to the 20% level. This is comparable to the precision typically obtained using end-of-mission Kepler data

With the wealth of information from large surveys and observational campaigns in the contemporary era, it is critical to properly exploit the data to constrain the parameters of Galactic chemical evolution models and quantify the associated uncertainties. We aim at constraining the two-infall chemical evolution models for the solar annulus using the measured chemical abundance ratios and seismically inferred age of stars in the APOKASC sample. In the revised two-infall chemical evolution models by Spitoni et al. (2019), a significant delay of $\sim 4.3$ Gyr has been invoked between the two episodes of gas accretion. In this work, we wish to test its robustness and statistically confirm/quantify the delay. For the first time, a Bayesian framework based on Markov Chain Monte Carlo methods has been used for fitting the two-infall chemical evolution models to the data. In addition to fitting the data for stars in the APOKASC sample, our best fit models also reproduce other important observational constraints of the chemical evolution of the disk: i) present day stellar surface density; ii) present day supernova and star formation rates; iii) the metallicity distribution function; and iv) solar abundance values. We found a significant delay between the two gas accretion episodes for various models explored with different values for the star formation efficiencies. The values for the delay lie in the range $4.5-5.5$ Gyr. The results suggest that the APOKASC sample carries the signature of delayed gas-rich merger, with the dilution as main process determining the shape of low-$\alpha$ stars in the abundance ratios space.

Galactic disc chemical evolution models generally ignore azimuthal surface density variation that can introduce chemical abundance azimuthal gradients. Recent observations, however, have revealed chemical abundance changes with azimuth in the gas and stellar components of both the Milky Way and external galaxies. To quantify the effects of spiral arm density fluctuations on the azimuthal variations of the oxygen and iron abundances in disc galaxies. We develop a new 2D galactic disc chemical evolution model, capable of following not just radial but also azimuthal inhomogeneities. The density fluctuations resulting from a Milky Way-like N-body disc formation simulation produce azimuthal variations in the oxygen abundance gradients of the order of 0.1 dex. Moreover, in agreement with the most recent observations in external galaxies, the azimuthal variations are more evident in the outer galactic regions. Using a simple analytical model, we show that the largest fluctuations with azimuth result near the spiral structure corotation resonance, where the relative speed between spiral and gaseous disc is the slowest. In conclusion we provided a new 2D chemical evolution model capable of following azimuthal density variations. Density fluctuations extracted from a Milky Way-like dynamical model lead to a scatter in the azimuthal variations of the oxygen abundance gradient in agreement with observations in external galaxies. We interpret the presence of azimuthal scatter at all radii by the presence of multiple spiral modes moving at different pattern speeds, as found in both observations and numerical simulations.

In this work, we study the formation and chemical evolution of the Galactic bulge with particular focus on the abundance pattern ([Mg/Fe] vs. [Fe/H]), metallicity and age distribution functions. We consider detailed chemical evolution models for the Galactic bulge and inner disc, with the aim of shedding light on the connection between these components and the origin of bulge stars. In particular, we first present a model assuming a fast and intense star formation, with the majority of bulge stars forming on a timescale less than 1 Gyr. Then we analyze the possibility of two distinct stellar populations in the bulge, as suggested by Gaia-ESO and APOGEE data. These two populations, one metal poor and the other metal rich, can have had two different origins: i) the metal rich formed after a stop of roughly 250 Myr in the star formation rate of the bulge, or ii) the metal rich population is made of stars formed in the inner disc and brought into the bulge by the early secular evolution of the bar. We also examine the case of multiple star bursts in the bulge with consequent formation of multiple populations, as suggested by studies of microlensed stars. After comparing model results and observations, we suggest that the most likely scenario is that there are two main stellar populations, both made mainly by old stars (> 10 Gyr), with the metal rich and younger one formed from inner thin disc stars, in agreement with kinematical arguments. However, on the basis of dynamical simulations, we cannot completely exclude that the second population formed after a stop in the star formation during the bulge evolution, so that all the stars formed in situ.

We present the first chemical evolution model for Enceladus, alias the Gaia Sausage, to investigate the star formation history of one of the most massive satellites accreted by the Milky Way during a major merger event. Our best chemical evolution model for Enceladus nicely fits the observed stellar [$\alpha$/Fe]-[Fe/H] chemical abundance trends, and reproduces the observed stellar metallicity distribution function, by assuming low star formation efficiency, fast infall time scale, and mild outflow intensity. We predict a median age for Enceladus stars $12.33^{+0.92}_{-1.36}$ Gyr, and - at the time of the merger with our Galaxy ($\approx10$ Gyr ago from Helmi et al.) - we predict for Enceladus a total stellar mass $M_{\star} \approx 5 \times 10^{9}\,\text{M}_{\odot}$. By looking at the predictions of our best model, we discuss that merger events between the Galaxy and systems like Enceladus may have inhibited the gas accretion onto the Galaxy disc at high redshifts, heating up the gas in the halo. This scenario could explain the extended period of quenching in the star formation activity of our Galaxy about 10 Gyr ago, which is predicted by Milky Way chemical evolution models, in order to reproduce the observed bimodality in [$\alpha$/Fe]-[Fe/H] between thick- and thin-disc stars.

Precise stellar ages from asteroseismology have become available and can help setting stronger constraints on the evolution of the Galactic disc components. Recently, asteroseismology has confirmed a clear age difference in the solar annulus between two distinct sequences in the [$\alpha$/Fe] versus [Fe/H] abundance ratios relation: the high-$\alpha$ and low-$\alpha$ stellar populations. We aim at reproducing these new data with chemical evolution models including different assumptions for the history and number of accretion events. We tested two different approaches: a revised version of the `two-infall' model where the high-$\alpha$ phase forms by a fast gas accretion episode and the low-$\alpha$ sequence follows later from a slower gas infall rate, and the parallel formation scenario where the two disc sequences form coevally and independently. The revised `two-infall' model including uncertainties in age and metallicity is capable of reproducing: i) the [$\alpha$/Fe] vs. [Fe/H] abundance relation at different Galactic epochs, ii) the age$-$metallicity relation and the time evolution [$\alpha$/Fe]; iii) the age distribution of the high-$\alpha$ and low-$\alpha$ stellar populations, iv) the metallicity distribution function. The parallel approach is not capable of properly reproduce the stellar age distribution, in particular at old ages. In conclusion, the best chemical evolution model is the revised `two-infall' one, where a consistent delay of $\sim$4.3 Gyr in the beginning of the second gas accretion episode is a crucial assumption to reproduce stellar abundances and ages.

We tested the implementation of different IMFs in our model for the chemical evolution of ellipticals, with the aim of reproducing the observed relations of [Fe/H] and [Mg/Fe] abundances with galaxy mass in a sample of early-type galaxies selected from the SPIDER-SDSS catalog. Abundances in the catalog were derived from averaged spectra, obtained by stacking individual spectra according to central velocity dispersion, as a proxy of galaxy mass. We tested initial mass functions already used in a previous work, as well as two new models, based on low-mass tapered ("bimodal") IMFs, where the IMF becomes either (1) bottom-heavy in more massive galaxies, or (2) is time-dependent, switching from top-heavy to bottom-heavy in the course of galactic evolution. We found that observations could only be reproduced by models assuming either a constant, Salpeter IMF, or a time-dependent distribution, as other IMFs failed. We further tested the models by calculating their M/L ratios. We conclude that a constant, time-independent bottom-heavy IMF does not reproduce the data, especially the increase of the $[\alpha/Fe]$ ratio with galactic stellar mass, whereas a variable IMF, switching from top to bottom-heavy, can match observations. For the latter models, the IMF switch always occurs at the earliest possible considered time, i.e. $t_{\text{switch}}= 0.1$ Gyr.

In this paper, we study the formation and chemical evolution of the Milky Way disc with particular focus on the abundance patterns ([$\alpha$/Fe] vs. [Fe/H]) at different Galactocentric distances, the present-time abundance gradients along the disc and the time evolution of abundance gradients. We consider the chemical evolution models for the Galactic disc developed by Grisoni et al. (2017) for the solar neighborhood, both the two-infall and the one-infall ones, and we extend our analysis to the other Galactocentric distances. In particular, we examine the processes which mainly influence the formation of the abundance gradients: the inside-out scenario, a variable star formation efficiency, and radial gas flows. We compare our model results with recent abundance patterns obtained along the Galactic disc from the APOGEE survey and with abundance gradients observed from Cepheids, open clusters, HII regions and PNe. We conclude that the inside-out scenario is a key ingredient, but cannot be the only one to explain abundance patterns at different Galactocentric distances and abundance gradients. Further ingredients, such as radial gas flows and variable star formation efficiency, are needed to reproduce the observed features in the thin disc. The evolution of abundance gradients with time is also shown, although firm conclusions cannot still be drawn.

In the light of the new observational data related to fluorine abundances in the solar neighborhood stars, we present here chemical evolution models testing different fluorine nucleosynthesis prescriptions with the aim to best fit those new data related to the abundance ratios [F/O] vs. [O/H] and [F/Fe] vs. [Fe/H]. The adopted chemical evolution models are: i) the classical "two-infall" model which follows the chemical evolution of halo-thick disk and thin disk phases, ii) and the "one-infall" model designed only for the thin disk evolution. We tested the effects on the predicted fluorine abundance ratios of different nucleosynthesis yield sources: AGB stars, Wolf-Rayet stars, Type II and Type Ia supernovae, and novae. We find that the fluorine production is dominated by AGB stars but the Wolf-Rayet stars are required to reproduce the trend of the observed data in the solar neighborhood by J\"onsson et al. (2017a) with our chemical evolution models. In particular, the best model both for the "two-infall" and "one-infall" cases requires an increase by a factor of two of the Wolf-Rayet yields given by Meynet & Arnould (2000). We also show that the novae, even if their yields are still uncertain, could help to better reproduce the secondary behavior of F in the [F/O] vs. [O/H] relation. The inclusion of the fluorine production by Wolf-Rayet stars seems to be essential to reproduce the observed ratio [F/O] vs [O/H] in the solar neighborhood by J\"onsson et al. (2017a). Moreover, the inclusion of novae helps substantially to reproduce the observed fluorine secondary behavior.

We study the chemical evolution of the thick and thin discs of the Galaxy by comparing detailed chemical evolution models with recent data from the AMBRE Project. The data suggest that the stars in the thick and thin discs form two distinct sequences with the thick disc stars showing higher [{\alpha}/Fe] ratios. We adopt two different approaches to model the evolution of thick and thin discs. In particular, we adopt: i) a two-infall approach where the thick disc forms fast and before the thin disc and by means of a fast gas accretion episode, whereas the thin disc forms by means of a second accretion episode on a longer timescale; ii) a parallel approach, where the two discs form in parallel but at different rates. By comparing our model results with the observed [Mg/Fe] vs. [Fe/H] and the metallicity distribution functions in the two Galactic components, we conclude that the parallel approach can account for a group of {\alpha}-enhanced metal rich stars present in the data, whereas the two-infall approach cannot explain these stars unless they are the result of stellar migration. In both approaches, the thick disc has formed on a timescale of accretion of 0.1 Gyr, whereas the thin disc formed on a timescale of 7 Gyr in the solar region. In the two-infall approach a gap in star formation between the thick and thin disc formation of several hundreds of Myr should be present, at variance with the parallel approach where no gap is present.

We test the hypothesis that the classical and ultra-faint dwarf spheroidal satellites of the our Galaxy have been the building blocks of the Galactic halo by comparing their [O/Fe] and [Ba/Fe] vs. [Fe/H] patterns with the ones observed in Galactic halo stars. The [O/Fe] ratio deviates substantially from the observed abundance ratios in the Galactic halo stars for [Fe/H] > -2 dex, while they overlap for lower metallicities. On the other hand, for the neutron capture elements, the discrepancy is extended at all the metallicities, suggesting that the majority of stars in the halo are likely to have been formed in situ. We present the results for a model considering the effects of an enriched gas stripped from dwarf satellites on the chemical evolution of the Galactic halo. We find that the resulting chemical abundances of the halo stars depend on the adopted infall time-scale, and the presence of a threshold in the gas for star formation.

We will discuss some highlights concerning the chemical evolution of our Galaxy, the Milky Way. First we will describe the main ingredients necessary to build a model for the chemical evolution of the Milky Way. Then we will illustrate some Milky Way models which includes detailed stellar nucleosynthesis and compute the evolution of a large number of chemical elements, including C, N, O, $\alpha$-elements, Fe and heavier. The main observables and in particular the chemical abundances in stars and gas will be considered. A comparison theory-observations will follow and finally some conclusions from this astroarchaeological approach will be derived.

We explore the hypothesis that the classical and ultra-faint dwarf spheroidal satellites of the Milky Way have been the building blocks of the Galactic halo by comparing their [O/Fe] and [Ba/Fe] versus [Fe/H] patterns with the ones observed in Galactic halo stars. Oxygen abundances deviate substantially from the observed abundances in the Galactic halo stars for [Fe/H] values larger than -2 dex, while they overlap for lower metallicities. On the other hand, for the [Ba/Fe] ratio the discrepancy is extended at all [Fe/H] values, suggesting that the majority of stars in the halo are likely to have been formed in situ. Therefore, we suggest that [Ba/Fe] ratios are a better diagnostic than [O/Fe] ratios. Moreover, we show the effects of an enriched infall of gas with the same chemical abundances as the matter ejected and/or stripped from dwarf satellites of the Milky Way on the chemical evolution of the Galactic halo. We find that the resulting chemical abundances of the halo stars depend on the assumed infall time scale, and the presence of a threshold in the gas for star formation.

The Galactic habitable zone is defined as the region with highly enough metallicity to form planetary systems in which Earth-like planets could be born and might be capable of sustaining life surviving to the destructive effects of nearby supernova explosion events. Galactic chemical evolution models can be useful tools for studying the galactic habitable zones in different systems. Our aim here is to find the Galactic habitable zone using chemical evolution models for the Milky Way disc, adopting the most recent prescriptions for the evolution of dust and for the probability of finding planetary systems around M and FGK stars. Moreover, for the first time, we will express those probabilities in terms of the dust-to-gas ratio of the ISM in the solar neighborhood as computed by detailed chemical evolution models. At a fixed Galactic time and Galactocentric distance we determine the number of M and FGK stars having Earths (but no gas giant planets) which survived supernova explosions, using the formalism of our Paper I. The probabilities of finding terrestrial planets but not gas giant planets around M stars deviate substantially from the ones around FGK stars for supersolar values of [Fe/H]. For both FGK and M stars the maximum number of stars hosting habitable planets is at 8 kpc from the Galactic Centre, if destructive effects by supernova explosions are taken into account. At the present time the total number of M stars with habitable planets are $\simeq$ 10 times the number of FGK stars. Moreover, we provide a sixth order polynomial fit (and a linear one but more approximated) for the relation found with chemical evolution models in the solar neighborhood between the [Fe/H] abundances and the dust-to-gas ratio.

Abridged: We used the fourth internal data release of the Gaia-ESO survey to characterize the bulge chemistry, spatial distribution, kinematics, and to compare it chemically with the thin and thick disks. The sample consist on ~2500 red clump stars in 11 bulge fields ($-10^\circ\leq l\leq+8^\circ$ and $-10^\circ\leq b\leq-4^\circ$), and a set of ~6300 disk stars selected for comparison. The bulge MDF is confirmed to be bimodal across the whole sampled area, with metal-poor stars dominating at high latitudes. The metal-rich stars exhibit bar-like kinematics and display a bimodality in their magnitude distribution, a feature which is tightly associated with the X-shape bulge. They overlap with the metal-rich end of the thin disk sequence in the [Mg/Fe] vs. [Fe/H] plane. Metal-poor bulge stars have a more isotropic hot kinematics and do not participate in the X-shape bulge. With similar Mg-enhancement levels, the position of the metal-poor bulge sequence "knee" is observed at [Fe/H]$_{knee}=-0.37\pm0.09$, being 0.06 dex higher than that of the thick disk. It suggests a higher SFR for the bulge than for the thick disk. Finally, we present a chemical evolution model that suitably fits the whole bulge sequence by assuming a fast ($<1$ Gyr) intense burst of stellar formation at early epochs. We associate metal-rich stars with the B/P bulge formed from the secular evolution of the early thin disk. On the other hand, the metal-poor subpopulation might be the product of an early prompt dissipative collapse dominated by massive stars. Nevertheless, our results do not allow us to firmly rule out the possibility that these stars come from the secular evolution of the early thick disk. This is the first time that an analysis of the bulge MDF and $\alpha$-abundances has been performed in a large area on the basis of a homogeneous, fully spectroscopic analysis of high-resolution, high S/N data.

Analytical models of chemical evolution, including inflow and outflow of gas, are important tools to study how the metal content in galaxies evolves as a function of time. In this work, we present new analytical solutions for the evolution of the gas mass, total mass and metallicity of a galactic system, when a decaying exponential infall rate of gas and galactic winds are assumed. We apply our model to characterize a sample of local star-forming and passive galaxies from the Sloan Digital Sky Survey data, with the aim of reproducing their observed mass-metallicity relation; in this way, we can derive how the two populations of star-forming and passive galaxies differ in their particular distribution of ages, formation time scales, infall masses and mass loading factors. We find that the local passive galaxies are on average older and assembled on shorter typical time-scales than the local star-forming ones; on the other hand, the larger mass star-forming galaxies show generally older ages and longer typical formation time-scales compared with the smaller mass star-forming galaxies. Finally, we conclude that the local star-forming galaxies experience stronger galactic winds than the passive galaxy population. We explore the effect of assuming different initial mass functions in our model, showing that to reproduce the observed mass-metallicity relation stronger winds are requested if the initial mass function is top-heavy. Finally, our analytical models predict the assumed sample of local galaxies to lie on a tight surface in the 3D space defined by stellar metallicity, star formation rate and stellar mass, thus mimicking the well-known "fundamental relation".

We present a theoretical method for solving the chemical evolution of galaxies, by assuming an instantaneous recycling approximation for chemical elements restored by massive stars and the Delay Time Distribution formalism for the delayed chemical enrichment by Type Ia Supernovae. The galaxy gas mass assembly history, together with the assumed stellar yields and initial mass function, represent the starting point of this method. We derive a simple and general equation which closely relates the Laplace transforms of the galaxy gas accretion history and star formation history, which can be used to simplify the problem of retrieving these quantities in the galaxy evolution models assuming a linear Schmidt-Kennicutt law. We find that - once the galaxy star formation history has been reconstructed from our assumptions - the differential equation for the evolution of the chemical element $X$ can be suitably solved with classical methods. We apply our model to reproduce the [O/Fe] and [Si/Fe] vs. [Fe/H] chemical abundance patterns as observed at the solar neighborhood, by assuming a decaying exponential infall rate of gas and different delay time distributions for Type Ia Supernovae; we also explore the effect of assuming a nonlinear Schmidt-Kennicutt law, with the index of the power law being $k=1.4$. Although approximate, we conclude that our model with the single degenerate scenario for Type Ia Supernovae provides the best agreement with the observed set of data. Our method can be used by other complementary galaxy stellar population synthesis models to predict also the chemical evolution of galaxies.

Asymptotic Giant Branch stars are known to produce `cosmic' fluorine but it is uncertain whether these stars are the main producers of fluorine in the Solar neighborhood or if any of the other proposed formation sites, type II supernovae and/or Wolf-Rayet stars, are more important. Recent articles have proposed both Asymptotic Giant Branch stars as well as type II supernovae as the dominant sources of fluorine in the Solar neighborhood. In this paper we set out to determine the fluorine abundance in a sample of 49 nearby, bright K-giants for which we previously have determined the stellar parameters as well as alpha abundances homogeneously from optical high-resolution spectra. The fluorine abundance is determined from a 2.3 $\mu$m HF molecular line observed with the spectrometer Phoenix. We compare the fluorine abundances with those of alpha elements mainly produced in type II supernovae and find that fluorine and the alpha-elements do not evolve in lock-step, ruling out type II supernovae as the dominating producers of fluorine in the Solar neighborhood. Furthermore, we find a secondary behavior of fluorine with respect to oxygen, which is another evidence against the type II supernovae playing a large role in the production of fluorine in the Solar neighborhood. This secondary behavior of fluorine will put new constraints on stellar models of the other two suggested production sites: Asymptotic Giant Branch stars and Wolf-Rayet stars.

According to the current cosmological cold dark matter paradigm, the Galactic halo could have been the result of the assemblage of smaller structures. Here we explore the hypothesis that the classical and ultra-faint dwarf spheroidal satellites of the Milky Way have been the building blocks of the Galactic halo by comparing their [$\alpha$/Fe] and [Ba/Fe] versus [Fe/H] patterns with the ones observed in Galactic halo stars. The $\alpha$ elements deviate substantially from the observed abundances in the Galactic halo stars for [Fe/H] values larger than -2 dex, while they overlap for lower metallicities. On the other hand, for the [Ba/Fe] ratio the discrepancy is extended at all [Fe/H] values, suggesting that the majority of stars in the halo are likely to have been formed in situ. Therefore, we suggest that [Ba/Fe] ratios are a better diagnostic than [$\alpha$/Fe] ratios. Moreover, for the first time we consider the effects of an enriched infall of gas with the same chemical abundances as the matter ejected and/or stripped from dwarf satellites of the Milky Way on the chemical evolution of the Galactic halo. We find that the resulting chemical abundances of the halo stars depend on the assumed infall time scale, and the presence of a threshold in the gas for star formation. In particular, in models with an infall timescale for the halo around 0.8 Gyr coupled with a threshold in the surface gas density for the star formation (4 $\mathrm{M}_{\odot}\,\mathrm{pc}^{-2}$), and the enriched infall from dwarf spheroidal satellites, the first halo stars formed show [Fe/H]$>$-2.4 dex. In this case, to explain [$\alpha$/Fe] data for stars with [Fe/H]$<$-2.4 dex we need stars formed in dSph systems.

In the last years, more and more interest has been devoted to analytical solutions, including inflow and outflow, to study the metallicity enrichment in galaxies. In this framework, we assume a star formation rate which follows a linear Schmidt law, and we present new analytical solutions for the evolution of the metallicity (Z) in galaxies. In particular, we take into account environmental effects including primordial and enriched gas infall, outflow, different star formation efficiencies, and galactic fountains. The enriched infall is included to take into account galaxy-galaxy interactions. Our main results can be summarized as: i) when a linear Schmidt law of star formation is assumed, the resulting time evolution of the metallicity Z is the same either for a closed-box model or for an outflow model. ii) The mass-metallicity relation for galaxies which suffer a chemically enriched infall, originating from another evolved galaxy with no pre-enriched gas, is shifted down in parallel at lower Z values, if compared the closed box model. iii) When a galaxy suffers at the same time a primordial infall and a chemically enriched one, the primordial infall always dominates the chemical evolution. iv) We present new solutions for the metallicity evolution in a galaxy which suffers galactic fountains and an enriched infall from another galaxy at the same time. The analytical solutions presented here can be very important to study the metallicity (oxygen), which is measured in high-redshift objects. These solutions can be very useful: a) in the context of cosmological semi-analytical models for galaxy formation and evolution, and b) for the study of compact groups of galaxies.

In the last years, stellar migration in galactic discs has been the subject of several investigations. However, its impact on the chemical evolution of the Milky Way still needs to be fully quantified. In this paper, we aim at imposing some constraints on the significance of this phenomenon by considering its influence on the chemical evolution of the Milky Way thin disc. We do not investigate the physical mechanisms underlying the migration of stars. Rather, we introduce a simple, heuristic treatment of stellar migration in a detailed chemical evolution model for the thin disc of the Milky Way, which already includes radial gas flows and reproduces several observational constraints for the solar vicinity and the whole Galactic disc. When stellar migration is implemented according to the results of chemo-dynamical simulations by Minchev et. al. (2013) and finite stellar velocities of 1 km s$^{-1}$ are taken into account, the high-metallicity tail of the metallicity distribution function of long-lived thin-disc stars is well reproduced. By exploring the velocity space, we find that the migrating stars must travel with velocities in the range 0.5 -2 km s$^{-1}$ to properly reproduce the high-metallicity tail of the metallicity distribution. We confirm previous findings by other authors that the observed spread in the age-metallicity relation of solar neighbourhood stars can be explained by the presence of stars which originated at different Galactocentric distances, and we conclude that the chemical properties of stars currently observed in the solar vicinity do suggest that stellar migration is present to some extent.

We present detailed chemical evolution models for the Milky Way and M31 in presence of radial gas flows. These models follow in detail the evolution of several chemical elements (H, He, CNO, $\alpha$ elements, Fe-peak elements) in space and time. The contribution of supernovae of different type to chemical enrichment is taken into account. We find that an inside-out formation of the disks coupled with radial gas inflows of variable speed can reproduce very well the observed abundance gradients in both galaxies. We also discuss the effects of other parameters, such as a threshold in the gas density for star formation and efficiency of star formation varying with galactic radius. Moreover, for the first time we compute the galactic habitable zone in our Galaxy and M31 in presence of radial gas flows. The main effect is to enhance the number of stars hosting a habitable planet with respect to the models without radial flow, in the region of maximum probability for this occurrence. In the Milky Way the maximum number of stars hosting habitable planets is at 8 kpc from the Galactic center, and the model with radial gas flows predicts a number of planets which is 38% larger than that predicted by the classical model.

We aim at reproducing the chemical evolution of the bulge of M31 by means of a detailed chemical evolution model, including radial gas flows coming from the disk. We study the impact of the initial mass function, the star formation rate and the time scale for bulge formation on the metallicity distribution function of stars. We compute several models of chemical evolution using the metallicity distribution of dwarf stars as an observational constraint for the bulge of M31. Then, by means of the model which best reproduces the metallicity distribution function, we predict the [X/Fe] vs. [Fe/H] relations for several chemical elements (O, Mg, Si, Ca, C, N). Our best model for the bulge of M31 is obtained by means of a robust statistical method and assumes a Salpeter initial mass function, a Schmidt-Kennicutt law for star formation with an exponent k=1.5, an efficiency of star formation of $\sim 15\pm 0.27\, Gyr^{-1}$, and an infall timescale of $\sim 0.10\pm 0.03$Gyr. Our results suggest that the bulge of M31 formed very quickly by means of an intense star formation rate and an initial mass function flatter than in the solar vicinity but similar to that inferred for the Milky Way bulge. The [$\alpha$/Fe] ratios in the stars of the bulge of M31 should be high for most of the [Fe/H] range, as is observed in the Milky Way bulge. These predictions await future data to be proven.

We review chemical evolution models developed for M31 as well as the abundance determinations available for this galaxy. Then we present a recent chemical evolution model for M31 including radial gas flows and galactic fountains along the disk, as well as a model for the bulge. Our models are predicting the evolution of the abundances of several chemical species such as H, He, C, N, O, Ne, Mg, Si, S, Ca and Fe. From comparison between model predictions and observations we can derive some constraints on the evolution of the disk and the bulge of M31. We reach the conclusions that Andromeda must have evolved faster than the Milky Way and inside-out, and that its bulge formed much faster than the disk on a timescale $\leq$ 0.5 Gyr. Finally, we present a study where we apply the model developed for the disk of M31 in order to study the probability of finding galactic habitable zones in this galaxy.

The galactic habitable zone is defined as the region with sufficient abundance of heavy elements to form planetary systems in which Earth-like planets could be born and might be capable of sustaining life, after surviving to close supernova explosion events. Galactic chemical evolution models can be useful for studying the galactic habitable zones in different systems. We apply detailed chemical evolution models including radial gas flows to study the galactic habitable zones in our Galaxy and M31. We compare the results to the relative galactic habitable zones found with "classical" (independent ring) models, where no gas inflows were included. For both the Milky Way and Andromeda, the main effect of the gas radial inflows is to enhance the number of stars hosting a habitable planet with respect to the "classical" model results, in the region of maximum probability for this occurrence, relative to the classical model results. These results are obtained by taking into account the supernova destruction processes. In particular, we find that in the Milky Way the maximum number of stars hosting habitable planets is at 8 kpc from the Galactic center, and the model with radial flows predicts a number which is 38% larger than what predicted by the classical model. For Andromeda we find that the maximum number of stars with habitable planets is at 16 kpc from the center and that in the case of radial flows this number is larger by 10 % relative to the stars predicted by the classical model.

We compute the abundance gradients along the disk of the Milky Way by means of the two-infall model: in particular, the gradients of oxygen and iron and their temporal evolution. First, we explore the effects of several physical processes which influence the formation and evolution of abundance gradients. They are: i) the inside-out formation of the disk, ii) a threshold in the gas density for star formation, iii) a variable star formation efficiency along the disk, iv) radial flows and their speed, and v) different total surface mass density (gas plus stars) distributions for the halo. We are able to reproduce at best the present day gradients of oxygen and iron if we assume an inside-out formation, no threshold gas density, a constant efficiency of star formation along the disk and radial gas flows. It is particularly important the choice of the velocity pattern for radial flows and the combination of this velocity pattern with the surface mass density distribution in the halo. Having selected the best model, we then explore the evolution of abundance gradients in time and find that the gradients in general steepen in time and that at redshift z~3 there is a gradient inversion in the inner regions of the disk, in the sense that at early epochs the oxygen abundance decreases toward the Galactic center. This effect, which has been observed, is naturally produced by our models if an inside-out formation of the disk and and a constant star formation efficiency are assumed. The inversion is due to the fact that in the inside-out formation a strong infall of primordial gas, contrasting chemical enrichment, is present in the innermost disk regions at early times. The gradient inversion remains also in the presence of radial flows, either with constant or variable speed in time, and this is a new result.

We compute the Type Ia supernova rates in typical elliptical galaxies by varying the progenitor models for Type Ia supernovae. To do that a formalism which takes into account the delay distribution function (DTD) of the explosion times and a given star formation history is adopted. Then the chemical evolution for ellipticals with baryonic initial masses $10^{10}$, $10^{11}$ and $10^{12} M_{\odot}$ is computed, and the mass of Fe produced by each galaxy is precisely estimated. We also compute the expected Fe mass ejected by ellipticals in typical galaxy clusters (e.g. Coma and Virgo), under different assumptions about Type Ia SN progenitors. As a last step, we compute the cosmic Type Ia SN rate in an unitary volume of the Universe by adopting several cosmic star formation rates and compare it with the available and recent observational data. Unfortunately, no firm conclusions can be derived only from the cosmic SNIa rate, neither on SNIa progenitors nor on the cosmic star formation rate. Finally, by analysing all our results together, and by taking into account previous chemical evolution results, we try to constrain the best Type Ia progenitor model. We conclude that the best progenitor models for Type Ia SNe are still the single degenerate model, the double degenerate wide model, and the empirical bimodal model. All these models require the existence of prompt Type Ia supernovae, exploding in the first 100 Myr since the beginning of star formation, although their fraction should not exceed 15-20% in order to fit chemical abundances in galaxies.

Galactic fountains and radial gas flows are very important ingredients in modeling the chemical evolution of galactic disks. Our aim here is to study the effects of galactic fountains and radial gas flows in the chemical evolution of the disk of M31. We adopt a ballistic method to study the effects of galactic fountains on the chemical enrichment of the M31 disk. We find that the landing coordinate for the fountains in M31 is no more than 1 kpc from the starting point, thus producing negligible effect on the chemical evolution of the disk. We find that the delay time in the enrichment process due to fountains is no longer than 100 Myr and this timescale also produces negligible effects on the results. Then, we compute the chemical evolution of the M31 disk with radial gas flows produced by the infall of extragalactic material and fountains. We find that a moderate inside-out formation of the disk coupled with radial flows of variable speed can very well reproduce the observed gradient. We discuss also the effects of other parameters such a threshold in the gas density for star formation and an efficiency of star formation varying with the galactic radius. We conclude that the most important physical processes in creating disk gradients are the inside-out formation and the radial gas flows. More data on abundance gradients both locally and at high redshift are necessary to confirm this conclusion.

Until recently the dynamical evolution of the interstellar medium (ISM) was simulated using collisional ionization equilibrium (CIE) conditions. However, the ISM is a dynamical system, in which the plasma is naturally driven out of equilibrium due to atomic and dynamic processes operating on different timescales. A step forward in the field comprises a multi-fluid approach taking into account the joint thermal and dynamical evolutions of the ISM gas.

The thick disk rotation--metallicity correlation, \partial V_\phi/\partial[Fe/H] =40\div 50 km s^{-1}dex^{-1} represents an important signature of the formation processes of the galactic disk. We use nondissipative numerical simulations to follow the evolution of a Milky Way (MW)-like disk to verify if secular dynamical processes can account for this correlation in the old thick disk stellar population. We followed the evolution of an ancient disk population represented by 10 million particles whose chemical abundances were assigned by assuming a cosmologically plausible radial metallicity gradient with lower metallicity in the inner regions, as expected for the 10-Gyr-old MW. Essentially, inner disk stars move towards the outer regions and populate layers located at higher |z|. A rotation--metallicity correlation appears, which well resembles the behaviour observed in our Galaxy at a galactocentric distance between 8 kpc and 10 kpc. In particular,we measure a correlation of \partial V_\phi/\partial[Fe/H]\simeq 60 km s^{-1}dex^{-1} for particles at 1.5 kpc < |z| < 2.0 kpc that persists up to 6 Gyr. Our pure N-body models can account for the V_\phi vs. [Fe/H] correlation observed in the thick disk of our Galaxy, suggesting that processes internal to the disk such as heating and radial migration play a role in the formation of this old stellar component. In this scenario, the positive rotation-metallicity correlation of the old thick disk population would represent the relic signature of an ancient "inverse" chemical (radial) gradient in the inner Galaxy, which resulted from accretion of primordial gas.

Using three-dimensional non-equilibrium ionization (NEI) hydrodynamical simulation of the interstellar medium (ISM), we study the electron density, $n_{e}$, in the Galactic disk and compare it with the values derived from dispersion measures towards pulsars with known distances located up to 200 pc on either side of the Galactic midplane. The simulation results, consistent with observations, can be summarized as follows: (i) the DMs in the simulated disk lie between the maximum and minimum observed values, (ii) the log <n_e> derived from lines of sight crossing the simulated disk follows a Gaussian distribution centered at \mu=-1.4 with a dispersion \sigma=0.21, thus, the Galactic midplane <n_e>=0.04\pm 0.01$ cm$^{-3}$, (iii) the highest electron concentration by mass (up to 80%) is in the thermally unstable regime (200<T<10^{3.9} K), (iv) the volume occupation fraction of the warm ionized medium is 4.9-6%, and (v) the electrons have a clumpy distribution along the lines of sight.

E(A+M)PEC traces the ionization structure, cooling and emission spectra of plasmas. It is written in OpenCL, runs in NVIDIA Graphics Processor Units and can be coupled to any HD or MHD code to follow the dynamical and thermal evolution of any plasma in, e.g., the interstellar medium (ISM).

We present the first to date three-dimensional high-resolution hydrodynamical simulation tracing the non-equilibrium ionization evolution (using the Eborae Atomic and Molecular Plasma Emission Code - E(A+M)PEC) of the Local Bubble and Loop I superbubbles embedded in a turbulent supernova-driven interstellar medium.

The majority of chemical evolution models assume that the Galactic disk forms by means of infall of gas and divide the disk into several independent rings without exchange of matter between them. However, if gas infall is important, radial gas flows should be taken into account as a dynamical consequence of infall. The aim of this paper is to test the effect of radial gas flows on detailed chemical evolution models (one-infall and two-infall) for the Milky Way disk with different prescriptions for the infall law and star formation rate. We found, that with a gas radial inflow of constant speed the metallicity gradient tends to steepen. Taking into account a constant time scale for the infall rate along the Galaxy disk and radial flows with a constant speed, we obtained a too flat gradient, at variance with data, implying that an inside-out formation and/or a variable gas flow speed are required. To reproduce the observed gradients the gas flow should increase in modulus with the galactocentric distance, both in the one-infall and two-infall models. However, the inside-out disk formation coupled with a threshold in the gas density (only in the two-infall model) for star formation and/or a variable efficiency of star formation with galactocentric distance can also reproduce the observed gradients without radial flows. We showed that the radial flows can be the most important process in reproducing abundance gradients but only with a variable gas speed. Finally, one should consider that uncertainties in the data concerning gradients prevent us to draw firm conclusions. Future more detailed data will help to ascertain whether the radial flows are a necessary ingredient in the formation and evolution of the Galactic disk and disks in general.

The existence of a mass-metallicity (MZ) relation in star forming galaxies at all redshift has been recently established. We aim at studying some possible physical mechanisms contributing to the MZ relation by adopting analytical solutions of chemical evolution models including infall and outflow. We explore the hypotheses of a variable galactic wind rate, infall rate and yield per stellar generation (i.e. a variation in the IMF), as possible causes for the MZ relation. By means of analytical models we compute the expected O abundance for galaxies of a given total baryonic mass and gas mass.The stellar mass is derived observationally and the gas mass is derived by inverting the Kennicutt law of star formation, once the star formation rate is known. Then we test how the parameters describing the outflow, infall and IMF should vary to reproduce the MZ relation, and we exclude the cases where such a variation leads to unrealistic situations. We find that a galactic wind rate increasing with decreasing galactic mass or a variable IMF are both viable solutions for the MZ relation. A variable infall rate instead is not acceptable. It is difficult to disentangle among the outflow and IMF solutions only by considering the MZ relation, and other observational constraints should be taken into account to select a specific solution. For example, a variable efficiency of star formation increasing with galactic mass can also reproduce the MZ relation and explain the downsizing in star formation suggested for ellipticals. The best solution could be a variable efficiency of star formation coupled with galactic winds, which are indeed observed in low mass galaxies.

The majority of galactic chemical evolution models assumes the instantaneous mixing approximation (IMA). This assumption is probably not realistic as indicated by the existence of chemical inhomogeneities, although current chemical evolution models of the Milky Way can reproduce the majority of the observational constraints under the IMA. The aim of this paper is to test whether relaxing this approximation in a detailed chemical evolution model can improve or worsen the agreement with observations. To do that, we investigated two possible causes for relaxing of the instantaneous mixing: i) the ``galactic fountain time delay effect'' and ii) the ``metal cooling time delay effect''. We found that the effect of galactic fountains is negligible if an average time delay of 0.1 Gyr, as suggested in a previous paper, is assumed. Longer time delays produce differences in the results but they are not realistic. We also found that the O abundance gradient in the disk is not affected by galactic fountains. The metal cooling time delays produce strong effects on the evolution of the chemical abundances only if we adopt stellar yields depending on metallicity. If instead, the yields computed for to the solar chemical composition are adopted, negligible effects are produced, as in the case of the galactic fountain delay. The relaxation of the IMA by means of the galactic fountain model, where the delay is considered only for massive stars and only in the disk, does not affect the chemical evolution results. The combination of metal dependent yields and time delay in the chemical enrichment from all stars starting from the halo phase, instead, produces results at variance with observations.

Our aim is to show how different hypotheses about Type Ia supernova progenitors can affect Galactic chemical evolution. We include different Type Ia SN progenitor models, identified by their distribution of time delays, in a very detailed chemical evolution model for the Milky Way which follows the evolution of several chemical species. We test the single degenerate and the double degenerate models for supernova Ia progenitors, as well as other more empirical models based on differences in the time delay distributions. We find that assuming the single degenerate or the double degenerate scenario produces negligible differences in the predicted [O/Fe] vs. [Fe/H] relation. On the other hand, assuming a percentage of prompt (exploding in the first 100 Myr) Type Ia supernovae of 50%, or that the maximum Type Ia rate is reached after 3-4 Gyr from the beginning of star formation, as suggested by several authors, produces more noticeable effects on the [O/Fe] trend. However, given the spread still existing in the observational data no model can be firmly excluded on the basis of only the [O/Fe] ratios. On the other hand, when the predictions of the different models are compared with the G-dwarf metallicity distribution, the scenarios with very few prompt Type Ia supernovae can be excluded. Models including the single degenerate or double degenerate scenario with a percentage of 10-13% of prompt Type Ia supernovae produce results in very good agreement with the observations. A fraction of prompt Type Ia supernovae larger than 30% worsens the agreement with observations and the same occurs if no prompt Type Ia supernovae are allowed. In particular, two empirical models for the Type Ia SN progenitors can be excluded: the one without prompt Type Ia supernovae and the one assuming delay time distribution going like t^{-0.5}.

(Abridged) The aim of this paper is to study the basic equations of the chemical evolution of galaxies with gas flows. We focus on models in which the outflow is differential, namely in which the heavy elements (or some of the heavy elements) can leave the parent galaxy more easily than other chemical species such as H and He. We study the chemical evolution of galaxies in the framework of simple models. This allows us to solve analytically the equations for the evolution of gas masses and metallicities. We find new analytical solutions for various cases in which the effects of winds and infall are taken into account. Differential galactic winds have the effect of reducing the global metallicity of a galaxy, with the amount of reduction increasing with the ejection efficiency of the metals. Abundance ratios are predicted to remain constant throughout the whole evolution of the galaxy, even in the presence of differential winds. One way to change them is by assuming differential winds with different ejection efficiencies for different elements. However, simple models apply only to elements produced on short timescales, namely all by Type II SNe, and therefore large differences in the ejection efficiencies of different metals are unlikely. Variations in abundance ratios such as [O/Fe] in galaxies, without including the Fe production by Type Ia supernovae, can in principle be obtained by assuming an unlikely different efficiency in the loss of O relative to Fe from Type II supernovae. Therefore, we conclude that it is not realistic to ignore Type Ia supernovae and that the delayed production of some chemical elements relative to others (time-delay model) remains the most plausible explanation for the evolution of alpha-elements relative to Fe.

The aim of this paper is to calculate the expansion law and chemical enrichment of a supershell powered by the energetic feedback of a typical Galactic OB association at various galactocentric radii. We study then the orbits of the fragments created when the supershell breaks out and we compare their kinetic and chemical properties with the available observations of high - and intermediate - velocity clouds. We use the Kompaneets (1960) approximation for the evolution of the superbubble driven by sequential supernova explosions and we compute the abundances of oxygen and iron residing in the thin cold supershell. We assume that supershells are fragmented by means of Rayleigh-Taylor instabilities and we follow the orbit of the clouds either ballistically or by means of a hybrid model considering viscous interaction between the clouds and the extra-planar gas.Given the self-similarity of the Kompaneets solutions, clouds are always formed ~ 448 pc above the plane. If the initial metallicity is solar, the pollution from dying stars of the OB association has a negligible effect on the chemical composition of the clouds. The maximum height reached by the clouds above the plane seldom exceeds 2 kpc and when averaging over different throwing angles, the landing coordinate differs from the throwing coordinate ~ 1 kpc at most. The range of heights and [O/Fe] ratios spun by our clouds suggest us that the high velocity clouds cannot have a Galactic origin, whereas intermediate velocity clouds have kinematic properties similar to our modeled clouds but overabundance observed for the [O/Fe] ratios which can be reproduced only with initial metallicities which are too low compared for those of the Galaxy disk.