Recovery of the star formation history of the LMC from the VISTA survey of the Magellanic system

The VISTA near infrared survey of the Magellanic System (VMC) will provide deep YJK(s) photometry reaching stars in the oldest turn-off point throughout the Magellanic Clouds (MCs). As part of the preparation for the survey, we aim to access the accuracy in the star formation history (SFH) that can be expected from VMC data, in particular for the Large Magellanic Cloud (LMC). To this aim, we first simulate VMC images containing not only the LMC stellar populations but also the foreground Milky Way (MW) stars and background galaxies. The simulations cover the whole range of density of LMC field stars. We then perform aperture photometry over these simulated images, access the expected levels of photometric errors and incompleteness, and apply the classical technique of SFH-recovery based on the reconstruction of colour-magnitude diagrams (CMD) via the minimisation of a chi-squared-like statistics. We verify that the foreground MW stars are accurately recovered by the minimisation algorithms, whereas the background galaxies can be largely eliminated from the CMD analysis due to their particular colours and morphologies. We then evaluate the expected errors in the recovered star formation rate as a function of stellar age, SFR(t), starting from models with a known age-metallicity relation (AMR). It turns out that, for a given sky area, the random errors for ages older than similar to 0.4 Gyr seem to be independent of the crowding. This can be explained by a counterbalancing effect between the loss of stars from a decrease in the completeness and the gain of stars from an increase in the stellar density. For a spatial resolution of similar to 0.1 deg(2), the random errors in SFR(t) will be below 20% for this wide range of ages. On the other hand, due to the lower stellar statistics for stars younger than similar to 0.4 Gyr, the outer LMC regions will require larger areas to achieve the same level of accuracy in the SFR( t). If we consider the AMR as unknown, the SFH-recovery algorithm is able to accurately recover the input AMR, at the price of an increase of random errors in the SFR(t) by a factor of about 2.5. Experiments of SFH-recovery performed for varying distance modulus and reddening indicate that these parameters can be determined with (relative) accuracies of Delta(m-M)(0) similar to 0.02 mag and Delta E(B-V) similar to 0.01 mag, for each individual field over the LMC. The propagation of these errors in the SFR(t) implies systematic errors below 30%. This level of accuracy in the SFR(t) can reveal significant imprints in the dynamical evolution of this unique and nearby stellar system, as well as possible signatures of the past interaction between the MCs and the MW.

DECISION TREE CLASSIFIERS FOR STAR/GALAXY SEPARATION

We study the star/galaxy classification efficiency of 13 different decision tree algorithms applied to photometric objects in the Sloan Digital Sky Survey Data Release Seven (SDSS-DR7). Each algorithm is defined by a set of parameters which, when varied, produce different final classification trees. We extensively explore the parameter space of each algorithm, using the set of 884,126 SDSS objects with spectroscopic data as the training set. The efficiency of star-galaxy separation is measured using the completeness function. We find that the Functional Tree algorithm (FT) yields the best results as measured by the mean completeness in two magnitude intervals: 14 <= r <= 21 (85.2%) and r >= 19 (82.1%). We compare the performance of the tree generated with the optimal FT configuration to the classifications provided by the SDSS parametric classifier, 2DPHOT, and Ball et al. We find that our FT classifier is comparable to or better in completeness over the full magnitude range 15 <= r <= 21, with much lower contamination than all but the Ball et al. classifier. At the faintest magnitudes (r > 19), our classifier is the only one that maintains high completeness (> 80%) while simultaneously achieving low contamination (similar to 2.5%). We also examine the SDSS parametric classifier (psfMag - modelMag) to see if the dividing line between stars and galaxies can be adjusted to improve the classifier. We find that currently stars in close pairs are often misclassified as galaxies, and suggest a new cut to improve the classifier. Finally, we apply our FT classifier to separate stars from galaxies in the full set of 69,545,326 SDSS photometric objects in the magnitude range 14 <= r <= 21.

The merging cluster Abell 1758 revisited: multi-wavelength observations and numerical simulations

