We present a detailed analysis of 108 helium-line (DB) white dwarfs based on model atmosphere fits to high signal-to-noise optical spectroscopy. We derive a mean mass of 0.67 Mo for our sample, with a dispersion of only 0.09 Mo. White dwarfs also showing hydrogen lines, the DBA stars, comprise 44% of our sample, and their mass distribution appears similar to that of DB stars. As in our previous investigation, we find no evidence for the existence of low-mass (M < 0.5 Mo) DB white dwarfs. We derive a luminosity function based on a subset of DB white dwarfs identified in the Palomar-Green survey. We show that 20% of all white dwarfs in the temperature range of interest are DB stars, although the fraction drops to half this value above Teff ~ 20,000 K. We also show that the persistence of DB stars with no hydrogen features at low temperatures is difficult to reconcile with a scenario involving accretion from the interstellar medium, often invoked to account for the observed hydrogen abundances in DBA stars. We present evidence for the existence of two different evolutionary channels that produce DB white dwarfs: the standard model where DA stars are transformed into DB stars through the convective dilution of a thin hydrogen layer, and a second channel where DB stars retain a helium-atmosphere throughout their evolution. We finally demonstrate that the instability strip of pulsating V777 Her white dwarfs contains no nonvariables, if the hydrogen content of these stars is properly accounted for.

We present new analytical distribution functions for anisotropic spherical galaxies. They have the density profiles of the gamma-models, which allow a wide range of central density slopes, and are widely used to fit elliptical galaxies and the bulges of spiral galaxies. Most of our models belong to two two-parameter families. One of these parameters is the slope gamma of the central density cusp. The other allows a wide range of varying radial and tangential anisotropies, at either small or large radii. We give analytical formulas for their distribution functions, velocity dispersions, and the manner in which energy and transverse velocity are distributed between orbits. We also give some of their observable properties, including line-of-sight velocity profiles which have been computed numerically. Our models can be used to provide a useful tool for creating initial conditions for N-body and Monte Carlo simulations.

We study bisymmetric modes of angular wavenumber 2 for flat stellar disks in potentials with smooth cores. Stars either all circulate in the same direction or a small fraction may counter-rotate. The bisymmetric modes are unstable unless there is a sufficiently large external halo or bulge. We find two modes for each disk: a more central fundamental mode and a more extensive and more spiral (trailing) secondary mode. The fundamental mode is particularly sensitive to the orbital population in the central part of the disk. Depending on that population, it varies from a small compact bar to a trailing spiral that is almost as wound as the secondary mode. All modes rotate too rapidly for there to be an inner Lindblad resonance. They transfer angular momentum from the central to the outer regions of the disk. Most of them release gravitational energy and convert it to kinetic energy, which also flows outwards through the disk.

The Jeans equations relate the second-order velocity moments to the density and potential of a stellar system. For general three-dimensional stellar systems, there are three equations and six independent moments. By assuming that the potential is triaxial and of separable Staeckel form, the mixed moments vanish in confocal ellipsoidal coordinates. Consequently, the three Jeans equations and three remaining non-vanishing moments form a closed system of three highly-symmetric coupled first-order partial differential equations in three variables. These equations were first derived by Lynden-Bell, over 40 years ago, but have resisted solution by standard methods. We present the general solution here. We consider the two-dimensional limiting cases first. We solve their Jeans equations by a new method which superposes singular solutions. The singular solutions, which are new, are standard Riemann-Green functions. The two-dimensional solutions are applied to non-axisymmetric discs, oblate and prolate spheroids, and also to the scale-free triaxial limit. We then extend the method of singular solutions to the triaxial case, and obtain a full solution. The general solution can be expressed in terms of complete (hyper)elliptic integrals which can be evaluated in a straightforward way, and provides the full set of second moments which can support a triaxial density distribution in a separable triaxial potential. (abridged)

