We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

We used the PMAS integral field spectrograph to obtain large sets of radial velocities in the central regions of three northern Galactic globular clusters: M3, M13, and M92. By applying the novel technique of crowded field 3D spectroscopy, we measured radial velocities for about 80 stars within the central ~ 10 arcsec of each cluster. These are by far the largest spectroscopic datasets obtained in the innermost parts of these clusters up to now. To obtain kinematical data across the whole extent of the clusters, we complement our data with measurements available in the literature. We combine our velocity measurements with surface brightness profiles to analyse the internal dynamics of each cluster using spherical Jeans models, and investigate whether our data provide evidence for an intermediate-mass black hole in any of the clusters. The surface brightness profiles reveal that all three clusters are consistent with a core profile, although shallow cusps cannot be excluded. We find that spherical Jeans models with a constant mass-to-light ratio provide a good overall representation of the kinematical data. A massive black hole is required in none of the three clusters to explain the observed kinematics. Our 1sigma (3sigma) upper limits are 5300 M_sun (12000 M_sun) for M3, 8600 M_sun (13000 M_sun) for M13, and 980 M_sun (2700 M_sun) for M92. A puzzling circumstance is the existence of several potential high velocity stars in M3 and M13, as their presence can account for the majority of the discrepancies that we find in our mass limits compared to M92.

We measure the shape of the velocity ellipsoid in two late-type spiral galaxies (Hubble types Sc and Scd) and combine these results with our previous analyses of six early-type spirals (Sa to Sbc) to probe the relation between galaxy morphology and the ratio of the vertical and radial dispersions. We confirm at much higher significance (99.9 percent) our prior detection of a tight correlation between these quantities. We explore the trends of the magnitude and shape of the velocity ellipsoid axes with galaxy properties (colour, gas surface mass density, and spiral arm structure). The observed relationships allow for an observational identification of the radial and vertical disk heating agents in external disk galaxies.

