Adaptive optics observations of the flattened nuclear star cluster in the nearby edge-on spiral galaxy NGC 4244 using the Gemini Near-Infrared Integral Field Spectrograph (NIFS) have revealed clear rotation. Using these kinematics plus 2MASS photometry we construct a series of axisymmetric two-component particle dynamical models with our improved version of NMAGIC, a flexible Chi^2-made-to-measure code. The models consist of a nuclear cluster disc embedded within a spheroidal particle population. We find a mass for the nuclear star cluster of M=1.6^+0.5_-0.2 x 10^7 M_sun within ~42.4 pc (2"). We also explore the presence of an intermediate mass black hole and show that models with a black hole as massive as M_bh = 5.0 x 10^5 M_sun are consistent with the available data. Regardless of whether a black hole is present or not, the nuclear cluster is vertically anisotropic (beta_z < 0), as was found with earlier two-integral models. We then use the models as initial conditions for N-body simulations. These simulations show that the nuclear star cluster is stable against non-axisymmetric perturbations. We also explore the effect of the nuclear cluster accreting star clusters at various inclinations. Accretion of a star cluster with mass 13% that of the nuclear cluster is already enough to destroy the vertical anisotropy, regardless of orbital inclination.

We create dynamical models of the massive elliptical galaxy, NGC 4649, using the N-body made-to-measure code, NMAGIC, and kinematic constraints from long-slit and planetary nebula (PN) data. We explore a range of potentials based on previous determinations from X-ray observations and a dynamical model fitting globular cluster (GC) velocities and a stellar density profile. The X-ray mass distributions are similar in the central region but have varying outer slopes, while the GC mass profile is higher in the central region and on the upper end of the range further out. Our models cannot differentiate between the potentials in the central region, and therefore if non-thermal pressures or multi-phase components are present in the hot gas, they must be smaller than previously inferred. In the halo, we find that the PN velocities are sensitive tracers of the mass, preferring a less massive halo than that derived from the GC mass profile, but similar to one of the mass distributions derived from X-rays. Our results show that the GCs may form a dynamically distinct system, and that the properties of the hot gas derived from X-rays in the outer halo have considerable uncertainties that need to be better understood. Estimating the mass in stars using photometric information and a stellar population mass-to-light ratio, we infer a dark matter mass fraction in NGC 4649 of ~0.39 at 1Re (10.5 kpc) and ~0.78 at 4Re. We find that the stellar orbits are isotropic to mildly radial in the central ~6 kpc depending on the potential assumed. Further out, the orbital structure becomes slightly more radial along R and more isotropic along z, regardless of the potential assumed. In the equatorial plane, azimuthal velocity dispersions dominate over meridional velocity dispersions, implying that meridional velocity anisotropy is the mechanism for flattening the stellar system.

We present new planetary nebulae (PNe) positions, radial velocities, and magnitudes for 6 early-type galaxies obtained with the Planetary Nebulae Spectrograph, their two-dimensional velocity and velocity dispersion fields. We extend this study to include an additional 10 early-type galaxies with PNe radial velocity measurements available from the literature, to obtain a broader description of the outer-halo kinematics in early-type galaxies. These data extend the information derived from stellar kinematics to typically up to ~8 Re. The combination of photometry, stellar and PNe kinematics shows: i) good agreement between the PNe number density and the stellar surface brightness in the region where the two data sets overlap; ii) good agreement between PNe and stellar kinematics; iii) that the mean rms velocity profiles fall into two groups: with of the galaxies characterized by slowly decreasing profiles and the remainder having steeply falling profiles; iv) a larger variety of velocity dispersion profiles; v) that twists and misalignments in the velocity fields are more frequent at large radii, including some fast rotators; vi) that outer haloes are characterised by more complex radial profiles of the specific angular momentum-related lambda_R parameter than observed within 1Re; vii) that many objects are more rotationally dominated at large radii than in their central parts; and viii) that the halo kinematics are correlated with other galaxy properties, such as total luminosity, isophotal shape, total stellar mass, V/sigma, and alpha parameter, with a clear separation between fast and slow rotators.

