Finding outlier light curves in catalogues of periodic variable stars

We present a methodology to discover outliers in catalogues of periodic light curves. We use a cross-correlation as the measure of 'similarity' between two individual light curves, and then classify light curves with lowest average 'similarity' as outliers. We performed the analysis on catalogues of periodic variable stars of known type from the MACHO and OGLE projects. This analysis was carried out in Fourier space and we established that our method correctly identifies light curves that do not belong to those catalogues as outliers. We show how an approximation to this method, carried out in real space, can scale to large data sets that will be available in the near future such as those anticipated from the Panoramic Survey Telescope & Rapid Response System (Pan-STARRS) and Large Synoptic Survey Telescope (LSST).     

Photometric redshifts for weak lensing tomography from space: the role of optical and near infrared photometry

We study in detail the photometric redshift requirements needed for tomographic weak gravitational lensing in order to measure accurately the dark energy equation of state. In particular, we examine how ground-based photometry  ( u , g , r , i , z , y )  can be complemented by space-based near-infrared (near-IR) photometry  ( J , H )  , e.g. onboard the planned DUNE satellite. Using realistic photometric redshift simulations and an artificial neural network photo- z method we evaluate the figure of merit for the dark energy parameters  ( w ₀, w _a )  . We consider a DUNE -like broad optical filter supplemented with ground-based multiband optical data from surveys like the Dark Energy Survey, Pan-STARRS and LSST. We show that the dark energy figure of merit would be improved by a factor of 1.3–1.7 if IR filters are added onboard DUNE . Furthermore we show that with IR data catastrophic photo- z outliers can be removed effectively. There is an interplay between the choice of filters, the magnitude limits and the removal of outliers. We draw attention to the dependence of the results on the galaxy formation scenarios encoded into the mock galaxies, e.g. the galaxy reddening. For example, very deep u -band data could be as effective as the IR. We also find that about  10⁵–10⁶  spectroscopic redshifts are needed for calibration of the full survey.     

The Durham/UKST Galaxy Redshift Survey – III. Large-scale structure via the two-point correlation function

We have investigated the statistical clustering properties of galaxies by calculating the two-point galaxy correlation function from the optically selected Durham/UKST Galaxy Redshift Survey. This survey is magnitude-limited to b _J∼17, contains ∼2500 galaxies sampled at a rate of one-in-three and surveys a ∼4×10⁶ ( h ⁻¹ Mpc)³ volume of space. We have empirically determined the optimal method of estimating the two-point correlation function from just such a magnitude-limited survey. Applying our methods to this survey, we find that our redshift-space results agree well with those from previous optical surveys. In particular, we confirm the previously claimed detections of large-scale power out to ∼40 h ⁻¹ Mpc scales. We compare with two common models of cosmological structure formation and find that our two-point correlation function has power significantly in excess of the standard cold dark matter model in the 10–30 h ⁻¹ Mpc region. We therefore support the observational results of the APM galaxy survey. Given that only the redshift-space clustering can be measured directly, we use standard modelling methods and indirectly estimate the real-space two-point correlation function from the projected two-point correlation function. We then invert this projected correlation function to obtain an estimate of the spatial two-point correlation function in real space. This correlation function in real space has a lower amplitude than that in redshift space, but a steeper slope.     

Wisps and knots in the central Crab nebula

