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.

Cross-entropy optimizer: a new tool to study precession in astrophysical jets

Evidence of jet precession in many galactic and extragalactic sources has been reported in the literature. Much of this evidence is based on studies of the kinematics of the jet knots, which depends on the correct identification of the components to determine their respective proper motions and position angles on the plane of the sky. Identification problems related to fitting procedures, as well as observations poorly sampled in time, may influence the follow-up of the components in time, which consequently might contribute to a misinterpretation of the data. In order to deal with these limitations, we introduce a very powerful statistical tool to analyse jet precession: the cross-entropy method for continuous multi-extremal optimization. Only based on the raw data of the jet components (right ascension and declination offsets from the core), the cross-entropy method searches for the precession model parameters that better represent the data. In this work we present a large number of tests to validate this technique, using synthetic precessing jets built from a given set of precession parameters. With the aim of recovering these parameters, we applied the cross-entropy method to our precession model, varying exhaustively the quantities associated with the method. Our results have shown that even in the most challenging tests, the cross-entropy method was able to find the correct parameters within a 1 per cent level. Even for a non-precessing jet, our optimization method could point out successfully the lack of precession.

PLANETARY NEBULAE AND THE CHEMICAL EVOLUTION OF THE MAGELLANIC CLOUDS

The determination of accurate chemical abundances of planetary nebulae (PN) in different galaxies allows us to obtain important constraints on chemical evolution models for these systems. We have a long-term program to derive abundances in the galaxies of the Local Group, particularly the Large and Small Magellanic Clouds. In this work, we present our new results on these objects and discuss their implications in view of recent abundance determinations in the literature. In particular, we obtain distance-independent correlations involving He, N, O, Ne, S, and Ar, and compare the results with data from our own Galaxy and other galaxies in the Local Group. As a result of our observational program, we have a large database of PN in the Galaxy and the Magellanic Clouds, so that we can obtain reliable constraints on the nucleosynthesis processes in the progenitor stars in galaxies of different metallicities.

MODELING VERY LONG BASELINE INTERFEROMETRIC IMAGES WITH THE CROSS-ENTROPY GLOBAL OPTIMIZATION TECHNIQUE

We present a new technique for obtaining model fittings to very long baseline interferometric images of astrophysical jets. The method minimizes a performance function proportional to the sum of the squared difference between the model and observed images. The model image is constructed by summing N(s) elliptical Gaussian sources characterized by six parameters: two-dimensional peak position, peak intensity, eccentricity, amplitude, and orientation angle of the major axis. We present results for the fitting of two main benchmark jets: the first constructed from three individual Gaussian sources, the second formed by five Gaussian sources. Both jets were analyzed by our cross-entropy technique in finite and infinite signal-to-noise regimes, the background noise chosen to mimic that found in interferometric radio maps. Those images were constructed to simulate most of the conditions encountered in interferometric images of active galactic nuclei. We show that the cross-entropy technique is capable of recovering the parameters of the sources with a similar accuracy to that obtained from the very traditional Astronomical Image Processing System Package task IMFIT when the image is relatively simple (e. g., few components). For more complex interferometric maps, our method displays superior performance in recovering the parameters of the jet components. Our methodology is also able to show quantitatively the number of individual components present in an image. An additional application of the cross-entropy technique to a real image of a BL Lac object is shown and discussed. Our results indicate that our cross-entropy model-fitting technique must be used in situations involving the analysis of complex emission regions having more than three sources, even though it is substantially slower than current model-fitting tasks (at least 10,000 times slower for a single processor, depending on the number of sources to be optimized). As in the case of any model fitting performed in the image plane, caution is required in analyzing images constructed from a poorly sampled (u, v) plane.

Global properties of `ordinary` early-type galaxies: photometry and spectroscopy of stars and globular clusters in NGC 4494

