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.

Analysis of total visual and CCD V-broadband observations of Comet C/1995 O1 (Hale-Bopp): 1995-2001

We developed a general method for determination of water production rates from groundbased visual observations and applied it to Comet Hale-Bopp. Our main objective is to extend the method to include total visual magnitude observations obtained with CCD detector and V filter in the analysis of total visual magnitudes. We compare the CCD V-broadband careful observations of Liller [Liller, W. Pre-perihelion CCD photometry of Comet 1995 01 (Hale-Bopp). Planet. Space Sci. 45, 1505-1513, 1997; Liller, W. CCD photometry of Comet C/1995 O1 (Hale-Bopp): 1995-2000. Int. Comet Quart. 23(3), 93-97, 2001] with the total visual magnitude observations from experienced international observers found in the International Comet Quarterly (ICQ) archive. A data set of similar to 400 CCD observations covering about the same 6 years time span of the similar to 12,000 ICQ selected total visual magnitude observations were used in the analysis. A least-square method applied to the water production rates, yields power laws as a function of the heliocentric distances for the pre- and post-perihelion phases. The average dimension of the nucleus as well as its effective active area is determined and compared with values published in the literature. (C) 2009 COSPAR. Published by Elsevier Ltd. All rights reserved.

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.

An analysis of the composite stellar population in M32

We obtained long-slit spectra of high signal-to-noise ratio of the galaxy M32 with the Gemini Multi-Object Spectrograph at the Gemini-North telescope. We analysed the integrated spectra by means of full spectral fitting in order to extract the mixture of stellar populations that best represents its composite nature. Three different galactic radii were analysed, from the nuclear region out to 2 arcmin from the centre. This allows us to compare, for the first time, the results of integrated light spectroscopy with those of resolved colour-magnitude diagrams from the literature. As a main result we propose that an ancient and an intermediate-age population co-exist in M32, and that the balance between these two populations change between the nucleus and outside one effective radius (1r(eff)) in the sense that the contribution from the intermediate population is larger at the nuclear region. We retrieve a smaller signal of a young population at all radii whose origin is unclear and may be a contamination from horizontal branch stars, such as the ones identified by Brown et al. in the nuclear region. We compare our metallicity distribution function for a region 1 to 2 arcmin from the centre to the one obtained with photometric data by Grillmair et al. Both distributions are broad, but our spectroscopically derived distribution has a significant component with [Z/Z(circle dot)] <= -1, which is not found by Grillmair et al.

Comparative study of productivity of the Rosetta target Comet 67P/Churyumov-Gerasimenko

In an attempt to evaluate correlations between several properties of comets we report the results of a cometary research involving a criterious analysis of gas and dust mass production rates in Comets 67P/Churyumov-Gerasimenko (main target of Rosetta Mission), 1P/Halley, Hyakutake (C/1996 B2), and 46P/Wirtanen and make a comparison between them. (C) 2009 COSPAR. Published by Elsevier Ltd. All rights reserved.

Local star formation triggered by supernova shocks in magnetized diffuse neutral clouds

In this work, considering the impact of a supernova remnant (SNR) with a neutral magnetized cloud we derived analytically a set of conditions that are favourable for driving gravitational instability in the cloud and thus star formation. Using these conditions, we have built diagrams of the SNR radius, R(SNR), versus the initial cloud density, n(c), that constrain a domain in the parameter space where star formation is allowed. This work is an extension to previous study performed without considering magnetic fields (Melioli et al. 2006, hereafter Paper I). The diagrams are also tested with fully three-dimensional MHD radiative cooling simulations involving a SNR and a self-gravitating cloud and we find that the numerical analysis is consistent with the results predicted by the diagrams. While the inclusion of a homogeneous magnetic field approximately perpendicular to the impact velocity of the SNR with an intensity similar to 1 mu G within the cloud results only a small shrinking of the star formation zone in the diagram relative to that without magnetic field, a larger magnetic field (similar to 10 mu G) causes a significant shrinking, as expected. Though derived from simple analytical considerations these diagrams provide a useful tool for identifying sites where star formation could be triggered by the impact of a supernova blast wave. Applications of them to a few regions of our own Galaxy (e.g. the large CO shell in the direction of Cassiopeia, and the Edge Cloud 2 in the direction of the Scorpious constellation) have revealed that star formation in those sites could have been triggered by shock waves from SNRs for specific values of the initial neutral cloud density and the SNR radius. Finally, we have evaluated the effective star formation efficiency for this sort of interaction and found that it is generally smaller than the observed values in our own Galaxy (SFE similar to 0.01-0.3). This result is consistent with previous work in the literature and also suggests that the mechanism presently investigated, though very powerful to drive structure formation, supersonic turbulence and eventually, local star formation, does not seem to be sufficient to drive global star formation in normal star-forming galaxies, not even when the magnetic field in the neutral clouds is neglected.

Multi-Conjugate Adaptive Optics VLT imaging of the distant old open cluster FSR 1415

We employ the recently installed near-infrared Multi-Conjugate Adaptive Optics demonstrator (MAD) to determine the basic properties of a newly identified, old and distant, Galactic open cluster (FSR 1415). The MAD facility remarkably approaches the diffraction limit, reaching a resolution of 0.07 arcsec (in K), that is also uniform in a field of similar to 1.8 arcmin in diameter. The MAD facility provides photometry that is 50 per cent complete at K similar to 19. This corresponds to about 2.5 mag below the cluster main-sequence turn-off. This high-quality data set allows us to derive an accurate heliocentric distance of 8.6 kpc, a metallicity close to solar and an age of similar to 2.5 Gyr. On the other hand, the deepness of the data allows us to reconstruct (completeness-corrected) mass functions (MFs) indicating a relatively massive cluster, with a flat core MF. The Very Large Telescope/MAD capabilities will therefore provide fundamental data for identifying/analysing other faint and distant open clusters in the Galaxy III and IV quadrants.