Using the Inamori Magellan Areal Camera and Spectrograph (IMACS) integral-field unit (IFU) on the 6.5m Magellan telescope, we have designed the first statistically significant investigation of the two-dimensional distribution and kinematics of ionized gas and stars in the central kiloparsec regions of a well-matched sample of Seyfert and inactive control galaxies selected from the Sloan Digital Sky Survey. The goals of the project are to use the fine spatial sampling (0.2 arcsec/pixel) and large wavelength coverage (4000-7000A) of the IMACS-IFU to search for dynamical triggers of nuclear activity in the central region where active galactic nucleus (AGN) activity and dynamical timescales become comparable, to identify and assess the impact of AGN-driven outflows on the host galaxy and to provide a definitive sample of local galaxy kinematics for comparison with future three-dimensional kinematic studies of high-redshift systems. In this paper, we provide the first detailed description of the procedure to reduce and calibrate data from the IMACS-IFU in `long mode' to obtain two-dimensional maps of the distribution and kinematics of ionized gas and stars. The sample selection criteria are presented, observing strategy described and resulting maps of the sample galaxies presented along with a description of the observed properties of each galaxy and the overall observed properties of the sample.

We have used the VIMOS integral field spectrograph to map the emission line properties in a sample of 24 star forming galaxies selected from the SDSS database. In this data paper we present and describe the sample, and explore some basic properties of SDSS galaxies with resolved emission line fields. We fit the Halpha+[NII] emission lines in each spectrum to derive maps of continuum, Halpha flux, velocity and velocity dispersion. The Halpha, Hbeta, [NII] and [OIII] emission lines are also fit in summed spectra for circular annuli of increasing radius. A simple mass model is used to estimate dynamical mass within 10 kpc, which compared to estimates of stellar mass shows that between 10 and 100% of total mass is in stars. We present plots showing the radial behaviour of EW[Halpha], u-i colour and emission line ratios. Although EW[Halpha] and u-i colour trace current or recent star formation, the radial profiles are often quite different. Whilst line ratios do vary with annular radius, radial gradients in galaxies with central line ratios typical of AGN or LINERS are mild, with a hard component of ionization required out to large radii. We use our VIMOS maps to quantify the fraction of Halpha emission contained within the SDSS fibre, taking the ratio of total Halpha flux to that of a simulated SDSS fibre. A comparison of the flux ratios to colour-based SDSS extrapolations shows a 175% dispersion in the ratio of estimated to actual corrections in normal star forming galaxies, with larger errors in galaxies containing AGN. We find a strong correlation between indicators of nuclear activity: galaxies with AGN-like line ratios and/or radio emission frequently show enhanced dispersion peaks in their cores, requiring non-thermal sources of heating. Altogether, about half of the galaxies in our sample show no evidence for nuclear activity or non-thermal heating.

The reduction of integral-field spectrograph (IFS) data is demanding work. Many repetitive operations are required in order to convert raw data into, typically a large number of, spectra. This effort can be markedly simplified through the use of a tool or pipeline, which is designed to complete many of the repetitive operations without human interaction. Here we present our semi-automatic data-reduction tool p3d that is designed to be used with fiber-fed IFSs. Important components of p3d include a novel algorithm for automatic finding and tracing of spectra on the detector, and two methods of optimal spectrum extraction in addition to standard aperture extraction. p3d also provides tools to combine several images, perform wavelength calibration and flat field data. p3d is at the moment configured for four IFSs. In order to evaluate its performance we have tested the different components of the tool. For these tests we used both simulated and observational data. We demonstrate that for three of the IFSs a correction for so-called cross-talk due to overlapping spectra on the detector is required. Without such a correction spectra will be inaccurate, in particular if there is a significant intensity gradient across the object. Our tests showed that p3d is able to produce accurate results. p3d is a highly general and freely available tool. It is easily extended to include improved algorithms, new visualization tools and support for additional instruments. The program code can be downloaded from the p3d-project web site http://p3d.sourceforge.net

We investigate the effect of dust on the observed rotation rate of a stellar bar as measured by application of the Tremaine & Weinberg (TW) method. As the method requires that the tracer satisfies the continuity equation, it has been applied largely to early-type barred galaxies. We show using numerical simulations of barred galaxies that dust attenuation factors typically found in these systems change the observed bar pattern speed by 20-40%. We also address the effect of star formation on the TW method and find that it does not change the results significantly. This suggests that applications of the TW method can be extended to include barred galaxies covering the full range of Hubble type.

MUSE, a 2nd generation VLT instrument, will become the world's largest integral field spectrograph. It will be an AO assisted instrument which, in a single exposure, covers the wavelength range from 465 to 930 nm with an average resolution of 3000 over a field of view of 1'x1' with 0.2'' spatial sampling. Both the complexity and the rate of the data are a challenge for the data processing of this instrument. We will give an overview of the data processing scheme that has been designed for MUSE. Specifically, we will use only a single resampling step from the raw data to the reduced data product. This allows us to improve data quality, accurately propagate variance, and minimize spreading of artifacts and correlated noise. This approach necessitates changes to the standard way in which reduction steps like wavelength calibration and sky subtraction are carried out, but can be expanded to include combination of multiple exposures.

We have used the VIMOS IFU to map the properties of the Seyfert 1.9 galaxy LEDA 135736. These maps reveal a number of interesting features including: an Extended Narrow Line Region detectable out to 9 kpc, an area of intense star formation located at a projected distance of 12 kpc from the centre, an elliptical companion galaxy, and kinematic features, aligned along the long-axis of the ENLR, that are consistent with radio jet-driven mass outflow. We propose that the ENLR results from extra-planar gas ionized by the AGN, and that the AGN in turn might be triggered by interaction with the companion galaxy, which can also explain the burst of star formation and morphological features. Only about two percent of the ENLR's kinetic energy is in the mass outflow. We infer from this that the bulk of mechanical energy imparted by the jet is used to heat this gas.

We investigate the effect of dust on the observed rotation rate of a stellar bar. The only direct way to measure this quantity relies on the Tremaine & Weinberg method which requires that the tracer satisfies the continuity equation. Thus it has been applied largely to early-type barred galaxies. We show using numerical simulations of barred galaxies that dust attenuation factors typically found in these systems change the observed bar pattern speed by 20-40 percent. We also address the effect of star formation on the TW method and find that it does not change the results significantly. The results presented here suggest that applications of the TW method can be extended to include barred galaxies covering the full range of Hubble type.

The Multi Unit Spectroscopic Explorer (MUSE) is a second-generation VLT panoramic integral-field spectrograph under preliminary design study. MUSE has a field of 1x1 arcmin**2 sampled at 0.2x0.2 arcsec**2 and is assisted by the VLT ground layer adaptive optics ESO facility using four laser guide stars. The simultaneous spectral range is 465-930 nm, at a resolution of R~3000. MUSE couples the discovery potential of a large imaging device to the measuring capabilities of a high-quality spectrograph, while taking advantage of the increased spatial resolution provided by adaptive optics. This makes MUSE a unique and tremendously powerful instrument for discovering and characterizing objects that lie beyond the reach of even the deepest imaging surveys. MUSE has also a high spatial resolution mode with 7.5x7.5 arcsec**2 field of view sampled at 25 milli-arcsec. In this mode MUSE should be able to obtain diffraction limited data-cubes in the 600-930 nm wavelength range. Although the MUSE design has been optimized for the study of galaxy formation and evolution, it has a wide range of possible applications; e.g. monitoring of outer planets atmosphere, environment of young stellar objects, super massive black holes and active nuclei in nearby galaxies or massive spectroscopic surveys of stellar fields in the Milky Way and nearby galaxies.

(Abridged) We study the nuclear star clusters in spiral galaxies of various Hubble types using spectra obtained with STIS on-board HST. We observed the nuclear clusters in 40 galaxies, selected from two previous HST/WFPC2 imaging surveys. The spectra provide a better separation of cluster light from underlying galaxy light than is possible with ground-based spectra. To infer the star formation history, metallicity and dust extinction, we fit weighted superpositions of single-age stellar population templates to the spectra. The luminosity-weighted age ranges from 10 Myrs to 10 Gyrs. The stellar populations of NCs are generally best fit as a mixture of populations of different ages. This indicates that NCs did not form in a single event, but instead they had additional star formation long after the oldest stars formed. On average, the sample clusters in late-type spirals have a younger luminosity-weighted mean age than those in early-type spirals (log(age/yr) = 8.37+/-0.25 vs. 9.23+/-0.21). The average cluster masses are smaller in late-type spirals than in early-type spirals (log(M/Msun) = 6.25+/-0.21 vs. 7.63+/-0.24), and exceed the masses typical of globular clusters. The cluster mass correlates strongly with both the Hubble type of the host galaxy and the luminosity of its bulge. The latter correlation has the same slope as the well-known correlation between supermassive black hole mass and bulge luminosity. The properties of both nuclear clusters and black holes are therefore intimately connected to the properties of the host galaxy.

We present a datacube covering the central 10 arcsec of the archetypal active galaxy NGC 1068 over a wavelength range 4200 to 5400 Angstrom. obtained during the commissioning of the integral field unit (IFU) of the Gemini Multiobject Spectrograph (GMOS) installed on the Gemini-North telescope. The datacube shows a complex emission line morphology in the [OIII] doublet and Hbeta line. To describe this structure phenomenologically we construct an atlas of velocity components derived from multiple Gaussian component fits to the emission lines. The atlas contains many features which cannot be readily associated with distinct physical structures. While some components are likely to be associated with the expected biconical outflow, others are suggestive of high velocity flows or disk-like structures. As a first step towards interpretation, we seek to identify the stellar disk using kinematical maps derived from the Mgb absorption line feature at 5170 Angstrom and make associations between this and gaseous components in the atlas of emission line components.

We review the science case for the Laser Guide Star system being built for the William Herschel Telescope (WHT) on La Palma. When used in combination with the NAOMI Adaptive Optics system and the OASIS visible-wavelength Integral Field Spectrograph, we demonstrate that there are substantial, exciting areas of astrophysical research in which the WHT can contribute.

High-velocity galactic outflows, driven by intense bursts of star formation and black hole accretion, are invoked by current theories of galaxy formation to terminate star formation in the most massive galaxies and to deposit heavy elements in the intergalactic medium. From existing observational evidence on high-redshift galaxies, it is unclear whether such outflows are localized to regions of intense star formation just a few kiloparsecs in extent, or whether they instead have a significant impact on the entire galaxy and its surroundings. Here we present two-dimensional spectroscopy of a star-forming galaxy at redshift z=3.09 (seen 11.5 Gyr ago, when the Universe was 20 per cent of its current age): its spatially extended Ly-alpha emission appears to be absorbed by HI in a foreground screen covering the entire galaxy, with a lateral extent of at least 100 kpc and remarkable velocity coherence. It was plausibly ejected from the galaxy during a starburst several 1E8 yr earlier and has subsequently swept up gas from the surrounding intergalactic medium and cooled. This demonstrates the galaxy-wide impact of high-redshift superwinds.

We investigate the central regions of 23 spiral galaxies using archival NICMOS imaging and STIS spectroscopy. The sample is taken from our program to determine the masses of central massive black holes (MBH) in 54 nearby spiral galaxies. Stars are likely to contribute significantly to any dynamical central mass concentration that we find in our MBH program and this paper is part of a series to investigate the nuclear properties of these galaxies. We use the Nuker law to fit surface brightness profiles, derived from the NICMOS images, to look for nuclear star clusters and find possible extended sources in 3 of the 23 galaxies studied (13 per cent). The fact that this fraction is lower than that inferred from optical Hubble Space Telescope studies is probably due to the greater spatial resolution of those studies. Using R-H and J-H colors and equivalent widths of H-alpha emission (from the STIS spectra) we investigate the nature of the stellar population with evolutionary models. Under the assumption of hot stars ionizing the gas, as opposed to a weak AGN, we find that there are young stellar populations (~10-20 Myr) however these data do not allow us to determine what percentage of the total nuclear stellar population they form. Also, in an attempt to find any unknown AGN we use [N II] and [S II] line flux ratios (relative to H-alpha) and find tentative evidence for weak AGN in NGC 1300 and NGC 4536.

We present Space Telescope Imaging Spectrograph emission-line spectra of the central regions of the spiral galaxies NGC 1300 and NGC 2748. From the derived kinematics of the nuclear gas we have found evidence for central supermassive black holes in both galaxies. The estimated mass of the black hole in NGC 1300 is 6.6 (+6.3, -3.2) x 10^7 solar masses and in NGC 2748 is 4.4 (+3.5, -3.6) x 10^7 solar masses (both at the 95% confidence level). These two black hole mass estimates contribute to the poorly sampled low-mass end of the nuclear black hole mass spectrum.

We present surface photometry for the central regions of a sample of 48 spiral galaxies (mostly unbarred and barred of types Sbc or Sc) observed with the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope. Surface brightness profiles were derived and modeled with a Nuker law. We also analyzed archival Wide Field Planetary Camera~2 images with a larger field of view, available for 18 galaxies in our sample. We modeled the extracted bulge surface brightness profiles with an exponential, a de Vaucouleurs or a Sersic profile. In agreement with previous studies, we find that bulges of Sbc galaxies fall into two categories: bulges well described by an exponential profile and those well described by an de Vaucouleurs profile. Only one galaxy requires the use of a more general Sersic profile to properly describe the bulge. Nuclear photometrically distinct components are found in ~55% of the galaxies. For those that we classify as star clusters based on their resolved extent we find absolute magnitudes that are brighter on average than those previously identified in spiral galaxies. This might be due to a bias in our sample toward star forming galaxies, combined with a trend for star forming galaxies to host brighter central clusters.

We present results from an HST study of the morphology and kinematics of NGC 6240. This merging galaxy with a double nucleus is one of the nearest and best-studied ultraluminous infrared galaxies. HST resolves both nuclei into seperate components. The distance between the northern and southern optical/near-infrared components is greater than that observed in radio and X-ray studies, arguing that even in K-band we may not be seeing all the way through the dust to the true nuclei. The ionized gas does not display rotation around either of the nuclei, or equilibrium motion in general. There is a strong velocity gradient between the nuclei, similar to what is seen in CO data. There is no such gradient in our stellar kinematics. The velocity dispersion of the gas is larger than expected for a cold disk. We also map and model the emission-line velocity field at an off-nuclear position where a steep velocity gradient was previously detected in ground-based data. Overall, the data indicate that line-of-sight projection effects, dust absorption, non-equilibrium merger dynamics, and the possible influence of a wind may be playing an important role in the observed kinematics. Chandra observations of hard X-rays have shown that both of the nuclei contain an Active Galactic Nucleus (AGN). The HST data show no clear sign of the two AGNs: neither continuum nor narrow-band imaging shows evidence for unresolved components in the nuclei, and there are no increased emission line widths or rapid rotation near the nuclei. This underscores the importance of X-ray data for identifying AGNs in highly dust-enshrouded environments.

The pattern speed in one of the fundamental parameters that determines the structure of barred galaxies. This quantity is usually derived from indirect methods or by employing model assumptions. The number of bar pattern speeds derived using the model independent Tremaine and Weinberg technique is still very limited. We present the results of model independent measurements of the bar pattern speed in four galaxies ranging in Hubble type from SB0 to SBbc. Three of the four galaxies in our sample are consistent with bars being fast rotators. The lack of slow bars is consistent with previous observations and suggests that barred galaxies do not have centrally concentrated dark matter halos. This contradicts simulations of cosmological structure formation and observations of the central mass concentration in nonbarred galaxies.

Current understanding of the secular evolution of galactic disks suggests that this process is dominated by two or more heating mechanisms, which increase the random motions of stars in the disk. In particular, the gravitational influence of giant molecular clouds and irregularities in the spiral potential have been proposed to explain the observed velocity dispersions in the solar neighborhood. Each of these mechanisms acts on different components of the stellar velocities, which affects the ratio of the vertical and radial components of the stellar velocity dispersion since the relative strengths of giant molecular clouds and spiral irregularities vary with Hubble type. A study of this ratio as function of Hubble type has the potential to provide strong constraints on disk heating mechanisms. We present major and minor axis stellar kinematics for four spiral galaxies of Hubble type from Sa to Sbc, and use the data to infer the ratios sigma_z/sigma_R in the galaxy disks. The results combined with two galaxies studied previously and with Milky Way data show that the ratio is generally in the range 0.5 - 0.8. There is a marginally significant trend of decreasing ratio with advancing Hubble type, consistent with the predictions of disk heating theories. However, the errors on individual measurements are large, and the absence of any trend is consistent with the data at the 1-sigma level. As a byproduct of our study, we find that three of the four galaxies in our sample have a central drop in their stellar line-of-sight velocity dispersion, a phenomenon that is increasingly observed in spiral galaxies. [ABRIDGED]

We present HST/STIS spectra of the Sbc spiral galaxy NGC 4041 which were used to map the velocity field of the gas in its nuclear region. We detect the presence of a compact (r~0.4" ~40 pc), high surface brightness, rotating nuclear disk co-spatial with a nuclear star cluster. The disk is characterized by a rotation curve with a peak to peak amplitude of ~40 km/s and is systematically blueshifted by 10 - 20 km/s with respect to the galaxy systemic velocity. With the standard assumption of constant mass-to-light ratio and with the nuclear disk inclination taken from the outer disk, we find that a dark point mass of 1(-0.7;+0.6) 10^7 Msun is needed to reproduce the observed rotation curve. However the observed blueshift suggests the possibility that the nuclear disk could be dynamically decoupled. Following this line of reasoning we relax the standard assumptions and find that the kinematical data can be accounted for by the stellar mass provided that either the central mass-to-light ratio is increased by a factor of ~2 or that the inclination is allowed to vary. This model results in a 3 sigma upper limit of 6 10^6 Msun on the mass of any nuclear black hole. Overall, our analysis only allows us to set an upper limit of 2 10^7 Msun on the mass of the nuclear black hole. If this upper limit is taken in conjunction with an estimated bulge B magnitude of -17.7 and with a central stellar velocity dispersion of 95 km/s, then these results are not inconsistent with both the MBH-Lsph and the MBH-sigma correlations. Constraints on black hole masses in spiral galaxies of types as late as Sbc are still very scarce and therefore the present result adds an important new datapoint to our understanding of black hole demography.

We analyze HST/STIS spectra (see Paper I) of the central region of the dense globular cluster M15. We infer the velocities of 64 individual stars, two-thirds of which have their velocity measured for the first time. This triples the number of stars with measured velocities in the central 1 arcsec of M15 and doubles the number in the central 2 arcsec. Combined with existing ground-based data we obtain the radial profiles of the projected kinematical quantities. The RMS velocity sigma_RMS rises to 14 km/s in the central few arcsec, somewhat higher than the values of 10-12 km/s inferred previously from ground-based data. To interpret the results we construct dynamical models based on the Jeans equation, which imply that M15 must have a central concentration of non-luminous material. If this is due to a single black hole, then its mass is M_BH = (3.9 +/- 2.2) x 10^3 solar masses. This is consistent with the relation between M_BH and sigma_RMS that has been established for galaxies. Also, the existence of intermediate-mass black holes in globular clusters is consistent with several scenarios for globular cluster evolution proposed in the literature. Therefore, these results may have important implications for our understanding of the evolution of globular clusters, the growth of black holes, the connection between globular cluster and galaxy formation, and the nature of the recently discovered `ultra-luminous' X-ray sources in nearby galaxies. Instead of a single black hole, M15 could have a central concentration of dark remnants (e.g., neutron stars) due to mass segregation. However, the best-fitting Fokker-Planck models that have previously been constructed for M15 do not predict a central mass concentration that is sufficient to explain the observed kinematics.[ABRIDGED]