We present a comprehensive analysis of the spatial, kinematic and chemical properties of stars and globular clusters (GCs) in the `ordinary` elliptical galaxy NGC 4494 using data from the Keck and Subaru telescopes. We derive galaxy surface brightness and colour profiles out to large galactocentric radii. We compare the latter to metallicities derived using the near-infrared Calcium Triplet. We obtain stellar kinematics out to similar to 3.5 effective radii. The latter appear flattened or elongated beyond similar to 1.8 effective radii in contrast to the relatively round photometric isophotes. In fact, NGC 4494 may be a flattened galaxy, possibly even an S0, seen at an inclination of similar to 45 degrees. We publish a catalogue of 431 GC candidates brighter than i(0) = 24 based on the photometry, of which 109 are confirmed spectroscopically and 54 have measured spectroscopic metallicities. We also report the discovery of three spectroscopically confirmed ultra-compact dwarfs around NGC 4494 with measured metallicities of -0.4 less than or similar to [Fe/H] less than or similar to -0.3. Based on their properties, we conclude that they are simply bright GCs. The metal-poor GCs are found to be rotating with similar amplitude as the galaxy stars, while the metal-rich GCs show marginal rotation. We supplement our analysis with available literature data and results. Using model predictions of galaxy formation, and a suite of merger simulations, we find that many of the observational properties of NGC 4494 may be explained by formation in a relatively recent gas-rich major merger. Complete studies of individual galaxies incorporating a range of observational avenues and methods such as the one presented here will be an invaluable tool for constraining the fine details of galaxy formation models, especially at large galactocentric radii.

Probing the 2D kinematic structure of early-type galaxies out to three effective radii

We detail an innovative new technique for measuring the two-dimensional (2D) velocity moments (rotation velocity, velocity dispersion and Gauss-Hermite coefficients h(3) and h(4)) of the stellar populations of galaxy haloes using spectra from Keck DEIMOS (Deep Imaging Multi-Object Spectrograph) multi-object spectroscopic observations. The data are used to reconstruct 2D rotation velocity maps. Here we present data for five nearby early-type galaxies to similar to three effective radii. We provide significant insights into the global kinematic structure of these galaxies, and challenge the accepted morphological classification in several cases. We show that between one and three effective radii the velocity dispersion declines very slowly, if at all, in all five galaxies. For the two galaxies with velocity dispersion profiles available from planetary nebulae data we find very good agreement with our stellar profiles. We find a variety of rotation profiles beyond one effective radius, i.e. rotation speed remaining constant, decreasing and increasing with radius. These results are of particular importance to studies which attempt to classify galaxies by their kinematic structure within one effective radius, such as the recent definition of fast- and slow-rotator classes by the Spectrographic Areal Unit for Research on Optical Nebulae project. Our data suggest that the rotator class may change when larger galactocentric radii are probed. This has important implications for dynamical modelling of early-type galaxies. The data from this study are available on-line.

The evolution of the mass-metallicity relation in SDSS galaxies uncovered by astropaleontology

We have obtained the mass-metallicity (M-Z) relation at different lookback times for the same set of galaxies from the Sloan Digital Sky Survey, using the stellar metallicities estimated with our spectral synthesis code STARLIGHT. We have found that this relation steepens and spans a wider range in both mass and metallicity at higher redshifts. We have modelled the time evolution of stellar metallicity with a closed-box chemical evolution model, for galaxies of different types and masses. Our results suggest that the M-Z relation for galaxies with present-day stellar masses down to 10(10) M(circle dot) is mainly driven by the history of star formation and not by inflows or outflows.

The influence of halo assembly on galaxies and galaxy groups

