The impact of a merger episode in the galactic disc white dwarf population

In this paper we analyse the consequences in the white dwarf population of a hypothetical merger episode in our Galactic disc. We have studied several different merging scenarios with our Monte Carlo simulator. For each one of these scenarios we have derived the main characteristics of the resulting white dwarf population and we have compared them with the available observational data, namely the white dwarf luminosity function and the kinematic properties of the white dwarf population. Our results indicate that very recent (less than ∼6 Gyr ago) and massive (∼16 per cent of the mass of our Galaxy) merger episodes are quite unlikely in view of the available kinematical properties of the disc white dwarf population. Smaller merger episodes (of the order of ∼4 per cent of the mass of our Galaxy) are, however, compatible with our current knowledge of those kinematical properties. Finally, we prove that the white dwarf luminosity function is quite insensitive to such a merger episode.     

Formation of ω Centauri from an ancient nucleated dwarf galaxy in the young Galactic disc

We first present a self-consistent dynamical model in which ω Cen is formed from an ancient nucleated dwarf galaxy merging with the first generation of the Galactic thin disc in a retrograde manner with respect to the Galactic rotation. Our numerical simulations demonstrate that during merging between the Galaxy and the ω Cen host dwarf with   M _B≃−14 mag  and its nucleus mass of  10⁷ M⊙  , the outer stellar envelope of the dwarf is nearly completely stripped, whereas the central nucleus can survive from the tidal stripping because of its compactness. The developed naked nucleus has a very bound retrograde orbit around the young Galactic disc, as observed for ω Cen, with apocentre and pericentre distances of ∼8 and ∼1 kpc, respectively. The Galactic tidal force can induce radial inflow of gas to the centre of the dwarf and consequently triggers moderately strong nuclear starbursts in a repetitive manner. This result implies that efficient nuclear chemical enrichment resulting from the later starbursts can be closely associated with the origin of the observed relatively young and metal-rich stars in ω Cen. Dynamical heating by the ω Cen host can transform the young thin disc into the thick disc during merging.     

The white dwarf luminosity function – II. The effect of the measurement errors and other biases

The disc white dwarf luminosity function is an important tool for studying the solar neighbourhood, since it allows the determination of several Galactic parameters, the most important one being the age of the Galactic disc. However, only the     method has been employed so far for observationally determining the white dwarf luminosity function, whereas for other kind of luminosity functions several other methods have been frequently used. Moreover, the procedures to determine the white dwarf luminosity function are not free of biases. These biases have two different origins: they can either be of statistical nature or a consequence of the measurement errors. In a previous paper we carried out an in-depth study of the first category of biases for several luminosity function estimators. In this paper we focus on the biases introduced by the measurement errors and on the effects of the degree of contamination of the input sample used to build the disc white dwarf luminosity function by different kinematical populations. To assess the extent of these biases we use a Monte Carlo simulator to generate a controlled synthetic population and analyse the behaviour of the disc white dwarf luminosity function for several assumptions about the magnitude of the measurement errors and for several degrees of contamination, comparing the performances of the most robust luminosity function estimators under such conditions.     

The field of NGC 6397 as a test for Galactic models

Taking advantage of recent Hubble Space Telescope ( HST ) data for field stars in the region of the Galactic globular cluster NGC 6397, we tested the predictions of several Galactic models with star counts reaching a largely unexplored range of magnitudes, down to V ∼26.5. After updating the input stellar ( V − I ) colours, we found that the two-component Bahcall–Soneira (B&S) model can be put into satisfactory agreement with observations for suitable choices of disc/spheroid luminosity functions (LFs). However, if one assumes the disc LF of Gould, Bahcall and Flynn together with the spheroid LF of Gould, Flynn and Bahcall, there is no way to reconcile the predicted and observed V -magnitude distribution. We also analysed the agreement between observed and predicted magnitude and colour distributions for two selected models with a thick disc component. Even in this case there are suitable combinations of model parameters and faint-magnitude LFs that can give a reasonable agreement with observational star counts in both magnitude and colour. However, the above-quoted combination of Gould et al. LFs again gives predictions in clear disagreement with observations.     

Something about the structure of the Galaxy

