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.     

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.     

Red giant stellar collisions in the Galactic Centre

We show that collisions with stellar-mass black holes can partially explain the absence of bright giant stars in the Galactic Centre, first noted by Genzel et al. We show that the missing objects are low-mass giants and asymptotic giant branch stars in the range  1–3 M_⊙  . Using detailed stellar evolution calculations, we find that to prevent these objects from evolving to become visible in the depleted K bands, we require that they suffer collisions on the red giant branch, and we calculate the fractional envelope mass losses required. Using a combination of smoothed particle hydrodynamic calculations, restricted three-body analysis and Monte Carlo simulations, we compute the expected collision rates between giants and black holes, and between giants and main-sequence stars in the Galactic Centre. We show that collisions can plausibly explain the missing giants in the  10.5 < K < 12  band. However, depleting the brighter  ( K < 10.5)  objects out to the required radius would require a large population of black hole impactors which would in turn deplete the  10.5 < K < 12  giants in a region much larger than is observed. We conclude that collisions with stellar-mass black holes cannot account for the depletion of the very brightest giants, and we use our results to place limits on the population of stellar-mass black holes in the Galactic Centre.     

Hipparcos and the age of the Galactic disc

We use the Hipparcos colour–magnitude diagram of field stars with Tycho colours to make a new minimum age estimate for the Galactic disc. The method is based on fits to the red envelope of subgiants in the Hipparcos colour–magnitude diagram with synthetic isochrones covering the range of disc metal abundance. The colours and luminosities of the isochrones as a function of abundance are checked using new techniques involving 'red-clump' stars in the giant branch region and on the main sequence using G and K dwarfs. We derive a minimum disc age of 8 Gyr, in good agreement with other methods.     

The structure of the Galactic bar

We present a deep near-infrared wide-angle photometric analysis of the structure of the inner Galactic bar and central disc. The presence of a triaxial structure at the centre of the Galaxy is confirmed, consistent with a bar inclined at  22°± 55  from the Sun—Galactic Centre line, extending to approximately 2.5 kpc from the Galactic Centre and with a rather small axis ratio. A feature at  ℓ=−98  not aligned with this triaxiality suggests the existence of a second structure in the inner Galaxy, a double triaxiality or an inner ring. We argue that this is likely to be the signature of the end of the Galactic bar, at approximately 2.5–3 kpc, which is circumscribed by an inner pseudo-ring. No thick dust lane preceding the bar is detected and a hole in the dust distribution of the disc inside the bar radius is inferred.     

The structure of the Galactic halo: SDSS versus SuperCOSMOS

The halo structure at high Galactic latitudes near both the north and south poles is studied using Sloan Digital Sky Survey (SDSS) and SuperCOSMOS data. For the south cap halo, the archive of the SuperCOSMOS photographic photometry sky survey is used. The coincident source rate between SuperCOSMOS data in B _J band from 16.5 to 20.5 mag and SDSS data is about 92 per cent, in a common sky area in the south. While that in the R _F band is about 85 per cent from 16.5 to 19.5 mag. Transformed to the SuperCOSMOS system and downgraded to the limiting magnitudes of SuperCOSMOS, the star counts in the North Galactic Cap from SDSS show up to an  16.9 ± 6.3  per cent  asymmetric ratio (defined as relative fluctuations over the rotational symmetry structure) in the B _J band, and up to  13.5 ± 6.7  per cent  asymmetric ratio in the R _F band. From SuperCOSMOS B _J and R _F bands, the structure of the Southern Galactic hemisphere does not show the same obvious asymmetric structures as the northern sky does in both the original and downgraded SDSS star counts. An axisymmetric halo model with n = 2.8 and q = 0.7 can fit the projected number density from SuperCOSMOS fairly well, with an average error of about 9.17 per cent. By careful analysis of the difference of star counts between the downgraded SDSS northern halo data and SuperCOSMOS southern halo data, it is shown that no asymmetry can be detected in the South Galactic Cap at the accuracy of SuperCOSMOS, and the Virgo overdensity is likely a foreign component in the Galactic halo.     

Detection of the old stellar component of the major Galactic bar