It has been reported that there is an error in the figure in Dull et al. (1997, D97) that shows the radial M/L profile in Fokker-Planck models of M15. We discuss how this modifies the interpretation of our kinematical data. These imply the existence of a dark and compact mass component near the center of M15, either a single black hole (BH) or a collection of dark remnants that have sunk to the cluster center due to mass segregation. We previously showed that the latter interpretation is in conflict with the D97 M/L profile, which supported the BH interpretation. We repeat our analysis here with the corrected D97 profile. Models without a BH are now found to be statistically acceptable (within 1-sigma), although inclusion of a BH still provides a marginally better fit. It does not necessarily follow that dark remnants are now the preferred interpretation of the data. The D97 models, as well as N-body models by Baumgardt et al. (2002), assume that all neutron stars are retained during cluster evolution. This conflicts with predictions of the neutron-star retention rate (typically below 10 per cent) based on pulsar kick velocities. The presence of a single BH therefore continues to be a viable interpretation of the data. The best fit BH mass with the corrected D97 M/L profile is M_BH = 1.7^{+2.7}_{-1.7} x 10^3 solar masses, and with a constant M/L it is M_BH = 3.2^{+2.2}_{-2.2} x 10^3 solar masses. A model that includes both neutron star escape and mass segregation would probably yield a value between these numbers. This agrees with the correlation between velocity dispersion and BH mass inferred for galaxies. However, with the presently available models and data it is neither uniquely implied nor ruled out that M15 has an intermediate-mass BH.