The fast-spinning Crab pulsar (∼30 turn s⁻¹), which powers the massive expansion and synchrotron emission of the entire Crab nebula, is surrounded by quasi-stationary features such as fibrous arc-like wisps and bright polar knots in the radial range of 2×10¹⁶≲ r ≲2×10¹⁷ cm, as revealed by high-resolution (∼0.1 arcsec) images from the Wide Field and Planetary Camera 2 (WFPC2) on board the Hubble Space Telescope ( HST ). The spin-down energy flux (∼5×10³⁸ erg s⁻¹) from the pulsar to the luminous outer nebula, which occupies the radial range 0.1≲ r ≲2 pc, is generally believed to be transported by a magnetized relativistic outflow of an electron–positron e^± pair plasma. It is then puzzling that mysterious structures like wisps and knots, although intrinsically dynamic in synchrotron emission, remain quasi-stationary on time-scales of a few days to a week in the relativistic pulsar wind. Here we demonstrate that, as a result of slightly inhomogeneous wind streams emanating from the rotating pulsar, fast magnetohydrodynamic (MHD) shock waves are expected to appear in the pulsar wind at relevant radial distances in the forms of wisps and knots. While forward fast MHD shocks move outward with a speed close to the speed of light c , reverse fast MHD shocks may appear quasi-stationary in space under appropriate conditions. In addition, Alfvénic fluctuations in the shocked magnetized pulsar wind can effectively scatter synchrotron beams from gyrating relativistic electrons and positrons.     

A circular statistical method for extracting rotation measures

We propose a new method for the extraction of rotation measures from spectral polarization data. The method is based on maximum likelihood analysis and takes into account the circular nature of the polarization data. The method is unbiased and statistically more efficient than the standard χ ² procedure. We also find that the method is computationally much more convenient than the standard χ ² procedure if the number of data points is very large. We find that for most sources the method gives results very close to the standard χ ² minimization procedure. We give results for all the cases for which the method gives significantly different results from χ ² minimization. We also make a λ ³ fit to the data in order to extract non-Faraday rotation behaviour for those sources for which a large number of data points are available.     

Fitting the spectrum of the X-ray background: the effects ofhigh-metallicity absorption

Recent work by Risaliti, Maiolino & Salvati suggests that more than half of all Seyfert 2 galaxies in the local Universe are Compton-thick ( N _H>10²⁴ cm⁻²). This has implications for AGN synthesis models for the X-ray background, the flexibility of which for the inclusion of large numbers of high- z type 2 sources we examine here. We highlight the importance of Compton down-scattering in determining the individual source spectra and the fit to the X-ray background spectrum, and demonstrate how parameter space 'opens up' considerably if a super-solar iron abundance is assumed for the absorbing material. This is illustrated with a model which satisfies the present constraints, but which predicts substantial numbers of type 2 sources at the faint flux levels soon to be probed for the first time by the Chandra and XMM missions. We demonstrate also how a strong negative K -correction facilitates the detection of sources with 10^∼24 N _H10²⁵ cm⁻² out to the highest redshifts at which they could plausibly exist.     

Radio observations of comet 9P/Tempel 1 with the Australia Telescope facilities during the Deep Impact encounter

