The spiral structure of the Galaxy revealed by CS sources and evidence for the 4:1 resonance

We present a map of the spiral structure of the Galaxy, as traced by molecular carbon monosulphide (CS) emission associated with IRAS sources which are believed to be compact H II regions. The CS line velocities are used to determine the kinematic distances of the sources in order to investigate their distribution in the galactic plane. This allows us to use 870 objects to trace the arms, a number larger than that of previous studies based on classical H II regions. The distance ambiguity of the kinematic distances, when it exists, is solved by different procedures, including the latitude distribution and an analysis of the longitude-velocity diagram. The study of the spiral structure is complemented with other tracers: open clusters, Cepheids, methanol masers and H II regions. The well-defined spiral arms are seen to be confined inside the corotation radius, as is often the case in spiral galaxies. We identify a square-shaped sub-structure in the CS map with that predicted by stellar orbits at the 4:1 resonance (four epicycle oscillations in one turn around the galactic centre). The sub-structure is found at the expected radius, based on the known pattern rotation speed and epicycle frequency curve. An inner arm presents an end with strong inwards curvature and intense star formation that we tentatively associate with the region where this arm surrounds the extremity of the bar, as seen in many barred galaxies. Finally, a new arm with concave curvature is found in the Sagitta to Cepheus region of the sky. The observed arms are interpreted in terms of perturbations similar to grooves in the gravitational potential of the disc, produced by crowding of stellar orbits.

Diagnostic tool to analyse colour-magnitude diagrams of poorly populated stellar concentrations

The dynamical processes that lead to open cluster disruption cause its mass to decrease. To investigate such processes from the observational point of view, it is important to identify open cluster remnants (OCRs), which are intrinsically poorly populated. Due to their nature, distinguishing them from field star fluctuations is still an unresolved issue. In this work, we developed a statistical diagnostic tool to distinguish poorly populated star concentrations from background field fluctuations. We use 2MASS photometry to explore one of the conditions required for a stellar group to be a physical group: to produce distinct sequences in a colour-magnitude diagram (CMD). We use automated tools to (i) derive the limiting radius; (ii) decontaminate the field and assign membership probabilities; (iii) fit isochrones; and (iv) compare object and field CMDs, considering the isochrone solution, in order to verify the similarity. If the object cannot be statistically considered as a field fluctuation, we derive its probable age, distance modulus, reddening and uncertainties in a self-consistent way. As a test, we apply the tool to open clusters and comparison fields. Finally, we study the OCR candidates DoDz 6, NGC 272, ESO 435 SC48 and ESO 325 SC15. The tool is optimized to treat these low-statistic objects and to separate the best OCR candidates for studies on kinematics and chemical composition. The study of the possible OCRs will certainly provide a deep understanding of OCR properties and constraints for theoretical models, including insights into the evolution of open clusters and dissolution rates.

Breaks in the radial oxygen abundance and corotation radius of three spiral galaxies

The correlation between the breaks in the metallicity distribution and the corotation radius of spiral galaxies has been already advocated in the past and is predicted by a chemodynamical model of our Galaxy that effectively introduces the role of spiral arms in the star formation rate. In this work, we present photometric and spectroscopic observations made with the Gemini Telescope for three of the best candidates of spiral galaxies to have the corotation inside the optical disc: IC 0167, NGC 1042 and NGC 6907. We observed the most intense and well-distributed H ii regions of these galaxies, deriving reliable galactocentric distances and oxygen abundances by applying different statistical methods. From these results, we confirm the presence of variations in the gradients of metallicity of these galaxies that are possibly correlated with the corotation resonance.

Morphological properties of superclusters of galaxies