We analyse a sample of 507 evolved (OH/IR) stars in the region (10°>ℓ>−45°), (| b |<3°). We derive average ages for subsets of this sample, and use those sets as beacons for the evolution of the Galaxy. In the bulge, the oldest OH/IR stars in the plane are 7.5 Gyr (1.3 M_⊙), and in the disc 2.7 Gyr (2.3 M_⊙). The vertical distribution of almost all AGB stars in the disc is found to be nearly exponential, with scaleheight increasing from 100 pc for stars ≲1 Gyr old to 500 pc for stars ≳5 Gyr old. There may be a small, disjunct population of OH/IR stars. The radial distribution of AGB stars is dictated by the metallicity gradient. Unequivocal morphological evidence is presented for the existence of a central bar, but parameters can be constrained only for a given spatial-density model. Using a variety of indicators, we identify the radii of the inner ultraharmonic (2.5 kpc) and corotation resonance (3.5 kpc). We show that the 3-kpc arm is likely to be an inner ring, as observed in other barred galaxies, by identifying a group of evolved stars that is connected to the 3-kpc H  i filament. Also, using several observed features, we argue that an inner-Lindblad resonance exists, at ∼1–1.5 kpc. The compositions of OH/IR populations within 1 kpc of the Galactic Centre give insight into the bar-driven evolution of the inner regions. We suggest that the bar is ∼8 Gyr old, relatively weak (SAB), and may be in a final stage of its existence.     

The APM survey for cool carbon stars in the Galactic halo – I

A by-product of the APM high-redshift quasar survey was the discovery of several distant (20–100 kpc) N-type carbon stars at high galactic latitude. Following on from this, we have started a systematic all-sky survey at galactic latitudes ⊢ b ⊢>30° to find further examples of these rare objects, and we report here on the results from the first season of follow-up spectroscopy. Faint, high-latitude carbon (FHLC) giants make excellent probes of the kinematic structure of the outer Galactic halo. Therefore, in addition to detailed spectrophotometry covering a wide wavelength range, we have obtained high-resolution (∼1 Å) spectra centred on the CN bands at ∼8000 Å, and have derived accurate (≲10 km s⁻¹) radial velocities for the known FHLC stars. From the initial phase of our survey covering ≈6500 deg², we find a surface density of faint N-type carbon stars in the halo of ≈1 per 200 deg², roughly a factor of 4 less than the surface density of CH-type carbon stars in the halo. Intermediate-age, N-type carbon stars seem unlikely to have formed in the halo in isolation from other star-forming regions, and one possibility that we are investigating is that they either arise from the disruption of tidally captured dwarf satellite galaxies or are a manifestation of the long-sought optical component of the Magellanic Stream.     

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.     

High-velocity stars from decay of small stellar systems

In this paper we present numerical results on the decay of small stellar systems under different initial conditions (multiplicity 3 ≤  N  ≤ 10, and various mass spectra, initial velocities and initial configurations). The numerical treatment uses the CHAIN1 code (Mikkola &38; Aarseth). Particular attention is paid to the distribution of high-velocity escapers: we define these as stars with velocity above 30 km s⁻¹. These numerical experiments confirm that small N -body systems are dynamically unstable and produce cascades of escapers in the process of their decay. It is shown that the fraction of stars that escape from small dense stellar systems with an escape velocity greater than 30 km s⁻¹ is ∼1 per cent for all systems treated here. This relatively small fraction must be considered in relation to the rate of star formation in the Galaxy in small groups: this could explain some moderately high-velocity stars observed in the Galactic disc and possibly some young stars with relatively high metallicity in the thick disc.     

High spatial resolution Galactic 3D extinction mapping with IPHAS

