The SuperCOSMOS Sky Survey – II. Image detection, parametrization, classification and photometry

In this, the second in a series of three papers concerning the SuperCOSMOS Sky Survey, we describe the methods for image detection, parametrization, classification and photometry. We demonstrate the internal and external accuracy of our object parameters. Using examples from the first release of data, the South Galactic Cap survey, we show that our image detection completeness is close to 100 per cent to within ∼1.5 mag of the nominal plate limits. We show that for the B _J survey data, the image classification is externally > 99 per cent reliable to B _J∼19.5 . Internally, the image classification is reliable at a level of > 90 per cent to B _J∼21 , R ∼19 . The photometric accuracy of our data is typically ∼0.3 mag with respect to external data for m >15 . Internally, the relative photometric accuracy in restricted position and magnitude ranges can be as accurate as ∼5 per cent for well-exposed stellar images. Colours are externally accurate to σ _{B − R , R − I} ∼0.07 at B _J∼16.5 , rising to σ _{B − R , R − I} ∼0.16 at B _J∼20 .     

The asymmetric structure of the Galactic halo

Using the stellar photometry catalogue based on the latest data release (DR4) of the Sloan Digital Sky Survey (SDSS), a study of the Galactic structure using star counts is carried out for selected areas of the sky. The sample areas are selected along a circle at a Galactic latitude of +60°, and 10 strips of high Galactic latitude along different longitudes. Direct statistics of the data show that the surface densities of ℓ from 180° to 360° are systematically higher than those of ℓ from 0° to 180°, defining a region of overdensity (in the direction of Virgo) and another one of underdensity (in the direction of Ursa Major) with respect to an axisymmetric model. It is shown by comparing the results from star counts in the ( g − r ) colour that the density deviations are due to an asymmetry of the stellar density in the halo. Theoretical models for the surface density profile are built and star counts are performed using a triaxial halo of which the parameters are constrained by observational data. Two possible reasons for the asymmetric structure are discussed.     

The complex iron line in NGC 7469 observed by BeppoSAX

We compare the large-scale galaxy clustering in the new Sloan Digital Sky Survey (SDSS) early data release (EDR) with the clustering in the APM Galaxy Survey. We cut out pixel maps (identical in size and shape) from the SDSS and APM data to allow a direct comparison of the clustering. Here we concentrate our analysis on an equatorial SDSS strip in the South Galactic Cap (EDR/SGC) of 166 deg², 25 wide and 65° long . Only galaxies with Petrosian magnitudes  16.8< g '<19.8  are included to match the surface density of the  17< b _J<20  APM pixel maps (median depth of ∼400  h ⁻¹ Mpc). Both the amplitude and the shape of the angular two-point function and variance turn out to be in very good agreement with the APM on scales smaller than 2° (or ≲15  h ⁻¹ Mpc). The three-point function and skewness are also in excellent agreement within a 90 per cent confidence level. On the one hand these results illustrate that the EDR data and SDSS software pipelines work well and are suitable to carry out analysis of large-scale clustering. On the other hand they confirm that large-scale clustering analysis is becoming 'repeatable' and therefore that our conclusions for structure formation models seem to stand on solid scientific grounds.     

Kinematics of SDSS subdwarfs: structure and substructure of the Milky Way halo

We construct a new sample of ∼1700 solar neighbourhood halo subdwarfs from the Sloan Digital Sky Survey (SDSS), selected using a reduced proper-motion diagram. Radial velocities come from the SDSS spectra and proper motions from the light-motion curve catalogue of Bramich et al. Using a photometric parallax relation to estimate distances gives us the full phase-space coordinates. Typical velocity errors are in the range  30–50 km s⁻¹  . This halo sample is one of the largest constructed to date and the disc contamination is at a level of ≲1 per cent. This enables us to calculate the halo velocity dispersion to excellent accuracy. We find that the velocity dispersion tensor is aligned in spherical polar coordinates and that  (σ _r , σ_φ, σ_θ) = (143 ± 2, 82 ± 2, 77 ± 2) km s⁻¹  . The stellar halo exhibits no net rotation, although the distribution of   v _φ  shows tentative evidence for asymmetry. The kinematics are consistent with a mildly flattened stellar density falling with distance like   r ^−3.75  .
Using the full phase-space coordinates, we look for signs of kinematic substructure in the stellar halo. We find evidence for four discrete overdensities localized in angular momentum and suggest that they may be possible accretion remnants. The most prominent is the solar neighbourhood stream previously identified by Helmi et al., but the remaining three are new. One of these overdensities is potentially associated with a group of four globular clusters (NGC 5466, NGC 6934, M2 and M13) and raises the possibility that these could have been accreted as part of a much larger progenitor.     