We present near-infrared colour–magnitude diagrams and star counts for a number of regions along the Galactic plane. It is shown that along the l =27°, b =0° line of sight there is a feature at 5.7±0.7 kpc with a density of stars at least a factor of 2 and probably more than a factor of 5 times that of the disc at the same position. This feature forms a distinct clump on an H versus J − H diagram and is seen at all longitudes from the bulge to about l =28°, but at no longitude greater than this. The distance to the feature at l =20° is about 0.5 kpc further than at l =27°, and by l =10° it has merged with, or has become, the bulge. Given that at l =27° and l =21° there is also a clustering of very young stars, the only component that can reasonably explain what is seen is a bar with half-length of around 4 kpc and a position angle of about 43°±7°.     

Large-amplitude variables near the Galactic Centre

We report here the results of a 4-yr K -band (2.2 μm) survey for large-amplitude variable stars in a     area centred on the Galactic Centre. A total of 409 likely long-period variables (LPVs) were detected, for which positions, amplitudes, average magnitudes and periods were obtained whenever possible. The surface density of LPVs is more than ten times greater than in the Sgr I Baade window at        
The limits of completeness arising from interstellar and circumstellar absorption are discussed. Most of the area suffers interstellar extinction of     The shorter-period LPVs are less luminous than the longer-period ones and may be slightly under-represented in the data. Extremely heavy extinction     which affects the probability of detecting variables, occurs in less than 25 per cent of the area.
Almost all of the LPVs are Miras or OH/IR stars, with periods ranging from 150 d to about 800 d. K -band counterparts have been found for 59 per cent of the 109 known OH sources in the field. The average period of the variables found is 427 d, while that of the OH/IR stars is 524 d. For comparison, the average period in the Sgr I window, which contains no known OH/IR stars, is 333 d and only two stars are detected with     The survey field also contains a number of long-period, large-amplitude variables that are not OH emitters.     

A near-infrared survey of Miras and the distance to the Galactic Centre

We report the results of a near-infrared survey for long-period variables in a field of view of 20× 30 arcmin² towards the Galactic Centre (GC). We have detected 1364 variables, of which 348 are identified with those reported in Glass et al. We present a catalogue and photometric measurements for the detected variables and discuss their nature. We also establish a method for the simultaneous estimation of distances and extinctions using the period–luminosity relations for the JHK _s bands. Our method is applicable to Miras with periods in the range 100–350 d and mean magnitudes available in two or more filter bands. While J band means are often unavailable for our objects because of the large extinction, we estimated distances and extinctions for 143 Miras whose H - and   K _s  -band mean magnitudes are obtained. We find that most are located at the same distance to within our accuracy. Assuming that the barycentre of these Miras corresponds to the GC, we estimate its distance modulus to be  14.58 ± 0.02 (stat.) ± 0.11 (syst.) mag  , corresponding to  8.24 ± 0.08 (stat.) ± 0.42 (syst.) kpc  . We have assumed the distance modulus to the Large Magellanic Cloud to be 18.45 mag, and the uncertainty in this quantity is included in the above systematic error. We also discuss the large and highly variable extinction. Its value ranges from 1.5 mag to larger than 4 mag in     except towards the thicker dark nebulae and it varies in a complicated way with the line of sight. We have identified mid-infrared counterparts in the Spitzer /IRAC catalogue of Ramírez et al. for most of our variables and find that they follow rather narrow period–luminosity relations in the 3.6–8.0 μm wavelength range.     

A new, complete sample of faint B stars in the Galactic halo

A new sample of 31 faint B and A0 stars is reported, 30 of which comprise a complete sample within the limits ( U − V )<0 and 10.0< B ≲18.0. The sample is based on low- and intermediate-resolution spectrophotometry of colour-excess objects selected in the US survey. Atmospheric parameters for the stars are derived through the use of synthetic colours, Balmer-line strengths, and model-atmosphere fitting. The atmospheric parameters and preliminary metallicity estimates indicate that most of the stars are distributed along the blue horizontal branch, with low metallicities ([Fe/H]∼−1.0) and with both the first and second Newell gaps present. However, nine of the B/A0 stars can be identified as candidate main-sequence stars, based on evidence of high metallicities ([Fe/H]∼0) and/or derived effective temperatures and surface gravities which place them close to the main-sequence relation. The completeness characteristics of the sample are discussed, and its surface density is compared to that of other recently isolated B-star samples. The sample exhibits a shallow integral number-count slope. This new sample will help provide increased statistical coverage of the B-star population in the Galactic halo through its relatively faint magnitude-completeness limits and its relatively red colour-completeness limit.