We present an algorithm (MEAD, for 'Mapping Extinction Against Distance') which will determine intrinsic  ( r '− i ')  colour, extinction, and distance for early-A to K4 stars extracted from the IPHAS   r '/ i '/Hα  photometric data base. These data can be binned up to map extinction in three dimensions across the northern Galactic plane. The large size of the IPHAS data base (∼200 million unique objects), the accuracy of the digital photometry it contains and its faint limiting magnitude  ( r '∼ 20)  allow extinction to be mapped with fine angular (∼10 arcmin) and distance (∼ 0.1 kpc) resolution to distances of up to 10 kpc, outside the solar circle. High reddening within the solar circle on occasion brings this range down to ∼2 kpc. The resolution achieved, both in angle and depth, greatly exceeds that of previous empirical 3D extinction maps, enabling the structure of the Galactic Plane to be studied in increased detail. MEAD accounts for the effect of the survey magnitude limits, photometric errors, unresolved interstellar medium (ISM) substructure and binarity. The impact of metallicity variations, within the range typical of the Galactic disc is small. The accuracy and reliability of MEAD are tested through the use of simulated photometry created with Monte Carlo sampling techniques. The success of this algorithm is demonstrated on a selection of fields and the results are compared to the literature.     

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.     

Stellar populations in the Canis Major overdensity

We performed a photometric multicolour survey of the core of the Canis Major overdensity at     , reaching   V ∼ 22  and covering  0.3 × 1.0  arcmin². The main aim is to unravel the complex mixture of stellar populations toward this Galactic direction, where in the recent past important signatures of an accretion event have been claimed to be detected. While our previous investigations were based on disjointed pointings aimed at revealing the large-scale structure of the third Galactic Quadrant, we now focus on a complete coverage of a smaller field centred on the Canis Major overdensity. A large wavelength baseline, in the UBVRI bands, allows us to build up a suite of colour–colour and colour–magnitude diagrams, providing a much better diagnostic tool to disentangle the stellar populations of the region. In fact, the simple use of one colour–magnitude diagram, widely employed in all the previous studies defending the existence of the Canis Major galaxy, does not allow one to separate the effects of the different parameters (reddening, age, metallicity and distance) involved in the interpretation of data, forcing to rely on heavy modelling. In agreement with our previous studies, in the same general region of the Milky Way, we recognize a young stellar population compatible with the expected structure and extension of the Local (Orion) and Outer (Norma–Cygnus) spiral arms in the Third Galactic Quadrant. Moreover, we interpret the conspicuous intermediate-age metal-poor population as belonging to the Galactic thick disc, distorted by the effect of strong disc warping at this latitude, and to the Galactic halo.     

Could the Galactic disc heating be due to globular cluster impacts?

So far, six mechanisms have been proposed to account for the Galactic disc heating. Of these, the most important appear to be a combination of scattering of stars by molecular clouds and by spiral arms. We study a further mechanism, namely the repeated disc impact of the original Galactic globular cluster population up to the present. We find that globular clusters could have contributed at most a small fraction of the current vertical energy of the disc, as they could heat the whole disc to  σ _z = 5.5 km s⁻¹  (c.f. the observed 18 and 39 km s⁻¹ for the thick and thin discs, respectively). We find that the rate of rise of disc heat (  α= 0.22  in  σ _z ∼ t ^α  with t being time) is close to that found for scattering by molecular clouds.     

The orbit and mass of the Sagittarius dwarf galaxy

Possible orbital histories of the Sgr dwarf galaxy are explored. A special-purpose N -body code is used to construct the first models of the Milky Way–Sgr dwarf system in which both the Milky Way and the Sgr dwarf are represented by full N -body systems and followed for a Hubble time. These models are used to calibrate a semi-analytic model of the Sgr dwarf's orbit that enables us to explore a wider parameter space than is accessible to the N -body models. We conclude that the extant data on the Sgr dwarf are compatible with a wide range of orbital histories. At one extreme the Sgr dwarf initially possesses ∼10¹¹ M_⊙ and starts from a Galactocentric distance R _D(0)≳200 kpc. At the other extreme the Sgr dwarf starts with ∼10⁹ M_⊙ and R _D(0)∼60 kpc, similar to its present apocentric distance. In all cases the Sgr dwarf is initially dark matter dominated and the current velocity dispersion of the Sgr dwarf's dark matter is tightly constrained to be 21±2 km s⁻¹. This number is probably compatible with the smaller measured dispersion of the Sgr dwarf's stars because of (i) the dynamical difference between dark and luminous matter, and (ii) velocity anisotropy.     

An Inventory of Gravitational Microlenses Toward the Galactic Bulge

Although microlensing experiments toward the Galactic bulge were originally initiated to check the feasibility of the experiments, they have now become an important tool which allows one to investigate Galactic matter composition. However, previous determination of the lens mass function was not based on the actual number of lenses even for the well-determined population of stellar lenses, but rather on arbitrarily assumed functional forms such as power laws, and thus the derived mass functions were subject to large uncertainties. In this paper, we take a different approach in which we first estimate the event rate distribution expected from observationally well-constrained populations of lenses, and then test other possible lens populations. By comparing the determined event rate distribution Γ( t _E) for various mass function models of lens populations with the observed distribution, we find that stars and white dwarfs explain just ∼ 50 per cent of the total observed events even including very faint stars just above the hydrogen-burning limit. Additionally, the expected time-scale distribution of events caused by these known populations of lenses deviates significantly from the observed distribution, especially in the short time-scale region. However, if the rest of the dynamical mass of the bulge (∼ 2.1 × 10¹⁰ M) is composed of brown dwarfs, the expected event rate distribution matches the observation well.     

The environment of formation as a second parameter for globular cluster classification

We perform an evolutionary multivariate analysis of a sample of 54 Galactic globular clusters with high-quality colour–magnitude diagrams and well-determined ages. The four parameters adopted for the analysis are: metallicity, age, maximum temperature on the horizontal branch and absolute V magnitude. Our cladistic analysis breaks the sample into three novel groups. An a posteriori kinematical analysis puts groups 1 and 2 in the halo, and group 3 in the thick disc. The halo and disc clusters separately follow a luminosity–metallicity relation of much weaker slope than galaxies. This property is used to propose a new criterion for distinguishing halo and disc clusters. A comparison of the distinct properties of the two halo groups with those of Galactic halo field stars indicates that the clusters of group 1 originated in the inner halo, while those of group 2 formed in the outer halo of the Galaxy. The inner halo clusters were presumably initially the most massive one, which allowed the formation of more strongly helium-enriched second generation stars, thus explaining the presence of Cepheids and of very hot horizontal-branch stars exclusively in this group. We thus conclude that the 'second parameter' is linked to the environment in which globular clusters form, the inner halo favouring the formation of the most massive clusters which subsequently become more strongly self-enriched than their counterparts of the galactic outer halo and disc.     

Three-dimensional evolution of worms and chimneys in the Galactic disc

Three-dimensional modelling of the flow of gas and plasma in a section of the Galaxy has been carried out to study the evolution and formation of Galactic chimneys and worms. It is found that clustered supernovae located on either side of the Galactic plane are sources for the formation of well-collimated chimneys, having widths of ∼     . The thick gas disc may have a role in the collimation of chimneys. Channel maps of disc gas, obtained from the simulations, show the presence of sheet-like structures running perpendicular to the Galactic plane and resembling worms. Worms are believed to result from the break-up of the shells and supershells. However, the simulations show that although some worms correlate well with the debris of broken shells/supershells, others do not. They are cold gas that has been accelerated in the disc and rise on to the thick gas disc.     

The birthrate of magnetars

Magnetars, neutron stars with ultrastrong magnetic fields  ( B ∼ 10¹⁴−10¹⁵G)  , manifest their exotic nature in the form of soft gamma-ray repeaters and anomalous X-ray pulsars. This study estimates the birthrate of magnetars to be ∼0.22 per century with a Galactic population comprising ∼17 objects. A population synthesis was carried out based on the five anomalous X-ray pulsars detected in the ROSAT All Sky Survey by comparing their number to that of massive OB stars in a well-defined volume. Additionally, the group of seven X-ray dim isolated neutron stars detected in the same survey were found to have a birthrate of ∼2 per century with a Galactic population of ∼22 000 objects.     

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.     

Searching for hidden Wolf–Rayet stars in the Galactic plane – 15 new Wolf–Rayet stars

We report the discovery of 15 previously unknown Wolf–Rayet (WR) stars found as part of an infrared (IR) broad-band study of candidate WR stars in the Galaxy. We have derived an empirically based selection algorithm which has selected ∼5000 WR candidate stars located within the Galactic plane drawn from the Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (mid-IR) and Two-Micron All-Sky Survey (near-IR) catalogues. Spectroscopic follow-up of 184 of these reveals 11 nitrogen-rich (WN) and four carbon-rich (WC) WR stars. Early WC subtypes are absent from our sample and none shows evidence for circumstellar dust emission. Of the candidates which are not WR stars, ∼120 displayed hydrogen emission-line features in their spectra. Spectral features suggest that the majority of these are in fact B supergiants/hypergiants, ∼40 of these are identified Be/B[e] candidates.
Here, we present the optical spectra for six of the newly detected WR stars, and the near-IR spectra for the remaining nine of our sample. With a WR yield rate of ∼7 per cent and a massive star detection rate of ∼65 per cent, initial results suggest that this method is one of the most successful means for locating evolved, massive stars in the Galaxy.