We present first results of a study of the halo kinematics for a sample of early type galaxies using planetary nebulae (PNe) as kinematical tracers. PNe allow to extend up to several effective radii (Re) the information from absorption line kinematics (confined to within 1 or 2 Re), providing valuable information and constraints for merger simulations and galaxy formation models. We find that the specific angular momentum per unit mass has a more complex radial dependence when the halo region is taken into account and that the halo velocity dispersion is related to the total galaxy luminosity, isophotal shape, and number of PNe per unit of luminosity

We present new Planetary Nebula Spectrograph observations of the ordinary elliptical galaxy NGC 4494, resulting in positions and velocities of 255 PNe out to 7 effective radii (25 kpc). We also present new wide-field surface photometry from MMT/Megacam, and long-slit stellar kinematics from VLT/FORS2. The spatial and kinematical distributions of the PNe agree with the field stars in the region of overlap. The mean rotation is relatively low, with a possible kinematic axis twist outside 1 Re. The velocity dispersion profile declines with radius, though not very steeply, down to ~70 km/s at the last data point. We have constructed spherical dynamical models of the system, including Jeans analyses with multi-component LCDM-motivated galaxies as well as logarithmic potentials. These models include special attention to orbital anisotropy, which we constrain using fourth-order velocity moments. Given several different sets of modelling methods and assumptions, we find consistent results for the mass profile within the radial range constrained by the data. Some dark matter (DM) is required by the data; our best-fit solution has a radially anisotropic stellar halo, a plausible stellar mass-to-light ratio, and a DM halo with an unexpectedly low central density. We find that this result does not substantially change with a flattened axisymmetric model. Taken together with other results for galaxy halo masses, we find suggestions for a puzzling pattern wherein most intermediate-luminosity galaxies have very low concentration halos, while some high-mass ellipticals have very high concentrations. We discuss some possible implications of these results for DM and galaxy formation.

Recent results from the Planetary Nebula Spectrograph (PN.S) survey have revealed a rapidly falling velocity dispersion profile in the nearby elliptical galaxy NGC 3379, casting doubts on whether this intermediate-luminosity galaxy has the kind of dark matter halo expected in LambdaCDM cosmology. We present a detailed dynamical study of this galaxy, combining long-slit spectroscopy, SAURON integral-field data, and PN.S velocities, reaching to more than seven effective radii (R_e). We construct spherical and axisymmetric dynamical models for these data with the flexible made-to-measure NMAGIC code, in a sequence of gravitational potentials with varying dark halo mass. We find that the data are consistent both with near-isotropic spherical systems dominated by the stellar mass, and with models in massive halos with strongly radially anisotropic outer parts (beta >~ 0.8 at 7R_e). Formal likelihood limits would exclude (at 1 sigma) the model with stars only, as well as halo models with v_circ(7R_e) >~ 250 km/s. A sequence of more realistic axisymmetric models of different inclinations and a small number of triaxial tests confirm the spherical results. All valid models fitting all the data are dynamically stable over Gyrs, including the most anisotropic ones. NGC 3379 may well have a dark matter halo as predicted by recent merger models within LambdaCDM cosmology, provided its outer envelope is strongly radially anisotropic. (abridged)

We present a dynamical study of NGC 4697, an almost edge-on, intermediate-luminosity, E4 elliptical galaxy, combining new surface brightness photometry, new as well as published long-slit absorption line kinematic data, and published planetary nebulae (PNe) velocity data. The combined kinematic data set extends out to ~= 5' ~= 4.5 R_e and allows us to probe the galaxy's outer halo. For the first time, we model such a dataset with the new and flexible Chi^2-made-to-measure particle code NMAGIC. We extend NMAGIC to include seeing effects, introduce an efficient scheme to estimate the mass-to-light ratio, and incorporate a maximum likelihood technique to account for discrete velocity measurements. For modelling the PNe kinematics we use line-of-sight velocities and velocity dispersions computed on two different spatial grids, and we also use the individual velocity measurements with the likelihood method, in order to make sure that our results are not biased by the way we treat the PNe measurements. We generate axisymmetric self-consistent models as well models including various dark matter halos. These models fit all the mean velocity and velocity dispersion data with Chi^2/N<1, both in the case with only luminous matter and in potentials including quite massive halos. The likelihood analysis together with the velocity histograms suggest that models with low density halos such that the circular velocity v_c ~< 200 km/s at R_e are not consistent with the data. A range of massive halos with v_c ~= 250 km/s at 5R_e fit the PN data best. To derive stronger results would require PN velocities at even larger radii. The best fitting models are slightly radially anisotropic; the anisotropy parameter beta ~= 0.3 at the center, increasing to beta ~= 0.5 at radii >~ 2R_e.