We studied superclusters of galaxies in a volume-limited sample extracted from the Sloan Digital Sky Survey Data Release 7 and from mock catalogues based on a semi-analytical model of galaxy evolution in the Millennium Simulation. A density field method was applied to a sample of galaxies brighter than M(r) = -21+5 log h(100) to identify superclusters, taking into account selection and boundary effects. In order to evaluate the influence of the threshold density, we have chosen two thresholds: the first maximizes the number of objects (D1) and the second constrains the maximum supercluster size to similar to 120 h(-1) Mpc (D2). We have performed a morphological analysis, using Minkowski Functionals, based on a parameter, which increases monotonically from filaments to pancakes. An anticorrelation was found between supercluster richness (and total luminosity or size) and the morphological parameter, indicating that filamentary structures tend to be richer, larger and more luminous than pancakes in both observed and mock catalogues. We have also used the mock samples to compare supercluster morphologies identified in position and velocity spaces, concluding that our morphological classification is not biased by the peculiar velocities. Monte Carlo simulations designed to investigate the reliability of our results with respect to random fluctuations show that these results are robust. Our analysis indicates that filaments and pancakes present different luminosity and size distributions.

Nuclear and extended spectra of NGC 1068-II. Near-infrared stellar population synthesis

We performed stellar population synthesis on the nuclear and extended regions of NGC 1068 by means of near-infrared spectroscopy to disentangle their spectral energy distribution components. This is the first time that such a technique is applied to the whole 0.8-2.4 mu m wavelength interval in this galaxy. NGC 1068 is one of the nearest and probably the most studied Seyfert 2 galaxy, becoming an excellent laboratory to study the interaction between black holes, the jets that they can produce and the medium in which they propagate. Our main result is that traces of young stellar population are found at similar to 100 pc south of the nucleus. The contribution of a power-law continuum in the centre is about 25 per cent, which is expected if the light is scattered from a Seyfert 1 nucleus. We find peaks in the contribution of the featureless continuum about 100-150 pc from the nucleus on both sides. They might be associated with regions where the jet encounters dense clouds. Further support to this scenario is given by the peaks of hot dust distribution found around these same regions and the H(2) emission-line profile, leading us to propose that the peaks might be associated to regions where stars are being formed. Hot dust also has an important contribution to the nuclear region, reinforcing the idea of the presence of a dense, circumnuclear torus in this galaxy. Cold dust appears mostly in the south direction, which supports the view that the south-west emission is behind the plane of the galaxy and is extinguished very likely by dust in the plane. Intermediate-age stellar population contributes significantly to the continuum, especially in the inner 200 pc.

Alternative diagnostic diagrams and the `forgotten` population of weak line galaxies in the SDSS

A numerous population of weak line galaxies (WLGs) is often left out of statistical studies on emission-line galaxies (ELGs) due to the absence of an adequate classification scheme, since classical diagnostic diagrams, such as [O iii]/H beta versus [N ii]/H alpha (the BPT diagram), require the measurement of at least four emission lines. This paper aims to remedy this situation by transposing the usual divisory lines between star-forming (SF) galaxies and active galactic nuclei (AGN) hosts and between Seyferts and LINERs to diagrams that are more economical in terms of line quality requirements. By doing this, we rescue from the classification limbo a substantial number of sources and modify the global census of ELGs. More specifically, (1) we use the Sloan Digital Sky Survey Data Release 7 to constitute a suitable sample of 280 000 ELGs, one-third of which are WLGs. (2) Galaxies with strong emission lines are classified using the widely applied criteria of Kewley et al., Kauffmann et al. and Stasinska et al. to distinguish SF galaxies and AGN hosts and Kewley et al. to distinguish Seyferts from LINERs. (3) We transpose these classification schemes to alternative diagrams keeping [N ii]/H alpha as a horizontal axis, but replacing H beta by a stronger line (H alpha or [O ii]), or substituting the ionization-level sensitive [O iii]/H beta ratio with the equivalent width of H alpha (W(H alpha)). Optimized equations for the transposed divisory lines are provided. (4) We show that nothing significant is lost in the translation, but that the new diagrams allow one to classify up to 50 per cent more ELGs. (5) Introducing WLGs in the census of galaxies in the local Universe increases the proportion of metal-rich SF galaxies and especially LINERs. In the course of this analysis, we were led to make the following points. (i) The Kewley et al. BPT line for galaxy classification is generally ill-used. (ii) Replacing [O iii]/H beta by W(H alpha) in the classification introduces a change in the philosophy of the distinction between LINERs and Seyferts, but not in its results. Because the W(H alpha) versus [N ii]/H alpha diagram can be applied to the largest sample of ELGs without loss of discriminating power between Seyferts and LINERs, we recommend its use in further studies. (iii) The dichotomy between Seyferts and LINERs is washed out by WLGs in the BPT plane, but it subsists in other diagnostic diagrams. This suggests that the right wing in the BPT diagram is indeed populated by at least two classes, tentatively identified with bona fide AGN and `retired` galaxies that have stopped forming stars and are ionized by their old stellar populations.

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.

