Galactic fountain flows in high-resolution hydrodynamical simulations of galaxy disks

Feedback from the most massive components of a young stellar cluster deeply affects the surrounding ISM driving an expanding over-pressured hot gas cavity in it. In spiral galaxies these structures may have sufficient energy to break the disk and eject large amount of material into the halo. The cycling of this gas, which eventually will fall back onto the disk, is known as galactic fountains. We aim at better understanding the dynamics of such fountain flow in a Galactic context, frame the problem in a more dynamic environment possibly learning about its connection and regulation to the local driving mechanism and understand its role as a metal diffusion channel. The interaction of the fountain with a hot corona is hereby analyzed, trying to understand the properties and evolution of the extraplanar material. We perform high resolution hydrodynamical simulations with the moving-mesh code AREPO to model the multi-phase ISM of a Milky Way type galaxy. A non-equilibrium chemical network is included to self consistently follow the evolution of the main coolants of the ISM. Spiral arm perturbations in the potential are considered so that large molecular gas structures are able to dynamically form here, self shielded from the interstellar radiation field. We model the effect of SN feedback from a new-born stellar cluster inside such a giant molecular cloud, as the driving force of the fountain. Passive Lagrangian tracer particles are used in conjunction to the SN energy deposition to model and study diffusion of freshly synthesized metals. We find that both interactions with hot coronal gas and local ISM properties and motions are equally important in shaping the fountain. We notice a bimodal morphology where most of the ejected gas is in a cold $10^4$ K clumpy state while the majority of the affected volume is occupied by a hot diffuse medium. While only about 20\% of the produced metals stay local, most of them quickly diffuse through this hot regime to great scales.

Analytic modelling of the galactic fountain in the Milky Way potential

Le galassie a spirale, come la Via Lattea, sono caratterizzate dalla presenza di gas freddo e formazione stellare e vengono perciò chiamate star-forming. Per creare nuove stelle è necessaria una sufficiente riserva di gas, la cui disponibilità governa l’evoluzione della galassia stessa. Finora, non è stato individuato con certezza un meccanismo che possa alimentare la formazione di nuove stelle nelle galassie star-forming. Una delle possibili sorgenti di tale gas è l’alone galattico caldo (corona galattica) il cui raffreddamento e successivo accrescimento possono essere stimolati dal processo di fontana galattica. L’esplosione di supernovae porta nubi di gas freddo in orbita al di sopra del disco stellare; queste nubi raggiungono altezze dell’ordine del kiloparsec, interagendo con la corona di gas caldo. Il moto delle nubi all’interno di un mezzo meno denso comporta l’instaurarsi dell’instabilità di Kelvin-Helmholtz, che ’strappa’ gas dalle nubi e causa la condensazione di materia coronale. Quest’ultima viene quindi accresciuta e, ricadendo sul disco, trasferisce nuovo materiale alla galassia e ne alimenta la formazione stellare. Lo scopo di questa tesi è derivare un modello analitico di fontana galattica che consenta di ottenere una formulazione analitica per il tempo orbitale, cioè il tempo richiesto alle nubi per ricadere sul disco galattico. Infatti, più u tempo le nubi impiegano per attraversare il materiale coronale caldo e ricadere sul disco, più materiale viene accresciuto durante l’orbita. Conoscendo i tempi orbitali sarebbe possibile calcolare il tasso di accrescimento legato al fenomeno di fontana e studiarne l’andamento con il raggio del disco. Questo modello potrebbe rivelarsi utile per lo studio dell’impatto della fontana nell’evoluzione globale del disco galattico e della chimica dell’intera galassia.

Galactic outflows and the pollution of the galactic environment by supernovae

We here explore the effects of the SN explosions into the halo of star-forming galaxies like the Milky Way. Successive randomly distributed and clustered SNe explosions cause the formation of hot superbubbles that drive either fountains or galactic winds above the galactic disk, depending on the amount and concentration of energy that is injected by the SNe. 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 up to about 5 kpc in the halo and thus allow for the formation of the so called intermediate-velocity-clouds (IVCs) which are often observed in the halos of disk galaxies. The high-velocity-clouds that are also observed but 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 to the halo with the gas in the fountains. Due to angular momentum losses to the halo, we find that the fountain material falls back to smaller radii and is not largely spread over the galactic disk. Instead, the SNe ejecta fall nearby the region where the fountain was produced, a result which is consistent with recent chemical models of the galaxy. The fall back material leads to the formation of new generations of molecular clouds and to supersonic turbulence feedback in the disk. (C) 2009 COSPAR. Published by Elsevier Ltd. All rights reserved.