In this paper, we study the variations of groups (galaxy properties according to the assembly history in Sloan Digital Sky Survey Data Release 6 (SDSS-DR6) selected groups. Using mock SDSS group catalogues, we find two suitable indicators of group formation time: (i) the isolation of the group, defined as the distance to the nearest neighbour ill terms of its virial radius and 00 the concentration. measured as the groups inner density calculated using the fifth nearest bright galaxy to the groups centre. Groups Within narrow ranges of Mass ill the mock catalO-Lie show increasing Ifl-OLIP alle With isolation and concentration. However, in the observational data the stellar age, as indicated by the spectral type, only shows a correlation with concentration. We study groups of similar mass and different assembly history. finding important differences ill their galaxy population. Particularly, ill high-mass SDSS groups. the number of members. mass-to-light ratios, red galaxy fractions and the magnitude difference between the brightest and second-brightest group galaxies, show different trends as a function of isolation and concentration, even when it is expected that the latter two quantities correlate with group age. Conversely. low-mass SDSS groups appear to be less sensitive to their assembly history. The correlations detected in the SDSS are not consistent with the trends measured in the mock catalogues. However, discrepancies can he explained in terms of the disagreement found in the a-e-isolation trends, suggesting that the model might be overestimating the effects of, environment, We discuss how the modelling of the cold gas ill satellite galaxies could be responsible for this problem. These results call be Used to improve our Understanding of the evolution of galaxies ill high-density environments.

Blowout bifurcation and spatial mode excitation in the bubbling transition to turbulence

The transition to turbulence (spatio-temporal chaos) in a wide class of spatially extended dynamical system is due to the loss of transversal stability of a chaotic attractor lying on a homogeneous manifold (in the Fourier phase space of the system) causing spatial mode excitation Since the latter manifests as intermittent spikes this has been called a bubbling transition We present numerical evidences that this transition occurs due to the so called blowout bifurcation whereby the attractor as a whole loses transversal stability and becomes a chaotic saddle We used a nonlinear three-wave interacting model with spatial diffusion as an example of this transition (C) 2010 Elsevier B V All rights reserved

Low-dimensional chaos and wave turbulence in plasmas

We investigated drift-wave turbulence in the plasma edge of a small tokamak by considering solutions of the Hasegawa-Mima equation involving three interacting modes in Fourier space. The resulting low-dimensional dynamics presented periodic as well as chaotic evolution of the Fourier-mode amplitudes, and we performed the control of chaotic behaviour through the application of a fourth resonant wave of small amplitude.

Generation of Vortices and Observation of Quantum Turbulence in an Oscillating Bose-Einstein Condensate

We report on the experimental observation of vortex formation and production of tangled vortex distribution in an atomic BEC of (87)Rb atoms submitted to an external oscillatory perturbation. The oscillatory perturbations start by exciting quadrupolar and scissors modes of the condensate. Then regular vortices are observed finally evolving to a vortex tangle configuration. The vortex tangle is a signature of the presence of a turbulent regime in the cloud. We also show that this turbulent cloud has suppression of the aspect ratio inversion typically observed in quantum degenerate bosonic gases during free expansion.

Transition to quantum turbulence in finite-size superfluids

A novel concept of quantum turbulence in finite size superfluids, such as trapped bosonic atoms, is discussed. We have used an atomic (87)Rb Bose-Einstein condensate (BEC) to study the emergence of this phenomenon. In our experiment, the transition to the quantum turbulent regime is characterized by a tangled vortex lines formation, controlled by the amplitude and time duration of the excitation produced by an external oscillating field. A simple model is suggested to account for the experimental observations. The transition from the non-turbulent to the turbulent regime is a rather gradual crossover. But it takes place in a sharp enough way, allowing for the definition of an effective critical line separating the regimes. Quantum turbulence emerging in a finite-size superfluid may be a new idea helpful for revealing important features associated to turbulence, a more general and broad phenomenon. [GRAPHICS] Amplitude versus elapsed time diagram of magnetically excited BEC superfluid, presenting the evolution from the non-turbulent regime, with well separated vortices, to the turbulent regimes, with tangled vortices (C) 2011 by Astro Ltd. Published exclusively by WILEY-VCH Verlag GmbH & Co. KGaA

Route to turbulence in a trapped Bose-Einstein condensate

We have studied a Bose-Einstein condensate of (87)Rb atoms under an oscillatory excitation. For a fixed frequency of excitation, we have explored how the values of amplitude and time of excitation must be combined in order to produce quantum turbulence in the condensate. Depending on the combination of these parameters different behaviors are observed in the sample. For the lowest values of time and amplitude of excitation, we observe a bending of the main axis of the cloud. Increasing the amplitude of excitation we observe an increasing number of vortices. The vortex state can evolve into the turbulent regime if the parameters of excitation are driven up to a certain set of combinations. If the value of the parameters of these combinations is exceeded, all vorticity disappears and the condensate enters into a different regime which we have identified as the granular phase. Our results are summarized in a diagram of amplitude versus time of excitation in which the different structures can be identified. We also present numerical simulations of the Gross-Pitaevskii equation which support our observations. (C) 2011 by Astro Ltd. Published exclusively by WILEY-VCH Verlag GmbH & Co. KGaA

Lie algebra methods for the applications to the statistical theory of turbulence

Approximate Lie symmetries of the Navier-Stokes equations are used for the applications to scaling phenomenon arising in turbulence. In particular, we show that the Lie symmetries of the Euler equations are inherited by the Navier-Stokes equations in the form of approximate symmetries that allows to involve the Reynolds number dependence into scaling laws. Moreover, the optimal systems of all finite-dimensional Lie subalgebras of the approximate symmetry transformations of the Navier-Stokes are constructed. We show how the scaling groups obtained can be used to introduce the Reynolds number dependence into scaling laws explicitly for stationary parallel turbulent shear flows. This is demonstrated in the framework of a new approach to derive scaling laws based on symmetry analysis [11]-[13].