Constraints on the mass and abundance of black holes in the Galactic halo: the high-mass limit

We establish constraints on the mass and abundance of black holes in the Galactic halo by determining their impact on globular clusters, which are conventionally considered to be little evolved. Using detailed Monte Carlo simulations and simple evolutionary models, we argue that black holes with masses M _bh≳(1–3)×10⁶ M_⊙ can comprise no more than a fraction f _bh≈0.17 of the total halo density at Galactocentric radius R ≈8 kpc. This bound arises from requiring stability of the cluster mass function. A more restrictive bound may be derived if we demand that the probability of destruction of any given, low-mass M _c≈(2.5–7.5)×10⁴ M_⊙] globular cluster not exceed 50 per cent; this bound is f _bh≲0.025–0.05 at R ≈8 kpc. This constraint improves those based on disc heating and dynamical friction arguments as well as current lensing results. At smaller radius the constraint on f _bh strengthens, while at larger radius an increased fraction of black holes is allowed.     

Inversion of stellar statistics equation for the Galactic bulge

A method based on Lucy's iterative algorithm is developed to invert the equation of stellar statistics for the Galactic bulge and is then applied to the K -band star counts from the Two-Micron Galactic Survey in a number of off-plane regions (10°>| b |>2°, | l |<15°). The top end of the K -band luminosity function is derived and the morphology of the stellar density function is fitted to triaxial ellipsoids, assuming a non-variable luminosity function within the bulge. The results, which have already been outlined by López-Corredoira et al., are shown in this paper with a full explanation of the steps of the inversion: the luminosity function shows a sharp decrease brighter than M _K =−8.0  mag when compared with the disc population; the bulge fits triaxial ellipsoids with the major axis in the Galactic plane at an angle with the line of sight to the Galactic centre of 12° in the first quadrant; the axial ratios are 1:0.54:0.33, and the distance of the Sun from the centre of the triaxial ellipsoid is 7860 pc. The major–minor axial ratio of the ellipsoids is found not to be constant, the best fit to the gradient being K _z =(8.4±1.7)×exp(− t /(2000±920) pc), where t is the distance along the major axis of the ellipsoid in parsecs. However, the interpretation of this is controversial. An eccentricity of the true density-ellipsoid gradient and a population gradient are two possible explanations. The best fit for the stellar density, for 1300 pc< t <3000 pc, is calculated for both cases, assuming an ellipsoidal distribution with constant axial ratios, and when K _z is allowed to vary. From these, the total number of bulge stars is ∼3×10¹⁰ or ∼4×10¹⁰, respectively.     

Deep galaxy counts, extragalactic background light and the stellar baryon budget

We assess the constraints imposed by the observed extragalactic background light (EBL) on the cosmic history of star formation and the stellar-mass density today. The logarithmic slope of the galaxy number–magnitude relation from the Southern Hubble Deep Field imaging survey is flatter than 0.4 in all seven UBVIJHK optical bandpasses, i.e. the light from resolved galaxies has converged from the UV to the near-IR. We find a lower limit to the surface brightness of the optical extragalactic sky of about 15 nW m⁻² sr⁻¹, comparable to the intensity of the far-IR background from COBE data. Assuming a Salpeter initial mass function with a lower cut-off consistent with observations of M subdwarf disc stars, we set a lower limit of Ω_g+s h ²>0.0013  I ₅₀ to the visible (processed gas + stars) mass density required to generate an EBL at a level of 50  I ₅₀ nW m⁻² sr⁻¹; our 'best-guess' value is Ω_g+s h ²≈0.0031  I ₅₀. Motivated by the recent microlensing results of the MACHO collaboration, we consider the possibility that massive dark haloes around spiral galaxies are composed of faint white dwarfs, and show that only a small fraction (≲5 per cent) of the nucleosynthetic baryons can be locked in the remnants of intermediate-mass stars forming at z _F≲5, as the bright early phases of such haloes would otherwise overproduce the observed EBL.     

Mass models of the Milky Way

