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.

COMPARISONS OF THE INTERSTELLAR MAGNETIC FIELD DIRECTIONS OBTAINED FROM THE IBEX RIBBON AND INTERSTELLAR POLARIZATIONS

Variations in the spatial configuration of the interstellar magnetic field (ISMF) near the Sun can be constrained by comparing the ISMF direction at the heliosphere found from the Interstellar Boundary Explorer (IBEX) spacecraft observations of a ""Ribbon"" of energetic neutral atoms (ENAs), with the ISMF direction derived from optical polarization data for stars within similar to 40 pc. Using interstellar polarization observations toward similar to 30 nearby stars within similar to 90 degrees of the heliosphere nose, we find that the best fits to the polarization position angles are obtained for a magnetic pole directed toward ecliptic coordinates of lambda, beta similar to 263 degrees, 37 degrees (or galactic coordinates of l, b similar to 38 degrees, 23 degrees), with uncertainties of +/- 35 degrees based on the broad minimum of the best fits and the range of data quality. This magnetic pole is 33 degrees from the magnetic pole that is defined by the center of the arc of the ENA Ribbon. The IBEX ENA ribbon is seen in sight lines that are perpendicular to the ISMF as it drapes over the heliosphere. The similarity of the polarization and Ribbon directions for the local ISMF suggests that the local field is coherent over scale sizes of tens of parsecs. The ISMF vector direction is nearly perpendicular to the flow of local interstellar material (ISM) through the local standard of rest, supporting a possible local ISM origin related to an evolved expanding magnetized shell. The local ISMF direction is found to have a curious geometry with respect to the cosmic microwave background dipole moment.

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 bubbles or the boiling pot?: an ecosystemic approach to culture, environment and quality of life

For the diagnosis and prognosis of the problems of quality of life, a multidisciplinary ecosystemic approach encompasses four dimensions of being-in-the-world, as donors and recipients: intimate, interactive, social and biophysical. Social, cultural and environmental vulnerabilities are understood and dealt with, in different circumstances of space and time, as the conjugated effect of all dimensions of being-in-the-world, as they induce the events (deficits and assets), cope with consequences (desired or undesired) and contribute for change. Instead of fragmented and reduced representations of reality, diagnosis and prognosis of cultural, educational, environmental and health problems considers the connections (assets) and ruptures (deficits) between the different dimensions, providing a planning model to develop and evaluate research, teaching programmes, public policies and field projects. The methodology is participatory, experiential and reflexive; heuristic-hermeneutic processes unveil cultural and epistemic paradigms that orient subject-object relationships; giving people the opportunity to reflect on their own realities, engage in new experiences and find new ways to live better in a better world. The proposal is a creative model for thought and practice, providing many opportunities for discussion, debate and development of holistic projects integrating different scientific domains (social sciences, psychology, education, philosophy, etc.)

PRECESSING JETS AND X-RAY BUBBLES FROM NGC 1275 (3C 84) IN THE PERSEUS GALAXY CLUSTER: A VIEW FROM THREE-DIMENSIONAL NUMERICAL SIMULATIONS

The Perseus galaxy cluster is known to present multiple and misaligned pairs of cavities seen in X-rays, as well as twisted kiloparsec-scale jets at radio wavelengths; both morphologies suggest that the active galactic nucleus (AGN) jet is subject to precession. In this work, we performed three-dimensional hydrodynamical simulations of the interaction between a precessing AGN jet and the warm intracluster medium plasma, whose dynamics are coupled to a Navarro-Frenk-White dark matter gravitational potential. The AGN jet inflates cavities that become buoyantly unstable and rise up out of the cluster core. We found that under certain circumstances precession can originate multiple pairs of bubbles. For the physical conditions in the Perseus cluster, multiple pairs of bubbles are obtained for a jet precession opening angle >40 degrees acting for at least three precession periods, reproducing both radio and X-ray maps well. Based on such conditions, assuming that the Bardeen-Peterson effect is dominant, we studied the evolution of the precession opening angle of this system. We were able to constrain the ratio between the accretion disk and the black hole angular momenta as 0.7-1.4. We were also able to constrain the present precession angle to 30 degrees-40 degrees, as well as the approximate age of the inflated bubbles to 100-150 Myr.