Context. Cluster properties can be more distinctly studied in pairs of clusters, where we expect the effects of interactions to be strong. Aims. We here discuss the properties of the double cluster Abell 1758 at a redshift z similar to 0.279. These clusters show strong evidence for merging. Methods. We analyse the optical properties of the North and South cluster of Abell 1758 based on deep imaging obtained with the Canada-France-Hawaii Telescope (CFHT) archive Megaprime/Megacam camera in the g' and r' bands, covering a total region of about 1.05 x 1.16 deg(2), or 16.1 x 17.6 Mpc(2). Our X-ray analysis is based on archive XMM-Newton images. Numerical simulations were performed using an N-body algorithm to treat the dark-matter component, a semi-analytical galaxy-formation model for the evolution of the galaxies and a grid-based hydrodynamic code with a parts per million (PPM) scheme for the dynamics of the intra-cluster medium. We computed galaxy luminosity functions (GLFs) and 2D temperature and metallicity maps of the X-ray gas, which we then compared to the results of our numerical simulations. Results. The GLFs of Abell 1758 North are well fit by Schechter functions in the g' and r' bands, but with a small excess of bright galaxies, particularly in the r' band; their faint-end slopes are similar in both bands. In contrast, the GLFs of Abell 1758 South are not well fit by Schechter functions: excesses of bright galaxies are seen in both bands; the faint-end of the GLF is not very well defined in g'. The GLF computed from our numerical simulations assuming a halo mass-luminosity relation agrees with those derived from the observations. From the X-ray analysis, the most striking features are structures in the metal distribution. We found two elongated regions of high metallicity in Abell 1758 North with two peaks towards the centre. In contrast, Abell 1758 South shows a deficit of metals in its central regions. Comparing observational results to those derived from numerical simulations, we could mimic the most prominent features present in the metallicity map and propose an explanation for the dynamical history of the cluster. We found in particular that in the metal-rich elongated regions of the North cluster, winds had been more efficient than ram-pressure stripping in transporting metal-enriched gas to the outskirts. Conclusions. We confirm the merging structure of the North and South clusters, both at optical and X-ray wavelengths.

Hubble parameter reconstruction from a principal component analysis: minimizing the bias

Aims. A model-independent reconstruction of the cosmic expansion rate is essential to a robust analysis of cosmological observations. Our goal is to demonstrate that current data are able to provide reasonable constraints on the behavior of the Hubble parameter with redshift, independently of any cosmological model or underlying gravity theory. Methods. Using type Ia supernova data, we show that it is possible to analytically calculate the Fisher matrix components in a Hubble parameter analysis without assumptions about the energy content of the Universe. We used a principal component analysis to reconstruct the Hubble parameter as a linear combination of the Fisher matrix eigenvectors (principal components). To suppress the bias introduced by the high redshift behavior of the components, we considered the value of the Hubble parameter at high redshift as a free parameter. We first tested our procedure using a mock sample of type Ia supernova observations, we then applied it to the real data compiled by the Sloan Digital Sky Survey (SDSS) group. Results. In the mock sample analysis, we demonstrate that it is possible to drastically suppress the bias introduced by the high redshift behavior of the principal components. Applying our procedure to the real data, we show that it allows us to determine the behavior of the Hubble parameter with reasonable uncertainty, without introducing any ad-hoc parameterizations. Beyond that, our reconstruction agrees with completely independent measurements of the Hubble parameter obtained from red-envelope galaxies.

The clusters Abell 222 and Abell 223: a multi-wavelength view

Context. The Abell 222 and 223 clusters are located at an average redshift z similar to 0.21 and are separated by 0.26 deg. Signatures of mergers have been previously found in these clusters, both in X-rays and at optical wavelengths, thus motivating our study. In X-rays, they are relatively bright, and Abell 223 shows a double structure. A filament has also been detected between the clusters both at optical and X-ray wavelengths. Aims. We analyse the optical properties of these two clusters based on deep imaging in two bands, derive their galaxy luminosity functions (GLFs) and correlate these properties with X-ray characteristics derived from XMM-Newton data. Methods. The optical part of our study is based on archive images obtained with the CFHT Megaprime/Megacam camera, covering a total region of about 1 deg(2), or 12.3 x 12.3 Mpc(2) at a redshift of 0.21. The X-ray analysis is based on archive XMM-Newton images. Results. The GLFs of Abell 222 in the g' and r' bands are well fit by a Schechter function; the GLF is steeper in r' than in g'. For Abell 223, the GLFs in both bands require a second component at bright magnitudes, added to a Schechter function; they are similar in both bands. The Serna & Gerbal method allows to separate well the two clusters. No obvious filamentary structures are detected at very large scales around the clusters, but a third cluster at the same redshift, Abell 209, is located at a projected distance of 19.2 Mpc. X-ray temperature and metallicity maps reveal that the temperature and metallicity of the X-ray gas are quite homogeneous in Abell 222, while they are very perturbed in Abell 223. Conclusions. The Abell 222/Abell 223 system is complex. The two clusters that form this structure present very different dynamical states. Abell 222 is a smaller, less massive and almost isothermal cluster. On the other hand, Abell 223 is more massive and has most probably been crossed by a subcluster on its way to the northeast. As a consequence, the temperature distribution is very inhomogeneous. Signs of recent interactions are also detected in the optical data where this cluster shows a ""perturbed"" GLF. In summary, the multiwavelength analyses of Abell 222 and Abell 223 are used to investigate the connection between the ICM and the cluster galaxy properties in an interacting system.

