ANALYSIS OF THE PICO DOS DIAS SURVEY HERBIG Ae/Be CANDIDATES

A large sample of Herbig Ae/Be (HAeBe) candidates, distributed in different Galactic regions south to declination +30 degrees, were identified by the Pico dos Dias Survey (a search for young stellar objects based on IRAS colors). Most of the candidates are nearby or associated with star-forming clouds, but several others are considered isolated objects. Aiming to verify the young nature of 93 HAeBe candidates, we searched for additional information that could be useful to confirm if they are pre-main-sequence (PMS) stars or evolved objects, which coincidentally show similar IRAS colors. By adopting a spectral index that is related to the amount of infrared excess and the shape of the spectral energy distribution, we have classified the sample according to three groups, which are analyzed on the basis of (1) circumstellar luminosity; (2) spatial distribution; (3) optical polarization; (4) near-infrared colors; (5) stellar parameters (mass, age, effective temperature); and (5) intensity of emission lines. Our analysis indicates that only 76% of the studied sample, mainly the group with intermediate to low levels of circumstellar emission, can be more confidently considered PMS stars. The nature of the remaining stars, which are in the other group that contains the highest levels of infrared excess, remains to be confirmed. They share the same characteristics of evolved objects, requiring complementary studies in order to correctly classify them. At least seven objects show characteristics typical of post-asymptotic giant branch or proto-planetary nebulae.

HR Del REMNANT ANATOMY USING TWO-DIMENSIONAL SPECTRAL DATA AND THREE-DIMENSIONAL PHOTOIONIZATION SHELL MODELS

The HR Del nova remnant was observed with the IFU-GMOS at Gemini North. The spatially resolved spectral data cube was used in the kinematic, morphological, and abundance analysis of the ejecta. The line maps show a very clumpy shell with two main symmetric structures. The first one is the outer part of the shell seen in H alpha, which forms two rings projected in the sky plane. These ring structures correspond to a closed hourglass shape, first proposed by Harman & O'Brien. The equatorial emission enhancement is caused by the superimposed hourglass structures in the line of sight. The second structure seen only in the [O III] and [N II] maps is located along the polar directions inside the hourglass structure. Abundance gradients between the polar caps and equatorial region were not found. However, the outer part of the shell seems to be less abundant in oxygen and nitrogen than the inner regions. Detailed 2.5-dimensional photoionization modeling of the three-dimensional shell was performed using the mass distribution inferred from the observations and the presence of mass clumps. The resulting model grids are used to constrain the physical properties of the shell as well as the central ionizing source. A sequence of three-dimensional clumpy models including a disk-shaped ionization source is able to reproduce the ionization gradients between polar and equatorial regions of the shell. Differences between shell axial ratios in different lines can also be explained by aspherical illumination. A total shell mass of 9 x 10(-4) M(circle dot) is derived from these models. We estimate that 50%-70% of the shell mass is contained in neutral clumps with density contrast up to a factor of 30.

Galaxy Disruption Caught in the Act

Direct evidence of stellar material from galaxy disruption in the intra-cluster medium (ICM) relies on challenging observations of individual stars, planetary nebulae and diffuse optical light. Here we show that the ultra-compact dwarf galaxies (UCDs) we have discovered in the Fornax Cluster are a new and easy-to-measure probe of disruption in the ICM. We present spectroscopic observations supporting the hypothesis that the UCDs are the remnant nuclei of tidally threshed dwarf galaxies. Deep optical imaging of the cluster has revealed a 43-kpc long arc of tidal debris, flanking a nucleated dwarf elliptical (dE,N) cluster member. We may be witnessing galaxy threshing in action.

Deep GMRT radio observations and a multi-wavelength study of the region around HESS J1858+020