The bubbles or the boiling pot? An ecosystemic approach to culture, environment and quality of life

For the diagnosis and prognosis of the problems of quality of life, a multidisciplinary ecosystemic approach encompasses four dimensions of being-in-the-world, as donors and recipients: intimate, interactive, social and biophysical. Social, cultural and environmental vulnerabilities are understood and dealt with, in different circumstances of space and time, as the conjugated effect of all dimensions of being-in-the-world, as they induce the events (deficits and assets), cope with consequences (desired or undesired) and contribute for change. Instead of fragmented and reduced representations of reality, diagnosis and prognosis of cultural, educational, environmental and health problems considers the connections (assets) and ruptures (deficits) between the different dimensions, providing a planning model to develop and evaluate research, teaching programmes, public policies and field projects. The methodology is participatory, experiential and reflexive; heuristic-hermeneutic processes unveil cultural and epistemic paradigms that orient subject-object relationships; giving people the opportunity to reflect on their own realities, engage in new experiences and find new ways to live better in a better world. The proposal is a creative model for thought and practice, providing many opportunities for discussion, debate and development of holistic projects integrating different scientific domains (social sciences, psychology, education, philosophy, etc.).

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.

H(2) infrared line emission from the ionized region of planetary nebulae

Context. The analysis and interpretation of the H(2) line emission from planetary nebulae have been done in the literature by assuming that the molecule survives only in regions where the hydrogen is neutral, as in photodissociation, neutral clumps, or shocked regions. However, there is strong observational and theoretical evidence that at least part of the H(2) emission is produced inside the ionized region of these objects. Aims. The aim of the present work is to calculate and analyze the infrared line emission of H(2) produced inside the ionized region of planetary nebulae using a one-dimensional photoionization code. Methods. The photoionization code Aangaba was improved in order to calculate the statistical population of the H(2) energy levels, as well as the intensity of the H(2) infrared emission lines in the physical conditions typical of planetary nebulae. A grid of models was obtained and the results then analyzed and compared with the observational data. Results. We show that the contribution of the ionized region to the H(2) line emission can be important, particularly in the case of nebulae with high-temperature central stars. This result explains why H(2) emission is more frequently observed in bipolar planetary nebulae (Gatley's rule), since this kind of object typically has hotter stars. Collisional excitation plays an important role in populating the rovibrational levels of the electronic ground state of H(2) molecules. Radiative mechanisms are also important, particularly for the upper vibrational levels. Formation pumping can have minor effects on the line intensities produced by de-excitation from very high rotational levels, especially in dense and dusty environments. We included the effect of the H(2) molecule on the thermal equilibrium of the gas, concluding that, in the ionized region, H(2) only contributes to the thermal equilibrium in the case of a very high temperature of the central star or a high dust-to-gas ratio, mainly through collisional de-excitation.

First stars XIV. Sulfur abundances in extremely metal-poor stars

Context. Precise S abundances are important in the study of the early chemical evolution of the Galaxy. In particular the site of the formation remains uncertain because, at low metallicity, the trend of this alpha-element versus [Fe/H] remains unclear. Moreover, although sulfur is not bound significantly in dust grains in the ISM, it seems to behave differently in DLAs and old metal-poor stars. Aims. We attempt a precise measurement of the S abundance in a sample of extremely metal-poor stars observed with the ESO VLT equipped with UVES, taking into account NLTE and 3D effects. Methods. The NLTE profiles of the lines of multiplet 1 of S I were computed with a version of the program MULTI, including opacity sources from ATLAS9 and based on a new model atom for S. These profiles were fitted to the observed spectra. Results. We find that sulfur in EMP stars behaves like the other alpha-elements, with [S/Fe] remaining approximately constant below [Fe/H] = -3. However, [S/Mg] seems to decrease slightly with increasing [Mg/H]. The overall abundance patterns of O, Na, Mg, Al, S, and K are most closely matched by the SN model yields by Heger & Woosley. The [S/Zn] ratio in EMP stars is solar, as also found in DLAs. We derive an upper limit to the sulfur abundance [S/Fe] < +0.5 for the ultra metal-poor star CS 22949-037. This, along with a previously reported measurement of zinc, argues against the conjecture that the light-element abundance pattern of this star (and by analogy, the hyper iron-poor stars HE 0107-5240 and HE 1327-2326) would be due to dust depletion.

First stars XIII. Two extremely metal-poor RR Lyrae stars

Context. The chemical composition of extremely metal-poor stars (EMP stars; [Fe/H] < similar to -3) is a unique tracer of early nucleosynthesis in the Galaxy. As such stars are rare, we wish to find classes of luminous stars which can be studied at high spectral resolution. Aims. We aim to determine the detailed chemical composition of the two EMP stars CS 30317-056 and CS 22881-039, originally thought to be red horizontal-branch (RHB) stars, and compare it to earlier results for EMP stars as well as to nucleosynthesis yields from various supernova (SN) models. In the analysis, we discovered that our targets are in fact the two most metal-poor RR Lyrae stars known. Methods. Our detailed abundance analysis, taking into account the variability of the stars, is based on VLT/UVES spectra (R similar or equal to 43 000) and 1D LTE OSMARCS model atmospheres and synthetic spectra. For comparison with SN models we also estimate NLTE corrections for a number of elements. Results. We derive LTE abundances for the 16 elements O, Na, Mg, Al, Si, S, Ca, Sc, Ti, Cr, Mn, Fe, Co, Ni, Sr and Ba, in good agreement with earlier values for EMP dwarf, giant and RHB stars. Li and C are not detected in either star. NLTE abundance corrections are newly calculated for O and Mg and taken from the literature for other elements. The resulting abundance pattern is best matched by model yields for supernova explosions with high energy and/or significant asphericity effects. Conclusions. Our results indicate that, except for Li and C, the surface composition of EMP RR Lyr stars is not significantly affected by mass loss, mixing or diffusion processes; hence, EMP RR Lyr stars should also be useful tracers of the chemical evolution of the early Galactic halo. The observed abundance ratios indicate that these stars were born from an ISM polluted by energetic, massive (25-40 M(circle dot)) and/or aspherical supernovae, but the NLTE corrections for Sc and certain other elements do play a role in the choice of model.

The chemical evolution of IC 10

Context. Dwarf irregular galaxies are relatively simple unevolved objects where it is easy to test models of galactic chemical evolution. Aims. We attempt to determine the star formation and gas accretion history of IC 10, a local dwarf irregular for which abundance, gas, and mass determinations are available. Methods. We apply detailed chemical evolution models to predict the evolution of several chemical elements (He, O, N, S) and compared our predictions with the observational data. We consider additional constraints such as the present-time gas fraction, the star formation rate (SFR), and the total estimated mass of IC 10. We assume a dark matter halo for this galaxy and study the development of a galactic wind. We consider different star formation regimes: bursting and continuous. We explore different wind situations: i) normal wind, where all the gas is lost at the same rate and ii) metal-enhanced wind, where metals produced by supernovae are preferentially lost. We study a case without wind. We vary the star formation efficiency (SFE), the wind efficiency, and the time scale of the gas infall, which are the most important parameters in our models. Results. We find that only models with metal-enhanced galactic winds can reproduce the properties of IC 10. The star formation must have proceeded in bursts rather than continuously and the bursts must have been less numerous than similar to 10 over the whole galactic lifetime. Finally, IC 10 must have formed by a slow process of gas accretion with a timescale of the order of 8 Gyr.