On the influence of non-thermal pressure on the mass determination of galaxy clusters

Aims. Given that in most cases just thermal pressure is taken into account in the hydrostatic equilibrium equation to estimate galaxy cluster mass, the main purpose of this paper is to consider the contribution of all three non-thermal components to total mass measurements. The non-thermal pressure is composed by cosmic rays, turbulence and magnetic pressures. Methods. To estimate the thermal pressure we used public XMM-Newton archival data of five Abell clusters to derive temperature and density profiles. To describe the magnetic pressure, we assume a radial distribution for the magnetic field, B(r) proportional to rho(alpha)(g). To seek generality we assume alpha within the range of 0.5 to 0.9, as indicated by observations and numerical simulations. Turbulent motions and bulk velocities add a turbulent pressure, which is considered using an estimate from numerical simulations. For this component, we assume an isotropic pressure, P(turb) = 1/3 rho(g)(sigma(2)(r) + sigma(2)(t)). We also consider the contribution of cosmic ray pressure, P(cr) proportional to r(-0.5). Thus, besides the gas (thermal) pressure, we include these three non-thermal components in the magnetohydrostatic equilibrium equation and compare the total mass estimates with the values obtained without them. Results. A consistent description for the non-thermal component could yield a variation in mass estimates that extends from 10% to similar to 30%. We verified that in the inner parts of cool core clusters the cosmic ray component is comparable to the magnetic pressure, while in non-cool core clusters the cosmic ray component is dominant. For cool core clusters the magnetic pressure is the dominant component, contributing more than 50% of the total mass variation due to non-thermal pressure components. However, for non-cool core clusters, the major influence comes from the cosmic ray pressure that accounts for more than 80% of the total mass variation due to non-thermal pressure effects. For our sample, the maximum influence of the turbulent component to the total mass variation can be almost 20%. Although all of the assumptions agree with previous works, it is important to notice that our results rely on the specific parametrization adopted in this work. We show that this analysis can be regarded as a starting point for a more detailed and refined exploration of the influence of non-thermal pressure in the intra-cluster medium (ICM).

Substructure and Dynamics of the Fornax Cluster

We present the first dynamical analysis of a galaxy cluster to include a large fraction of dwarf galaxies. Our sample of 108 Fornax Cluster members measured with the UK Schmidt Telescope FLAIR-II spectrograph contains 55 dwarf galaxies (15.5 > b(j) > 18.0 or -16 > M-B > -13.5). H alpha emission shows that of the dwarfs are star forming, twice the fraction implied by morphological classifications. The total sample has a mean velocity of 1493 +/- 36 kms s(-1) and a velocity dispersion of 374 +/- 26 km s(-1). The dwarf galaxies form a distinct population: their velocity dispersion (429 +/- 41 km s(-1)) is larger than that of the giants () at the 98% confidence level. This suggests that the dwarf population is dominated by infalling objects whereas the giants are virialized. The Fornax system has two components, the main Fornax Cluster centered on NGC 1399 with cz = 1478 km s(-1) and sigma (cz) = 370 km s(-1) and a subcluster centered 3 degrees to the southwest including NGC 1316 with cz = 1583 km s(-1) and sigma (cz) = 377 km s(-1). This partition is preferred over a single cluster at the 99% confidence level. The subcluster, a site of intense star formation, is bound to Fornax and probably infalling toward the cluster core for the first time. We discuss the implications of this substructure for distance estimates of the Fornax Cluster. We determine the cluster mass profile using the method of Diaferio, which does not assume a virialized sample. The mass within a projected radius of 1.4 Mpc is (7 +/- 2) x 10(13) M-., and the mass-to-light ratio is 300 +/- 100 M-./L-.. The mass is consistent with values derived from the projected mass virial estimator and X-ray measurements at smaller radii.

X-ray properties of the Parkes sample of flat-spectrum radio sources: dust in radio-loud quasars?