SURFACE ALFVEN WAVE DAMPING IN A THREE-DIMENSIONAL SIMULATION OF THE SOLAR WIND

Here we investigate the contribution of surface Alfven wave damping to the heating of the solar wind in minima conditions. These waves are present in the regions of strong inhomogeneities in density or magnetic field (e.g., the border between open and closed magnetic field lines). Using a three-dimensional (3D) magnetohydrodynamics (MHD) model, we calculate the surface Alfven wave damping contribution between 1 and 4 R(circle dot) (solar radii), the region of interest for both acceleration and coronal heating. We consider waves with frequencies lower than those that are damped in the chromosphere and on the order of those dominating the heliosphere: 3 x 10(-6) to 10(-1) Hz. In the region between open and closed field lines, within a few R(circle dot) of the surface, no other major source of damping has been suggested for the low frequency waves we consider here. This work is the first to study surface Alfven waves in a 3D environment without assuming a priori a geometry of field lines or magnetic and density profiles. We demonstrate that projection effects from the plane of the sky to 3D are significant in the calculation of field line expansion. We determine that waves with frequencies >2.8 x 10(-4) Hz are damped between 1 and 4 R(circle dot). In quiet-Sun regions, surface Alfven waves are damped at further distances compared to active regions, thus carrying additional wave energy into the corona. We compare the surface Alfven wave contribution to the heating by a variable polytropic index and find it as an order of magnitude larger than needed for quiet-Sun regions. For active regions, the contribution to the heating is 20%. As it has been argued that a variable gamma acts as turbulence, our results indicate that surface Alfven wave damping is comparable to turbulence in the lower corona. This damping mechanism should be included self-consistently as an energy driver for the wind in global MHD models.

Further probing the nature of FSR 1767

With Two-Micron All-Sky Survey (2MASS) photometry and proper motions, Bonatto et al. suggested that FSR 1767 is a globular cluster (GC), while with J and K NTT/SOFI photometry Froebrich, Meusinger & Scholz concluded that it is not a star cluster. In this study, we combine previous and new evidence that are consistent with a GC. For instance, we show that the horizontal branch (HB) and red giant branch (RGB) stars, besides sharing a common proper motion, have radial density profiles that consistently follow the King`s law independently. Reddening maps around FSR 1767 are built using the bulge RGB as reference and also Schlegel`s extinction values to study local absorptions. Both approaches provide similar maps and show that FSR 1767 is not located in a dust window, which otherwise might have produced the stellar overdensity. Besides, neighbouring regions of similar reddening as FSR 1767 do not present the blue HB stars that are a conspicuous feature in the colour-magnitude diagram of FSR 1767. We report the presence of a compact group of stars located in the central parts of FSR 1767. It appears to be a detached post-collapse core, similar to those of other nearby low-luminosity GCs projected towards the bulge. We note that while the NTT/SOFI photometry of the star cluster FSR 1716 matches perfectly that from 2MASS, it shows a considerable offset for FSR 1767. We discuss the possible reasons why both photometries differ. We confirm our previous structural and photometric fundamental parameters for FSR 1767, which are consistent with a GC.