We present radio observations of comet 9P/Tempel 1 associated with the Deep Impact spacecraft collision of 2005 July 4. Weak 18-cm OH emission was detected with the Parkes 64-m telescope, in data averaged over July 4–6, at a level of  12 ± 3 mJy km s⁻¹  , corresponding to OH production rate  2.8 × 10²⁸  molecules s⁻¹ (Despois et al. inversion model, or  1.0 × 10²⁸ s⁻¹  for the Schleicher & A'Hearn model). We did not detect the HCN 1–0 line with the Mopra 22-m telescope over the period July 2–6. The 3σ limit of 0.06 K km s⁻¹ for HCN on July 4 after the impact gives the limit to the HCN production rate of  <1.8 × 10²⁵ s⁻¹  . We did not detect the HCN 1–0 line, 6.7 GHz CH₃OH line or 3.4-mm continuum with the Australia Telescope Compact Array (ATCA) on July 4, giving further limits on any small-scale structure due to an outburst. The 3σ limit on HCN emission of 2.5 K km s⁻¹ from the ATCA around impact corresponds to limit < 4 × 10²⁹ HCN molecules released by the impact.     

Upper limits on the hot Jupiter fraction in the field of NGC 7789

We describe a method of estimating the abundance of short-period extra-solar planets based on the results of a photometric survey for planetary transits. We apply the method to a 21-night survey with the 2.5-m Isaac Newton Telescope of ∼32 000 stars in a ∼0.5 × 0.5 deg² field including the open cluster NGC 7789. From the colour–magnitude diagram, we estimate the mass and radius of each star by comparison with the cluster main sequence. We search for injected synthetic transits throughout the light curve of each star in order to determine their recovery rate, and thus calculate the expected number of transit detections and false alarms in the survey. We take proper account of the photometric accuracy, time sampling of the observations and criteria (signal-to-noise ratio and number of transits) adopted for transit detection. Assuming that none of the transit candidates found in the survey will be confirmed as real planets, we place conservative upper limits on the abundance of planets as a function of planet radius, orbital period and spectral type.     

A search for propylene oxide and glycine in Sagittarius B2 (LMH) and Orion

We have used the Mopra Telescope to search for glycine and the simple chiral molecule propylene oxide in the Sgr B2 (LMH) and Orion KL, in the 3-mm band. We have not detected either species, but have been able to put sensitive upper limits on the abundances of both molecules. The 3σ upper limits derived for glycine conformer I are  3.7 × 10¹⁴ cm⁻²  in both Orion-KL and Sgr B2 (LMH), comparable to the reported detections of conformer I by Kuan et al. However, as our values are 3σ upper limits rather than detections we conclude that this weighs against confirming the detection of Kuan et al. We find upper limits for the glycine II column density of  7.7 × 10¹² cm⁻²  in both Orion-KL and Sgr B2 (LMH), in agreement with the results of Combes et al. The results presented here show that glycine conformer II is not present in the extended gas at the levels detected by Kuan et al. for conformer I. Our ATCA results have ruled out the detection of glycine (both conformers I and II) in the compact hot core of the LMH at the levels reported, so we conclude that it is unlikely that Kuan et al. have detected glycine in either Sgr B2 or Orion-KL. We find upper limits for propylene oxide abundance of  3.0 × 10¹⁴ cm⁻²  in Orion-KL and  6.7 × 10¹⁴ cm⁻²  in Sgr B2 (LMH). We have detected fourteen features in Sgr B2 and four features in Orion-KL which have not previously been reported in the interstellar medium, but have not been able to plausibly assign these transitions to any carrier.     

Deep echelle spectrophotometry of S 311, a Galactic H ii region located outside the solar circle

We present echelle spectrophotometry of the Galactic H  ii region S 311. The data have been taken with the Very Large Telescope Ultraviolet-Visual Echelle Spectrograph in the 3100–10 400 Å range. We have measured the intensities of 263 emission lines; 178 are permitted lines of H⁰, D⁰ (deuterium), He⁰, C⁰, C⁺, N⁰, N⁺, O⁰, O⁺, S⁺, Si⁰, Si⁺, Ar⁰ and Fe⁰; some of them are produced by recombination and others mainly by fluorescence. Physical conditions have been derived using different continuum- and line-intensity ratios. We have derived He⁺, C⁺⁺ and O⁺⁺ ionic abundances from pure recombination lines as well as abundances from collisionally excited lines for a large number of ions of different elements. We have obtained consistent estimations of t ² applying different methods. We have found that the temperature fluctuations paradigm is consistent with the T _e(He  i ) versus T _e(H  i ) relation for H  ii regions, in contrast with what has been found for planetary nebulae. We report the detection of deuterium Balmer lines up to Dδ in the blue wings of the hydrogen lines, whose excitation mechanism seems to be continuum fluorescence.     

Chandra unveils a binary active galactic nucleus in Mrk 463

We analyse Chandra , XMM–Newton and Hubble Space Telescope ( HST ) data of the double-nucleus Ultraluminous Infrared Galaxy (ULIRG), Mrk 463. The Chandra detection of two luminous  ( L _2–10 keV= 1.5 × 10⁴³  and  3.8 × 10⁴² erg cm⁻² s⁻¹)  , unresolved nuclei in Mrk 463 indicates that this galaxy hosts a binary active galactic nucleus (AGN), with a projected separation of ≃3.8 kpc (  3.83 ± 0.01  arcsec). While the East nucleus was already known to be a type 2 Seyfert (and this is further confirmed by our Chandra detection of a neutral iron line), this is the first unambiguous evidence in favour of the AGN nature of the West nucleus. Mrk 463 is therefore the clearest case so far for a binary AGN, after NGC 6240.     

Near-infrared/optical counterparts of hotspots in radio galaxies

We present new high spatial resolution Very Large Telescope (VLT) and Very Large Array (VLA) observations of a sample of nine low-power (   P _1.4 GHz≤ 10²⁵  W Hz⁻¹) radio hotspots. Infrared/optical emission is definitely detected in four of the nine observed objects, resulting in a detection rate of at least 45 per cent. This emission is interpreted as synchrotron radiation from the electrons accelerated in the hotspots. The integrated spectra of these hotspots reveal typical break frequencies between 10⁵ and 10⁶ GHz, two orders of magnitude higher than typically found in high-power hotspots. This supports the idea that in low-power hotspots with their relatively low magnetic field strengths, electrons emit most of their energy at higher frequencies. A simple spectral ageing analysis would imply that the emitting electrons have been injected into the hotspot volume less than  ∼10³  years ago. We discuss possible scenarios to explain the lack of older electrons in the hotspot region. In particular, the extended morphology of the near-infrared/optical emission would suggest that efficient re-acceleration mechanisms rejuvenate the electron populations.     

CO abundances in a protostellar cloud: freeze-out and desorption in the envelope and outflow of L483

CO isotopes are able to probe the different components in protostellar clouds. These components, core, envelope and outflow have distinct physical conditions, and sometimes more than one component contributes to the observed line profile. In this study, we determine how CO isotope abundances are altered by the physical conditions in the different components. We use a 3D molecular line transport code to simulate the emission of four CO isotopomers, ¹²CO   J = 2 → 1, ¹³CO J = 2 → 1  , C¹⁸O   J = 2 → 1  and C¹⁷O   J = 2 → 1  from the Class 0/1 object L483, which contains a cold quiescent core, an infalling envelope and a clear outflow. Our models replicate James Clerk Maxwell Telescope (JCMT) line observations with the inclusion of freeze-out, a density profile and infall. Our model profiles of ¹²CO and ¹³CO have a large linewidth due to a high-velocity jet. These profiles replicate the process of more abundant material being susceptible to a jet. C¹⁸O and C¹⁷O do not display such a large linewidth as they trace denser quiescent material deep in the cloud.     

Two serendipitous low-mass LMC clusters discovered with HST1

We present V and I photometry of two open clusters in the LMC down to V ∼26. The clusters were imaged with the Wide Field and Planetary Camera 2 (WFPC2) on board the Hubble Space Telescope ( HST ), as part of the Medium Deep Survey Key Project. Both are low-luminosity ( M_V ∼−3.5), low-mass ( M ∼10³ M⊙) systems. The chance discovery of these two clusters in two parallel WFPC2 fields suggests a significant incompleteness in the LMC cluster census near the bar. One of the clusters is roughly elliptical and compact, with a steep light profile, a central surface brightness μ _V (0)∼20.2 mag arcsec⁻², a half-light radius r _hl∼0.9 pc (total visual major diameter D ∼3 pc) and an estimated mass M ∼1500 M⊙. From the colour–magnitude diagram and isochrone fits we estimate its age as τ∼(2–5)×10⁸ yr. Its mass function has a fitted slope of Γ=Δlogφ( M )/Δlog M =−1.8±0.7 in the range probed (0.9≲ M /M⊙≲4.5). The other cluster is more irregular and sparse, having shallower density and surface brightness profiles. We obtain Γ=−1.2±0.4, and estimate its mass as M ∼400 M⊙. A derived upper limit for its age is τ≲5×10⁸ yr. Both clusters have mass functions with slopes similar to that of R136, a massive LMC cluster, for which HST results indicate Γ∼−1.2. They also seem to be relaxed in their cores and well contained in their tidal radii.     

The European Large Area ISO Survey – IV. The preliminary 90-μm luminosity function

We present the luminosity function of 90-μm-selected galaxies from the European Large Area ISO Survey (ELAIS), extending to z =0.3. Their luminosities are in the range 10⁹65⁻² L /L_⊙<10¹², i.e. non-ultraluminous. From our sample of 37 reliably detected galaxies in the ELAIS S1 region from the Efstathiou et al. S ₉₀100 mJy data base, we have found optical, 15-μm or 1.4-GHz identifications for 24 (65 per cent). We have obtained 2dF and UK Schmidt FLAIR spectroscopy of 89 per cent of identifications to rigid multivariate flux limits. We construct a luminosity function assuming that (i) our spectroscopic subset is an unbiased sparse sample, and (ii) there are no galaxies that would not be represented in our spectroscopic sample at any redshift. We argue that we can be confident of both assumptions. We find that the luminosity function is well described by the local 100-μm luminosity function of Rowan-Robinson, Helou & Walker. Assuming this local normalization, we derive luminosity evolution of (1+ z )^2.45±0.85 (95 per cent confidence). We argue that star formation dominates the bolometric luminosities of these galaxies, and we derive comoving star formation rates in broad agreement with the Flores et al. and Rowan-Robinson et al. mid-infrared-based estimates.     

Interstellar CO towards X Persei (HD 24534) – II. Two-component model

The observations made by the Goddard High Resolution Spectrograph (GHRS) aboard the Hubble Space Telescope ( HST ) of molecular CO in absorbing gas towards X Persei are reported. The two-component statistical equilibrium model incorporating radiative excitation of CO by line emission at the same velocity that originates in nearby molecular clouds has been used to reproduce high-resolution GHRS spectra. Earlier analysis indicates that the cloud has a complex structure and at least a two-component model should be used to obtain accurate results. The spectra obtained from the International Ultraviolet Explorer ( IUE ) were used to complement GHRS data and constrain the space of possible solutions. The new oscillator strengths recommended by Eidelsberg et al. for A–X bands have been used. The results show that one of the components may be attributed to the Perseus OB2 molecular cloud, and the other component to an extension of the Taurus dark cloud. The total CO column density N (CO)=(1.0±0.2)×10¹⁶ cm⁻² has been determined. According to the results about 85 per cent of the observed CO belongs to an extension of the Taurus dark cloud. The CO radiation that originates in nearby molecular clouds may be the dominant excitation mechanism of the observed CO. The early results of ¹³CO line analysis indicate a ¹³CO/¹²CO ratio of about 40.     

H 0 from an orientation-unbiased sample of Sunyaev–Zel'dovich and X-ray clusters

We have observed the Sunyaev–Zel'dovich (SZ) effect in a sample of five moderate-redshift clusters with the Ryle Telescope, and used them in conjunction with X-ray imaging and spectral data from ROSAT and ASCA to measure the Hubble constant. This sample was chosen with a strict X-ray flux limit using both the Bright Cluster Sample and the Northern ROSAT All-Sky Survey (RASS) cluster catalogues to be well above the surface brightness limit of the RASS, and hence to be unbiased with respect to the orientation of the cluster. This controls a major potential systematic effect in the SZ/X-ray method of measuring H ₀. Taking the weighted geometric mean of the results and including the main sources of error, namely the noise in the SZ measurement, the uncertainty in the X-ray temperatures and the unknown ellipticity and substructure of the clusters, we find   H ₀= 59⁺¹⁰₋₉ (random)⁺⁸₋₇(systematic) km s⁻¹ Mpc⁻¹  assuming a standard cold dark matter model with  Ω_M= 1.0, Ω_Λ= 0.0  or   H ₀= 66⁺¹¹₋₁₀ ⁺⁹₋₈ km  s⁻¹ Mpc⁻¹  if  Ω_M= 0.3, Ω_Λ= 0.7  .     

On the formation of massive galaxies: a simultaneous study of number density, size and intrinsic colour evolution in GOODS

The evolution of number density, size and intrinsic colour is determined for a volume-limited sample of visually classified early-type galaxies selected from the Hubble Space Telescope /Advanced Camera for Surveys images of the Great Observatories Origins Deep Survey (GOODS) North and South fields (version 2). The sample comprises 457 galaxies over 320 arcmin² with stellar masses above  3 × 10¹⁰ M_⊙  in the redshift range  0.4 < z < 1.2  . Our data allow a simultaneous study of number density, intrinsic colour distribution and size. We find that the most massive systems  (≳3 × 10¹¹ M_⊙)  do not show any appreciable change in comoving number density or size in our data. Furthermore, when including the results from 2dF galaxy redshift survey, we find that the number density of massive early-type galaxies is consistent with no evolution between   z = 1.2  and 0, i.e. over an epoch spanning more than half of the current age of the Universe. We find large discrepancies between the predictions of semi-analytic models. Massive galaxies show very homogeneous intrinsic colour distributions, with nearly flat radial colour gradients, but with a significant negative correlation between stellar mass and colour gradient, such that red cores appear predominantly in massive galaxies. The distribution of half-light radii – when compared to   z ∼ 0  and   z > 1  samples – is compatible with the predictions of semi-analytic models relating size evolution to the amount of dissipation during major mergers.     

GeV emission from neutron-rich internal shocks of some long γ-ray bursts

In the neutron-rich internal shocks model for γ-ray bursts (GRBs), the Lorentz factors (LFs) of ion shells are variable, and so are the LFs of accompanying neutron shells. For slow neutron shells with a typical LF of approximate tens, the typical β-decay radius is  ∼10¹⁴–10¹⁵ cm  . As GRBs last long enough  [ T ₉₀ > 14(1 + z ) s]  , one earlier but slower ejected neutron shell will be swept successively by later ejected ion shells in the range  ∼10¹³–10¹⁵ cm  , where slow neutrons have decayed significantly. Part of the thermal energy released in the interaction will be given to the electrons. These accelerated electrons will mainly be cooled by the prompt soft γ-rays and give rise to GeV emission. This kind of GeV emission is particularly important for some very long GRBs and is detectable for the upcoming satellite Gamma-Ray Large Area Space Telescope (GLAST).     

The Millennium Galaxy Catalogue: 16 ≤BMGC < 24 galaxy counts and the calibration of the local galaxy luminosity function

The Millennium Galaxy Catalogue (MGC) is a 37.5 deg², medium-deep, B -band imaging survey along the celestial equator, taken with the Wide Field Camera on the Isaac Newton Telescope. The survey region is contained within the regions of both the Two Degree Field Galaxy Redshift Survey (2dFGRS) and the Sloan Digital Sky Survey Early Data Release (SDSS-EDR). The survey has a uniform isophotal detection limit of 26 mag arcsec⁻² and it provides a robust, well-defined catalogue of stars and galaxies in the range  16 ≤ B _MGC < 24 mag  .
Here we describe the survey strategy, the photometric and astrometric calibration, source detection and analysis, and present the galaxy number counts that connect the bright and faint galaxy populations within a single survey. We argue that these counts represent the state of the art and use them to constrain the normalizations (φ*) of a number of recent estimates of the local galaxy luminosity function. We find that the 2dFGRS, SDSS Commissioning Data (CD), ESO Slice Project, Century Survey, Durham/UKST, Mt Stromlo/APM, SSRS2 and NOG luminosity functions require a revision of their published φ* values by factors of  1.05 ± 0.05, 0.76 ± 0.10, 1.02 ± 0.22, 1.02 ± 0.16, 1.16 ± 0.28, 1.75 ± 0.37, 1.40 ± 0.26  and  1.01 ± 0.39  , respectively. After renormalizing the galaxy luminosity functions we find a mean local b _J luminosity density of     . ¹