We investigate the X-ray properties of the Parkes sample of Bat-spectrum radio sources using data from the ROSAT All-Sky Survey and archival pointed PSPC observations. In total, 163 of the 323 sources are detected. For the remaining 160 sources, 2 sigma upper limits to the X-ray flux are derived. We present power-law photon indices in the 0.1-2.4 keV energy band for 115 sources, which were determined either with a hardness ratio technique or from direct fits to pointed PSPC data if a sufficient number of photons were available. The average photon index is <Gamma > = 1.95(-0.12)(+0.13) for flat-spectrum radio-loud quasars, <Gamma > = 1.70(-0.24)(+0.23) for galaxies, and <Gamma > = 2.40(-0.31)(+0.12) for BL Lac objects. The soft X-ray photon index is correlated with redshift and with radio spectral index in the sense that sources at high redshift and/or with flat (or inverted) radio spectra have flatter X-ray spectra on average. The results are in accord with orientation-dependent unification schemes for radio-loud active galactic nuclei. Webster et al. discovered many sources with unusually red optical continua among the quasars of this sample, and interpreted this result in terms of extinction by dust. Although the X-ray spectra in general do not show excess absorption, we find that low-redshift optically red quasars have significantly lower soft X-ray luminosities on average than objects with blue optical continua. The difference disappears for higher redshifts, as is expected for intrinsic absorption by cold gas associated with the dust. In addition, the scatter in log(f(x)/f(o)) is consistent with the observed optical extinction, contrary to previous claims based on optically or X-ray selected samples. Although alternative explanations for the red optical continua cannot be excluded with the present X-ray data, we note that the observed X-ray properties are consistent with the idea that dust plays an important role in some of the radio-loud quasars with red optical continua.

An extragalactic HI cloud with no optical counterpart?

We report the discovery, from the H I Parkes All-Sky Survey (HIPASS), of an isolated cloud of neutral hydrogen, which we believe to be extragalactic. The H I mass of the cloud (HIPASS J1712-64) is very low, 1.7 x 10(7) M-circle dot, using an estimated distance of similar to 3.2 Mpc. Most significantly, we have found no optical companion to this object to very faint limits [mu(B) similar to 27 mag arcsec(-2)]. HIPASS J1712-64 appears to be a binary system similar to, but much less massive than, H I 1225 + 01 (the Virgo H. I cloud) and has a size of at least 15 kpc. The mean velocity dispersion measured with the Australia Telescope Compact Array (ATCA) is only 4 km s(-1) for the main component and, because of the weak or nonexistent star formation, possibly reflects the thermal line width (T < 2000 K) rather than bulk motion or turbulence. The peak column density for HIPASS J1712-64, from the combined Parkes and ATCA data, is only 3.5 x 1019 cm(-2), which is estimated to be a factor of 2 below the critical threshold for star formation. Apart from its significantly higher velocity, the properties of HIPASS J1712-64 are similar to the recently recognized class of compact high-velocity clouds. We therefore consider the evidence for a Local Group or Galactic origin, although a more plausible alternative is that HIPASS J1712-64 was ejected from the interacting Magellanic Cloud-Galaxy system at perigalacticon similar to 2 x 10(8) yr ago.

Quasars from a complete spectroscopic survey

With the advent of multi-fibre spectrographs such as the 'Two-Degree Field' (2dF) instrument at the Angle-Australian Telescope, quasar surveys that are free of any preselection of candidates and any biases this implies have become possible for the first time. The first of these is that which is being undertaken as part of the Fornax Spectroscopic Survey, a survey of the area around the Fornax Cluster of galaxies, and aims to obtain the spectra of all objects in the magnitude range 16.5 < b(j) < 19.7. To date, 3679 objects in the central pi -deg(2) area have been successfully identified from their spectral characteristics. Of these, 71 are found to be quasars, 61 with redshifts 0.3 < z < 2.2 and 10 with redshifts z > 2.2. Using this complete quasar sample, a new determination of quasar number counts is made, enabling an independent check of existing quasars surveys. Cumulative counts per square degree at a magnitude limit of b(j) < 19.5 are found to be 11.5 +/- 2.2 for 0.3 < z < 2.2, 2.22 +/- 0.93 for z > 2.2 and 13.7 +/- 3.1 for z > 0.3. Given the likely detection of extra quasars in the Fornax survey, we make a more detailed examination of existing quasar selection techniques. First, looking at the use of a stellar criterion, four of the 71 quasars are 'non-stellar' on the basis of the automated plate measuring facility (APM) b(j) classification, however inspection shows all are consistent with stellar, but misclassified due to image confusion. Examining the ultraviolet excess and multicolour selection techniques, for the selection criteria investigated, ultraviolet excess would find 69 +/- 6 per cent of our 0.3 < z < 2.2 quasars and only 50(-18)(+14), per cent of our z > 2.2 quasars, while the completeness level for multicolour selection is found to be 90(-4)(+3) per cent for 0.3 < z < 2.2 quasars and 80(-12)(+14) per cent for z > 2.2 quasars. The extra quasars detected by our all-object survey thus have unusually red star-like colours, and this appears to be a result of the continuum shape rather than any emission features. An intrinsic dust extinction model may, at least partly, account for the red colours.