We investigate the manner in which lenticular galaxies are formed by studying their stellar kinematics: an S0 formed from a fading spiral galaxy should display similar cold outer disc kinematics to its progenitor, while an S0 formed in a minor merger should be more dominated by random motions. In a pilot study to attempt to distinguish between these scenarios, we have measured the planetary nebula (PN) kinematics of the nearby S0 system NGC 1023. Using the Planetary Nebula Spectrograph, we have detected and measured the line-of-sight velocities of 204 candidate PNe in the field of this galaxy. Out to intermediate radii, the system displays the kinematics of a normal rotationally-supported disc system. After correction of its rotational velocities for asymmetric drift, the galaxy lies just below the spiral galaxy Tully-Fisher relation, as one would expect for a fading system. However, at larger radii the kinematics undergo a gradual but major transition to random motion with little rotation. This transition does not seem to reflect a change in the viewing geometry or the presence of a distinct halo component, since the number counts of PNe follow the same simple exponential decline as the stellar continuum with the same projected disc ellipticity out to large radii. The galaxy's small companion, NGC 1023A, does not seem to be large enough to have caused the observed modification either. This combination of properties would seem to indicate a complex evolutionary history in either the transition to form an S0 or in the past life of the spiral galaxy from which the S0 formed. More data sets of this type from both spirals and S0s are needed in order to definitively determine the relationship between these types of system.

We examine the dark matter properties of nearby early-type galaxies using planetary nebulae (PNe) as mass probes. We have designed a specialised instrument, the Planetary Nebula Spectrograph (PN.S) operating at the William Herschel telescope, with the purpose of measuring PN velocities with best efficiency. The primary scientific objective of this custom-built instrument is the study of the PN kinematics in 12 ordinary round galaxies. Preliminary results showing a dearth of dark matter in ordinary galaxies (Romanowsky et al. 2003) are now confirmed by the first complete PN.S datasets. On the other hand early-type galaxies with a "regular" dark matter content are starting to be observed among the brighter PN.S target sample, thus confirming a correlation between the global dark-to-luminous mass virial ratio (f_DM=M_DM/M_star) and the galaxy luminosity and mass.

Red clump giants in the Galactic bulge are approximate standard candles and hence they can be used as distance indicators. We compute the proper motion dispersions of RCG stars in the Galactic bulge using the proper motion catalogue from the second phase of the Optical Gravitational Microlensing Experiment (OGLE-II, Sumi et al. 2004) for 45 fields. The proper motion dispersions are measured to a few per cent accuracy due to the large number of stars in the fields. The observational sample is comprised of 577736 stars. These observed data are compared to a state-of-the-art particle simulation of the Galactic bulge region. The predictions are in rough agreement with observations, but appear to be too anisotropic in the velocity ellipsoid. We note that there is significant field-to-field variation in the observed proper motion dispersions. This could either be a real feature, or due to some unknown systematic effect.

We describe a made-to-measure algorithm for constructing N-particle models of stellar systems from observational data (Chi-Squared-M2M), extending earlier ideas by Syer and Tremaine. The algorithm properly accounts for observational errors, is flexible, and can be applied to various systems and geometries. We implement this algorithm in a parallel code NMAGIC and carry out a sequence of tests to illustrate its power and performance: (i) We reconstruct an isotropic Hernquist model from density moments and projected kinematics and recover the correct differential energy distribution and intrinsic kinematics. (ii) We build a self-consistent oblate three-integral maximum rotator model and compare how the distribution function is recovered from integral field and slit kinematic data. (iii) We create a non-rotating and a figure rotating triaxial stellar particle model, reproduce the projected kinematics of the figure rotating system by a non-rotating system of the same intrinsic shape, and illustrate the signature of pattern rotation in this model. From these tests we comment on the dependence of the results from Chi-Squared-M2M on the initial model, the geometry, and the amount of available data.