Hydrodynamical simulations of Galactic fountains - II. Evolution of multiple fountains

The ejection of gas out of the disc in late-type galaxies is related to star formation and is mainly due to the explosion of Type II supernovae (SN II). In a previous paper, we considered the evolution of a single Galactic fountain, that is, a fountain powered by a single SN cluster. Using three-dimensional 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 centre and the presence of a hot gaseous halo. As a natural followup, this paper investigates the dynamical evolution of multiple generations of fountains generated by similar to 100 OB associations. 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. Most of the results of the previous paper have been confirmed, like for example the formation of intermediate velocity clouds above the disc by the multiple fountains. Also, this 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 disc chemical abundance. The multiple fountain simulations also allowed us to consistently calculate the feedback of the star formation on the halo gas. We found that the hot gas gains about 10 per cent of all the SN II energy produced in the disc. Thus, the SN feedback more than compensate for the halo radiative losses and allow for a quasi steady-state disc-halo circulation to exist. 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 haloes around the disc 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.

The gaseous halo of the Milky Way

In the last decade, sensitive observations have revealed that disc galaxies are surrounded by multiphase gaseous halos produced by the circulation of gas from the discs to the environment and vice-versa. This Thesis is a study of the gaseous halo of the Milky Way carried out via the modelling of the HI emission and the available absorption-line data. We fitted simple kinematical models to the HI LAB Survey and found that the Galaxy has a massive (~3x10^8 Mo) HI halo extending a few kiloparsecs above the plane. This layer rotates more slowly than the disc and shows a global inflow motion, a kinematics similar to that observed in the HI halos of nearby galaxies. We built a dynamical model of the galactic fountain to reproduce the properties of this layer. In this model, fountain clouds are ejected from the disc by SN feedback and - as suggested by hydrodynamical simulations - triggers the cooling of coronal gas, which is entrained by the cloud wakes and accretes onto the disc when the clouds fall back. For a proper choice of the parameters, the model reproduces well the HI data and predicts an accretion of coronal gas onto the disc at a rate of 2 Mo/yr. We extended this model to the warm-hot component of the halo, showing that most of the ion absorption features observed towards background sources are consistent with being produced in the turbulent wakes that lag behind the fountain clouds. Specifically, the column densities, positions, and velocities of the absorbers are well reproduced by our model. Finally, we studied the gas content of galaxies extracted from a cosmological N-body+SPH simulation, and found that an HI halo with the forementioned properties is not observed, probably due ti the relatively low resolution of the simulations.

Hydrodynamical simulations of Galactic fountains - I. Evolution of single fountains

The ejection of the gas out of the disc 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 on the redistribution of the freshly delivered metals over the disc. 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 supernova explosions generated by a single OB association are considered. A typical fountain powered by 100 Type II supernovae may eject material up to similar to 2 kpc which than collapses back mostly in the form of dense, cold clouds and filaments. The majority of the gas lifted up by the fountains falls back on the disc remaining within a radial distance Delta R = 0.5 kpc from the place where the fountain originated. This localized circulation of disc 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.

A SEARCH FOR FIELD HORIZONTAL-BRANCH STARS NEAR THE SOUTH GALACTIC POLE