The Las Campanas/AAT rich cluster survey - I. Precision and reliability of the photometric catalogue

The Las Campanas Observatory and Anglo-Australian Telescope Rich Cluster Survey (LARCS) is a panoramic imaging and spectroscopic survey of an X-ray luminosity-selected sample of 21 clusters of galaxies at 0.07 < z < 0.16. Charge-coupled device (CCD) imaging was obtained in B and R of typically 2 degrees wide regions centred on the 21 clusters, and the galaxy sample selected from the imaging is being used for an on-going spectroscopic survey of the clusters with the 2dF spectrograph on the Anglo-Australian Telescope. This paper presents the reduction of the imaging data and the photometric analysis used in the survey. Based on an overlapping area of 12.3 deg(2) we compare the CCD-based LARCS catalogue with the photographic-based galaxy catalogue used for the input to the 2dF Galaxy Redshift Survey (2dFGRS) from the APM, to the completeness of the GRS/APM catalogue, b(J) = 19.45. This comparison confirms the reliability of the photometry across our mosaics and between the clusters in our survey. This comparison also provides useful information concerning the properties of the GRS/APM. The stellar contamination in the GRS/APM galaxy catalogue is confirmed as around 5-10 per cent, as originally estimated. However, using the superior sensitivity and spatial resolution in the LARCS survey evidence is found for four distinct populations of galaxies that are systematically omitted from the GRS/APM catalogue. The characteristics of the 'missing' galaxy populations are described, reasons for their absence examined and the impact they will have on the conclusions drawn from the 2dF Galaxy Redshift Survey are discussed.

Using the 6dF galaxy redshift survey to detect gravitationally-lensed quasars

It is possible to detect gravitationally-lensed quasars spectroscopically if the spectra obtained during galaxy surveys are searched for the presence of quasar emission lines. The up-coming 6 degree Field (6dF) redshift survey on the United Kingdom Schmidt Telescope will involve obtaining similar to 10(5) spectra of near-infrared selected galaxies to a magnitude limit of K = 13. Applying previously developed techniques implies that at least one lens should be discovered in the 6dF survey, but that as many as ten could be found if quasars typically have B-J - K similar or equal to 8. In this model there could be up to fifty lensed quasars in the, sample, but most of them could only be detected by infrared spectroscopy.

The Las Campanas/AAT Rich Cluster Survey - II. The environmental dependence of galaxy colours in clusters at z similar to 0.1

We present a photometric investigation of the variation in galaxy colour with environment in 11 X-ray-luminous clusters at 0.07 less than or equal to z less than or equal to 0.16 taken from the Las Campanas/AAT Rich Cluster Survey. We study the properties of the galaxy populations in individual clusters, and take advantage of the homogeneity of the sample to combine the clusters together to investigate weaker trends in the composite sample. We find that modal colours of galaxies lying on the colour-magnitude relation in the clusters become bluer by d(B - R)/dr(p) = -0.022 +/- 0.004 from the cluster core out to a projected radius of r(p) = 6 Mpc, further out in radius than any previous study. We also examine the variation in modal galaxy colour with local galaxy density, 2, for galaxies lying close to the colour-magnitude relation, and find that the median colour shifts bluewards by d(B - R)/d log(10)(Sigma) = -0.076 +/- 0.009 with decreasing local density across three orders of magnitude. We show that the position of the red envelope of galaxies in the colour-magnitude relation does not vary as a function of projected radius or density within the clusters, suggesting that the change in the modal colour results from an increasing fraction of bluer galaxies within the colour-magnitude relation, rather than a change in the colours of the whole population. We show that this shift in the colour-magnitude relations with projected radius and local density is greater than that expected from the changing morphological mix based on the local morphology-density relation. We therefore conclude that we are seeing a real change in the properties of galaxies on the colour-magnitude relation in the outskirts of clusters. The simplest interpretation of this result (and similar constraints in local clusters) is that an increasing fraction of galaxies in the lower density regions at large radii within clusters exhibit signatures of star formation in the recent past, signatures which are not seen in the evolved galaxies in the highest density regions.