Context. There are a number of very high energy sources in the Galaxy that remain unidentified. Multi-wavelength and variability studies, and catalogue searches, are powerful tools to identify the physical counterpart, given the uncertainty in the source location and extension. Aims. This work carries out a thorough multi-wavelength study of the unidentified, very high energy source HESS J1858+020 and its environs. Methods. We have performed Giant Metrewave Radio Telescope observations at 610 MHz and 1.4 GHz to obtain a deep, low-frequency radio image of the region surrounding HESS J1858+020. We analysed archival radio, infrared, and X-ray data as well. This observational information, combined with molecular data, catalogue sources, and a nearby Fermi gamma-ray detection of unidentified origin, are combined to explore possible counterparts to the very high energy source. Results. We provide with a deep radio image of a supernova remnant that might be related to the GeV and TeV emission in the region. We confirm the presence of an H ii region next to the supernova remnant and coincident with molecular emission. A potential region of star formation is also identified. We identify several radio and X-ray sources in the surroundings. Some of these sources are known planetary nebulae, whereas others may be non-thermal extended emitters and embedded young stellar objects. Three old, background Galactic pulsars also neighbour HESS J1858+020 along the line of sight. Conclusions. The region surrounding HESS J1858+020 is rich in molecular structures and non-thermal objects that may potentially be linked to this unidentified very high energy source. In particular, a supernova remnant interacting with nearby molecular clouds may be a good candidate, but a star forming region, or a non-thermal radio source of yet unclear nature, may also be behind the gamma-ray source. The neighbouring pulsars, despite being old and distant, cannot be discarded as candidates. Further observational studies are needed, however, to narrow the search for a counterpart to the HESS source.

The Coordinated Radio and Infrared Survey for High-mass Star Formation. II. Source Catalog

The CORNISH project is the highest resolution radio continuum survey of the Galactic plane to date. It is the 5 GHz radio continuum part of a series of multi-wavelength surveys that focus on the northern GLIMPSE region (10° < l < 65°), observed by the Spitzer satellite in the mid-infrared. Observations with the Very Large Array in B and BnA configurations have yielded a 1.''5 resolution Stokes I map with a root mean square noise level better than 0.4 mJy beam 1. Here we describe the data-processing methods and data characteristics, and present a new, uniform catalog of compact radio emission. This includes an implementation of automatic deconvolution that provides much more reliable imaging than standard CLEANing. A rigorous investigation of the noise characteristics and reliability of source detection has been carried out. We show that the survey is optimized to detect emission on size scales up to 14'' and for unresolved sources the catalog is more than 90% complete at a flux density of 3.9 mJy. We have detected 3062 sources above a 7σ detection limit and present their ensemble properties. The catalog is highly reliable away from regions containing poorly sampled extended emission, which comprise less than 2% of the survey area. Imaging problems have been mitigated by down-weighting the shortest spacings and potential artifacts flagged via a rigorous manual inspection with reference to the Spitzer infrared data. We present images of the most common source types found: H II regions, planetary nebulae, and radio galaxies. The CORNISH data and catalog are available online at http://cornish.leeds.ac.uk.

CHEMICAL EVOLUTION OF THE GALACTIC BULGE: SINGLE AND DOUBLE INFALL MODELS

Recent work, has produced a wealth of data concerning the chemical evolution of the Galactic bulge, both for stars and nebulae. Present theoretical models generally adopt it limited range of such constraints, frequenfly using it single chemical element (usually iron), which is not enough to describe it unambiguously. In this work, we take into account contraints involving,9 Many chemical elements as possible, basically obtained from bulge nebulae and stars. Our main goal is to show that different scenarios can describe, at least partially the abundance distribution and several dishuice-independent correlations for these objects . Three classes of models were developed. The first is it one-zone, single-infall model, the. Second is it one-zone, double-infall model and the third is a multizone, double-infall model. We show that a one-zone model with it single infall episode is able to reproduce some of the observational data, but the best results tire achieved using it multizone, double-infall model.

Modelling the warm H(2) infrared emission of the Helix nebula cometary knots