Searches for field horizontal-branch (FHB) stars in the halo of the Galaxy in the past have been carried out by several techniques, such as objective-prism surveys and visual or infrared photometric surveys. By choosing adequate color criteria, it is possible to improve the efficiency of identifying bona fide FHB stars among the other objects that exhibit similar characteristics, such as main-sequence A-stars, blue stragglers, subdwarfs, etc. In this work, we report the results of a spectroscopic survey carried out near the south Galactic pole intended to validate FHB stars originally selected from the HK objective-prism survey of Beers and colleagues, based on near-infrared color indices. A comparison between the stellar spectra obtained in this survey with theoretical stellar atmosphere models allows us to determine T(eff), log g, and [Fe/H] for 13 stars in the sample. Stellar temperatures were calculated from measured (B-V)(o), when this measurement was available (16 stars). The color index criteria adopted in this work are shown to correctly classify 30% of the sample as FHB, 25% as non-FHB (main-sequence stars and subdwarfes), whereas 40% could not be distinguished between FHB and main-sequence stars. We compare the efficacy of different color criteria in the literature intended to select FHB stars, and discuss the use of the Mg II 4481 line to estimate the metallicity.

Abundances in the Galactic bulge: results from planetary nebulae and giant stars

Context. Our understanding of the chemical evolution (CE) of the Galactic bulge requires the determination of abundances in large samples of giant stars and planetary nebulae (PNe). Studies based on high resolution spectroscopy of giant stars in several fields of the Galactic bulge obtained with very large telescopes have allowed important progress. Aims. We discuss PNe abundances in the Galactic bulge and compare these results with those presented in the literature for giant stars. Methods. We present the largest, high-quality data-set available for PNe in the direction of the Galactic bulge (inner-disk/bulge). For comparison purposes, we also consider a sample of PNe in the Large Magellanic Cloud (LMC). We derive the element abundances in a consistent way for all the PNe studied. By comparing the abundances for the bulge, inner-disk, and LMC, we identify elements that have not been modified during the evolution of the PN progenitor and can be used to trace the bulge chemical enrichment history. We then compare the PN abundances with abundances of bulge field giant. Results. At the metallicity of the bulge, we find that the abundances of O and Ne are close to the values for the interstellar medium at the time of the PN progenitor formation, and hence these elements can be used as tracers of the bulge CE, in the same way as S and Ar, which are not expected to be affected by nucleosynthetic processes during the evolution of the PN progenitors. The PN oxygen abundance distribution is shifted to lower values by 0.3 dex with respect to the distribution given by giants. A similar shift appears to occur for Ne and S. We discuss possible reasons for this PNe-giant discrepancy and conclude that this is probably due to systematic errors in the abundance derivations in either giants or PNe (or both). We issue an important warning concerning the use of absolute abundances in CE studies.

NGC 7097: THE ACTIVE GALACTIC NUCLEUS AND ITS MIRROR, REVEALED BY PRINCIPAL COMPONENT ANALYSIS TOMOGRAPHY

Three-dimensional spectroscopy techniques are becoming more and more popular, producing an increasing number of large data cubes. The challenge of extracting information from these cubes requires the development of new techniques for data processing and analysis. We apply the recently developed technique of principal component analysis (PCA) tomography to a data cube from the center of the elliptical galaxy NGC 7097 and show that this technique is effective in decomposing the data into physically interpretable information. We find that the first five principal components of our data are associated with distinct physical characteristics. In particular, we detect a low-ionization nuclear-emitting region (LINER) with a weak broad component in the Balmer lines. Two images of the LINER are present in our data, one seen through a disk of gas and dust, and the other after scattering by free electrons and/or dust particles in the ionization cone. Furthermore, we extract the spectrum of the LINER, decontaminated from stellar and extended nebular emission, using only the technique of PCA tomography. We anticipate that the scattered image has polarized light due to its scattered nature.

The stellar content of obscured Galactic giant H II regions. VI. W51A

