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.     

Microlensing in dark matter haloes

Using eight dark matter haloes extracted from fully self-consistent cosmological N -body simulations, we perform microlensing experiments. A hypothetical observer is placed at a distance of 8.5 kpc from the centre of the halo measuring optical depths, event durations and event rates towards the direction of the Large Magellanic Cloud. We simulate 1600 microlensing experiments for each halo. Assuming that the whole halo consists of massive astronomical compact halo objects (MACHOs),   f = 1.0  , and a single MACHO mass is   m _M= 1.0 M_⊙  , the simulations yield mean values of  τ= 4.7^+5.0_−2.2× 10⁻⁷  and  Γ= 1.6^+1.3_−0.6× 10⁻⁶  events star⁻¹ yr⁻¹. We find that triaxiality and substructure can have major effects on the measured values so that τ and Γ values of up to three times the mean can be found. If we fit our values of τ and Γ to the MACHO collaboration observations, we find   f = 0.23^+0.15_−0.13  and   m _M= 0.44^+0.24_−0.16  . Five out of the eight haloes under investigation produce f and m _M values mainly concentrated within these bounds.     

The kinematics and zero-point of the log P–〈MK〉 relation for Galactic RR Lyrae variables via statistical parallax

We use accurate absolute proper motions and Two-Micron All-Sky Survey   K_s   -band apparent magnitudes for 364 Galactic RR Lyrae variables to determine the kinematical parameters of the Galactic RR Lyrae population and constrain the zero-point of the   K_s   -band period–luminosity relation for these stars via statistical parallax. We find the mean velocities of the halo- and thick-disc RR Lyrae populations in the solar neighbourhood to be  [ U ₀(Halo), V ₀(Halo), W ₀(Halo)]= (−12 ± 10, −217 ± 9, −6 ± 6) km s⁻¹  and  [ U ₀(Disc), V ₀(Disc), W ₀(Disc)]= (−15 ± 7, −44 ± 7, −25 ± 5) km s⁻¹  , respectively, and the corresponding components of the velocity-dispersion ellipsoids,  [σ V_R (Halo), σ V _θ(Halo), σ W (Halo)]= (167 ± 9, 86 ± 6, 78 ± 5) km s⁻¹  and  [σ V_R (Disc), σ V _θ(Disc), σ W (Disc)]= (55 ± 7, 44 ± 6, 30 ± 4) km s⁻¹  , respectively. The fraction of thick-disc stars is estimated at  0.25 ± 0.03  . The corrected infrared period–luminosity relation is     , implying a Large Magellanic Cloud (LMC) distance modulus of  18.27 ± 0.08  and a solar Galactocentric distance of  7.58 ± 0.40 kpc  . Our results suggest no or slightly prograde rotation for the population of halo RR Lyraes in the Milky Way.     

Parallax microlensing events in the OGLE II data base toward the Galactic bulge

We present a systematic search for parallax microlensing events among a total of 512 microlensing candidates in the OGLE II data base for the  1997–1999  seasons. We fit each microlensing candidate with both the standard microlensing model and a parallax model that accounts for the Earth's motion around the Sun. We then search for the parallax signature by comparing the χ ² of the standard and parallax models. For the events which show a significant improvement, we further use the 'duration' of the event and the signal-to-noise ratio as criteria to separate true parallax events from other noisy microlensing events. We have discovered one convincing new candidate, sc33_4505, and seven other marginal cases. The convincing candidate (sc33_4505) is caused by a slow-moving, and likely low-mass, object, similar to other known parallax events. We found that irregular sampling and gaps between observing seasons hamper the recovery of parallax events. We have also searched for long-duration events that do not show parallax signatures. The lack of parallax effects in a microlensing event puts a lower limit on the Einstein radius projected on to the observer plane, which in turn imposes a lower limit on the lens mass divided by the relative lens–source parallax. Most of the constraints are however quite weak.     

Microlensing of tidal debris on the Magellanic great circle

Increasing evidence suggests that the Galactic halo is lumpy on kpc scales as a result of the accretion of at least a dozen small galaxies [Large and Small Magellanic Clouds (LMC/SMC), Sgr, Fornax, etc.]. Faint stars in such lumpy structures can significantly microlense a background star with an optical depth of 10⁻⁷–10⁻⁶, which is comparable to the observed value to the LMC. The observed microlensing events towards the LMC can be explained by a tidal debris tail from the progenitor of the Magellanic Clouds and Magellanic Stream. The LMC stars can either lense stars in the debris tail a few kpc behind the LMC, or be lensed by stars in the part of the debris tail in front of the LMC. The models are consistent with an elementary particle dominated Galactic halo without massive compact halo objects (MACHOs). They also differ from Sahu's LMC-self-lensing model by predicting a higher optical depth and event rate and lower concentration of events to the LMC centre.     