A parametrized model of the mass distribution within the Milky Way is fitted to the available observational constraints. The most important single parameter is the ratio of the scalelength R _d* of the stellar disc to R ₀. The disc and bulge dominate v _c( R ) at R ≲ R ₀ only for R _d,*/ R ₀≲0.3. Since the only knowledge we have of the halo derives from studies like the present one, we allow it to contribute to the density at all radii. When allowed this freedom, however, the halo causes changes in assumptions relating to R  ≪  R ₀ to affect profoundly the structure of the best-fitting model at R  ≫  R ₀. For example, changing the disc slightly from an exponential surface-density profile significantly changes the form of v _c( R ) at R  ≫  R ₀, where the disc makes a negligible contribution to v _c. Moreover, minor changes in the constraints can cause the halo to develop a deep hole at its centre that is not physically plausible. These problems call into question the proposition that flat rotation curves arise because galaxies have physically distinct haloes rather than outwards-increasing mass-to-light ratios.   The mass distribution of the Galaxy and the relative importance of its various components will remain very uncertain until more observational data can be used to constrain mass models. Data that constrain the Galactic force field at z ≳ R and at R  >  R ₀ are especially important.     

Contribution to the search for binaries among Am stars – VIII. New spectroscopic orbits of eight systems and statistical study of a sample of 91 Am stars

Red clump giant (RCG) stars can be used as distance indicators to trace the mass distribution of the Galactic bar. We use RCG stars from 44 bulge fields from the OGLE-II microlensing collaboration data base to constrain analytic triaxial models for the Galactic bar. We find the bar major-axis is oriented at an angle of 24°–27° to the Sun–Galactic Centre line-of-sight. The ratio of semimajor and semiminor bar axis scalelengths in the Galactic plane   x ₀, y ₀  , and vertical bar scalelength z ₀, is   x ₀ :  y ₀ :  z ₀= 10 : 3.5 : 2.6  , suggesting a slightly more prolate bar structure than the working model of Gerhard which gives the scalelength ratios as   x ₀ :  y ₀ :  z ₀= 10 : 4 : 3  .     

The 6dF Galaxy Survey: final redshift release (DR3) and southern large-scale structures

We report the final redshift release of the 6dF Galaxy Survey (6dFGS), a combined redshift and peculiar velocity survey over the southern sky  (| b | > 10°)  . Its 136 304 spectra have yielded 110 256 new extragalactic redshifts and a new catalogue of 125 071 galaxies making near-complete samples with  ( K , H , J , r _F, b _J) ≤ (12.65, 12.95, 13.75, 15.60, 16.75)  . The median redshift of the survey is 0.053. Survey data, including images, spectra, photometry and redshifts, are available through an online data base. We describe changes to the information in the data base since earlier interim data releases. Future releases will include velocity dispersions, distances and peculiar velocities for the brightest early-type galaxies, comprising about 10 per cent of the sample. Here we provide redshift maps of the southern local Universe with   z ≤ 0.1  , showing nearby large-scale structures in hitherto unseen detail. A number of regions known previously to have a paucity of galaxies are confirmed as significantly underdense regions. The URL of the 6dFGS data base is http://www-wfau.roe.ac.uk/6dFGS .     

Near-infrared star counts as a test of a Galactic bar structure

We present survey data in the narrow-band L filter (nb L ), taken at UKIRT, for a total area of 277 arcmin², roughly equally divided between four regions at zero Galactic latitude and longitudes ±4.°3 and ±2.°3. The 80 per cent completeness level for these observations is at roughly magnitude 11.0. This magnitude limit, owing to the low coefficient for interstellar extinction at this wavelength ( A _{nb L} =0.047 A_V ), allows us to observe bulge giants. We match the nb L magnitudes with DENIS survey K magnitudes, and find 95 per cent of nb L sources are matched to K sources. Constructing colour–magnitude diagrams, we deredden the magnitudes and find evidence for a longitude-dependent asymmetry in the source counts. We find that there are ∼15 per cent and ∼5 per cent more sources at the negative longitude than at the corresponding positive longitude for the fields at ±4.°3 and ±2.°3, respectively. This is compared with the predictions of some Galactic bar models. We find an asymmetry in the expected sense, which favours gas dynamical models and the recent deconvolution of surface photometry data over earlier treatments of photometric data.     

On the intrinsic bias in detecting caustic crossings between Galactic halo and self-lensing events in the Magellanic Clouds

In this paper, we investigate the intrinsic bias in detecting caustic crossings between the Galactic halo and self-lensing gravitational microlensing events in the Magellanic Clouds. For this, we determine the region for optimal caustic-crossing detection in the parameter space of the physical binary separations, ℓ, and the total binary lens mass, M , and find that the optimal regions for both populations of events are similar to each other. In particular, if the Galactic halo is composed of lenses with the claimed average mass of 〈 M 〉∼0.5 M_⊙, the optimal binary separation range of Galactic halo events of 3.5 au≲ℓ≲14 au matches well with that of a Magellanic Cloud self-lensing event caused by a binary lens with a total mass of M ∼1 M_⊙; well within the mass range of the most probable lens population of stars in the Magellanic Clouds. Therefore, our computation implies that if the binary fractions and the distributions of binary separations of the two populations of lenses are not significantly different from each other, there is no strong detection bias against Galactic halo caustic-crossing events.     