TOWARD UNDERSTANDING THE B[e] PHENOMENON. III. PROPERTIES OF THE OPTICAL COUNTERPART OF IRAS 00470+6429

FS CMa type stars are a group of Galactic objects with the B[e] phenomenon. They exhibit strong emission-line spectra and infrared excesses, which are most likely due to recently formed circumstellar dust. The group content and identification criteria were described in the first two papers of the series. In this paper we report our spectroscopic and photometric observations of the optical counterpart of IRAS 00470+6429 obtained in 2003-2008. The optical spectrum is dominated by emission lines, most of which have P Cyg type profiles. We detected significant brightness variations, which may include a regular component, and variable spectral line profiles in both shape and position. The presence of a weak Li I 6708 angstrom line in the spectrum suggests that the object is most likely a binary system with a B2-B3 spectral-type primary companion of a luminosity log L/L(circle dot) = 3.9 +/- 0.3 and a late-type secondary companion. We estimate a distance toward the object to be 2.0 +/- 0.3 kpc from the Sun.

Searching for s-Process-Enhanced Metal-Poor Stars

We present preliminary results for the estimation of barium [Ba/Fe], and strontium [Sr/Fe], abundances ratios using medium-resolution spectra (1-2 angstrom). We established a calibration between the abundance ratios and line indices for Ba and Sr, using multiple regression and artificial neural network techniques. A comparison between the two techniques (showing the advantage of the latter), as well as a discussion of future work, is presented.

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.

Simplified quartessence cosmology

We propose a new class of accelerating world models unifying the cosmological dark sector (dark matter and dark energy). All the models are described by a simplified version of the Chaplygin gas quartessence cosmology. It is found that even for Omega(k) not equal 0, this quartessence scenario depends only on a pair of parameters which can severely be constrained by the cosmological tests. As an example we perform a joint analysis involving the latest SNe type la data and the recent Sloan Digital Sky Survey measurement of baryon acoustic oscillations. In our analysis we have considered the SNe type la Union sample compiled by Kowalski et al. [M. Kowalski et al., Astrophys. J. 686 (2008) 749, arXiv:0804.4142]. At 95.4% (c.l.), we find for BAD + Union sample, alpha = 0.81(-0.04)(+0.04) and Omega(Q4) = 1.15(-0.17)(+0.16) The best-fit for this simplified quartessence scenario is a spatially closed Universe and its reduced chi(2) is exactly the same of the flat concordance model (Lambda CDM). (C) 2009 Elsevier B.V. All rights reserved.

A simplified approach for Chaplygin-type cosmologies

A new class of accelerating cosmological models driven by a one-parameter version of the general Chaplygin-type equation of state is proposed. The simplified version is naturally obtained from causality considerations with basis on the adiabatic sound speed vs plus the observed accelerating stage of the universe. We show that very stringent constraints on the unique free parameter a describing the simplified Chaplygin model can be obtained from a joint analysis involving the latest SNe type la data and the recent Sloan Digital Sky Survey measurement of baryon acoustic oscillations (BAO). In our analysis we have considered separately the SNe type la gold sample measured by [A.G. Riess et al.. Astrophys. J. 607 (2004) 665] and the supernova legacy survey (SNLS) from [P. Astier et al., Astron. Astrophys. 447 (2006) 31]. At 95.4% (c.l.), we find for BAO + gold sample, 0.91 <= alpha <= 1.0 and Omega(m) = 0.28(-0.048)(+0.043) while BAO + SNLS analysis provides 0.94 <= alpha <= 1.0 and Omega(m) = 0.27(-0.045)(+0.048). (C) 2008 Elsevier B.V. All rights reserved.