Late superhumps and the stream–disc impact in IY UMa

We present two simple dynamical models for Sagittarius based on N -body simulations of the progressive disruption of a satellite galaxy orbiting for 12.5 Gyr within a realistic Galactic potential. In both models the satellite initially has observable properties similar to those of current outlying dwarfs; in one case it is purely stellar while in the other it is embedded in an extended massive halo. The purely stellar progenitor is a King model with a total velocity dispersion of 18.9 km s⁻¹, a core radius of 0.44 kpc and a tidal radius of 3 kpc. The initial stellar distribution in the other case follows a King profile with the same core radius, a slightly larger total velocity dispersion and similar extent. Both these models are consistent with all published data on the current Sagittarius system, they match not only the observed properties of the main body of Sagittarius, but also those reported for unbound debris at larger distances.     

The growth of central and satellite galaxies in cosmological smoothed particle hydrodynamics simulations

We examine the accretion and merger histories of central and satellite galaxies in a smoothed particle hydrodynamics (SPH) cosmological simulation that resolves galaxies down to  7 × 10⁹ M_⊙  . Most friends-of-friends haloes in the simulation have a distinct central galaxy, typically 2–5 times more massive than the most massive satellite. As expected, satellites have systematically higher assembly redshifts than central galaxies of the same baryonic mass, and satellites in more massive haloes form earlier. However, contrary to the simplest expectations, satellite galaxies continue to accrete gas and convert it to stars; the gas accretion declines steadily over a period of 0.5–1 Gyr after the satellite halo merges with a larger parent halo. Satellites in a cluster mass halo eventually begin to lose baryonic mass. Typically, satellites in our simulation are 0.1–0.2 mag bluer than in models that assume no gas accretion on to satellites after a halo merger. Since   z = 1  , 27 per cent of central galaxies (above  3 × 10¹⁰ M_⊙  ) and 22 per cent of present-day satellite galaxies have merged with a smaller system above a 1:4 mass ratio; about half of the satellite mergers occurred after the galaxy became a satellite and half before. In effect, satellite galaxies can remain 'central' objects of halo substructures, with continuing accretion and mergers, making the transition in assembly histories and physical properties a gradual one. Implementing such a gradual transformation in semi-analytic models would improve their agreement with observed colour distributions of satellite galaxies in groups and with the observed colour dependence of galaxy clustering.     

Kinematics of massive stars in the Small Magellanic Cloud

We present radial velocities for 2045 stars in the Small Magellanic Cloud (SMC), obtained from the 2dF survey by Evans et al. The great majority of these stars are of OBA type, tracing the dynamics of the young stellar population. Dividing the sample into ad hoc 'bar' and 'wing' samples (north and south, respectively, of the line:  δ=−77°50'+[4α]'  , where α is the right ascension in minutes of time) we find that the velocities in the SMC bar show a gradient of 26.3 ± 1.6 km s⁻¹ deg⁻¹ at a position angle of 126°± 4°. The derived gradient in the bar is robust to the adopted line of demarcation between the two samples. The largest redshifts are found in the SMC wing, in which the velocity distribution appears distinct from that in the bar, most probably a consequence of the interaction between the Magellanic Clouds that is predicted to have occurred 0.2 Gyr ago. The mean velocity for all stars in the sample is +172.0 ± 0.2 km s⁻¹ (redshifted by ∼20 km s⁻¹ when compared to published results for older populations), with a velocity dispersion of 30 km s⁻¹.     

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.     

Spectroscopy of the optical Einstein ring 0047 – 2808