Analysis of the thin layer of Galactic warm ionized gas in the range

We present an analysis of the thin layer of Galactic warm ionized gas at an angular resolution ∼10 arcmin. This is carried out using radio continuum data at 1.4, 2.7 and 5 GHz in the coordinate region     . For this purpose, we evaluate the zero level of the 2.7- and 5-GHz surveys using auxiliary data at 2.3 GHz and 408 MHz. The derived zero-level corrections are   T _zero(2.7 GHz) = 0.15 ± 0.06 K  and   T _zero(5 GHz) = 0.1 ± 0.05 K  . We separate the thermal (free–free) and non-thermal (synchrotron) component by means of a spectral analysis performed adopting an antenna temperature spectral index −2.1 for the free–free emission, a realistic spatial distribution of indices for the synchrotron radiation and by fitting, pixel-by-pixel, the Galactic spectral index. We find that at 5 GHz, for  | b | = 0°  , the fraction of thermal emission reaches a maximum value of 82 per cent, while at 1.4 GHz, the corresponding value is 68 per cent. In addition, for the thermal emission, the analysis indicates a dominant contribution of the diffuse component relative to the source component associated with discrete H  ii regions.     

Maximum-likelihood method for estimating the mass and period distributions of extrasolar planets

We investigate the distribution of mass M and orbital period P of extrasolar planets, taking account of selection effects caused by the limited velocity precision and duration of existing surveys. We fit the data on 72 planets to a power-law distribution of the form  d n = CM ^−α P ^−β(d M / M )(d P / P )  , and find  α= 0.11 ± 0.10  ,  β=−0.27 ± 0.06  for   M ≲ 10 M _J  , where   M _J  is the mass of Jupiter. The correlation coefficient between these two exponents is −0.31, indicating that uncertainties in the two distributions are coupled. We estimate that 4 per cent of solar-type stars have companions in the range  1 M _J < M < 10 M _J  ,  2_d < P < 10 yr  .     

Distant field blue horizontal branch stars – III. Identification of a probable outer halo stream at Galactocentric distance r= 70 kpc