We present K-band spectra of newly born OB stars in the obscured Galactic giant H II region W51A and approximate to 0.8 '' angular resolution images in the J, H, and K(S)-bands. Four objects have been spectroscopically classified as O-type stars. The mean spectroscopic parallax of the four stars gives a distance of 2.0 +/- 0.3 kpc (error in the mean), significantly smaller than the radio recombination line kinematic value of 5.5 kpc or the values derived from maser proper motion observations (6-8 kpc). The number of Lyman continuum photons from the contribution of all massive stars (NLyc approximate to 1.5 x 10(50) s(-1)) is in good agreement with that inferred from radio recombination lines (NLyc = 1.3 x 10(50) s(-1)) after accounting for the smaller distance derived here. We present analysis of archival high angular resolution images (NAOS CONICA at VLT and T-ReCS at Gemini) of the compact region W51 IRS 2. The K(S)-band images resolve the infrared source IRS 2 indicating that it is a very young compact H II region. Sources IRS 2E was resolved into compact cluster (within 660 AU of projected distance) of three objects, but one of them is just bright extended emission. W51d1 and W51d2 were identified with compact clusters of three objects (maybe four in the case of W51d1) each one. Although IRS 2E is the brightest source in the K-band and at 12.6 mu m, it is not clearly associated with a radio continuum source. Our spectrum of IRS 2E shows, similar to previous work, strong emission in Br gamma and He I, as well as three forbidden emission lines of Fe III and emission lines of molecular hydrogen (H(2)) marking it as a massive young stellar object.

Chemical similarities between Galactic bulge and local thick disk red giant stars

Context. The evolution of the Milky Way bulge and its relationship with the other Galactic populations is still poorly understood. The bulge has been suggested to be either a merger-driven classical bulge or the product of a dynamical instability of the inner disk. Aims. To probe the star formation history, the initial mass function and stellar nucleosynthesis of the bulge, we performed an elemental abundance analysis of bulge red giant stars. We also completed an identical study of local thin disk, thick disk and halo giants to establish the chemical differences and similarities between the various populations. Methods. High-resolution infrared spectra of 19 bulge giants and 49 comparison giants in the solar neighborhood were acquired with Gemini/Phoenix. All stars have similar stellar parameters but cover a broad range in metallicity. A standard 1D local thermodynamic equilibrium analysis yielded the abundances of C, N, O and Fe. A homogeneous and differential analysis of the bulge, halo, thin disk and thick disk stars ensured that systematic errors were minimized. Results. We confirm the well-established differences for [O/Fe] (at a given metallicity) between the local thin and thick disks. For the elements investigated, we find no chemical distinction between the bulge and the local thick disk, which is in contrast to previous studies relying on literature values for disk dwarf stars in the solar neighborhood. Conclusions. Our findings suggest that the bulge and local thick disk experienced similar, but not necessarily shared, chemical evolution histories. We argue that their formation timescales, star formation rates and initial mass functions were similar.

New constraints on the chemical evolution of the solar neighbourhood and Galactic disc(s) Improved astrophysical parameters for the Geneva-Copenhagen Survey

We present a re-analysis of the Geneva-Copenhagen survey, which benefits from the infrared flux method to improve the accuracy of the derived stellar effective temperatures and uses the latter to build a consistent and improved metallicity scale. Metallicities are calibrated on high-resolution spectroscopy and checked against four open clusters and a moving group, showing excellent consistency. The new temperature and metallicity scales provide a better match to theoretical isochrones, which are used for a Bayesian analysis of stellar ages. With respect to previous analyses, our stars are on average 100 K hotter and 0.1 dex more metal rich, which shift the peak of the metallicity distribution function around the solar value. From Stromgren photometry we are able to derive for the first time a proxy for [alpha/Fe] abundances, which enables us to perform a tentative dissection of the chemical thin and thick disc. We find evidence for the latter being composed of an old, mildly but systematically alpha-enhanced population that extends to super solar metallicities, in agreement with spectroscopic studies. Our revision offers the largest existing kinematically unbiased sample of the solar neighbourhood that contains full information on kinematics, metallicities, and ages and thus provides better constraints on the physical processes relevant in the build-up of the Milky Way disc, enabling a better understanding of the Sun in a Galactic context.

Chemical evolution of the Galactic bulge as traced by microlensed dwarf and subgiant stars IV. Two bulge populations