We present optical and near-infrared spectroscopic observations of the optical Einstein ring 0047 – 2808. We detect both [O III ] lines λλ4959, 5007 near ∼ 2.3 μm, confirming the redshift of the lensed source as z  = 3.595. The Lyα line is redshifted relative to the [O III ] line by 140 ± 20 km s⁻¹. Similar velocity shifts have been seen in nearby starburst galaxies. The [O III ] line is very narrow, 130 km s⁻¹ FWHM. If the ring is the image of the centre of a galaxy, the one-dimensional stellar velocity dispersion σ = 55 km s⁻¹ is considerably smaller than the value predicted by Baugh et al. for the somewhat brighter Lyman-break galaxies. The Lyα line is significantly broader than the [O III ] line, probably due to resonant scattering. The stellar central velocity dispersion of the early-type deflector galaxy at z  = 0.485 is 250 ± 30 km s⁻¹. This value is in good agreement both with the value predicted from the radius of the Einstein ring (and a singular isothermal sphere model for the deflector), and with the value estimated from the D _n −σ relation.     

The pattern speed of the OH/IR stars in the Milky Way

We show how the continuity equation can be used to determine pattern speeds in the Milky Way Galaxy (MWG). This method, first discussed by Tremaine & Weinberg in the context of external galaxies, requires projected positions, ( l , b ), and line-of-sight velocities for a spatially complete sample of relaxed tracers. If the local standard of rest (LSR) has a zero velocity in the radial direction ( u _LSR), then the quantity that is measured is  Δ V ≡Ω_p R ₀- V _LSR  , where Ω_p is the pattern speed of the non-axisymmetric feature, R ₀ is the distance of the Sun from the Galactic centre and V _LSR is the tangential motion of the LSR, including the circular velocity. We use simple models to assess the reliability of the method for measuring a single, constant pattern speed of either a bar or spiral in the inner MWG. We then apply the method to the OH/IR stars in the ATCA/VLA OH 1612-MHz survey of Sevenster et al., finding  Δ V =252±41 km s^-1,  if   u _LSR=0  . Assuming further that   R ₀=8 kpc  and   V _LSR=220 km s^-1,  this gives  Ω_p=59±5 km s^-1 kpc^-1  with a possible systematic error of perhaps 10 km s⁻¹ kpc⁻¹. The non-axisymmetric feature for which we measure this pattern speed must be in the disc of the MWG.     

The host galaxy of a narrow-line Seyfert 1 galaxy, RE J1034+396, with X-ray quasi-periodic oscillations

Using simple stellar population synthesis, we model the bulge stellar contribution in the optical spectrum of a narrow-line Seyfert 1 galaxy, RE J1034+396. We find that its bulge stellar velocity dispersion is  67.7 ± 8 km s⁻¹  . The supermassive black hole (SMBH) mass is about  (1–4) × 10⁶ M_⊙  if it follows the well-known   M _BH–σ_*  relation found in quiescent galaxies. We also derive the SMBH mass from the Hβ second moment, which is consistent with that from its bulge stellar velocity dispersion. The SMBH mass of (1–4)  × 10⁶ M_⊙  implies that the X-ray quasi-periodic oscillation (QPO) of RE J1034+396 can be scaled to a high-frequency QPO at 27–108 Hz found in Galactic black hole binaries with a  10-M_⊙  black hole. With the mass distribution in different age stellar populations, we find that the mean specific star formation rate (SSFR) over the past 0.1 Gyr is  0.0163 ± 0.0011  Gyr⁻¹, the stellar mass in the logarithm is  10.155 ± 0.06  in units of solar mass and the current star formation rate is  0.23 ± 0.016 M_⊙ yr⁻¹  . For RE J1034+396, there is no relation between the Eddington ratio and the SSFR as suggested by Chen et al., despite a larger scatter in their relation. We also suggest that about 7.0 per cent of the total Hα luminosity and 50 per cent of the total [O  ii ] luminosity come from the star formation process.     

Is Ursa Major II the progenitor of the Orphan Stream?