Molecular hydrogen emission is commonly observed in planetary nebulae. Images taken in infrared H(2) emission lines show that at least part of the molecular emission is produced inside the ionized region. In the best studied case, the Helix nebula, the H(2) emission is produced inside cometary knots (CKs), comet-shaped structures believed to be clumps of dense neutral gas embedded within the ionized gas. Most of the H(2) emission of the CKs seems to be produced in a thin layer between the ionized diffuse gas and the neutral material of the knot, in a mini-photodissociation region (mini-PDR). However, PDR models published so far cannot fully explain all the characteristics of the H(2) emission of the CKs. In this work, we use the photoionization code AANGABA to study the H(2) emission of the CKs, particularly that produced in the interface H(+)/H(0) of the knot, where a significant fraction of the H(2) 1-0 S(1) emission seems to be produced. Our results show that the production of molecular hydrogen in such a region may explain several characteristics of the observed emission, particularly the high excitation temperature of the H(2) infrared lines. We find that the temperature derived from H(2) observations, even of a single knot, will depend very strongly on the observed transitions, with much higher temperatures derived from excited levels. We also proposed that the separation between the H alpha and [N II] peak emission observed in the images of CKs may be an effect of the distance of the knot from the star, since for knots farther from the central star the [N II] line is produced closer to the border of the CK than H alpha.

VVV DR1: The first data release of the Milky Way bulge and southern plane from the near-infrared ESO public survey VISTA variables in the Via Lactea

Context. The ESO public survey VISTA variables in the Via Lactea (VVV) started in 2010. VVV targets 562 sq. deg in the Galactic bulge and an adjacent plane region and is expected to run for about five years. Aims. We describe the progress of the survey observations in the first observing season, the observing strategy, and quality of the data obtained. Methods. The observations are carried out on the 4-m VISTA telescope in the ZYJHK(s) filters. In addition to the multi-band imaging the variability monitoring campaign in the K-s filter has started. Data reduction is carried out using the pipeline at the Cambridge Astronomical Survey Unit. The photometric and astrometric calibration is performed via the numerous 2MASS sources observed in each pointing. Results. The first data release contains the aperture photometry and astrometric catalogues for 348 individual pointings in the ZYJHK(s) filters taken in the 2010 observing season. The typical image quality is similar to 0 ''.9-1 ''.0. The stringent photometric and image quality requirements of the survey are satisfied in 100% of the JHK(s) images in the disk area and 90% of the JHK(s) images in the bulge area. The completeness in the Z and Y images is 84% in the disk, and 40% in the bulge. The first season catalogues contain 1.28 x 10(8) stellar sources in the bulge and 1.68 x 10(8) in the disk area detected in at least one of the photometric bands. The combined, multi-band catalogues contain more than 1.63 x 10(8) stellar sources. About 10% of these are double detections because of overlapping adjacent pointings. These overlapping multiple detections are used to characterise the quality of the data. The images in the JHK(s) bands extend typically similar to 4 mag deeper than 2MASS. The magnitude limit and photometric quality depend strongly on crowding in the inner Galactic regions. The astrometry for K-s = 15-18 mag has rms similar to 35-175 mas. Conclusions. The VVV Survey data products offer a unique dataset to map the stellar populations in the Galactic bulge and the adjacent plane and provide an exciting new tool for the study of the structure, content, and star-formation history of our Galaxy, as well as for investigations of the newly discovered star clusters, star-forming regions in the disk, high proper motion stars, asteroids, planetary nebulae, and other interesting objects.

Diffuse optical light in galaxy clusters. I. Abell 3888

We are undertaking a program to measure the characteristics of the intracluster light ( ICL; total flux, profile, color, and substructure) in a sample of 10 galaxy clusters with a range of cluster mass, morphology, and redshift. We present here the methods and results for the first cluster in that sample, A3888. We have identified an ICL component in A3888 in V and r that contains 13% +/- 5% of the total cluster light and extends to 700 h(70)(-1) kpc (similar to 0.3r(200)) from the center of the cluster. The ICL color in our smallest radial bin is V - r 0.3 +/- 0.1, similar to the central cluster elliptical galaxies. The ICL is redder than the galaxies at 400 h(70)(-1) kpc < r < 700 h(70)(-1) kpc, although the uncertainty in any one radial bin is high. Based on a comparison of V - r color with simple stellar models, the ICL contains a component that formed more than 7 Gyr ago ( at z less than 1) with a high-metallicity ( 1.0 Z(circle dot) < Z(ICL) less than or similar to 2.5 Z(circle dot)) and a more centralized component that contains stars formed within the past 5 Gyr ( at z similar to 1). The profile of the ICL can be roughly fitted by a shallow exponential in the outer regions and a steeper exponential in the central region. We also find a concentration of diffuse light around a small group of galaxies 1.4 h(70)(-1) Mpc from the center of the cluster. In addition, we find three low surface brightness features near the cluster center that are blue ( V - r 0.0) and contain a total flux of 0.1M*. Based on these observations and X-ray and galaxy morphology, we suggest that this cluster is entering a phase of significant merging of galaxy groups in the core, whereupon we expect the ICL fraction to grow significantly with the formation of a cD galaxy, as well as the infall of groups.