We present FOcal Reducer/low dispersion Spectrograph-1 spectra (from the European Southern Observatory's Very Large Telescope) of a sample of 34 faint  20.0 < g * < 21.1  A-type stars selected from the Sloan Digital Sky Survey Early Data Release, with the goal of measuring the velocity dispersion of blue horizontal branch (BHB) stars in the remote Galactic halo,   R ∼ 80 kpc  . We show that colour selection with  1.08 < u *− g * < 1.40  and  −0.2 < g *− r * < −0.04  minimizes contamination of the sample by less luminous blue stragglers. In classifying the stars we confine our attention to the 20 stars with spectra of signal-to-noise ratio >15 Å⁻¹. Classification produces a sample of eight BHB stars at distances  65–102 kpc  from the Sun (mean 80 kpc), which represents the most distant sample of Galactic stars with measured radial velocities. The dispersion of the measured radial component of the velocity with respect to the centre of the Galaxy is  58 ± 15 km s⁻¹  . This value is anomalously low in comparison with measured values for stars at smaller distances, as well as for satellites at similar distances. Seeking an explanation for the low measured velocity dispersion, further analysis reveals that six of the eight remote BHB stars are plausibly associated with a single orbit. Three previously known outer halo carbon stars also appear to belong to this stream. The velocity dispersion of all nine stars relative to the orbit is only  15 ± 4 km s⁻¹  . Further observations along the orbit are required to trace the full extent of this structure on the sky.     

The shapes of galaxies in the Sloan Digital Sky Survey

We determine the underlying shapes of spiral and elliptical galaxies in the Sloan Digital Sky Survey Data Release 6 (SDSS DR6) from the observed distribution of projected galaxy shapes, taking into account the effects of dust extinction and reddening. We assume that the underlying shapes of spirals and ellipticals are well approximated by triaxial ellipsoids. The elliptical galaxy data are consistent with oblate spheroids, with a correlation between luminosity and ellipticity: the mean values of minor to middle axis ratios are 0.41 ± 0.03 for   M _r ≈−18  ellipticals and 0.76 ± 0.04 for   M _r ≈−22.5  ellipticals. Ellipticals show almost no dependence of axial ratio on galaxy colour, implying a negligible dust optical depth.
There is a strong variation of spiral galaxy shapes with colour indicating the presence of dust. The intrinsic shapes of spiral galaxies in the SDSS DR6 are consistent with flat discs with a mean and dispersion of thickness to diameter ratio of (21 ± 2) per cent, and a face-on ellipticity, e , of  ln( e ) =−2.33 ± 0.79  . Not including the effects of dust in the model leads to discs that are systematically rounder by up to 60 per cent. More luminous spiral galaxies tend to have thicker and rounder discs than lower luminosity spirals. Both elliptical and spiral galaxies tend to be rounder for larger galaxies.
The marginalized value of the edge-on r -band dust extinction E ₀ in spiral galaxies is   E ₀≃ 0.45  mag for galaxies of median colours, increasing to   E ₀= 1  mag for   g − r > 0.9  and   E ₀= 1.9  for the luminous and most compact galaxies, with half-light radii  <2  h ⁻¹ kpc  .     

The Durham/UKST Galaxy Redshift Survey — V. The catalogue

We present the radial velocities and blue, optical magnitudes for all of the galaxies within the Durham/UKST Galaxy Redshift Survey. This catalogue consists of ∼2500 galaxy redshifts to a limiting apparent magnitude of B _J⋍17 mag, covering a ∼1500-deg² area around the South Galactic Pole. The galaxies in this survey were selected from the Edinburgh/Durham Southern Galaxy Catalogue and were sampled, in order of apparent magnitude, at a rate of one galaxy in every three. The spectroscopy was performed at the 1.2-m UK Schmidt Telescope in Australia using the FLAIR multi-object spectrograph. We show that our radial velocity measurements made with this instrument have an empirical accuracy of ±150 km s⁻¹. The observational techniques and data reduction procedures used in the construction of this survey are also discussed. This survey demonstrates that the UKST can be used to make a three-dimensional map of the large-scale galaxy distribution, via a redshift survey to b _J⋍17 mag, over a wide area of the sky.     

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     . ¹     

Extensive serendipitous X-ray coverage of a flare star with ROSAT

We report the serendipitous discovery of a flare star observed with the ROSAT X-ray observatory. From optical spectra, which show strong and variable emission lines of the hydrogen Balmer series and neutral helium, we classify this object as a M3.0Ve star, and estimate a distance of 52 pc from published photometry. Owing to the close proximity of the star (13.6 arcmin) to the calibration source and RS CVn binary AR Lacertae, long-term X-ray coverage is available in the ROSAT archive (∼50 h spanning 6.5 yr). Two large flare events occurred early in the mission (1990 June–July), and the end of a third flare was detected in 1996 June. One flare, observed with the Position Sensitive Proportional Counter (PSPC), had a peak luminosity L _X=1.1×10³⁰ erg s⁻¹, an e-folding rise time of 2.2 h and a decay time of 7 h. This decay time is one of the longest detected on a dMe star, providing evidence for the possibility of additional heating during the decay phase. A large High Resolution Imager (HRI) flare (peak L _X=2.9×10³⁰ erg s⁻¹) is also studied. The 'background' X-ray emission is also variable – evidence for low-level flaring or microflaring. We find that 59 per cent of the HRI counts and 68 per cent of the PSPC counts are caused by flares. At least 41 per cent of the HRI exposure time and 47 per cent of the PSPC are affected by detectable flare enhancement.     

The violent past of Cygnus X-2

Cygnus X-2 appears to be the descendant of an intermediate-mass X-ray binary (IMXB). Using Mazzitelli's stellar code we compute detailed evolutionary sequences for the system and find that its prehistory is sensitive to stellar input parameters, in particular the amount of core overshooting during the main-sequence phase. With standard assumptions for convective overshooting a case B mass transfer starting with a 3.5-M_⊙ donor star is the most likely evolutionary solution for Cygnus X-2. This makes the currently observed state rather short-lived, of order 3 Myr, and requires a formation rate > 10⁻⁷–10⁻⁶ yr⁻¹ of such systems in the Galaxy. Our calculations show that neutron star IMXBs with initially more massive donors (≳4 M_⊙) encounter a delayed dynamical instability; they are unlikely to survive this rapid mass transfer phase. We determine limits for the age and initial parameters of Cygnus X-2 and calculate possible dynamical orbits of the system in a realistic Galactic potential, given its observed radial velocity. We find trajectories which are consistent with a progenitor binary on a circular orbit in the Galactic plane inside the solar circle that received a kick velocity ≤200 km s⁻¹ at the birth of the neutron star. The simulations suggest that about 7 per cent of IMXBs receiving an arbitrary kick velocity from a standard kick velocity spectrum would end up in an orbit similar to Cygnus X-2, while about 10 per cent of them reach yet larger Galactocentric distances.