Prominent in the 'Field of Streams'– the Sloan Digital Sky Survey map of substructure in the Galactic halo – is an 'Orphan Stream' without obvious progenitor. In this numerical study, we show a possible connection between the newly found dwarf satellite Ursa Major II (UMa II) and the Orphan Stream. We provide numerical simulations of the disruption of UMa II that match the observational data on the position, distance and morphology of the Orphan Stream. We predict the radial velocity of UMa II as −100 km s⁻¹, as well as the existence of strong velocity gradients along the Orphan Stream. The velocity dispersion of UMa II is expected to be high, though this can be caused both by a high dark matter content or by the presence of unbound stars in a disrupted remnant. However, the existence of a gradient in the mean radial velocity across UMa II provides a clear-cut distinction between these possibilities. The simulations support the idea that some of the anomalous, young halo globular clusters like Palomar 1 or Arp 2 or Ruprecht 106 may be physically associated with the Orphan Stream.     

The discovery of Segue 2: a prototype of the population of satellites of satellites

We announce the discovery of a new Milky Way satellite Segue 2 found in the data of the Sloan Extension for Galactic Understanding and Exploration (SEGUE). We followed this up with deeper imaging and spectroscopy on the Multiple Mirror Telescope (MMT). From this, we derive a luminosity of   M _v =−2.5  , a half-light radius of 34 pc and a systemic velocity of  ∼−40 km s⁻¹  . Our data also provide evidence for a stream around Segue 2 at a similar heliocentric velocity, and the SEGUE data show that it is also present in neighbouring fields. We resolve the velocity dispersion of Segue 2 as 3.4 km s⁻¹ and the possible stream as  ∼7 km s⁻¹  . This object shows points of comparison with other recent discoveries, Segue 1, Boo II and Coma. We speculate that all four objects may be representatives of a population of satellites of satellites – survivors of accretion events that destroyed their larger but less dense parents. They are likely to have formed at redshifts   z > 10  and are good candidates for fossils of the reionization epoch.     

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.     

Cosmological origin of the lowest metallicity halo stars

We explore the predictions of the standard hierarchical clustering scenario of galaxy formation, regarding the numbers and metallicities of PopIII stars that are likely to be found within our Galaxy today. By PopIII we refer to stars formed at large redshift ( z >4), with low metallicities ([ Z /Z_⊙]<−2.5) and in small systems (total mass ≲ 2×10⁸ M_⊙) that are extremely sensitive to stellar feedback, and which through a prescribed merging history end up becoming part of the Milky Way today. An analytic, extended Press–Schechter formalism is used to obtain the mass functions of haloes which will host PopIII stars at a given redshift, and which will end up in Milky Way sized systems today. Each of these is modelled as a mini-galaxy, with a detailed treatment of the dark halo structure, angular momentum distribution, final gas temperature and disc instabilities, all of which determine the fraction of the baryons that are subject to star formation. The use of new primordial metallicity stellar evolutionary models allows us to trace the history of the stars formed, and give accurate estimates of their expected numbers today and their location in L /L_⊙ versus T /K Hertzsprung–Russell (HR) diagrams. A first comparison with observational data suggests that the initial mass function (IMF) of the first stars was increasingly high-mass weighted towards high redshifts, levelling off at z ≳9 at a characteristic stellar mass scale m _s=10–15 M_⊙.     

Dark matter in dwarf spheroidals – II. Observations and modelling of Draco

We present stellar radial velocity data for the Draco dwarf spheroidal (dSph) galaxy obtained using the AF2/WYFFOS instrument combination on the William Herschel Telescope. Our data set consists of 186 member stars, 159 of which have good quality velocities, extending to a magnitude   V ≈19.5  with a mean velocity precision of ≈2 km s⁻¹. As this survey is based on a high-precision photometric target list, it contains many more Draco members at large radii. For the first time, this allows a robust determination of the radial behaviour of the velocity dispersion in a dSph.
We find statistically strong evidence of a rising velocity dispersion consistent with a dark matter halo that has a gently rising rotation curve. There is a <2 σ signature of rotation about the long axis, inconsistent with tidal disruption as the source of the rising dispersion. By comparing our data set with earlier velocities, we find that Draco probably has a binary distribution and fraction comparable to those in the solar neighbourhood.
We apply a novel maximum likelihood algorithm and fit the velocity data to a two parameter spherical model with an adjustable dark matter content and velocity anisotropy. Draco is best fit by a weakly tangentially anisotropic distribution of stellar orbits in a dark matter halo with a very slowly rising rotation law  ( v _circ∝ r ^0.17)  . We are able to rule out both a mass-follows-light distribution and an extended halo with a harmonic core at the 2.5 to 3 σ significance level, depending on the details of our assumptions about Draco's stellar binary population. Our modelling lends support to the idea that the dark matter in dwarf spheroidals is distributed in the form of massive, nearly isothermal haloes.     

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.     