A new general catalogue of nebulae and clusters of stars, being the catalogue of the late Sir John F.W. Herschel.

Mode of access: Internet.

Tropospheric planetary wave dynamics and mixture modeling: Two preferred regimes and a regime shift

Investigation of preferred structures of planetary wave dynamics is addressed using multivariate Gaussian mixture models. The number of components in the mixture is obtained using order statistics of the mixing proportions, hence avoiding previous difficulties related to sample sizes and independence issues. The method is first applied to a few low-order stochastic dynamical systems and data from a general circulation model. The method is next applied to winter daily 500-hPa heights from 1949 to 2003 over the Northern Hemisphere. A spatial clustering algorithm is first applied to the leading two principal components (PCs) and shows significant clustering. The clustering is particularly robust for the first half of the record and less for the second half. The mixture model is then used to identify the clusters. Two highly significant extratropical planetary-scale preferred structures are obtained within the first two to four EOF state space. The first pattern shows a Pacific-North American (PNA) pattern and a negative North Atlantic Oscillation (NAO), and the second pattern is nearly opposite to the first one. It is also observed that some subspaces show multivariate Gaussianity, compatible with linearity, whereas others show multivariate non-Gaussianity. The same analysis is also applied to two subperiods, before and after 1978, and shows a similar regime behavior, with a slight stronger support for the first subperiod. In addition a significant regime shift is also observed between the two periods as well as a change in the shape of the distribution. The patterns associated with the regime shifts reflect essentially a PNA pattern and an NAO pattern consistent with the observed global warming effect on climate and the observed shift in sea surface temperature around the mid-1970s.

Climate entropy budget of the HadCM3 atmosphere-ocean general circulation model and of FAMOUS, its low resolution version

The entropy budget is calculated of the coupled atmosphere–ocean general circulation model HadCM3. Estimates of the different entropy sources and sinks of the climate system are obtained directly from the diabatic heating terms, and an approximate estimate of the planetary entropy production is also provided. The rate of material entropy production of the climate system is found to be ∼50 mW m−2 K−1, a value intermediate in the range 30–70 mW m−2 K−1 previously reported from different models. The largest part of this is due to sensible and latent heat transport (∼38 mW m−2 K−1). Another 13 mW m−2 K−1 is due to dissipation of kinetic energy in the atmosphere by friction and Reynolds stresses. Numerical entropy production in the atmosphere dynamical core is found to be about 0.7 mW m−2 K−1. The material entropy production within the ocean due to turbulent mixing is ∼1 mW m−2 K−1, a very small contribution to the material entropy production of the climate system. The rate of change of entropy of the model climate system is about 1 mW m−2 K−1 or less, which is comparable with the typical size of the fluctuations of the entropy sources due to interannual variability, and a more accurate closure of the budget than achieved by previous analyses. Results are similar for FAMOUS, which has a lower spatial resolution but similar formulation to HadCM3, while more substantial differences are found with respect to other models, suggesting that the formulation of the model has an important influence on the climate entropy budget. Since this is the first diagnosis of the entropy budget in a climate model of the type and complexity used for projection of twenty-first century climate change, it would be valuable if similar analyses were carried out for other such models.

Constraints on wave drag parameterization schemes for simulating the quasi-biennial oscillation. Part II: combined effects of gravity waves and equatorial planetary waves