The Jeans equations relate the second-order velocity moments to the density and potential of a stellar system. For general three-dimensional stellar systems, there are three equations, but these are not very helpful, as they contain six independent moments. By assuming that the potential is triaxial and of separable Staeckel form, the mixed moments vanish in confocal ellipsoidal coordinates. The three Jeans equations and three remaining non-vanishing moments form a closed system of three highly-symmetric coupled first-order partial differential equations in three variables. They were first derived by Lynden-Bell in 1960, but have resisted solution by standard methods. Here we present the general solution by superposition of singular solutions.

A method is developed to evaluate the magnifications of the images of galaxies with lensing potentials stratified on similar concentric ellipses. A simple contour integral is provided which enables the sums of the magnifications of even parity or odd parity or the central image to be easily calculated. The sums for pairs of images vary considerably with source position, while the signed sums can be remarkably uniform inside the tangential caustic in the absence of naked cusps. For a family of models in which the potential is a power-law of the elliptic radius, the number of visible images is found as a function of flattening, external shear and core radius. The magnification of the central image depends on the core radius and the slope of the potential. For typical source and lens redshifts, the missing central image leads to strong constraints; the mass distribution in the lensing galaxy must be nearly cusped, and the cusp must be isothermal or stronger. This is in accord with the cuspy cores seen in high resolution photometry of nearby, massive, early-type galaxies, which typically have the surface density falling like distance^{-1.3} outside a break radius of a few hundred parsecs. Cuspy cores by themselves can provide an explanation of the missing central images. Dark matter at large radii may alter the slope of the projected density; provided the slope remains isothermal or steeper and the break radius remains small, then the central image remains unobservable. The sensitivity of the radio maps must be increased fifty-fold to find the central images in abundance.

The multiple images of lensed quasars provide evidence on the mass distribution of the lensing galaxy. The lensing invariants are constructed from the positions of the images, their parities and their fluxes. They depend only on the structure of the lensing potential. The simplest is the magnification invariant, which is the sum of the signed magnifications of the images. Higher order invariants are the sums of products of the signed magnifications with positive or negative powers of the position coordinates of the images. We consider the case of the four and five image systems produced by elliptical power-law galaxies. This paper provides simple contour integrals for evaluating all their lensing invariants. For practical evaluation, this offers considerable advantages over the algebraic methods used previously. To illustrate the power of the contour integral method, we calculate full expansions in the source offset for all lensing invariants in the presence of arbitrary external shear. As an example, we use the elliptical power-law galaxies to fit to the data on the four images of the Einstein Cross (G2237+030). The lensing invariants play a role by reducing the dimensionality of the parameter space in which the chi-squared minimisation proceeds with consequent gains in accuracy and speed.

The contour integral method of Hunter & Qian is applied to axisymmetric galaxy models in which the distribution function (DF) is of the form f=f(E,L_z), where E and L_z are the classical integrals of motion in an axisymmetric potential. A practical way to construct the unique even part of the DF for such systems is presented. It is applied to models, both oblate and prolate, in which the mass density is stratified on similar concentric spheroids. The spheroids with scale-free densities are discussed in detail. These provide useful approximations to the behaviour of more realistic models in the limit of small and large radii. The self-consistent case is treated, as well as the case in which there are additional contributions to the potential from a central black hole or dark halo. The range of axis ratios and density profile slopes is determined for which spheroidal power--law cusps with a central black hole have a physical two--integral DF. More generally, the f(E,L_z) DFs are discussed for a set of spheroidal `(alpha,beta)-models', characterized by a power-law density cusp with slope alpha at small radii, and a power-law density fall-off with slope [alpha + 2 beta] at large radii. As an application, the DF is constructed for the (alpha,beta) model with a 1.8 x 10^6 solar mass black hole used by van der Marel et al. to interpret their high spatial resolution spectroscopic data for M32. The results confirm that the model fits the observed line-of-sight velocity profiles remarkably well. The model is used to calculate predictions for future spectroscopic observations with the HST.