In this series of two papers we describe a project with STIS on HST to measure the line-of-sight velocities of stars in the central few arcsec of the dense globular cluster M15. This first paper focuses on the observations and reduction of the data. We `scanned' the central region of M15 spectroscopically by consecutively placing an 0.1 arcsec slit at 18 adjacent positions. Careful and complicated data reduction and analysis were required. The spectral pixel size exceeds the velocity dispersion of M15, which puts the project at the limit of what is feasible with STIS. We applied corrections for the following effects: drifts in the STIS wavelength scale during an HST orbit; the orbital velocity component of HST along the line-of-sight to the cluster, and its variations during the HST orbit; and the apparent wavelength shift that is perceived for a star that is not centered in the slit. To make these corrections we present a detailed study of the HST/STIS point-spread function and line-spread function. We also present a catalog with positions and U, B and V magnitudes of 31,983 stars in M15, created from existing HST/WFPC2 data. We extract 19,200 one-dimensional spectra from the STIS data, each for a different position in M15, and each with a velocity scale accurate to better than 2.5 km/s. We develop an algorithm that co-adds the spectra for individual apertures and use it to extract spectra of individual stars with minimum blending and maximum S/N. In Paper II we use these spectra to extract the line-of-sight velocities of individual stars. The velocities constrain the central structure, dynamics and mass distribution of the cluster, and provide evidence for the existence of an intermediate-mass black hole in the center of the cluster.[ABRIDGED]

We have observed planetary nebulae (PNe) in several early-type galaxies using new techniques on 4- to 8-meter-class telescopes. We obtain the first large data sets (~ 100 velocities each) of PN kinematics in galaxies at > 15 Mpc, and present some preliminary dynamical results.

WFPC2 images and STIS spectroscopic observations are presented of the double nucleus in the merger system NGC 6240. We find that: (a) the kinematics of the ionized gas is similar to that of the molecular gas, despite a different morphology; (b) the gaseous and stellar kinematics are quite different, suggesting an early merger stage; (c) neither the gaseous nor the stellar kinematics show an obvious sign of the supermassive black hole believed to be responsible for the X-ray emission of NGC 6240; and (d) the steep off-nuclear velocity gradient is not due to a 10E11 solar mass black hole, in contrast to earlier suggestions.

We describe an on-going HST program aimed at determining the relationship between the nuclear black hole mass and bulge mass in spiral galaxies. We have selected a volume limited sample of 54 nearby spiral galaxies for which we already have ground based emission line rotation curves, CCD surface photometry and radio maps. We are now obtaining HST/STIS longslit observations of each of the galaxies in the sample in order to determine the nuclear Halpha rotation curve at high (~0.1") spatial resolution. We will use these data to measure the unresolved dark mass concentration at the nucleus of each object. Here we show the first results from observations of objects in the sample.

By making spectral absorption-line observations of disk galaxies at intermediate inclinations, we have been able to determine the amplitude of their constituent stars' random motions in three dimensions. This full measure of the shape of the velocity ellipsoid is a useful diagnostic for determining the "heating" mechanism responsible for creating the random motions. The analysis implies that the main heating process varies with galaxy type, with minor mergers dominating the heating in early-type disk galaxies, and spiral density waves the most important mechanism in late-type systems.

Various processes have been proposed to explain how galaxy discs acquire their thickness. A simple diagnostic for ascertaining this ``heating'' mechanism is provided by the ratio of the vertical to radial velocity dispersion components. In a previous paper we have developed a technique for measuring this ratio, and demonstrated its viability on the Sb system NGC 488. Here we present follow-up observations of the morphologically similar Sab galaxy NGC 2985, still only the second galaxy for which this ratio has been determined outside of the solar neighbourhood. The result is consistent with simple disc heating models which predict ratios of $\sigma_z / \sigma_R$ less than one

The planetary nebula populations of relatively nearby galaxies can be easily observed and provide both a distance estimate and a tool with which dynamical information can be obtained. Usually the requisite radial velocities are obtained by multi-object spectroscopy once the planetary nebulae have been located by direct imaging. Here we report on a technique for measuring planetary nebula kinematics using the double-beam ISIS spectrograph at the William Herschel Telescope in a novel slitless mode, which enables the detection and radial velocity measurements to be combined into a single step. The results on our first target, the Sab galaxy NGC 4736, allow the velocity dispersion of the stellar population in a disk galaxy to be traced out to four scale lengths for the first time and are consistent with a simple isothermal sheet model.

Theories of stellar orbit diffusion in disk galaxies predict different rates of increase of the velocity dispersions parallel and perpendicular to the disk plane, and it is therefore of interest to measure the different velocity dispersion components in galactic disks of different types. We show that it is possible to extract the three components of the velocity ellipsoid in an intermediate-inclination disk galaxy from measured line-of-sight velocity dispersions on the major and minor axes. On applying the method to observations of the Sb galaxy NGC 488, we find evidence for a higher ratio of vertical to radial dispersion in NGC 488 than in the solar neighbourhood of the Milky Way (the only other place where this quantity has ever been measured). The difference is qualitatively consistent with the notion that spiral structure has been relatively less important in the dynamical evolution of the disk of NGC 488 than molecular clouds.

Using a spectrum obtained under moderate (of order 1 arcsecond) seeing, we show that the double nucleus in M31 produces a strong kinematic signature even though the individual components are not spatially resolved. The signature consists of a significant asymmetric wing in the stellar velocity distribution close to the center of the system. The properties of the second nucleus derived from this analysis agree closely with those measured from high-spatial resolution Hubble Space Telescope images. Even Space Telescope only has sufficient resolution to study the structure of very nearby galactic nuclei photometrically; this spectroscopic approach offers a tool for detecting structure such as multiple nuclei in a wider sample of galaxy cores.