SDSS absolute magnitudes for thin-disc stars based on trigonometric parallaxes

We present a new luminosity–colour relation based on trigonometric parallaxes for thin-disc main-sequence stars in Sloan Digital Sky Survey (SDSS) photometry. We matched stars from the newly reduced Hipparcos catalogue with the ones taken from Two-Micron All-Sky Survey (2MASS) All-Sky Catalogue of Point Sources, and applied a series of constraints, i.e. relative parallax errors  (σ_π/π≤ 0.05)  , metallicity  (−0.30 ≤[M/H]≤ 0.20 dex)  , age  (0 ≤ t ≤ 10 Gyr)  and surface gravity  (log  g > 4)  , and obtained a sample of thin-disc main-sequence stars. Then, we used our previous transformation equations ( Bilir et al. 2008a ) between SDSS and 2MASS photometries and calibrated the   M_g   absolute magnitudes to the  ( g − r )₀  and  ( r − i )₀  colours. The transformation formulae between 2MASS and SDSS photometries along with the absolute magnitude calibration provide space densities for bright stars which saturate the SDSS magnitudes.     

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  .     

Characterizing the pulsations of the ZZ Ceti star KUV 02464+3239

We present the results on period search and modelling of the cool DAV star KUV 02464+3239. Our observations resolved the multiperiodic pulsational behaviour of the star. In agreement with its position near the red edge of the DAV instability strip, it shows large amplitude, long-period pulsation modes, and has a strongly non-sinusoidal light curve. We determined six frequencies as normal modes and revealed remarkable short-term amplitude variations. A rigorous test was performed for the possible source of amplitude variation: beating of modes, effect of noise, unresolved frequencies or rotational triplets. Among the best-fitting models resulting from a grid search, we selected three that gave   l = 1  solutions for the largest amplitude modes. These models had masses of 0.645, 0.650 and  0.680  M _⊙  . The three 'favoured' models have   M _H  between  2.5 × 10⁻⁵ and 6.3 × 10⁻⁶  M _*  and give 14.2–14.8 mas seismological parallax. The  0.645  M _⊙  (11 400 K) model also matches the spectroscopic  log  g   and   T _eff  within 1σ. We investigated the possibility of mode trapping and concluded that while it can explain high amplitude modes, it is not required.     

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.     

Constraining quasar host halo masses with the strength of nearby Lyα forest absorption

Using cosmological hydrodynamic simulations, we measure the mean transmitted flux in the Lyα forest for quasar sightlines that pass near a foreground quasar. We find that the trend of absorption with pixel quasar separation distance can be fitted using a simple power-law form including the usual correlation function parameters r ₀ and γ, so that     . From the simulations, we find the relation between r ₀ and quasar host mass, and formulate this as a way to estimate quasar host dark matter halo masses, quantifying uncertainties due to cosmological and IGM parameters, and redshift errors. With this method, we examine data for ∼9000 quasars from the Sloan Digital Sky Survey (SDSS) Data Release 5, assuming that the effect of ionizing radiation from quasars (the so-called transverse proximity effect) is unimportant (no evidence for it is seen in the data). We find that the best-fitting host halo mass for SDSS quasars with mean redshift z = 3 and absolute G -band magnitude −27.5 is  log  M /M_⊙= 12.68^+0.81_−0.67  . We also use the Lyman-Break Galaxy (LBG) and Lyα forest data of Adelberger et al. in a similar fashion to constrain the halo mass of LBGs to be  log₁₀  M /M_⊙= 11.41^+0.54_−0.59  , a factor of ∼20 lower than the bright quasars. In addition, we study the redshift distortions of the Lyα forest around quasars, using the simulations. We use the quadrupole to monopole ratio of the quasar Lyα forest correlation function as a measure of the squashing effect. We find its dependence on halo mass difficult to measure, but find that it may be useful for constraining cosmic geometry.     

NGC 2587: a sparse open cluster projected on to a populous star field