Testing fundamental physics with distant star clusters: theoretical models for pressure-supported stellar systems

We investigate the mean velocity dispersion and the velocity dispersion profile of stellar systems in modified Newtonian dynamics (MOND), using the N -body code n-mody , which is a particle-mesh-based code with a numerical MOND potential solver developed by Ciotti, Londrillo & Nipoti. We have calculated mean velocity dispersions for stellar systems following Plummer density distributions with masses in the range of 10⁴ to  10⁹ M_⊙  and which are either isolated or immersed in an external field. Our integrations reproduce previous analytic estimates for stellar velocities in systems in the deep MOND regime  ( a _i, a _e≪ a ₀)  , where the motion of stars is either dominated by internal accelerations  ( a _i≫ a _e)  or constant external accelerations  ( a _e≫ a _i)  . In addition, we derive for the first time analytic formulae for the line-of-sight velocity dispersion in the intermediate regime  ( a _i∼ a _e∼ a ₀)  . This allows for a much-improved comparison of MOND with observed velocity dispersions of stellar systems. We finally derive the velocity dispersion of the globular cluster Pal 14 as one of the outer Milky Way halo globular clusters that have recently been proposed as a differentiator between Newtonian and MONDian dynamics.     

A model of supernova feedback in galaxy formation

A model of supernova feedback during disc galaxy formation is developed. The model incorporates infall of cooling gas from a halo, and outflow of hot gas from a multiphase interstellar medium (ISM). The star formation rate is determined by balancing the energy dissipated in collisions between cold gas clouds with that supplied by supernovae in a disc marginally unstable to axisymmetric instabilities. Hot gas is created by thermal evaporation of cold gas clouds in supernova remnants, and criteria are derived to estimate the characteristic temperature and density of the hot component and hence the net mass outflow rate. A number of refinements of the model are investigated, including a simple model of a galactic fountain, the response of the cold component to the pressure of the hot gas, pressure-induced star formation and chemical evolution. The main conclusion of this paper is that low rates of star formation can expel a large fraction of the gas from a dwarf galaxy. For example, a galaxy with circular speed 50 km s¹ can expel 6080 per cent of its gas over a time-scale of 1 Gyr, with a star formation rate that never exceeds 0.1 M yr¹. Effective feedback can therefore take place in a quiescent mode and does not require strong bursts of star formation. Even a large galaxy, such as the Milky Way, might have lost as much as 20 per cent of its mass in a supernova-driven wind. The models developed here suggest that dwarf galaxies at high redshifts will have low average star formation rates and may contain extended gaseous discs of largely unprocessed gas. Such extended gaseous discs might explain the numbers, metallicities and metallicity dispersions of damped Lyman systems.     

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.     

New constraints on a triaxial model of the Galaxy