TURBULENCE AND THE FORMATION OF FILAMENTS, LOOPS, AND SHOCK FRONTS IN NGC 1275

NGC 1275, the central galaxy in the Perseus cluster, is the host of gigantic hot bipolar bubbles inflated by active galactic nucleus (AGN) jets observed in the radio as Perseus A. It presents a spectacular H alpha-emitting nebulosity surrounding NGC 1275, with loops and filaments of gas extending to over 50 kpc. The origin of the filaments is still unknown, but probably correlates with the mechanism responsible for the giant buoyant bubbles. We present 2.5 and three-dimensional magnetohydrodynamical (MHD) simulations of the central region of the cluster in which turbulent energy, possibly triggered by star formation and supernovae (SNe) explosions, is introduced. The simulations reveal that the turbulence injected by massive stars could be responsible for the nearly isotropic distribution of filaments and loops that drag magnetic fields upward as indicated by recent observations. Weak shell-like shock fronts propagating into the intracluster medium (ICM) with velocities of 100-500 km s(-1) are found, also resembling the observations. The isotropic outflow momentum of the turbulence slows the infall of the ICM, thus limiting further starburst activity in NGC 1275. As the turbulence is subsonic over most of the simulated volume, the turbulent kinetic energy is not efficiently converted into heat and additional heating is required to suppress the cooling flow at the core of the cluster. Simulations combining the MHD turbulence with the AGN outflow can reproduce the temperature radial profile observed around NGC 1275. While the AGN mechanism is the main heating source, the SNe are crucial to isotropize the energy distribution.