We present CCD photometry in the Johnson U , B and V and Kron–Cousins I passbands for the open cluster NGC 2587. The sample consists of 4406 stars reaching down to   V ∼ 21.0  . We developed a new method to clean statistically the colour–magnitude diagrams. NGC 2587 appears to be a sparse, relatively bright open cluster, with a few tens of members projected on to a populous star field. The comparatively bright F7/8 II type star HD 70927, located close to the cluster centre, seems not to be a member. Our analysis suggests that NGC 2587 is slightly younger than the Hyades and probably of solar metallicity. A cluster radius of roughly 8 arcmin was estimated from the radial stellar density profile. From 18 probable cluster members with measured proper motions, we derive the following mean values for NGC 2587:  μ_α=−4.3 ± 3.6 mas yr⁻¹  and  μ_δ=−2.5 ± 3.4 mas yr⁻¹  . Adopting the theoretical metal content   Z = 0.02  , which provides the best global fit, we derive a cluster age of  500⁺⁶⁰₋₅₀  . Simultaneously, colour excesses   E ( B − V ) = 0.10  and   E ( V − I ) = 0.15  and an apparent distance modulus of   V − M_V = 12.50  are obtained. The interstellar extinction in the cluster direction is found to follow the normal law. NGC 2587 is located at a distance of (2.70 ± 0.70) kpc from the Sun and ∼9.8 kpc from the Galactic centre.     

Satellite kinematics – II. The halo mass–luminosity relation of central galaxies in SDSS

The kinematics of satellite galaxies reflect the masses of the extended dark matter haloes in which they orbit, and thus shed light on the mass–luminosity relation (MLR) of their corresponding central galaxies. In this paper, we select a large sample of centrals and satellites from the Sloan Digital Sky Survey and measure the kinematics (velocity dispersions) of the satellite galaxies as a function of the r -band luminosity of the central galaxies. Using the analytical framework presented in More, van den Bosch & Cacciato, we use these data to infer both the mean and the scatter of the MLR of central galaxies, carefully taking account of selection effects and biases introduced by the stacking procedure. As expected, brighter centrals on average reside in more massive haloes. In addition, we find that the scatter in halo masses for centrals of a given luminosity,  σ_{log  M}   , also increases with increasing luminosity. As we demonstrate, this is consistent with  σ_{log  L}   , which reflects the scatter in the conditional probability function   P ( L _c| M )  , being independent of halo mass. Our analysis of the satellite kinematics yields  σ_{log  L} = 0.16  ±  0.04  , in excellent agreement with constraints from clustering and group catalogues, and with predictions from a semi-analytical model of galaxy formation. We thus conclude that the amount of stochasticity in galaxy formation, which is characterized by  σ_{log  L}   , is well constrained, independent of halo mass and in a good agreement with current models of galaxy formation.     

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.     

Submillimetre continuum images of the NGC 2024 star-forming ridge

We present 450- and 800-μm images, made with the James Clerk Maxwell Telescope, of the NGC 2024 molecular ridge. The seven previously known compact cores, FIR1–7, have been detected, and FIR5 has been resolved into a compact object and an associated extended source to the east. The estimated masses of the dense cores vary between 1.6 and 5.1 M⊙ per 14-arcsec beam, assuming a dust temperature of 30 K and a dust opacity of κ_800 μm = 0.002 m² kg⁻¹. A spectral index map made from the 450- and 800-μm images shows spatial variations, with the spectral index, α ( F _ν ∝ ν^α), being systematically lower towards the dense cores. We interpret this as evidence for a lower value of the frequency dependence of the dust opacity, β, towards the denser cores relative to the surrounding molecular material. This may indicate that grain growth is occurring in the cores, prior to planetesimal formation. By comparing the high-resolution 450-μm image with interferometer maps of the integrated CS(2–1) emission, the previously reported discrepancy between dust continuum emission and molecular line emission is found to be very localized. Depletion and temperature variations are discussed as possible explanations.     

The black hole mass–stellar velocity dispersion correlation: bulges versus pseudo-bulges

We investigate the correlation between the supermassive black holes (SMBHs) mass ( M _bh) and the stellar velocity dispersion  (σ_*)  in two types of host galaxies: the early-type bulges (disc galaxies with classical bulges or elliptical galaxies) and pseudo-bulges. In the form  log ( M _bh/M_⊙) =α+β log (σ_*/200 km s⁻¹)  , the best-fitting results for the 39 early-type bulges are the slope  β= 4.06 ± 0.28  and the normalization  α= 8.28 ± 0.05  ; the best-fitting results for the nine pseudo-bulges are  β= 4.5 ± 1.3  and  α= 7.50 ± 0.18  . Both relations have intrinsic scatter in  log  M _bh  of ≲0.27 dex. The   M _bh–σ_*  relation for pseudo-bulges is different from the relation in the early-type bulges over the 3σ significance level. The contrasting relations indicate the formation and growth histories of SMBHs depend on their host type. The discrepancy between the slope of the   M _bh–σ_*  relations using different definition of velocity dispersion vanishes in our sample, a uniform slope will constrain the coevolution theories of the SMBHs and their host galaxies more effectively. We also find the slope for the 'core' elliptical galaxies at the high-mass range of the relation appears steeper  (β≃ 5–6)  , which may be the imprint of their origin of dissipationless mergers.     