Based on high-resolution (R approximate to 42 000 to 48 000) and high signal-to-noise (S/N approximate to 50 to 150) spectra obtained with UVES/VLT, we present detailed elemental abundances (O, Na, Mg, Al, Si, Ca, Ti, Cr, Fe, Ni, Zn, Y, and Ba) and stellar ages for 12 new microlensed dwarf and subgiant stars in the Galactic bulge. Including previous microlensing events, the sample of homogeneously analysed bulge dwarfs has now grown to 26. The analysis is based on equivalent width measurements and standard 1-D LTE MARCS model stellar atmospheres. We also present NLTE Li abundances based on line synthesis of the (7)Li line at 670.8 nm. The results from the 26 microlensed dwarf and subgiant stars show that the bulge metallicity distribution (MDF) is double-peaked; one peak at [Fe/H] approximate to -0.6 and one at [Fe/H] approximate to +0.3, and with a dearth of stars around solar metallicity. This is in contrast to the MDF derived from red giants in Baade's window, which peaks at this exact value. A simple significance test shows that it is extremely unlikely to have such a gap in the microlensed dwarf star MDF if the dwarf stars are drawn from the giant star MDF. To resolve this issue we discuss several possibilities, but we can not settle on a conclusive solution for the observed differences. We further find that the metal-poor bulge dwarf stars are predominantly old with ages greater than 10 Gyr, while the metal-rich bulge dwarf stars show a wide range of ages. The metal-poor bulge sample is very similar to the Galactic thick disk in terms of average metallicity, elemental abundance trends, and stellar ages. Speculatively, the metal-rich bulge population might be the manifestation of the inner thin disk. If so, the two bulge populations could support the recent findings, based on kinematics, that there are no signatures of a classical bulge and that the Milky Way is a pure-disk galaxy. Also, recent claims of a flat IMF in the bulge based on the MDF of giant stars may have to be revised based on the MDF and abundance trends probed by our microlensed dwarf stars.

Discovery of VVV CL001 A low-mass globular cluster next to UKS 1 in the direction of the Galactic bulge

Context. It is not known how many globular clusters may remain undetected towards the Galactic bulge. Aims. One of the aims of the VISTA Variables in the Via Lactea (VVV) Survey is to accurately measure the physical parameters of the known globular clusters in the inner regions of the Milky Way and search for new ones, hidden in regions of large extinction. Methods. From deep near-infrared images, we derive deep JHK(S)-band photometry of a region surrounding the known globular cluster UKS 1 and reveal a new low-mass globular cluster candidate that we name VVV CL001. Results. We use the horizontal-branch red clump to measure E(B-V) similar to 2.2 mag, (m - M)(0) = 16.01 mag, and D = 15.9 kpc for the globular cluster UKS 1. On the basis of near-infrared colour-magnitude diagrams, we also find that VVV CL001 has E(B-V) similar to 2.0, and that it is at least as metal-poor as UKS 1, although its distance remains uncertain. Conclusions. Our finding confirms the previous projection that the central region of the Milky Way harbours more globular clusters. VVV CL001 and UKS 1 are good candidates for a physical cluster binary, but follow-up observations are needed to decide if they are located at the same distance and have similar radial velocities.

Probing the Galactic thick disc vertical properties and interfaces

Aims. This work investigates the properties (metallicity and kinematics) and interfaces of the Galactic thick disc as a function of height above the Galactic plane. The main aim is to study the thick disc in a place where it is the main component of the sample. Methods. We take advantage of former astrometric work in two fields of several square degrees in which accurate proper motions were measured down to V-magnitudes of 18.5 in two directions, one near the north galactic pole and the other at a galactic latitude of 46 degrees and galactic longitude near 0 degrees. Spectroscopic observations have been acquired in these two fields for a total of about 400 stars down to magnitude 18.0, at spectral resolutions of 3.5 to 6.25 angstrom. The spectra have been analysed with the code ETOILE, comparing the target stellar spectra with a grid of 1400 reference stellar spectra. This comparison allowed us to derive the parameters effective temperature, gravity, [Fe/H] and absolute magnitude for each target star. Results. The Metallicity Distribution Function (MDF) of the thin-thick-disc-halo system is derived for several height intervals between 0 and 5 kpc above the Galactic plane. The MDFs show a decrease of the ratio of the thin to thick disc stars between the first and second kilo-parsec. This is consistent with the classical modelling of the vertical density profile of the disc with 2 populations with different scale heights. A vertical metallicity gradient, partial derivative[Fe/H]/partial derivative z = -0.068 +/- 0.009 dex kpc(-1), is observed in the thick disc. It is discussed in terms of scenarios of formation of the thick disc.