Star formation history of Canis Major R1 I. Wide-Field X-ray study of the young stellar population

Aims. The CMa R1 star-forming region contains several compact clusters as well as many young early-B stars. It is associated with a well-known bright rimmed nebula, the nature of which is unclear (fossil HII region or supernova remnant). To help elucidate the nature of the nebula, our goal was to reconstruct the star-formation history of the CMa R1 region, including the previously unknown older, fainter low-mass stellar population, using X-rays. Methods. We analyzed images obtained with the ROSAT satellite, covering similar to 5 sq. deg. Complementary VRI photometry was performed with the Gemini South telescope. Colour-magnitude and colour-colour diagrams were used in conjunction with pre-main sequence evolutionary tracks to derive the masses and ages of the X-ray sources. Results. The ROSAT images show two distinct clusters. One is associated with the known optical clusters near Z CMa, to which similar to 40 members are added. The other, which we name the ""GU CMa"" cluster, is new, and contains similar to 60 members. The ROSAT sources are young stars with masses down to M(star) similar to 0.5 M(circle dot), and ages up to 10 Myr. The mass functions of the two clusters are similar, but the GU CMa cluster is older than the cluster around Z CMa by at least a few Myr. Also, the GU CMa cluster is away from any molecular cloud, implying that star formation must have ceased; on the contrary (as already known), star formation is very active in the Z CMa region.

Ultra-high-energy cosmic rays from centaurus A: jet interaction with gaseous shells

Ultra-high-energy cosmic rays (UHECRs), with energies above similar to 6 x 10(19) eV, seem to show a weak correlation with the distribution of matter relatively near to us in the universe. It has earlier been proposed that UHECRs could be accelerated in either the nucleus or the outer lobes of the nearby radio galaxy Cen A. We show that UHECR production at a spatially intermediate location about 15 kpc northeast from the nucleus, where the jet emerging from the nucleus is observed to strike a large star-forming shell of gas, is a plausible alternative. A relativistic jet is capable of accelerating lower energy heavy seed cosmic rays (CRs) to UHECRs on timescales comparable to the time it takes the jet to pierce the large gaseous cloud. In this model, many CRs arising from a starburst, with a composition enhanced in heavy elements near the knee region around PeV, are boosted to ultra-high energies by the relativistic shock of a newly oriented jet. This model matches the overall spectrum shown by the Auger data and also makes a prediction for the chemical composition as a function of particle energy. We thus predict an observable anisotropy in the composition at high energy in the sense that lighter nuclei should preferentially be seen toward the general direction of Cen A. Taking into consideration the magnetic field models for the Galactic disk and a Galactic magnetic wind, this scenario may resolve the discrepancy between HiRes and Auger results concerning the chemical composition of UHECRs.

EVOLUTION AND LIFETIME OF TRANSIENT CLUMPS IN THE TURBULENT INTERSTELLAR MEDIUM

We study the evolution of dense clumps and provide an argument that the existence of the clumps is not limited by their crossing times. We claim that the lifetimes of the clumps are determined by turbulent motions on a larger scale, and we predict the correlation of clump lifetime with column density. We use numerical simulations to successfully test this relation. In addition, we study the morphological asymmetry and the magnetization of the clumps as functions of their masses.