This study examines the effect of combining equatorial planetary wave drag and gravity wave drag in a one-dimensional zonal mean model of the quasi-biennial oscillation (QBO). Several different combinations of planetary wave and gravity wave drag schemes are considered in the investigations, with the aim being to assess which aspects of the different schemes affect the nature of the modeled QBO. Results show that it is possible to generate a realistic-looking QBO with various combinations of drag from the two types of waves, but there are some constraints on the wave input spectra and amplitudes. For example, if the phase speeds of the gravity waves in the input spectrum are large relative to those of the equatorial planetary waves, critical level absorption of the equatorial planetary waves may occur. The resulting mean-wind oscillation, in that case, is driven almost exclusively by the gravity wave drag, with only a small contribution from the planetary waves at low levels. With an appropriate choice of wave input parameters, it is possible to obtain a QBO with a realistic period and to which both types of waves contribute. This is the regime in which the terrestrial QBO appears to reside. There may also be constraints on the initial strength of the wind shear, and these are similar to the constraints that apply when gravity wave drag is used without any planetary wave drag. In recent years, it has been observed that, in order to simulate the QBO accurately, general circulation models require parameterized gravity wave drag, in addition to the drag from resolved planetary-scale waves, and that even if the planetary wave amplitudes are incorrect, the gravity wave drag can be adjusted to compensate. This study provides a basis for knowing that such a compensation is possible.

Nonlinear symmetric stability of planetary atmospheres

The energy–Casimir method is applied to the problem of symmetric stability in the context of a compressible, hydrostatic planetary atmosphere with a general equation of state. Formal stability criteria for symmetric disturbances to a zonally symmetric baroclinic flow are obtained. In the special case of a perfect gas the results of Stevens (1983) are recovered. Finite-amplitude stability conditions are also obtained that provide an upper bound on a certain positive-definite measure of disturbance amplitude.

SIMULATIONS OF WINDS OF WEAK-LINED T TAURI STARS. II. THE EFFECTS OF A TILTED MAGNETOSPHERE AND PLANETARY INTERACTIONS

Based on our previous work, we investigate here the effects on the wind and magnetospheric structures of weak-lined T Tauri stars due to a misalignment between the axis of rotation of the star and its magnetic dipole moment vector. In such a configuration, the system loses the axisymmetry presented in the aligned case, requiring a fully three-dimensional (3D) approach. We perform 3D numerical magnetohydrodynamic simulations of stellar winds and study the effects caused by different model parameters, namely the misalignment angle theta(t), the stellar period of rotation, the plasma-beta, and the heating index.. Our simulations take into account the interplay between the wind and the stellar magnetic field during the time evolution. The system reaches a periodic behavior with the same rotational period of the star. We show that the magnetic field lines present an oscillatory pattern. Furthermore, we obtain that by increasing theta(t), the wind velocity increases, especially in the case of strong magnetic field and relatively rapid stellar rotation. Our 3D, time-dependent wind models allow us to study the interaction of a magnetized wind with a magnetized extrasolar planet. Such interaction gives rise to reconnection, generating electrons that propagate along the planet`s magnetic field lines and produce electron cyclotron radiation at radio wavelengths. The power released in the interaction depends on the planet`s magnetic field intensity, its orbital radius, and on the stellar wind local characteristics. We find that a close-in Jupiter-like planet orbiting at 0.05 AU presents a radio power that is similar to 5 orders of magnitude larger than the one observed in Jupiter, which suggests that the stellar wind from a young star has the potential to generate strong planetary radio emission that could be detected in the near future with LOFAR. This radio power varies according to the phase of rotation of the star. For three selected simulations, we find a variation of the radio power of a factor 1.3-3.7, depending on theta(t). Moreover, we extend the investigation done in Vidotto et al. and analyze whether winds from misaligned stellar magnetospheres could cause a significant effect on planetary migration. Compared to the aligned case, we show that the timescale tau(w) for an appreciable radial motion of the planet is shorter for larger misalignment angles. While for the aligned case tau(w) similar or equal to 100 Myr, for a stellar magnetosphere tilted by theta(t) = 30 degrees, tau(w) ranges from similar to 40 to 70 Myr for a planet located at a radius of 0.05 AU. Further reduction on tau(w) might occur for even larger misalignment angles and/or different wind parameters.