CCD photometry of variable stars in the Magellanic Clouds — VI. The eclipsing binary HV 982 in the Large Magellanic Cloud

Using improved techniques, high-quality CCD uvbyVI photometry has been obtained for the eclipsing binary HV 982 in the Large Magellanic Cloud (LMC). International Ultraviolet Explorer ultraviolet spectrophotometry was also obtained. These data have been analysed using the Wilson–Devinney synthetic light-curve code and Kurucz low-metallicity model atmospheres as well as the EBOP code. The system is detached and the orbit is eccentric. Apsidal motion is detected with apsidal period 205 ± 7 yr. The effective temperatures of the components are found via flux fitting to be T _eff,1 = 28 000 ± 5000 K and T _eff,2 = 27 200 ± 5000 K. The large errors result from uncertainties over the appropriate interstellar extinction correction. The system plausibly comprises two ∼ 8 M stars of radius 6–7 R separated by ∼ 30 R. For pedagogical and historical interest, the near simultaneity of the eclipse minima at different wavelengths is used to constrain the constancy of the speed of light with wavelength and the mass of the photon, yielding m _γ < 10⁻⁴¹ kg. Because of the great distance to HV 982, this limit is some 10² times smaller than previously achieved with eclipse timings, but it is nevertheless 10 orders of magnitude less stringent than that which is provided by satellite measurements of planetary magnetic fields.     

The discovery of a low-mass, pre-main-sequence stellar association around γ Velorum

We report the serendipitous discovery of a population of low-mass, pre-main-sequence (PMS) stars in the direction of the Wolf–Rayet/O-star binary system γ ²  Vel and the Vela OB2 association. We argue that γ ²  Vel and the low-mass stars are truly associated and approximately coeval, and that both are at distances between 360 and 490 pc, disagreeing at the 2 σ level with the recent Hipparcos parallax of γ ²  Vel, but consistent with older distance estimates. Our results clearly have implications for the physical parameters of the γ ²  Vel system, but also offer an exciting opportunity to investigate the influence of high-mass stars on the mass function and circumstellar disc lifetimes of their lower mass PMS siblings.     

Constraints on Jupiters from observations of Galactic bulge microlensing events during 2000

We present observations of eight Galactic bulge microlensing events taken with the 1.0-m Jacobus Kapteyn Telescope (JKT) on La Palma during 2000 June and July. The JKT observing schedule was optimized using a prioritizing algorithm to automatically update the target list. For most of these events we have sampled the light curves at times where no information was available from the OGLE alert team. We assume a point-source point-lens (PSPL) model and perform a maximum likelihood fit to both our data and the OGLE data to constrain the event parameters of the fit. We then refit the data assuming a binary lens and proceed to calculate the probability of detecting planets with mass ratio   q = 10⁻³  . We have seen no clear signatures of planetary deviations on any of the eight events and we quantify constraints on the presence of planetary companions to the lensing stars. For two well-observed events, 2000BUL31 and 2000BUL33, our detection probabilities peak at ∼30 and ∼20 per cent respectively for   q = 10⁻³  and   a ∼ R _E  for a  Δχ²  threshold value of 60.     

Mass-to-light ratio gradients in early-type galaxy haloes

Owing to the fact that the near future should see a rapidly expanding set of probes of the halo masses of individual early-type galaxies, we introduce a convenient parameter for characterizing the halo masses from both observational and theoretical results:  ∇_ℓϒ  , the logarithmic radial gradient of the mass-to-light ratio. Using halo density profiles from Λ-cold dark matter (CDM) simulations, we derive predictions for this gradient for various galaxy luminosities and star formation efficiencies  ε_SF  . As a pilot study, we assemble the available  ∇_ℓϒ  data from kinematics in early-type galaxies – representing the first unbiased study of halo masses in a wide range of early-type galaxy luminosities – and find a correlation between luminosity and  ∇_ℓϒ  , such that the brightest galaxies appear the most dark-matter dominated. We find that the gradients in most of the brightest galaxies may fit in well with the ΛCDM predictions, but that there is also a population of fainter galaxies whose gradients are so low as to imply an unreasonably high star formation efficiency  ε_SF > 1  . This difficulty is eased if dark haloes are not assumed to have the standard ΛCDM profiles, but lower central concentrations.