We determine the most likely values of the free parameters of an N -body model for the Galaxy developed by Fux via a discrete–discrete comparison with the positions on the sky and line-of-sight velocities of an unbiased, homogeneous sample of OH/IR stars. Via Monte Carlo simulation, we find the plausibility of the best-fitting models, as well as the errors on the determined values. The parameters that are constrained best by these projected data are the total mass of the model and the viewing angle of the central bar, although the distribution of the latter has multiple maxima. The other two free parameters, the size of the bar and the (azimuthal) velocity of the Sun, are less well-constrained. The best model has a viewing angle of ∼ 44°, a semimajor axis of 2.5 kpc (corotation radius 4.5 kpc, pattern speed 46 km s⁻¹ kpc⁻¹), a bar mass of 1.7×10¹⁰ M_⊙ and a tangential velocity of the local standard of rest of 171 km s⁻¹. We argue that the lower values that are commonly found from stellar data for the viewing angle (∼25°) arise when too few coordinates are available, when the longitude range is too narrow or when low latitudes are excluded from the fit. The new constraints on the viewing angle of the Galactic bar from stellar line-of-sight velocities decrease further the ability of the distribution of the bar to account for the observed microlensing optical depth toward Baade's window: our model reproduces only half the observed value. The signal of triaxiality diminishes quickly with increasing latitude, fading within approximately 1 scaleheight (≲3°). This suggests that Baade's window is not a very appropriate region in which to sample bar properties.     

A mystery solved: the mass ratio of the dwarf nova EM Cygni

We have discovered that the spectrum of the well-known dwarf nova EM Cyg is contaminated by light from a K25V star (in addition to the K-type mass donor star). The K25V star contributes approximately 16 per cent of the light from the system and if not taken into account has a considerable effect upon radial velocity measurements of the mass donor star. We obtain a new radial velocity amplitude for the mass donor star of K ₂=202±3 km s¹, compared with the value of K ₂=135±3 km s¹ obtained in Stover, Robinson & Nather's classic study of EM Cyg. The revised value of the amplitude, combined with a measurement of rotational broadening of the mass donor, v  sin  i =140±6 km s¹, leads to a new mass ratio of q M ₂ M ₁=0.88±0.05. This solves a long-standing problem with EM Cyg, because Stover et al.'s measurements indicated a mass ratio q >1, a value that should have led to dynamically unstable mass transfer for the secondary mass deduced by Stover et al. The revised value of the mass ratio, combined with the orbital inclination i =67±2°, leads to masses of 0.99±0.12 M and 1.12±0.08 M for the mass donor and white dwarf respectively. The mass donor is evolved, because it has a later spectral type (K3) than its mass would imply.
We discuss whether the K star could be physically associated with EM Cyg or not, and present the results of the spectroscopic study.     

Tidal evolution of discy dwarf galaxies in the Milky Way potential: the formation of dwarf spheroidals

We conduct high-resolution collisionless N -body simulations to investigate the tidal evolution of dwarf galaxies on an eccentric orbit in the Milky Way (MW) potential. The dwarfs originally consist of a low surface brightness stellar disc embedded in a cosmologically motivated dark matter halo. During 10 Gyr of dynamical evolution and after five pericentre passages, the dwarfs suffer substantial mass loss and their stellar component undergoes a major morphological transformation from a disc to a bar and finally to a spheroid. The bar is preserved for most of the time as the angular momentum is transferred outside the galaxy. A dwarf spheroidal (dSph) galaxy is formed via gradual shortening of the bar. This work thus provides a comprehensive quantitative explanation of a potentially crucial morphological transformation mechanism for dwarf galaxies that operates in groups as well as in clusters. We compare three cases with different initial inclinations of the disc and find that the evolution is fastest when the disc is coplanar with the orbit. Despite the strong tidal perturbations and mass loss, the dwarfs remain dark matter dominated. For most of the time, the one-dimensional stellar velocity dispersion, σ, follows the maximum circular velocity, V _max, and they are both good tracers of the bound mass. Specifically, we find that   M _bound∝ V ^3.5_max  and     in agreement with earlier studies based on pure dark matter simulations. The latter relation is based on directly measuring the stellar kinematics of the simulated dwarf, and may thus be reliably used to map the observed stellar velocity dispersions of dSphs to halo circular velocities when addressing the missing satellites problem.