Metallicity gradients at large galactocentric radii using the near-infrared Calcium triplet

We describe a new spectroscopic technique for measuring radial metallicity gradients out to large galactocentric radii. We use the DEep Imaging Multi-Object Spectrograph (DEIMOS) on the Keck telescope and the galaxy spectrum extraction technique of Proctor et al. We also make use of the metallicity sensitive near-infrared Calcium ii triplet (CaT) features together with single stellar population models to obtain metallicities. Our technique is applied as a pilot study to a sample of three relatively nearby (< 30 Mpc) intermediate-mass to massive early-type galaxies. Results are compared with previous literature inner region values and generally show good agreement. We also include a comparison with profiles from dissipational disc-disc major merger simulations. Based on our new extended metallicity gradients combined with other observational evidence and theoretical predictions, we discuss possible formation scenarios for the galaxies in our sample. The limitations of our new technique are also discussed.

Abundances in the Galactic bulge

The metallicity distribution and abundance ratios of the Galactic bulge are reviewed. Issues raised by recent work of different groups, in particular the high metallicity end, the overabundance of alpha-elements in the bulge relative to the thick disc and the measurement of giants versus dwarfs, are discussed. Abundances in the old moderately metal-poor bulge globular clusters are described.

THE VORONOI TESSELLATION CLUSTER FINDER IN 2+1 DIMENSIONS

We present a detailed description of the Voronoi Tessellation (VT) cluster finder algorithm in 2+1 dimensions, which improves on past implementations of this technique. The need for cluster finder algorithms able to produce reliable cluster catalogs up to redshift 1 or beyond and down to 10(13.5) solar masses is paramount especially in light of upcoming surveys aiming at cosmological constraints from galaxy cluster number counts. We build the VT in photometric redshift shells and use the two-point correlation function of the galaxies in the field to both determine the density threshold for detection of cluster candidates and to establish their significance. This allows us to detect clusters in a self-consistent way without any assumptions about their astrophysical properties. We apply the VT to mock catalogs which extend to redshift 1.4 reproducing the ACDM cosmology and the clustering properties observed in the Sloan Digital Sky Survey data. An objective estimate of the cluster selection function in terms of the completeness and purity as a function of mass and redshift is as important as having a reliable cluster finder. We measure these quantities by matching the VT cluster catalog with the mock truth table. We show that the VT can produce a cluster catalog with completeness and purity > 80% for the redshift range up to similar to 1 and mass range down to similar to 10(13.5) solar masses.

Segregation effects according to the evolutionary stage of galaxy groups

We study segregation phenomena in 57 groups selected from the 2dF Percolation-Inferred Galaxy Groups (2PIGG) catalogue of galaxy groups. The sample corresponds to those systems located in areas of at least 80 per cent redshift coverage out to 10 times the radius of the groups. The dynamical state of the galaxy systems was determined after studying their velocity distributions. We have used the Anderson-Darling test to distinguish relaxed and non-relaxed systems. This analysis indicates that 84 per cent of groups have galaxy velocities consistent with the normal distribution, while 16 per cent of them have more complex underlying distributions. Properties of the member galaxies are investigated taking into account this classification. Our results indicate that galaxies in Gaussian groups are significantly more evolved than galaxies in non-relaxed systems out to distances of similar to 4R(200), presenting significantly redder (B - R) colours. We also find evidence that galaxies with M(R) <= -21.5 in Gaussian groups are closer to the condition of energy equipartition.

CDM accelerating cosmology as an alternative to Lambda CDM model

A new accelerating cosmology driven only by baryons plus cold dark matter (CDM) is proposed in the framework of general relativity. In this scenario the present accelerating stage of the Universe is powered by the negative pressure describing the gravitationally-induced particle production of cold dark matter particles. This kind of scenario has only one free parameter and the differential equation governing the evolution of the scale factor is exactly the same of the Lambda CDM model. For a spatially flat Universe, as predicted by inflation (Omega(dm) + Omega(baryon) = 1), it is found that the effectively observed matter density parameter is Omega(meff) = 1 - alpha, where alpha is the constant parameter specifying the CDM particle creation rate. The supernovae test based on the Union data (2008) requires alpha similar to 0.71 so that Omega(meff) similar to 0.29 as independently derived from weak gravitational lensing, the large scale structure and other complementary observations.

Observational constraints on the dark energy and dark matter mutual coupling

We examine different phenomenological interaction models for Dark Energy and Dark Matter by performing statistical joint analysis with observational data arising from the 182 Gold type la supernova samples, the shift parameter of the Cosmic Microwave Background given by the three-year Wilkinson Microwave Anisotropy Probe observations, the baryon acoustic oscillation measurement from the Sloan Digital Sky Survey and age estimates of 35 galaxies. Including the time-dependent observable, we add sensitivity of measurement and give complementary results for the fitting. The compatibility among three different data sets seem to imply that the coupling between dark energy and dark matter is a small positive value, which satisfies the requirement to solve the coincidence problem and the second law of thermodynamics, being compatible with previous estimates. (c) 2008 Elsevier B.V. All rights reserved.