Bars inhabit the majority of local-Universe disk galaxies and may be important drivers of galaxy evolution through the redistribution of gas and angular momentum within disks. We investigate the star formation and gas properties of bars in galaxies spanning a wide range of masses, environments, and star formation rates using the MaNGA galaxy survey. Using a robustly-defined sample of 684 barred galaxies, we find that fractional (or scaled) bar length correlates with the host's offset from the star-formation main sequence. Considering the morphology of the H$\alpha$ emission we separate barred galaxies into different categories, including barred, ringed, and central configurations, together with H$\alpha$ detected at the ends of a bar. We find that only low-mass galaxies host star formation along their bars, and that this is located predominantly at the leading edge of the bar itself. Our results are supported by recent simulations of massive galaxies, which show that the position of star formation within a bar is regulated by a combination of shear forces, turbulence and gas flows. We conclude that the physical properties of a bar are mostly governed by the existing stellar mass of the host galaxy, but that they also play an important role in the galaxy's ongoing star formation.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

We perform a "fossil record" analysis for ~800 low-redshift spiral galaxies, using STARLIGHT applied to integral field spectroscopic observations from the SDSS-IV MaNGA survey to obtain fully spatially-resolved high-resolution star formation histories (SFHs). From the SFHs, we are able to build maps indicating the present-day distribution of stellar populations of different ages in each galaxy. We find small negative mean age gradients in most spiral galaxies, especially at high stellar mass, which reflects the formation times of stellar populations at different galactocentric radii. We show that the youngest (<10^{8.5} years) populations exhibit significantly more extended distributions than the oldest (>10^{9.5} years), again with a strong dependence on stellar mass. By interpreting the radial profiles of "time slices" as indicative of the size of the galaxy at the time those populations had formed, we are able to trace the simultaneous growth in mass and size of the spiral galaxies over the last 10 Gyr. Despite finding that the evolution of the measured light-weighted radius is consistent with inside-out growth in the majority of spiral galaxies, the evolution of an equivalent mass-weighted radius has changed little over the same time period. Since radial migration effects are likely to be small, we conclude that the growth of disks in spiral galaxies has occurred predominantly through an inside-out mode (with the effect greatest in high-mass galaxies), but this has not had anywhere near as much impact on the distribution of mass within spiral galaxies.

Dust attenuation in star-forming spiral galaxies affects stars and gas in different ways due to local variations in dust geometry. We present spatially resolved measurements of dust attenuation for a sample of 232 such star-forming spiral galaxies, derived from spectra acquired by the SDSS-IV MaNGA survey. The dust attenuation affecting the stellar populations of these galaxies (obtained using full spectrum stellar population fitting methods) is compared with the dust attenuation in the gas (derived from the Balmer decrement). Both of these attenuation measures increase for local regions of galaxies with higher star formation rates; the dust attenuation affecting the stellar populations increases more so than the dust attenuation in the gas, causing the ratio of the dust attenuation affecting the stellar populations to the dust attenuation in the gas to decrease for local regions of galaxies with higher star formation rate densities. No systematic difference is discernible in any of these dust attenuation quantities between the spiral arm and inter-arm regions of the galaxies. While both the dust attenuation in the gas and the dust attenuation affecting the stellar populations decrease with galactocentric radius, the ratio of the two quantities does not vary with radius. This ratio does, however, decrease systematically as the stellar mass of the galaxy increases. Analysis of the radial profiles of the two dust attenuation measures suggests that there is a disproportionately high concentration of birth clouds (incorporating gas, young stars and clumpy dust) nearer to the centres of star-forming spiral galaxies.

We present the first attempt to select AGN using long-term NIR variability. By analysing the K-band light curves of all the galaxies in the UKIDSS Ultra Deep Survey, the deepest NIR survey over ~1 sq degree, we have isolated 393 variable AGN candidates. A comparison to other selection techniques shows that only half of the variable sources are also selected using either deep Chandra X-ray imaging or IRAC colour selection, suggesting that using NIR variability can locate AGN that are missed by more standard selection techniques. In particular, we find that long-term NIR variability identifies AGN at low luminosities and in host galaxies with low stellar masses, many of which appear relatively X-ray quiet.

Despite numerous efforts, it is still unclear whether lenticular galaxies (S0s) evolve from spirals whose star formation was suppressed, or formed trough mergers or disk instabilities. In this paper we present a pilot study of 21 S0 galaxies in extreme environments (field and cluster), and compare their spatially-resolved kinematics and global stellar populations. Our aim is to identify whether there are different mechanisms that form S0s in different environments. Our results show that the kinematics of S0 galaxies in field and cluster are, indeed, different. Lenticulars in the cluster are more rotationally supported, suggesting that they are formed through processes that involve the rapid consumption or removal of gas (e.g. starvation, ram pressure stripping). In contrast, S0s in the field are more pressure supported, suggesting that minor mergers served mostly to shape their kinematic properties. These results are independent of total mass, luminosity, or disk-to-bulge ratio. On the other hand, the mass-weighted age, metallicity, and star formation time-scale of the galaxies correlate more with mass than with environment, in agreement with known relations from previous work such as the one between mass and metallicity. Overall, our results re-enforce the idea that there are multiple mechanisms that produce S0s, and that both mass $and$ environment play key roles. A larger sample is highly desirable to confirm or refute the results and the interpretation of this pilot study.

Estimating the bar pattern speed (\Om{}) is one of the main challenges faced in understanding the role of stellar bars in galaxy dynamical evolution. This work aims to characterise different uncertainty sources affecting the Tremaine Weinberg (TW)-method to study the correlation between bar and galaxies physical parameters. We use a sample of 15 MaNGA SDSS-IV galaxies and 3 CALIFA galaxies from \cite{Aguerri2015}. We studied the errors related with (i) galaxy centre determination, (ii) disc position angle (PA) emphasising the difficulties triggered by outer non-axisymmetric structures besides the bar, (iii) the slits length and (iv) the spatial resolution. In average, the PA uncertainties range $\sim 15 \%$, the slit length $\sim 9 \%$ and the centring error $\sim 5 \%$. Reducing the spatial resolution increases the sensitivity to the PA error. Through Monte Carlo simulations, we estimate the probability distribution of the \R{} bar speed parameter. The present sample is composed of 7 slow, 4 fast and 7 ultrafast bars, with no trend with morphological types. Although uncertainties and low sample numbers may mask potential correlations between physical properties, we present a discussion of them: We observe an anti-correlation of \Om{} with the bar length and the stellar mass, suggesting that massive galaxies tend to host longer and slower bars. We also observe a correlation of the molecular gas fraction with \R{}, and a weak anti-correlation with \Om{}, suggesting that bars rotate slower in gaseous discs. Confirmation of such trends awaits future studies.

The Milky Way has been described as an anaemic spiral, but is its star formation rate (SFR) unusually low when compared to its peers? To answer this question, we define a sample of Milky Way Analogues (MWAs) based on stringent cuts on the best literature estimates of non-transient structural features for the Milky Way. This selection yields only 176 galaxies from the whole of the SDSS DR7 spectroscopic sample which have morphological classifications in GZ2, from which we infer SFRs from two separate indicators. The mean SFRs found are $\log(\rm{SFR}_{SED}/\rm{M}_{\odot}~\rm{yr}^{-1})=0.53$ with a standard deviation of 0.23 dex from SED fits, and $\log(\rm{SFR}_{W4}/\rm{M}_{\odot}~\rm{yr}^{-1})=0.68$ with a standard deviation of 0.41 dex from a mid-infrared calibration. The most recent estimate for the Milky Way's star formation rate of $\log(\rm{SFR}_{MW}/\rm{M}_{\odot}~\rm{yr}^{-1})=0.22$ fits well within 2$\sigma$ of these values, where $\sigma$ is the standard deviation of each of the SFR indicator distributions. We infer that the Milky Way, while being a galaxy with a somewhat low SFR, is not unusual when compared to similar galaxies.

We present the initial results of a census of 684 barred galaxies in the MaNGA galaxy survey. This large sample contains galaxies with a wide range of physical properties, and we attempt to link bar properties to key observables for the whole galaxy. We find the length of the bar, when normalised for galaxy size, is correlated with the distance of the galaxy from the star formation main sequence, with more passive galaxies hosting larger-scale bars. Ionised gas is observed along the bars of low-mass galaxies only, and these galaxies are generally star-forming and host short bars. Higher-mass galaxies do not contain H{\alpha} emission along their bars, however, but are more likely to host rings or H{\alpha} at the centre and ends of the bar. Our results suggest that different physical processes are at play in the formation and evolution of bars in low- and high-mass galaxies.

Spectra of galaxies contain a wealth of information about the stellar populations from which they are made. With integral-field unit (IFU) surveys, such data can be used to map out stellar population properties across the face of a galaxy, allowing one to go beyond simple radial profiles and study details of non-axisymmetric structure. To-date, however, such studies have been limited by the quality of available data and the power of spectral analysis tools. We now take the next step and study the barred spiral galaxy MCG+07-28-064 from observations obtained as part of the SDSS-IV MaNGA project. We find that we can decompose this galaxy into "time slices," which reveal the varying contributions that stars of differing ages make to its bar and spiral structure, offering new insight into the evolution of these features. We find evidence for the ongoing growth of the bar, including the most recent star formation on its leading edge, and for the underlying density wave responsible for spiral structure. This pilot study indicates that there is a wealth of untapped information on the spatial distribution of SFHs available in the current generation of IFU galaxy surveys.

Bars in galaxies are thought to stimulate both inflow of material and radial mixing along them. Observational evidence for this mixing has been inconclusive so far however, limiting the evaluation of the impact of bars on galaxy evolution. We now use results from the MaNGA integral field spectroscopic survey to characterise radial stellar age and metallicity gradients along the bar and outside the bar in 128 strongly barred galaxies. We find that age and metallicity gradients are flatter in the barred regions of almost all barred galaxies when compared to corresponding disk regions at the same radii. Our results re-emphasize the key fact that by azimuthally averaging integral field spectroscopic data one loses important information from non-axisymmetric galaxy components such as bars and spiral arms. We interpret our results as observational evidence that bars are radially mixing material in galaxies of all stellar masses, and for all bar morphologies and evolutionary stages.

The exact nature of the arms of spiral galaxies is still an open question. It has been widely assumed that spiral arms in galaxies with two distinct symmetrical arms are the products of density waves that propagate around the disk, with the spiral arms being visibly enhanced by the star formation that is triggered as the passing wave compresses gas in the galaxy disk. Such a persistent wave would propagate with an approximately constant angular speed, its pattern speed Op. The quasi-stationary density wave theory can be tested by measuring this quantity and showing that it does not vary with radius in the galaxy. Unfortunately, this measurement is difficult because Op is only indirectly connected to observables such as the stellar rotation speed. Here, we use the detailed information on stellar populations of the grand-design spiral galaxy UGC 3825, extracted from spectral mapping, to measure the offset between young stars of a known age and the spiral arm in which they formed, allowing the first direct measure of Op at a range of radii. The offset in this galaxy is found to be as expected for a pattern speed that varies little with radius, indicating consistency with a quasi-stationary density wave, and lending credence to this new method.

By applying spectroscopic decomposition methods to a sample of MaNGA early-type galaxies, we separate out spatially and kinematically distinct stellar populations, allowing us to explore the similarities and differences between galaxy bulges and discs, and how they affect the global properties of the galaxy. We find that the components have interesting variations in their stellar populations, and display different kinematics. Bulges tend to be consistently more metal rich than their disc counterparts, and while the ages of both components are comparable, there is an interesting tail of younger, more metal poor discs. Bulges and discs follow their own distinct kinematic relationships, both on the plane of the stellar spin parameter, lambda_R, and ellipticity, and in the relation between stellar mass and specific angular momentum, j, with the location of the galaxy as a whole on these planes being determined by how much bulge and disc it contains. As a check of the physical significance of the kinematic decompositions, we also dynamically model the individual galaxy components within the global potential of the galaxy. The resulting components exhibit kinematic parameters consistent with those from the spectroscopic decomposition, and though the dynamical modelling suffers from some degeneracies, the bulges and discs display systematically different intrinsic dynamical properties. This work demonstrates the value in considering the individual components of galaxies rather than treating them as a single entity, which neglects information that may be crucial in understanding where, when and how galaxies evolve into the systems we see today.

Deriving a simple, analytic galaxy star formation history (SFH) using observational data is a complex task without the proper tool to hand. We therefore present SNITCH, an open source code written in Python, developed to quickly (~2 minutes) infer the parameters describing an analytic SFH model from the emission and absorption features of a galaxy spectrum dominated by star formation gas ionisation. SNITCH uses the Flexible Stellar Population Synthesis models of Conroy et al. (2009), the MaNGA Data Analysis Pipeline and a Markov Chain Monte Carlo method in order to infer three parameters (time of quenching, rate of quenching and model metallicity) which best describe an exponentially declining quenching history. This code was written for use on the MaNGA spectral data cubes but is customisable by a user so that it can be used for any scenario where a galaxy spectrum has been obtained, and adapted to infer a user-defined analytic SFH model for specific science cases. Herein we outline the rigorous testing applied to SNITCH and show that it is both accurate and precise at deriving the SFH of a galaxy spectra. The tests suggest that SNITCH is sensitive to the most recent epoch of star formation but can also trace the quenching of star formation even if the true decline does not occur at an exponential rate. With the use of both an analytical SFH and only five spectral features, we advocate that this code be used as a comparative tool across a large population of spectra, either for integral field unit data cubes or across a population of galaxy spectra.

The gravitational potential well of an $M>10^{13}$ $M_\odot$ galaxy will lead to a gravitational redshift differential of order 1 km/s between the galaxy core and its outskirts. Current surveys of massive galaxies with spatially resolved spectroscopy have reached a size which makes feasible attempts to detect gravitational redshifts within galaxies. We use spectra from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) experiment to attempt a measurement of the averaged stellar redshift profile of large elliptical galaxies. We find that systematic effects (possibly related to charge transfer or wavelength calibration accuracy) make the standard MaNGA data pipeline unsuitable for measuring the relevant sub km/s wavelength shifts. We therefore develop a cross-correlation technique to mitigate these effects, but caution that we are working beyond the design accuracy of the MaNGA experiment. With a sample of $272$ galaxies in halos with $\log (M/M_\odot)>13$, we attempt a measurement of the gravitational redshift profile, achieving $1 \sigma$ errors of size $\sim 0.5$ km/s, but are unable to make a significant detection of the signal. Even without a detection, our measurement can be used to limit the dark matter mass within the half light radius of elliptical galaxies to $1.2$ times the stellar mass, at the 68% confidence level. We also perform weighting scheme tests and split sample tests, and address target selection issues and other relativistic effects, including the transverse Doppler effect and relativistic beaming of stars. Future detections could lead to new constraints on the galaxy mass distribution that are different from kinematic and lensing determinations and open a window on galaxy properties and tests of gravity.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

The MaNGA project has obtained IFU data for several thousand nearby galaxies, including barred galaxies. With the two dimensional spectral and kinematic information provided by IFUs, we can measure the pattern speed of a barred galaxy, which determines the bar dynamics. We apply the non-parametric method proposed by Tremaine \& Weinberg to estimate the bar pattern speed for 53 barred galaxies, making this the largest sample studied so far in this way. Our sample is selected from the MaNGA first public data release as part of SDSS Data Release 13 according mainly to the axis ratio and position angle difference between the bar and disc, while kinematic data is from the later SDSS Data Release 14. We have used both the photometric position angle from the photometric image and the kinematic position angle from the stellar velocity map to derive the pattern speed. Combining three independent bar length measurements and the circular velocity from Jeans Anisotropic modelling (JAM), we also determine the dimensionless ratio $\cal{R}$ of the corotation radius to the bar length. We find that the galaxy's position angle is the main uncertainty in determining the bar pattern speed. The kinematic position angle leads to fewer ultrafast bars than the photometric position angle, and this could be due to the method of measuring the kinematic position angle. We study the dependence of $\cal{R}$ values on galaxy properties such as the dark matter fraction from JAM modelling and the stellar age and metallicity from stellar population synthesis (SPS). A positive correlation between the bar length and bar strength is found: the longer the bar, the stronger the bar. However, no other significant correlations are found. This may result from errors in deriving the $\cal{R}$ values or from the complex formation and slowdown processes of galactic bars.

Gas stripping of spiral galaxies or mergers are thought to be the formation mechanisms of lenticular galaxies. In order to determine the conditions in which each scenario dominates, we derive stellar populations of both the bulge and disk regions of 279 lenticular galaxies in the MaNGA survey. We find a clear bimodality in stellar age and metallicity within the population of S0s and this is strongly correlated with stellar mass. Old and metal-rich bulges and disks belong to massive galaxies, and young and metal-poor bulges and disks are hosted by low-mass galaxies. From this we conclude that the bulges and disks are co-evolving. When the bulge and disk stellar ages are compared, we find that the bulge is almost always older than the disk for massive galaxies ($\textrm{M}_{\star} > 10^{10}~\textrm{M}_{\odot}$). The opposite is true for lower mass galaxies. We conclude that we see two separate populations of lenticular galaxies. The old, massive, and metal-rich population possess bulges that are predominantly older than their disks, which we speculate may have been caused by morphological or inside-out quenching. In contrast, the less massive and more metal-poor population have bulges with more recent star formation than their disks. We postulate they may be undergoing bulge rejuvenation (or disk fading), or compaction. Environment doesn't play a distinct role in the properties of either population. Our findings give weight to the notion that while the faded spiral scenario likely formed low-mass S0s, other processes, such as mergers, may be responsible for high-mass S0s.

In the hierarchical two-phase formation scenario, the extended halos of early type galaxies (ETGs) are expected to have different physical properties from those of the galaxies' central regions. This work aims at characterizing the kinematic properties of ETG halos using planetary nebulae (PNe) as tracers, which allow us to overcome the limitations of absorption line spectroscopy of continuum at low surface brightness. We present two-dimensional velocity and velocity dispersion fields for 33 ETGs, including both fast (FRs) and slow rotators (SRs), making this the largest kinematic survey to-date of extragalactic PNe. The velocity fields extend out to a median 5.6 effective radii (Re), with a range [3Re-13Re] for the PN.S ETGs. We complemented the PN kinematics with absorption line data from the literature. We find that ETGs typically show a kinematic transition between inner regions and halo. Estimated transition radii in units of Re anticorrelate with stellar mass. SRs have increased but still modest rotational support at large radii, while most of the FRs show a decrease in rotation, due to the fading of the stellar disk in the outer, more slowly rotating spheroid. 30% of the FRs are dominated by rotation also at large radii. Most ETGs have flat or slightly falling halo velocity dispersion profiles, but 15% of the sample have steeply falling profiles. 40% of the FRs show kinematic twists, misalignments, or rotation along two axes, indicating that they turn from oblate at the center to triaxial in the halo, consistently with the distribution of isophote twists from extended photometry. ETGs have more diverse kinematic properties in their halos than in the central regions. The observed kinematic transition to the halo and its dependence on stellar mass is consistent with LambdaCDM simulations and supports a two-phase formation scenario. [abridged]

Recovering the origins of lenticular galaxies can shed light on the understanding of galaxy formation and evolution, since they present properties that can be found in both elliptical and spiral galaxies. In this work we study the kinematics of the globular cluster (GC) systems of three lenticular galaxies located in low density environments (NGC 2768, NGC 3115 and NGC 7457), and compare them with the kinematics of planetary nebulae (PNe). The PNe and GC data come from the Planetary Nebulae Spectrograph and the SLUGGS Surveys. Through photometric spheroid-disc decomposition and PNe kinematics we find the probability for a given GC to belong to either the spheroid or the disc of its host galaxy or be rejected from the model. We find that there is no correlation between the components that the GCs are likely to belong to and their colours. Particularly, for NGC 2768 we find that its red GCs display rotation preferentially at inner radii (Re<1). In the case of the GC system of NGC 3115 we find a group of GCs with similar kinematics that are not likely to belong to either its spheroid nor disc. For NGC 7457 we find that 70% of its GCs are likely to belong to the disc. Overall, our results suggest that these galaxies assembled into S0s through different evolutionary paths. Mergers seem to have been very important for NGC 2768 and NGC 3115 while NGC 7457 is more likely to have experienced secular evolution.

Information on the star-formation histories of cD galaxies and their extended stellar haloes lie in their spectra. Therefore, to determine whether these structures evolved together or through a two-phase formation, we need to spectroscopically separate the light from each component. We present a pilot study to use BUDDI to fit and extract the spectra of the cD galaxy NGC 3311 and its halo in an Integral Field Spectroscopy datacube, and carry out a simple stellar populations analysis to study their star-formation histories. Using MUSE data, we were able to isolate the light of the galaxy and its halo throughout the datacube, giving spectra representing purely the light from each of these structures. The stellar populations analysis of the two components indicates that, in this case, the bulk of the stars in both the halo and the central galaxy are very old, but the halo is more metal poor and less $\alpha$-enriched than the galaxy. This result is consistent with the halo forming through the accretion of much smaller satellite galaxies with more extended star formation. It is noteworthy that the apparent gradients in age and metallicity indicators across the galaxy are entirely consistent with the radially-varying contributions of galaxy and halo components, which individually display no gradients. The success of this study is promising for its application to a larger sample of cD galaxies that are currently being observed by IFU surveys.

M49 (NGC 4472) is the dominant galaxy in subcluster B of the Virgo Cluster, and a benchmark for studying the build-up of the extended halos of brightest group galaxies in the outskirts of galaxy clusters. We investigate the kinematics in the outer halo of M49, look for substructures, and describe the transition to the surrounding intra-group light. As kinematic tracers we use planetary nebulae (PNe), combining kinematics from the extended Planetary Nebula Spectrograph (PN.S) early-type galaxy survey with our recent deep photometric sample. We study the position-velocity-plane for bright and faint PN populations out to 95 kpc radius, and employ a multi-Gaussian model for the velocity distribution to identify stellar populations with distinct kinematics and histories. We report the detection of stellar-kinematic substructure associated with the interaction of M49 with the dwarf irregular galaxy VCC 1249. We find two kinematically distinct PN populations associated with the main M49 halo and the extended intra-group light (IGL). These have velocity dispersions $\sigma_{halo}\simeq 170$ km/s and $\sigma_{IGL} \simeq400$ km/s at 10-80 kpc radii. The overall luminosity profile and velocity dispersion at $\sim80$ kpc are consistent with a flat circular velocity curve extrapolated from X-ray observations. The dispersion of the PNe associated with the IGL joins onto that of the satellite galaxies in subcluster B at $\sim100$ kpc radius. This is the first time that the transition from halo to IGL is observed based on the velocities of individual stars. Therefore the halo of M49, consisting of at least three distinct components, has undergone an extended accretion history within its parent group potential. The blue colours of the IGL component are consistent with a population of stars formed in low-mass galaxies at redshift $\sim0.5$ that has since evolved passively, as suggested by other data.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

We investigate radial gradients in the recent star formation history (SFH) using 1917 galaxies with $0.01<z<0.14$ and integral-field spectroscopy from the ongoing MaNGA survey. For each galaxy, we obtain two-dimensional maps and radial profiles for three spectroscopically-measured parameters that are sensitive to the recent SFH: D$_n$(4000) (the 4000\AA\ break), EW(H$\delta_A$) (equivalent width of the H$\delta$ absorption line), and EW(H$\alpha$) (equivalent width of the H$\alpha$ emission line). We find the majority of the spaxels in these galaxies are consistent with models with continuously declining star formation rate, indicating that starbursts occur rarely in local galaxies. We classify the galaxies into three classes: fully star-forming (SF), partly quenched (PQ) and totally quenched (TQ), according to the fraction of quenched area within 1.5 times the effective radius, $f_Q(1.5R_e)$. We find that galaxies less massive than $10^{10}$M$_{\odot}$ present at most weak radial gradients in all the diagnostic parameters. In contrast, galaxies with stellar mass above $10^{10}$M$_{\odot}$ present significant gradients in all the three diagnostic parameters if classified as SF or PQ, but show weak gradients in D$_n$(4000) and EW(H$\delta_A$) and no gradients in EW(H$\alpha$) if in the TQ class. This implies the existence of a critical stellar mass, at $\sim10^{10}$M$_{\odot}$, above which the star formation in a galaxy gets shutdown from the inside out. Galaxies tend to evolve synchronously from inner to outer regions before their mass reaches the critical value. We have further divided the sample at fixed mass by both bulge-to-total luminosity ratio and morphological type, finding that our conclusions hold regardless of these factors: it appears that the presence of a central dense object is not a driving parameter, but rather a byproduct of the star formation cessation process.

The decomposition of the rotation curve of galaxies into contribution from the disc and dark halo remains uncertain and depends on the adopted mass to light ratio (M/L) of the disc. Given the vertical velocity dispersion of stars and disc scale height, the disc surface mass density and hence the M/L can be estimated. We address a conceptual problem with previous measurements of the scale height and dispersion. When using this method, the dispersion and scale height must refer to the same population of stars. The scale height is obtained from near-IR studies of edge-on galaxies and is weighted towards older kinematically hotter stars, whereas the dispersion obtained from integrated light in the optical bands includes stars of all ages. We aim to extract the dispersion for the hotter stars, so that it can then be used with the correct scale height to obtain the disc surface mass density. We use a sample of planetary nebulae (PNe) as dynamical tracers in the face-on galaxy NGC 628. We extract two different dispersions from its velocity histogram -- representing the older and younger PNe. We also present complementary stellar absorption spectra in the inner regions of this galaxy and use a direct pixel fitting technique to extract the two components. Our analysis concludes that previous studies, which do not take account of the young disc, underestimate the disc surface mass density by a factor of ~ 2. This is sufficient to make a maximal disc for NGC 628 appear like a submaximal disc.

We study 379 central and 159 satellite early-type galaxies with two-dimensional kinematics from the integral-field survey Mapping Nearby Galaxies at APO (MaNGA) to determine how their angular momentum content depends on stellar and halo mass. Using the Yang et. al. (2007) group catalog, we identify central and satellite galaxies in groups with halo masses in the range 10^12.5 h^-1 M_sun < M_200b < 10^15 h^-1 M_sun. As in previous work, we see a sharp dependence on stellar mass, in the sense that ~ 70% of galaxies with stellar mass M_* > 10^11 h^-2 M_sun tend to have very little rotation, while nearly all galaxies at lower mass show some net rotation. The ~ 30% of high-mass galaxies that have significant rotation do not stand out in other galaxy properties except for a higher incidence of ionized gas emission. Our data are consistent with recent simulation results suggesting that major merging and gas accretion have more impact on the rotational support of lower-mass galaxies. When carefully matching the stellar mass distributions, we find no residual differences in angular momentum content between satellite and central galaxies at the 20\% level. Similarly, at fixed mass, galaxies have consistent rotation properties across a wide range of halo mass. However, we find that errors in classification of centrals and satellites with group finders systematically lowers differences between satellite and central galaxies at a level that is comparable to current measurement uncertainties. To improve constraints, the impact of group finding methods will have to be forward modeled via mock catalogs.

Do the theorised different formation mechanisms of fast and slow rotators produce an observable difference in their star formation histories? To study this we identify quenching slow rotators in the MaNGA sample by selecting those which lie below the star forming sequence and identify a sample of quenching fast rotators which were matched in stellar mass. This results in a total sample of 194 kinematically classified galaxies, which is agnostic to visual morphology. We use u-r and NUV-u colours from SDSS and GALEX and an existing inference package, STARPY, to conduct a first look at the onset time and exponentially declining rate of quenching of these galaxies. An Anderson-Darling test on the distribution of the inferred quenching rates across the two kinematic populations reveals they are statistically distinguishable ($3.2\sigma$). We find that fast rotators quench at a much wider range of rates than slow rotators, consistent with a wide variety of physical processes such as secular evolution, minor mergers, gas accretion and environmentally driven mechanisms. Quenching is more likely to occur at rapid rates ($\tau \lesssim 1~\rm{Gyr}$) for slow rotators, in agreement with theories suggesting slow rotators are formed in dynamically fast processes, such as major mergers. Interestingly, we also find that a subset of the fast rotators quench at these same rapid rates as the bulk of the slow rotator sample. We therefore discuss how the total gas mass of a merger, rather than the merger mass ratio, may decide a galaxy's ultimate kinematic fate.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

The "kinematic" morphology-density relation for early-type galaxies posits that those galaxies with low angular momentum are preferentially found in the highest-density regions of the universe. We use a large sample of galaxy groups with halo masses 10^12.5 < M_halo < 10^14.5 M_sun/h observed with the Mapping Nearby Galaxies at APO (MaNGA) survey to examine whether there is a correlation between local environment and rotational support that is independent of stellar mass. We find no compelling evidence for a relationship between the angular momentum content of early-type galaxies and either local overdensity or radial position within the group at fixed stellar mass.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

To ascertain whether photometric decompositions of galaxies into bulges and disks are astrophysically meaningful, we have developed a new technique to decompose spectral data cubes into separate bulge and disk components, subject only to the constraint that they reproduce the conventional photometric decomposition. These decompositions allow us to study the kinematic and stellar population properties of the individual components and how they vary with position, in order to assess their plausibility as discrete elements, and to start to reconstruct their distinct formation histories. An initial application of this method to CALIFA integral field unit observations of three isolated S0 galaxies confirms that in regions where both bulge and disc contribute significantly to the flux they can be physically and robustly decomposed into a rotating dispersion-dominated bulge component, and a rotating low-dispersion disc component. Analysis of the resulting stellar populations shows that the bulges of these galaxies have a range of ages relative to their discs, indicating that a variety of processes are necessary to describe their evolution. This simple test case indicates the broad potential for extracting from spectral data cubes the full spectral data of a wide variety of individual galaxy components, and for using such decompositions to understand the interplay between these various structures, and hence how such systems formed.

We study the internal radial gradients of stellar population properties within $1.5\;R_{\rm e}$ and analyse the impact of galaxy environment. We use a representative sample of 721 galaxies with masses ranging between $10^{9}\;M_{\odot}$ to $10^{11.5}\;M_{\odot}$ from the SDSS-IV survey MaNGA. We split this sample by morphology into early-type and late-type galaxies. Using the full spectral fitting code FIREFLY, we derive the light and mass-weighted stellar population properties age and metallicity, and calculate the gradients of these properties. We use three independent methods to quantify galaxy environment, namely the $N^{th}$ nearest neighbour, the tidal strength parameter $Q$ and distinguish between central and satellite galaxies. In our analysis, we find that early-type galaxies generally exhibit shallow light-weighted age gradients in agreement with the literature and mass-weighted median age gradients tend to be slightly positive. Late-type galaxies, instead, have negative light-weighted age gradients. We detect negative metallicity gradients in both early and late-type galaxies that correlate with galaxy mass, with the gradients being steeper and the correlation with mass being stronger in late-types. We find, however, that stellar population gradients, for both morphological classifications, have no significant correlation with galaxy environment for all three characterisations of environment. Our results suggest that galaxy mass is the main driver of stellar population gradients in both early and late-type galaxies, and any environmental dependence, if present at all, must be very subtle.

We study the internal gradients of stellar population properties within $1.5\;R_{\rm e}$ for a representative sample of 721 galaxies with stellar masses ranging between $10^{9}\;M_{\odot}$ to $10^{11.5}\;M_{\odot}$ from the SDSS-IV MaNGA IFU survey. Through the use of our full spectral fitting code FIREFLY, we derive light and mass-weighted stellar population properties and their radial gradients, as well as full star formation and metal enrichment histories. We also quanfify the impact that different stellar population models and full spectral fitting routines have on the derived stellar population properties, and the radial gradient measurements. In our analysis, we find that age gradients tend to be shallow for both early-type and late-type galaxies. {\em Mass-weighted} age gradients of early-types are positive ($\sim 0.09\; {\rm dex}/R_{\rm e}$) pointing to "outside-in" progression of star formation, while late-type galaxies have negative {\em light-weighted} age gradients ($\sim -0.11\; {\rm dex}/R_{\rm e}$), suggesting an "inside-out" formation of discs. We detect negative metallicity gradients in both early and late-type galaxies, but these are significantly steeper in late-types, suggesting that radial dependence of chemical enrichment processes and the effect of gas inflow and metal transport are far more pronounced in discs. Metallicity gradients of both morphological classes correlate with galaxy mass, with negative metallicity gradients becoming steeper with increasing galaxy mass. The correlation with mass is stronger for late-type galaxies, with a slope of $d(\nabla [Z/H])/d(\log M)\sim -0.2\pm 0.05\;$, compared to $d(\nabla [Z/H])/d(\log M)\sim -0.05\pm 0.05\;$ for early-types. This result suggests that the merger history plays a relatively small role in shaping metallicity gradients of galaxies.

M40 has always been something of an oddity in the Messier Catalogue, since it is just a pair of stars rather than an extended object. Doubts have also been expressed as to whether it is even a physical binary. Here, we note that the Hipparcos/Gaia parallaxes for these two stars remove all doubt: Messier 40 comprises a pair of entirely unrelated stars.

We study the properties of 66 galaxies with kinematically misaligned gas and stars from MaNGA survey. The fraction of kinematically misaligned galaxies varies with galaxy physical parameters, i.e. M*, SFR and sSFR. According to their sSFR, we further classify these 66 galaxies into three categories, 10 star-forming, 26 "Green Valley" and 30 quiescent ones. The properties of different types of kinematically misaligned galaxies are different in that the star-forming ones have positive gradient in D4000 and higher gas-phase metallicity, while the green valley/quiescent ones have negative D4000 gradients and lower gas-phase metallicity on average. There is evidence that all types of the kinematically misaligned galaxies tend to live in more isolated environment. Based on all these observational results, we propose a scenario for the formation of star forming galaxies with kinematically misaligned gas and stars - the progenitor accretes misaligned gas from a gas-rich dwarf or cosmic web, the cancellation of angular momentum from gas-gas collisions between the pre-existing gas and the accreted gas largely accelerates gas inflow, leading to fast centrally-concentrated star-formation. The higher metallicity is due to enrichment from this star formation. For the kinematically misaligned green valley and quiescent galaxies, they might be formed through gas-poor progenitors accreting kinematically misaligned gas from satellites which are smaller in mass.

Galaxies grow through both internal and external processes. In about 10% of nearby red galaxies with little star formation, gas and stars are counter-rotating, demonstrating the importance of external gas acquisition in these galaxies. However, systematic studies of such phenomena in blue, star-forming galaxies are rare, leaving uncertain the role of external gas acquisition in driving evolution of blue galaxies. Based on new measurements with integral field spectroscopy of a large representative galaxy sample, we find an appreciable fraction of counter-rotators among blue galaxies (9 out of 489 galaxies). The central regions of blue counter-rotators show younger stellar populations and more intense, ongoing star formation than their outer parts, indicating ongoing growth of the central regions. The result offers observational evidence that the acquisition of external gas in blue galaxies is possible; the interaction with pre-existing gas funnels the gas into nuclear regions (< 1 kpc) to form new stars.

With the availability of large integral-field unit (IFU) spectral surveys of nearby galaxies, there is now the potential to extract spectral information from across the bulges and discs of galaxies in a systematic way. This information can address questions such as how these components built up with time, how galaxies evolve and whether their evolution depends on other properties of the galaxy such as its mass or environment. We present BUDDI, a new approach to fit the two-dimensional light profiles of galaxies as a function of wavelength to extract the spectral properties of these galaxies' discs and bulges. The fitting is carried out using GalfitM, a modified form of Galfit which can fit multi-waveband images simultaneously. The benefit of this technique over traditional multi-waveband fits is that the stellar populations of each component can be constrained using knowledge over the whole image and spectrum available. The decomposition has been developed using commissioning data from the SDSS-IV Mapping Nearby Galaxies at APO (MaNGA) survey with redshifts z<0.14 and coverage of at least 1.5 effective radii for a spatial resolution of 2.5 arcsec FWHM and field of view of >22 arcsec, but can be applied to any IFU data of a nearby galaxy with similar or better spatial resolution and coverage. We present an overview of the fitting process, the results from our tests, and we finish with example stellar population analyses of early-type galaxies from the MaNGA survey to give an indication of the scientific potential of applying bulge-disc decomposition to IFU data.

The nature of warm, ionized gas outside of galaxies may illuminate several key galaxy evolutionary processes. A serendipitous observation by the MaNGA survey has revealed a large, asymmetric H$\alpha$ complex with no optical counterpart that extends $\approx8"$ ($\approx6.3$ kpc) beyond the effective radius of a dusty, starbursting galaxy. This H$\alpha$ extension is approximately three times the effective radius of the host galaxy and displays a tail-like morphology. We analyze its gas-phase metallicities, gaseous kinematics, and emission-line ratios, and discuss whether this H$\alpha$ extension could be diffuse ionized gas, a gas accretion event, or something else. We find that this warm, ionized gas structure is most consistent with gas accretion through recycled wind material, which could be an important process that regulates the low-mass end of the galaxy stellar mass function.

The MaNGA Survey (Mapping Nearby Galaxies at Apache Point Observatory) is one of three core programs in the Sloan Digital Sky Survey IV. It is obtaining integral field spectroscopy (IFS) for 10K nearby galaxies at a spectral resolution of R~2000 from 3,622-10,354A. The design of the survey is driven by a set of science requirements on the precision of estimates of the following properties: star formation rate surface density, gas metallicity, stellar population age, metallicity, and abundance ratio, and their gradients; stellar and gas kinematics; and enclosed gravitational mass as a function of radius. We describe how these science requirements set the depth of the observations and dictate sample selection. The majority of targeted galaxies are selected to ensure uniform spatial coverage in units of effective radius (Re) while maximizing spatial resolution. About 2/3 of the sample is covered out to 1.5Re (Primary sample), and 1/3 of the sample is covered to 2.5Re (Secondary sample). We describe the survey execution with details that would be useful in the design of similar future surveys. We also present statistics on the achieved data quality, specifically, the point spread function, sampling uniformity, spectral resolution, sky subtraction, and flux calibration. For our Primary sample, the median r-band signal-to-noise ratio is ~73 per 1.4A pixel for spectra stacked between 1-1.5 Re. Measurements of various galaxy properties from the first year data show that we are meeting or exceeding the defined requirements for the majority of our science goals.

Globular clusters (GCs) can be considered discrete, long-lived, dynamical tracers that retain crucial information about the assembly history of their parent galaxy. In this paper, we present a new catalogue of GC velocities and colours for the lenticular galaxy NGC 1023, we study their kinematics and spatial distribution, in comparison with the underlying stellar kinematics and surface brightness profile, and we test a new method for studying GC properties. Specifically, we decompose the galaxy light into its spheroid (assumed to represent the bulge + halo components) and disk components and use it to assign to each GC a probability of belonging to one of the two components. Then we model the galaxy kinematics, assuming a disk and spheroidal component, using planetary nebulae (PNe) and integrated stellar light. We use this kinematic model and the probability previously obtained from the photometry to recalculate for each GC its likelihood of being associated with the disk, the spheroid, or neither. We find that the reddest GCs are likely to be associated with the disk, as found for faint fuzzies in this same galaxy, suggesting that the disk of this S0 galaxy originated at z ~ 2. The majority of blue GCs are found likely to be associated with the spheroidal (hot) component. The method also allows us to identify objects that are unlikely to be in equilibrium with the system. In NGC1023 some of the rejected GCs form a substructure in phase space that is connected with NGC 1023 companion galaxy.

This article is the second in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our aim is to quantify the scatter, systematic bias and completeness of cluster masses derived from a diverse set of 25 galaxy-based methods using two contrasting mock galaxy catalogues based on a sophisticated halo occupation model and a semi-analytic model. Analysing 968 clusters, we find a wide range in the RMS errors in log M200c delivered by the different methods (0.18 to 1.08 dex, i.e., a factor of ~1.5 to 12), with abundance matching and richness methods providing the best results, irrespective of the input model assumptions. In addition, certain methods produce a significant number of catastrophic cases where the mass is under- or over-estimated by a factor greater than 10. Given the steeply falling high-mass end of the cluster mass function, we recommend that richness or abundance matching-based methods are used in conjunction with these methods as a sanity check for studies selecting high mass clusters. We see a stronger correlation of the recovered to input number of galaxies for both catalogues in comparison with the group/cluster mass, however, this does not guarantee that the correct member galaxies are being selected. We do not observe significantly higher scatter for either mock galaxy catalogues. Our results have implications for cosmological analyses that utilise the masses, richnesses, or abundances of clusters, which have different uncertainties when different methods are used.

We present an analysis of the data produced by the MaNGA prototype run (P-MaNGA), aiming to test how the radial gradients in recent star formation histories, as indicated by the 4000AA-break (D4000), Hdelta absorption (EW(Hd_A)) and Halpha emission (EW(Ha)) indices, can be useful for understanding disk growth and star formation cessation in local galaxies. We classify 12 galaxies observed on two P-MaNGA plates as either centrally quiescent (CQ) or centrally star-forming (CSF), according to whether D4000 measured in the central spaxel of each datacube exceeds 1.6. For each galaxy we generate both 2D maps and radial profiles of D4000, EW(Hd_A) and EW(Ha). We find that CSF galaxies generally show very weak or no radial variation in these diagnostics. In contrast, CQ galaxies present significant radial gradients, in the sense that D4000 decreases, while both EW(Hd_A) and EW(Ha) increase from the galactic center outward. The outer regions of the galaxies show greater scatter on diagrams relating the three parameters than their central parts. In particular, the clear separation between centrally-measured quiescent and star-forming galaxies in these diagnostic planes is largely filled in by the outer parts of galaxies whose global colors place them in the green valley, supporting the idea that the green valley represents a transition between blue-cloud and red-sequence phases, at least in our small sample. These results are consistent with a picture in which the cessation of star formation propagates from the center of a galaxy outwards as it moves to the red sequence.

We present an overview of a new integral field spectroscopic survey called MaNGA (Mapping Nearby Galaxies at Apache Point Observatory), one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) that began on 2014 July 1. MaNGA will investigate the internal kinematic structure and composition of gas and stars in an unprecedented sample of 10,000 nearby galaxies. We summarize essential characteristics of the instrument and survey design in the context of MaNGA's key science goals and present prototype observations to demonstrate MaNGA's scientific potential. MaNGA employs dithered observations with 17 fiber-bundle integral field units that vary in diameter from 12" (19 fibers) to 32" (127 fibers). Two dual-channel spectrographs provide simultaneous wavelength coverage over 3600-10300 A at R~2000. With a typical integration time of 3 hr, MaNGA reaches a target r-band signal-to-noise ratio of 4-8 (per A, per 2" fiber) at 23 AB mag per sq. arcsec, which is typical for the outskirts of MaNGA galaxies. Targets are selected with stellar mass greater than 1e9 Msun using SDSS-I redshifts and i-band luminosity to achieve uniform radial coverage in terms of the effective radius, an approximately flat distribution in stellar mass, and a sample spanning a wide range of environments. Analysis of our prototype observations demonstrates MaNGA's ability to probe gas ionization, shed light on recent star formation and quenching, enable dynamical modeling, decompose constituent components, and map the composition of stellar populations. MaNGA's spatially resolved spectra will enable an unprecedented study of the astrophysics of nearby galaxies in the coming 6 yr.

By studying the individual star-formation histories of the bulges and discs of lenticular (S0) galaxies, it is possible to build up a sequence of events that leads to the cessation of star formation and the consequent transformation from the progenitor spiral. In order to separate the bulge and disc stellar populations, we spectroscopically decomposed long-slit spectra of Virgo Cluster S0s into bulge and disc components. Analysis of the decomposed spectra shows that the most recent star formation activity in these galaxies occurred within the bulge regions, having been fuelled by residual gas from the disc. These results point towards a scenario where the star formation in the discs of spiral galaxies are quenched, followed by a final episode of star formation in the central regions from the gas that has been funnelled inwards through the disc.

We study the dynamics of faint stellar substructures around the Umbrella Galaxy, NGC 4651, which hosts a dramatic system of streams and shells formed through the tidal disruption of a nucleated dwarf elliptical galaxy. We elucidate the basic characteristics of the system (colours, luminosities, stellar masses) using multi-band Subaru/Suprime-Cam images. The implied stellar mass-ratio of the ongoing merger event is about 1:50. We identify candidate kinematic tracers (globular clusters, planetary nebulae, H ii regions), and follow up a subset with Keck/DEIMOS spectroscopy to obtain velocities. We find that 15 of the tracers are likely associated with halo substructures, including the probable stream progenitor nucleus. These objects delineate a kinematically cold feature in position-velocity phase space. We model the stream using single test-particle orbits, plus a rescaled pre-existing N-body simulation. We infer a very eccentric orbit with a period of roughly 0.35 Gyr and turning points at approximately 2-4 and 40 kpc, implying a recent passage of the satellite through the disc, which may have provoked the visible disturbances in the host galaxy. This work confirms that the kinematics of low surface brightness substructures can be recovered and modeled using discrete tracers - a breakthrough that opens up a fresh avenue for unraveling the detailed physics of minor merging.

The individual star formation histories of bulges and discs of lenticular (S0) galaxies can provide information on the processes involved in the quenching of their star formation and subsequent transformation from spirals. In order to study this transformation in dense environments, we have decomposed long-slit spectroscopic observations of a sample of 21 S0s from the Virgo Cluster to produce one-dimensional spectra representing purely the bulge and disc light for each galaxy. Analysis of the Lick indices within these spectra reveals that the bulges contain consistently younger and more metal-rich stellar populations than their surrounding discs, implying that the final episode of star formation within S0s occurs in their central regions. Analysis of the $\alpha$-element abundances in these components further presents a picture in which the final episode of star formation in the bulge is fueled using gas that has previously been chemically enriched in the disc, indicating the sequence of events in the transformation of these galaxies. Systems in which star formation in the disk was spread over a longer period contain bulges in which the final episode of star formation occurred more recently, as one might expect for an approximately coeval population in which the transformation from spiral to S0 occurred at different times. With data of this quality and the new analysis method deployed here, we can begin to describe this process in a quantitative manner for the first time.

This paper is the first in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our primary aim is to test the performance of these cluster mass estimation techniques on a diverse set of models that will increase in complexity. We begin by providing participating methods with data from a simple model that delivers idealised clusters, enabling us to quantify the underlying scatter intrinsic to these mass estimation techniques. The mock catalogue is based on a Halo Occupation Distribution (HOD) model that assumes spherical Navarro, Frenk and White (NFW) haloes truncated at R_200, with no substructure nor colour segregation, and with isotropic, isothermal Maxwellian velocities. We find that, above 10^14 M_solar, recovered cluster masses are correlated with the true underlying cluster mass with an intrinsic scatter of typically a factor of two. Below 10^14 M_solar, the scatter rises as the number of member galaxies drops and rapidly approaches an order of magnitude. We find that richness-based methods deliver the lowest scatter, but it is not clear whether such accuracy may simply be the result of using an over-simplistic model to populate the galaxies in their haloes. Even when given the true cluster membership, large scatter is observed for the majority non-richness-based approaches, suggesting that mass reconstruction with a low number of dynamical tracers is inherently problematic.

We use integral field spectroscopy of 13 disk galaxies in the cluster AC114 at z ~ 0.31 in an attempt to disentangle the physical processes responsible for the transformation of spiral galaxies in clusters. Our sample is selected to display a dominant young stellar population, as indicated by strong Hdelta absorption lines in their integrated spectra. Most of our galaxies lack the [OII] emission line, and hence ongoing star formation. They therefore possess `k+a' spectra, indicative of a recent truncation of star formation, possibly preceded by a starburst. Disky `k+a' galaxies are a promising candidate for the intermediate stage of the transformation from star-forming spiral galaxies to passive S0s. Our observations allow us to study the spatial distributions and the kinematics of the different stellar populations within the galaxies. We used three different indicators to evaluate the presence of a young population: the equivalent width of Hdelta, the luminosity-weighted fraction of A stars, and the fraction of the galaxy light attributable to simple stellar populations with ages between 0.5 and 1.5 Gyr. We find a mixture of behaviours, but are able to show that in most of galaxies the last episode of star-formation occured in an extended disk, similar to preceeding generations of stars, though somewhat more centrally concentrated. We thus exclude nuclear starbursts and violent gravitational interactions as causes of the star formation truncation. Gentler mechanisms, such as ram-pressure stripping or weak galaxy-galaxy interactions, appear to be responsible for ending star-formation in these intermediate-redshift cluster disk galaxies.

Faint fuzzies are metal-rich apparently-old star clusters with unusually large radii (7-15 pc), found mostly in S0 galaxies, whose source remain obscure. To identify their origins, we compare planetary nebulae and neutral hydrogen with faint fuzzy positions and line-of-sight velocities in NGC1023. In this way, we rule out scenarios in which these objects are associated with an on-going merger or with a spheroid population in NGC1023. Their kinematics are indistinguishable from the stellar disk population in this galaxy, and we conclude that faint fuzzies are most likely just remnant open clusters. Their observed association with S0s then simply reflects the difficulty of identifying such objects in later-type disk galaxies.

Lenticular galaxies have long been thought of as evolved spirals, but the processes involved to quench the star formation are still unclear. By studying the individual star formation histories of the bulges and disks of lenticulars, it is possible to look for clues to the processes that triggered their transformation from spirals. To accomplish this feat, we present a new method for spectroscopic bulge-disk decomposition, in which a long-slit spectrum is decomposed into two one-dimensional spectra representing purely the bulge and disk light. We present preliminary results from applying this method to lenticular galaxies in the Virgo and Fornax Clusters, in which we show that the most recent star formation activity in these galaxies occurred within the bulges. We also find that the bulges are in general more Fe-enriched than the disks of the same galaxy, and that this enrichment grows stronger as the age of the bulge becomes younger. These results point towards a scenario where the star formation in the disks of spiral galaxies are quenched, followed by a burst of star formation in the central regions from the gas that has been funnelled inwards through the disk.

The stellar kinematics of the spheroids and discs of S0 galaxies contain clues to their formation histories. Unfortunately, it is difficult to disentangle the two components and to recover their stellar kinematics in the faint outer parts of the galaxies using conventional absorption line spectroscopy. This paper therefore presents the stellar kinematics of six S0 galaxies derived from observations of planetary nebulae (PNe), obtained using the Planetary Nebula Spectrograph. To separate the kinematics of the two components, we use a maximum-likelihood method that combines the discrete kinematic data with a photometric component decomposition. The results of this analysis reveal that: the discs of S0 galaxies are rotationally supported; however, the amount of random motion in these discs is systematically higher than in comparable spiral galaxies; and the S0s lie around one magnitude below the Tully--Fisher relation for spiral galaxies, while their spheroids lie nearly one magnitude above the Faber--Jackson relation for ellipticals. All of these findings are consistent with a scenario in which spirals are converted into S0s through a process of mild harassment or "pestering," with their discs somewhat heated and their spheroid somewhat enhanced by the conversion process. In such a scenario, one might expect the properties of S0s to depend on environment. We do not see such an effect in this fairly small sample, although any differences would be diluted by the fact that the current location does not necessarily reflect the environment in which the transformation occurred. Similar observations of larger samples probing a broader range of environments, coupled with more detailed modelling of the transformation process to match the wide range of parameters that we have shown can now be measured, should take us from these first steps to the definitive answer as to how S0 galaxies form.

The origins of S0 galaxies remain obscure, with various mechanisms proposed for their formation, likely depending on environment. These mechanisms would imprint different signatures in the galaxies' stellar kinematics out to large radii, offering a method for distinguishing between them. We aim to study a sample of six S0 galaxies from a range of environments, and use planetary nebulae (PNe) as tracers of their stellar populations out to very large radii, to determine their kinematics in order to understand their origins. Using a special-purpose instrument, the Planetary Nebula Spectrograph, we observe and extract PNe catalogues for these six systems*. We show that the PNe have the same spatial distribution as the starlight, that the numbers of them are consistent with what would be expected in a comparable old stellar population in elliptical galaxies, and that their kinematics join smoothly onto those derived at smaller radii from conventional spectroscopy. The high-quality kinematic observations presented here form an excellent set for studying the detailed kinematics of S0 galaxies, in order to unravel their formation histories. We find that PNe are good tracers of stellar kinematics in these systems. We show that the recovered kinematics are largely dominated by rotational motion, although with significant random velocities in most cases.

In order to try and understand its origins, we present high-quality long-slit spectral observations of the counter-rotating stellar discs in the strange S0 galaxy NGC 4550. We kinematically decompose the spectra into two counter-rotating stellar components (plus a gaseous component), in order to study both their kinematics and their populations. The derived kinematics largely confirm what was known previously about the stellar discs, but trace them to larger radii with smaller errors; the fitted gaseous component allows us to trace the hydrogen emission lines for the first time, which are found to follow the same rather strange kinematics previously seen in the [OIII] line. Analysis of the populations of the two separate stellar components shows that the secondary disc has a significantly younger mean age than the primary disc, consistent with later star formation from the associated gaseous material. In addition, the secondary disc is somewhat brighter, also consistent with such additional star formation. However, these measurements cannot be self-consistently modelled by a scenario in which extra stars have been added to initially-identical counter-rotating stellar discs, which rules out Evans & Collett's (1994) elegant "separatrix-crossing" model for the formation of such massive counter-rotating discs from a single galaxy, leaving some form of unusual gas accretion history as the most likely formation mechanism.

There are only a few tracers available to probe the kinematics of individual early-type galaxies beyond one effective radius. Here we directly compare a sample of planetary nebulae (PNe), globular clusters (GCs) and galaxy starlight velocities out to ~4 effective radii, in the S0 galaxy NGC 2768. Using a bulge-to-disk decomposition of a K-band image we assign PNe and starlight to either the disk or the bulge. We show that the bulge PNe and bulge starlight follow the same radial density distribution as the red subpopulation of GCs, whereas the disk PNe and disk starlight are distinct components. We find good kinematic agreement between the three tracers to several effective radii (and with stellar data in the inner regions). Further support for the distinct nature of the two galaxy components come from our kinematic analysis. After separating the tracers into bulge and disk components we find the bulge to be a slowly rotating pressure-supported system, whereas the disk reveals a rapidly rising rotation curve with a declining velocity dispersion profile. The resulting V/sigma ratio for the disk resembles that of a spiral galaxy and hints at an origin for NGC 2768 as a transformed late-type galaxy. A two-component kinematic analysis for a sample of S0s will help to elucidate the nature of this class of galaxy.

Galaxy structure and morphology is nearly always studied using the light originating from stars, however ideally one is interested in measuring structure using the stellar mass distribution. Not only does stellar mass trace out the underlying distribution of matter, it also minimises the effects of star formation and dust on the appearance and structure of a galaxy. We present in this paper a study of the stellar mass distributions and structures of galaxies at z<1 as found within the GOODS fields. We use pixel by pixel K-corrections to construct stellar mass and mass-to-light ratio maps of 560 galaxies of known morphology at magnitudes z_{850}<24. We measure structural and size parameters using these stellar mass maps, as well as on ACS BViz band imaging. This includes investigating the structural CAS-Gini-M_{20} parameters and half-light radius for each galaxy. We compare structural parameters and half-light radii in the ACS z_{850}-band and stellar mass maps, finding no systematic bias introduced by measuring galaxy sizes in z_{850}. We furthermore investigate relations between structural parameters in the ACS BViz bands and stellar mass maps, and compare our result to previous morphological studies. Combinations of various parameters in stellar mass generally reveal clear separations between early and late type morphologies, but cannot easily distinguish between star formation and dynamically disturbed systems. We also show that while ellipticals and early-type spirals have fairly constant CAS values at z<1 we find a tendency for late-type spiral and peculiar morphological types to have a higher A(M_{*}) at higher redshift. We argue that this, and the large fraction of peculiars that appear spiral-like in stellar mass maps, are possible evidence for either an active bulge formation in some late-type disks at z<1 or the presence of minor merger events.

A new method for spectroscopic bulge-disc decomposition is presented, in which the spatial light profile in a two-dimensional spectrum is decomposed wavelength-by-wavelength into bulge and disc components, allowing separate one-dimensional spectra for each component to be constructed. This method has been applied to observations of a sample of nine S0s in the Fornax Cluster in order to obtain clean high-quality spectra of their individual bulge and disc components. So far this decomposition has only been fully successful when applied to galaxies with clean light profiles, consequently limiting the number of galaxies that could be separated into bulge and disc components. Lick index stellar population analysis of the component spectra reveals that in those galaxies where the bulge and disc could be distinguished, the bulges have systematically higher metallicities and younger stellar populations than the discs. This correlation is consistent with a picture in which S0 formation comprises the shutting down of star formation in the disc accompanied by a final burst of star formation in the bulge. The variation in spatial-fit parameters with wavelength also allows us to measure approximate colour gradients in the individual components. Such gradients were detected separately in both bulges and discs, in the sense that redder light is systematically more centrally concentrated in all components. However, a search for radial variations in the absorption line strengths determined for the individual components revealed that they are absent from the vast majority of S0 discs and bulges. The absence of gradients in line indices for most galaxies implies that the colour gradient cannot be attributed to age or metallicity variations, and is therefore most likely associated with varying degrees of obscuration by dust.

A candidate diffuse stellar substructure was previously reported in the halo of the nearby dwarf starburst galaxy NGC 4449 by Karachentsev et al. We map and analyze this feature using a unique combination of deep integrated-light images from the Black Bird 0.5-meter telescope, and high-resolution wide-field images from the 8-meter Subaru telescope, which resolve the nebulosity into a stream of red giant branch stars, and confirm its physical association with NGC 4449. The properties of the stream imply a massive dwarf spheroidal progenitor, which after complete disruption will deposit an amount of stellar mass that is comparable to the existing stellar halo of the main galaxy. The ratio between luminosity or stellar-mass between the two galaxies is ~1:50, while the indirectly measured dynamical mass-ratio, when including dark matter, may be ~1:10-1:5. This system may thus represent a "stealth" merger, where an infalling satellite galaxy is nearly undetectable by conventional means, yet has a substantial dynamical influence on its host galaxy. This singular discovery also suggests that satellite accretion can play a significant role in building up the stellar halos of low-mass galaxies, and possibly in triggering their starbursts.

The direct detection of dark matter on Earth depends crucially on its density and its velocity distribution on a milliparsec scale. Conventional N-body simulations are unable to access this scale, making the development of other approaches necessary. In this paper, we apply the method developed in Fantin et al. 2008 to a cosmologically-based merger tree, transforming it into a useful instrument to reproduce and analyse the merger history of a Milky Way-like system. The aim of the model is to investigate the implications of any ultra-fine structure for the current and next generation of directional dark matter detectors. We find that the velocity distribution of a Milky Way-like Galaxy is almost smooth, due to the overlap of many streams of particles generated by multiple mergers. Only the merger of a 10^10 Msun analyse can generate significant features in the ultra-local velocity distribution, detectable at the resolution attainable by current experiments.

We present absorption-line index gradients for a sample of S0 galaxies in the Fornax Cluster. The sample has been selected to span a wide range in galaxy mass, and the deep VLT-FORS2 spectroscopy allows us to explore the stellar populations all the way to the outer disk-dominated regions of these galaxies. We find that globally, in both bulges and disks, star formation ceased earliest in the most massive systems, as a further manifestation of downsizing. However, within many galaxies, we find an age gradient which indicates that star formation ended first in the outermost regions. Metallicity gradients, when detected, are always negative such that the galaxy centres are more metal-rich. This finding fits with a picture in which star formation continued in the central regions, with enriched material, after it had stopped in the outskirts. Age and metallicity gradients are correlated, suggesting that large differences in star formation history between the inner and outer parts of S0 galaxies yield large differences in their chemical enrichment. In agreement with previous results, we conclude that the radial variations in the stellar populations of S0 galaxies are compatible with the hypothesis that these galaxies are the descendants of spiral galaxies whose star formation has ceased. With the addition of radial gradient information, we are able to show that this shutdown of star formation occurred from the outside inward, with the later star formation in the central regions offering a plausible mechanism for enhancing the bulge light in these systems, as the transformation to more bulge-dominated S0 galaxies requires.

Boxy and peanut-shaped bulges are seen in about half of edge-on disc galaxies. Comparisons of the photometry and major-axis gas and stellar kinematics of these bulges to simulations of bar formation and evolution indicate that they are bars viewed in projection. If the properties of boxy bulges can be entirely explained by assuming they are bars, then this may imply that their hosts are pure disc galaxies with no classical bulge. A handful of these bulges, including that of the Milky Way, have been observed to rotate cylindrically, i.e. with a mean stellar velocity independent of height above the disc. In order to assess whether such behaviour is ubiquitous in boxy bulges, and whether a pure disc interpretation is consistent with their stellar populations, we have analysed the stellar kinematics and populations of the boxy or peanut-shaped bulges in a sample of five edge-on galaxies. We placed slits along the major axis of each galaxy and at three offset but parallel positions to build up spatial coverage. The boxy bulge of NGC3390 rotates perfectly cylindrically within the spatial extent and uncertainties of the data. This is consistent with the metallicity and alpha-element enhancement of the bulge, which are the same as in the disk. This galaxy is thus a pure disc galaxy. The boxy bulge of ESO311-G012 also rotates very close to cylindrically. The boxy bulge of NGC1381 is neither clearly cylindrically nor non-cylindrically rotating, but it has a negative vertical metallicity gradient and is alpha-enhanced with respect to its disc, suggesting a composite bulge comprised of a classical bulge and bar (and possibly a discy pseudobulge) [abridged] Even this relatively small sample is sufficient to demonstrate that boxy bulges display a range of rotational and population properties, indicating that they do not form a homogeneous class of object.

To investigate the origins of S0 galaxies, we present a new method of analyzing their stellar kinematics from discrete tracers such as planetary nebulae. This method involves binning the data in the radial direction so as to extract the most general possible non-parametric kinematic profiles, and using a maximum likelihood fit within each bin in order to make full use of the information in the discrete kinematic tracers. Both disk and spheroid kinematic components are fitted, with a two-dimensional decomposition of imaging data used to attribute to each tracer a probability of membership in the separate components. Likelihood clipping also allows us to identify objects whose properties are not consistent with the adopted model, rendering the technique robust against contaminants and able to identify additional kinematic features. The method is first tested on an N-body simulated galaxy to assess possible sources of systematic error associated with the structural and kinematic decomposition, which are found to be small. It is then applied to the S0 system NGC~1023, for which a planetary nebula catalogue has already been released and analyzed by (Noordermeer et al., 2008). The correct inclusion of the spheroidal component allows us to show that, contrary to previous claims, the stellar kinematics of this galaxy are indistinguishable from those of a normal spiral galaxy, indicating that it may have evolved directly from such a system via gas stripping or secular evolution. The method also successfully identifies a population of outliers whose kinematics are different from those of the main galaxy; these objects can be identified with a stellar stream associated with the companion galaxy NGC~1023A.

As part of our current programme to test LCDM predictions for dark matter (DM) haloes using extended kinematical observations of early-type galaxies, we present a dynamical analysis of the bright elliptical galaxy NGC 4374 (M84) based on ~450 Planetary Nebulae (PNe) velocities from the PN.Spectrograph, along with extended long-slit stellar kinematics. This is the first such analysis of a galaxy from our survey with a radially constant velocity dispersion profile. We find that the spatial and kinematical distributions of the PNe agree with the field stars in the region of overlap. The velocity kurtosis is consistent with zero at almost all radii. We construct a series of Jeans models, fitting both velocity dispersion and kurtosis to help break the mass-anisotropy degeneracy. Our mass models include DM halos either with shallow cores or with central cusps as predicted by cosmological simulations - along with the novel introduction in this context of adiabatic halo contraction from baryon infall. Both classes of models confirm a very massive dark halo around NGC 4374, demonstrating that PN kinematics data are well able to detect such haloes when present. Considering the default cosmological mass model, we confirm earlier suggestions that bright galaxies tend to have halo concentrations higher than LCDM predictions, but this is found to be solved if either a Salpeter IMF or adiabatic contraction with a Kroupa IMF is assumed. Thus for the first time a case is found where the PN dynamics may well be consistent with a standard dark matter halo. A cored halo can also fit the data, and prefers a stellar mass consistent with a Salpeter IMF. The less dramatic dark matter content found in lower-luminosity "ordinary" ellipticals suggests a bimodality in the halo properties which may be produced by divergent baryonic effects during their assembly histories.

This paper introduces the idea that the general mixing inequality obeyed by evolving stellar phase densities may place useful constraints on the possible history of the over-all galaxy population. We construct simple models for the full stellar phase space distributions of galaxies' disk and spheroidal components, and reproduce the well-known result that the maximum phase density of an elliptical galaxy is too high to be produced collisionlessly from a disk system, although we also show that the inclusion of a bulge component in the disk removes this evolutionary impediment. In order to draw more general conclusions about the evolution of the galaxy population, we use the Millennium Galaxy Catalogue to construct a model of the entire phase density distribution of stars in a representative sample of the local Universe. In such a composite population, we show that the mixing inequality rules out some evolutionary paths that are not prohibited by consideration of the maximum phase density alone, and thus show that the massive ellipticals in this population could not have formed purely from collisionless mergers of a low mass galaxy population like that found in the local Universe. Although the violation of the mixing inequality is in this case quite minor, and hence avoidable with a modest amount of non-collisionless star formation in the merger process, it does confirm the potential of this approach. The future measurement of stellar phase densities at higher redshift will allow this potential to be fully exploited, offering a new way to look at the possible pathways for galaxy evolution, and to learn about the environment of star formation through the way that this phase space becomes populated over time.

This paper briefly reviews some of the exciting studies of the local Universe that will be enabled with the European Extremely Large Telescope (E-ELT). As illustrative examples, it summarizes a few of the scientific goals that have been set for this instrument for studies of young starburst clusters, evolved stars, the Galactic centre, Galactic structure, nucleo-chronometry, the interpretation of the "Spite Plateau," and the properties of resolved stellar populations in external galaxies. It finishes on a note of warning that we really need to be even more innovative and adventurous if we are to justify the cost of the E-ELT.

Galaxies are faint enough when one observes just their light distributions, but in studying their full dynamical structure the stars are spread over the six dimensions of phase space rather than just the three spatial dimensions, making their densities very low indeed. This low signal is unfortunate, as stellar dynamics hold important clues to these systems' life histories, and the issue is compounded by the fact that the most interesting information comes from the faintest outer parts of galaxies, where dynamical timescales (and hence memories of past history) are longest. To extract this information, we have constructed a special-purpose instrument, the Planetary Nebula Spectrograph, which observes planetary nebulae as kinematic tracers of the stellar population, and allows one to study the stellar dynamics of galaxies down to extremely low surface brightnesses. Here, we present results from this instrument that illustrate how it can uncover the nature of low surface-brightness features such as thick disks by studying their kinematics, and trace faint kinematic populations that are photometrically undetectable.

Current theories of spiral and bar structure predict a variety of pattern speed behaviors, calling for detailed, direct measurement of the radial variation of pattern speeds. Our recently developed Radial Tremaine-Weinberg (TWR) method allows this goal to be achieved for the first time. Here we present TWR spiral pattern speed estimates for M101, IC 342, NGC 3938 and NGC 3344 in order to investigate whether spiral structure is steady or winding, whether spirals are described by multiple pattern speeds, and the relation between bar and spiral speeds. Where possible, we interpret our pattern speeds estimates according to the resonance radii associated with each (established with the disk angular rotation), and compare these to previous determinations. By analyzing the high-quality HI and CO data cubes available for these galaxies, we show that it is possible to determine directly multiple pattern speeds within these systems, and hence identify the characteristic signatures of the processes that drive the spiral structure. Even this small sample of galaxies reveals a surprisingly complex taxonomy, with the first direct evidence for the presence of resonant coupling of multiple patterns found in some systems, and the measurement of a simple single pattern speed in others. Overall, this study demonstrates that we are now in a position to uncover more of the apparently complex physics that lies behind spiral structure.

The procedure that is currently employed to allocate time on telescopes is horribly onerous on those unfortunate astronomers who serve on the committees that administer the process, and is in danger of complete collapse as the number of applications steadily increases. Here, an alternative is presented, whereby the task is distributed around the astronomical community, with a suitable mechanism design established to steer the outcome toward awarding this precious resource to those projects where there is a consensus across the community that the science is most exciting and innovative.

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.

Various laboratory-based experiments are underway attempting to detect dark matter directly. The event rates and detailed signals expected in these experiments depend on the dark matter phase space distribution on sub-milliparsec scales. These scales are many orders of magnitude smaller than those that can be resolved by conventional N-body simulations, so one cannot hope to use such tools to investigate the effect of mergers in the history of the Milky Way on the detailed phase-space structure probed by the current experiments. In this paper we present an alternative approach to investigating the results of such mergers, by studying a simplified model for a merger of a sub-halo with a larger parent halo. With an appropriate choice of parent halo potential, the evolution of material from the sub-halo can be expressed analytically in action-angle variables, so it is possible to obtain its entire orbit history very rapidly without numerical integration. Furthermore by evolving backwards in time, we can obtain arbitrarily-high spatial resolution for the current velocity distribution at a fixed point. Although this model cannot provide a detailed quantitative comparison with the Milky Way, its properties are sufficiently generic that it offers qualitative insight into the expected structure arising from a merger at a resolution that cannot be approached with full numerical simulations. Preliminary results indicate that the velocity-space distribution of dark matter particles remains characterized by discrete and well-defined peaks over an extended period of time, both for single and multi-merging systems, in contrast to the simple smooth velocity distributions sometimes assumed in predicting laboratory experiment detection rates.

The grand-design spiral galaxy M51 has long been a crucial target for theories of spiral structure. Studies of this iconic spiral can address the question of whether strong spiral structure is transient (e.g. interaction-driven) or long-lasting. As a clue to the origin of the structure in M51, we investigate evidence for radial variation in the spiral pattern speed using the radial Tremaine-Weinberg (TWR) method. We implement the method on CO observations tracing the ISM-dominant molecular component. Results from the method's numerical implementation--combined with regularization, which smooths intrinsically noisy solutions--indicate two distinct patterns speeds inside 4 kpc at our derived major axis PA=170 deg., both ending at corotation and both significantly higher than the conventionally adopted global value. Inspection of the rotation curve suggests that the pattern speed interior to 2 kpc lacks an ILR, consistent with the leading structure seen in HST near-IR observations. We also find tentative evidence for a lower pattern speed between 4 and 5.3 kpc measured by extending the regularized zone. As with the original TW method, uncertainty in major axis position angle (PA) is the largest source of error in the calculation; in this study, where \delta PA=+/-5 deg. a ~20% error is introduced to the parameters of the speeds at PA=170 deg. Accessory to this standard uncertainty, solutions with PA=175 deg. (also admitted by the data) exhibit only one pattern speed inside 4 kpc, and we consider this circumstance under the semblance of a radially varying PA.

In order to test the basic equations believed to dictate the dynamics of disk galaxies, we present and analyze deep two-dimensional spectral data obtained using the PPAK integral field unit for the early-type spiral systems NGC 2273, NGC 2985, NGC 3898 and NGC 5533. We describe the care needed to obtain and process such data to a point where reliable kinematic measurements can be obtained from these observations, and a new more optimal method for deriving the rotational motion and velocity dispersions in such disk systems. The data from NGC 2273 and NGC 2985 show systematic variations in velocity dispersion with azimuth, as one would expect if the shapes of their velocity ellipsoids are significantly anisotropic, while the hotter disks in NGC 3898 and NGC 5533 appear to have fairly isotropic velocity dispersions. Correcting the rotational motion for asymmetric drift using the derived velocity dispersions reproduces the rotation curves inferred from emission lines reasonably well, implying that this correction is quite robust, and that the use of the asymmetric drift equation is valid. NGC 2985 is sufficiently close to face on for the data, combined with the asymmetric drift equation, to determine all three components of the velocity ellipsoid. The principal axes of this velocity ellipsoid are found to be in the ratio sigma_z:sigma_phi:sigma_R ~ 0.7:0.7:1, which shows unequivocally that this disk distribution function respects a third integral of motion. The ratio is also consistent with the predictions of epicyclic theory, giving some confidence in the application of this approximation to even fairly early-type disk galaxies.

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)

Using high signal-to-noise ratio VLT/FORS2 long-slit spectroscopy, we have studied the properties of the central stellar populations and dynamics of a sample of S0 galaxies in the Fornax Cluster. The central absorption-line indices in these galaxies correlate well with the central velocity dispersions (Sigma0) in accordance with what previous studies found for elliptical galaxies. However, contrary to what it is usually assumed for cluster ellipticals, the observed correlations seem to be driven by systematic age and alpha-element abundance variations, and not changes in overall metallicity. We also found that the observed scatter in the Index-Sigma0 relations can be partially explained by the rotationally-supported nature of these systems. Indeed, even tighter correlations exist between the line indices and the maximum circular velocity of the galaxies. This study suggests that the dynamical mass is the physical property driving these correlations, and for S0 galaxies such masses have to be estimated assuming a large degree of rotational support. The observed trends imply that the most massive S0s have the shortest star-formation timescales and the oldest stellar populations.

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.

At the intersection of galactic dynamics, evolution and global structure, issues such as the relation between bars and spirals and the persistence of spiral patterns can be addressed through the characterization of the angular speeds of the patterns and their possible radial variation. The Radial Tremaine-Weinberg (TWR) Method, a generalized version of the Tremaine-Weinberg method for observationally determining a single, constant pattern speed, allows the pattern speed to vary arbitrarily with radius. Here, we perform tests of the TWR method with regularization on several simulated galaxy data sets. The regularization is employed as a means of smoothing intrinsically noisy solutions, as well as for testing model solutions of different radial dependence (e.g. constant, linear or quadratic). We test these facilities in studies of individual simulations, and demonstrate successful measurement of both bar and spiral pattern speeds in a single disk, secondary bar pattern speeds, and spiral winding (in the first application of a TW calculation to a spiral simulation). We also explore the major sources of error in the calculation and find uncertainty in the major axis position angle most dominant. In all cases, the method is able to extract pattern speed solutions where discernible patterns exist to within 20% of the known values, suggesting that the TWR method should be a valuable tool in the area of galactic dynamics. For utility, we also discuss the caveats in, and compile a prescription for, applications to real galaxies.

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.

We present a study of pixel Colour Magnitude Diagrams (pCMDs) for a sample of 69 nearby galaxies chosen to span a wide range of Hubble types. Our goal is to determine how useful a pixel approach is for studying galaxies according to their stellar light distributions and content. The galaxy images were analysed on a pixel-by-pixel basis to reveal the structure of the individual pCMDs. We find that the average surface brightness (or projected mass density) in each pixel varies according to galaxy type. Early-type galaxies exihibit a clear ``prime sequence'' and some pCMDs of face-on spirals reveal ``inverse-L'' structures. We find that the colour dispersion at a given magnitude is found to be approximately constant in early-type galaxies but this quantity varies in the mid and late-types. We investigate individual galaxies and find that the pCMDs can be used to pick out morphological features. We discuss the discovery of ``Red Hooks'' in the pCMDs of six early-type galaxies and two spirals and postulate their origins. We develop quantitative methods to characterise the pCMDs, including measures of the blue-to-red light ratio and colour distributions of each galaxy and we organise these by morphological type. We compare the colours of the pixels in each galaxy with the stellar population models of Bruzual & Charlot (2003) to calculate star formation histories for each galaxy type and compare these to the stellar mass within each pixel. Maps of pixel stellar mass and mass-to-light ratio are compared to galaxy images. We apply the pCMD technique to three galaxies in the Hubble Ultra Deep Field to test the usefulness of the analysis at high redshift. We propose that these results can be used as part of a new system of automated classification of galaxies that can be applied at high redshift.

We test the theory that lenticular (S0) galaxies form from spirals whose star formation has been shut down. We use the globular cluster specific frequency S_N, defined as the number of globular clusters normalised by the galaxy luminosity as a diagnostic. NTT/EMMI long-slit spectroscopic observations of 11 S0 galaxies at z < 0.006 are used to measure the absorption-line indices, Hdelta, Hgamma, Mgb, Fe5270 and Fe5335 within the central r_e/8. By inverting single-stellar population models, luminosity-weighted mean ages, metallicities and alpha-element abundance ratios are derived. We estimate the amount of fading a galaxy has undergone by comparing each galaxy's S_N with its deviation from the mean spiral S_N. Galaxies with higher S_N have older stellar populations. Moreover, we find that the zero-point and amount of fading is consistent with a scenario where lenticulars are formed by the quenching of star formation in spiral galaxies. Our data also rule out any formation method for S0s which creates a large number of new globular clusters. We confirm that previous results showing a relationship between S_N and color are driven by the S_N - Age relation. Five galaxies show detectable Hbeta, [OIII], Halpha or [NII] emission lines. However, only in the two youngest galaxies is this emission unambiguously from star formation. Our results are consistent with the theory that S0 galaxies are formed when gas in normal spirals is removed, possibly as a result of a change in environment. The on-going star formation in the youngest galaxies hints that the timescale of quenching is ~< 1 Gyr. We speculate, therefore, that the truncation of star formation is a rather gentle process unlikely to involve a rapid burst of star formation.

We present results from Planetary Nebula Spectrograph (PN.S) observations of the elliptical galaxy NGC 3379 and a description of the data reduction pipeline. We detected 214 planetary nebulae of which 191 are ascribed to NGC 3379, and 23 to the companion galaxy NGC 3384. Comparison with data from the literature show that the PN.S velocities have an internal error of <20km/s and a possible offset of similar magnitude. We present the results of kinematic modelling and show that the PN kinematics are consistent with absorption-line data in the region where they overlap. The resulting combined kinematic data set, running from the center of NGC 3379 out to more than seven effective radii (Reff), reveals a mean rotation velocity that is small compared to the random velocities, and a dispersion profile that declines rapidly with radius. From a series of Jeans dynamical models we find the B-band mass-to-light ratio inside 5 Reff to be 8 to 12 in solar units, and the dark matter fraction inside this radius to be less than 40%. We compare these and other results of dynamical analysis with those of dark-matter-dominated merger simulations, finding that significant discrepancies remain, reiterating the question of whether NGC 3379 has the kind of dark matter halo that the current LambdaCDM paradigm requires.

Based on FORS2-VLT long-slit spectroscopy, the analysis of the central absorption line indices of 9 S0 galaxies in the Fornax Cluster is presented. Central indices correlate with central velocity dispersions as observed in ellipticals. However, the stellar population properties of these S0s indicates that the observed trends are produced by relative differences in age and alpha-element abundances and not in metallicity ([Fe/H]) as previous studies have found in elliptical galaxies. The observed scatter in the line indices vs. velocity dispersion relations can be partially explained by the rotationally-supported nature of many of these systems. The presence of tighter line indices vs. maximum (circular) rotational velocity relations confirms this statement. It was also confirmed that the dynamical mass is the driving physical property of all these correlations and in our Fornax S0s it has to be estimated assuming rotational support.

We present a study of the local B and K-band Tully-Fisher Relation (TFR) between absolute magnitude and maximum circular speed in S0 galaxies. To make this study, we have combined kinematic data, including a new high-quality spectral data set from the Fornax Cluster, with homogeneous photometry from the RC3 and 2MASS catalogues, to construct the largest sample of S0 galaxies ever used in a study of the TFR. Independent of environment, S0 galaxies are found to lie systematically below the TFR for nearby spirals in both optical and infrared bands. This offset can be crudely interpreted as arising from the luminosity evolution of spiral galaxies that have faded since ceasing star formation. However, we also find a large scatter in the TFR. We show that most of this scatter is intrinsic, not due to the observational uncertainties. The presence of such a large scatter means that the population of S0 galaxies cannot have formed exclusively by the above simple fading mechanism after all transforming at a single epoch. To better understand the complexity of the transformation mechanism, we have searched for correlations between the offset from the TFR and other properties of the galaxies such as their structural properties, central velocity dispersions and ages (as estimated from line indices). For the Fornax Cluster data, the offset from the TFR relates with the estimated age of the stars in the individual galaxies, in the sense and of the magnitude expected if S0 galaxies had passively faded since being converted from spirals. This correlation implies that a significant part of the scatter in the TFR arises from the different times at which galaxies began their transformation.

We test the hypothesis that S0 galaxies are the descendants of fading spirals whose star formation has been shut down, by using the properties of their globular cluster (GC) systems. We estimate the amount by which the GC specific frequency (number of GCs per unit V-band luminosity) is enhanced in S0s relative to spirals. If the transformation hypothesis is correct, and no GCs are created or destroyed in the process, then this difference provides a measure of the degree to which the S0's V-band luminosity has faded relative to that of its spiral progenitor. We also explore whether the degree to which the GC specific frequency is enhanced in S0s correlates with the colour of the stellar population, as also predicted by this hypothesis in which galaxies become redder as they fade. We find that, on average, the GC specific frequency is a factor ~3 larger for S0s than for spirals, which can be interpreted as meaning that passively-evolving S0s have faded by about a factor of 3 from their spiral progenitors. This value fits remarkably well with the predictions of stellar population synthesis calculations, and the offset between the S0 and spiral Tully-Fisher relations. We also find that the global colours of S0 galaxies are strongly correlated with their GC specific frequencies: the redder the stellar population of an S0, the larger its specific frequency, as we might expect if we are catching different S0s at different stages of passively fading and reddening. Comparison to the predictions of stellar population synthesis models show that this explanation works quantitatively as well as qualitatively. These tests strongly support the hypothesis that S0 galaxies were once normal spirals, whose star formation was cut off, presumably due to a change of environment.

We have obtained long-slit spectroscopy for a sample of 9 S0 galaxies in the Fornax Cluster using the FORS2 spectrograph at the 8.2m ESO VLT. From these data, we have extracted the kinematic parameters, comprising the mean velocity, velocity dispersion and higher-moment H3 and H4 coefficients, as a function of position along the major axes of these galaxies. Comparison with published kinematics indicates that earlier data are often limited by their lower signal-to-noise ratio and relatively poor spectral resolution. The greater depth and higher dispersion of the new data mean that we reach well beyond the bulges of these systems, probing their disk kinematics in some detail for the first time. Qualitative inspection of the results for individual galaxies shows that they are not entirely simple systems, perhaps indicating a turbulent past. Nonetheless, we are able to derive reliable circular velocities for most of these systems, which points the way toward a study of their Tully-Fisher relation. This study, along with an analysis of the stellar populations of these systems out to large galactocentric distances, will form the bases of future papers exploiting these new high-quality data, hopefully shedding new light on the evolutionary history of these systems.

We present radial velocities for a sample of 723 planetary nebulae (PNe) in the disk and bulge of M31, measured using the WYFFOS fibre spectrograph on the William Herschel telescope. Velocities are determined using the [OIII] 5007 Angstrom emission line. Rotation and velocity dispersion are measured to a radius of 50 arcminutes (11.5 kpc), the first stellar rotation curve and velocity dispersion profile for M31 to such a radius. Our kinematics are consistent with rotational support at radii well beyond the bulge effective radius of 1.4kpc, although our data beyond a radius of 5kpc are limited. We present tentative evidence for kinematic substructure in the bulge of M31 to be studied fully in a later work. This paper is part of an ongoing project to constrain the total mass, mass distribution and velocity anisotropy of the disk, bulge and halo of M31.

We present a catalogue of positions, magnitudes and velocities for 3300 emission-line objects found by the Planetary Nebula Spectrograph in a survey of the Andromeda Galaxy, M31. Of these objects, 2615 are found likely to be planetary nebulae (PNe) associated with M31. Initial results from this survey include: the likely non-existence of Andromeda VIII; a universal PN luminosity function, with the exception of a small amount of obscuration, and a small offset in normalization between bulge and disk components; very faint kinematically-selected photometry implying no cut-off in the disk to beyond 4 scalelengths and no halo population in excess of the bulge out to 10 effective bulge radii; disk kinematics that show significant dispersion and asymmetric drift out to large radii, consistent with a warm flaring disk; and no sign of any variation in kinematics with PN luminosity, suggesting that PNe arise from a fairly uniform population of old stars.

The lifetime of the structure in grand design spiral galaxies is observationally ill-determined, but is essentially set by how accurately the pattern's rotation can be characterized by a single angular pattern speed. This paper derives a generalized version of the Tremaine-Weinberg method for observationally determining pattern speeds, in which the pattern speed is allowed to vary arbitrarily with radius. The departures of the derived pattern speed from a constant then provides a simple metric of the lifetime of the spiral structure. Application of this method to CO observations of NGC 1068 reveal that the pattern speed of the spiral structure in this galaxy varies rapidly with radius, and that the lifetime of the spiral structure is correspondingly very short. If this result turns out to be common in grand-design spiral galaxies, then these features will have to be viewed as highly transient phenomena.

We have obtained Integral Field Spectroscopy for 23 spiral bulges using INTEGRAL on the William Herschel Telescope and SPIRAL on the Anglo-Australian Telescope. This is the first 2D survey directed solely at the bulges of spiral galaxies. Eleven galaxies of the sample do not have previous measurements of the stellar velocity dispersion (sigma*). These data are designed to complement our Space Telescope Imaging Spectrograph program for estimating black hole masses in the range 10^6-10^8M_sun using gas kinematics from nucleated disks. These observations will serve to derive the stellar dynamical bulge properties using the traditional Mgb and CaII triplets. We use both Cross Correlation and Maximum Penalized Likelihood to determine projected sigma* in these systems and present radial velocity fields, major axis rotation curves, curves of growth and sigma* fields. Using the Cross Correlation to extract the low order 2D stellar dynamics we generally see coherent radial rotation and irregular velocity dispersion fields suggesting that sigma* is a non-trivial parameter to estimate.

We present XMM-Newton observations of the wake-radiogalaxy system 4C34.16, which shows a cool and dense wake trailing behind 4C34.16's host galaxy. A comparison with numerical simulations is enlightening, as they demonstrate that the wake is produced mainly by ram pressure stripping during the galactic motion though the surrounding cluster. The mass of the wake is a substantial fraction of the mass of an elliptical galaxy's X-ray halo. This observational fact supports a wake formation scenario similar to the one demonstrated numerically by Acreman et al (2003): the host galaxy of 4C34.16 has fallen into its cluster, and is currently crossing its central regions. A substantial fraction of its X-ray halo has been stripped by ram pressure, and remains behind to form the galaxy wake.

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.

This paper presents a brief review of the evidence for dark matter in the Universe on the scales of galaxies. In the interests of critically and objectively testing the dark matter paradigm on these scales, this evidence is weighed against that from the only other game in town, modified Newtonian dynamics. The verdict is not as clear cut as one might have hoped.

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

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.

Using the Planetary Nebula Spectrograph, we have observed and measured the velocities for some 2764 PNe in the disk and halo of the Andromeda galaxy. Preliminary analysis using a basic ring model shows a rotation curve in good agreement with that obtained from HI data out to ~20kpc. Some substructure has also been detected within the velocity field, which can be modeled as the continuation of the tidal-remnant known as the Southern Stream, as it passes through Andromeda's disk.

This paper looks for any correlation between the internal dynamics of elliptical galaxies and the relatively mild nuclear activity found in many such systems. We show that there is such a relation in the sense that the active ellipticals tend to be significantly less rotationally supported than their inactive cousins. The correlation can partly be related to the galaxies' luminosities: the brightest galaxies tend to be more active and less rotationally supported. However, even at lower luminosities the active and inactive galaxies seem to have systematically different dynamics. This variation suggests that there are significant large-scale structural differences between active and inactive elliptical galaxies, and hence that the existence of both types of system cannot just be the result of random sporadic nuclear activity.

The Planetary Nebula Spectrograph is a dedicated instrument for measuring radial velocity of individual Planetary Nebulae (PNe) in galaxies. This new instrument is providing crucial data with which to probe the structure of dark halos in the outskirts of elliptical galaxies in particular, which are traditionally lacking of easy interpretable kinematical tracers at large distance from the center. Preliminary results on a sample of intermediate luminosity galaxies have shown little dark matter within 5 ~ R_eff implying halos either not as massive or not as centrally concentrated as CDM predicts (Romanowsky et al. 2003). We briefly discuss whether this is consistent with a systematic trend of the dark matter content with the luminosity as observed in an extended sample of early-type galaxies.

We present a possible orbit for the Southern Stream of stars in M31, which connects it to the Northern Spur. Support for this model comes from the dynamics of planetary nebulae (PNe) in the disk of M31: analysis of a new sample of 2611 PNe obtained using the Planetary Nebula Spectrograph reveals ~20 objects whose kinematics are inconsistent with the normal components of the galaxy, but which lie at the right positions and velocities to connect the two photometric features via this orbit. The satellite galaxy M32 is coincident with the stream both in position and velocity, adding weight to the hypothesis that the stream comprises its tidal debris.

A new instrument is providing crucial data with which to probe the structure of dark halos in elliptical galaxies.

Planetary nebulae (PNe) may be the most promising tracers in the halos of early-type galaxies. We have used multi-object spectrographs on the WHT and the VLT, and the new Planetary Nebula Spectrograph on the WHT, to obtain hundreds of PN velocities in a small sample of nearby galaxies. These ellipticals show weak halo rotation, which may be consistent with ab initio models of galaxy formation, but not with more detailed major merger simulations. The galaxies near L* show evidence of a universal declining velocity dispersion profile, and dynamical models indicate the presence of little dark matter within 5 Reff -- implying halos either not as massive or not as centrally concentrated as CDM predicts.

We have found a correlation betwen the M/L global gradients and the structural parameters of the luminous components of a sample of 19 early-type galaxies. Such a correlation supports the hypothesis that there is a connection between the dark matter content and the evolution of the baryonic component in such systems.

This paper reviews the available information on the central density distribution and shape of the Milky Way's halo. At present, there is no strong evidence that the Milky Way's halo properties conflict with the predictions of cold dark matter (CDM): a primordial central power law cusp can be accommodated by the observations, and the current constraints on flattening are also consistent with the predictions of the theory. If you want to pick a fight with CDM, then the Milky Way is probably not the place to do it.

The kinematics of the outer parts of three intermediate-luminosity elliptical galaxies have been studied using the Planetary Nebula Spectrograph. The galaxies' velocity dispersion profiles are found to decline with radius; dynamical modeling of the data indicates the presence of little if any dark matter in these galaxies' halos. This surprising result conflicts with findings in other galaxy types, and poses a challenge to current galaxy formation theories.

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.

Like the majority of spiral galaxies, the Milky Way contains a central non-axisymmetric bar component. Our position in the Galactic plane renders it rather hard to see, but also allows us to make measurements of the bar that are completely unobtainable for any other system. This paper reviews the evidence for a bar that can be gleaned from the many extensive surveys of both gas and stars in the Milky Way. We introduce some simplified models to show how the basic properties of the bar can be inferred in a reasonably robust manner despite our unfavorable location, and how the complex geometry can be used to our advantage to obtain a unique three-dimensional view of the bar. The emerging picture of the Galactic bar is also placed in the broader context of current attempts to understand how such structures form and evolve in spiral galaxies.

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.

Planetary nebulae are now well established as probes of galaxy dynamics and as standard candles in distance determinations. Motivated by the need to improve the efficiency of planetary nebulae searches and the speed with which their radial velocities are determined, a dedicated instrument - the Planetary Nebulae Spectrograph or PN.S - has been designed and commissioned at the 4.2m William Herschel Telescope. The high optical efficiency of the spectrograph results in the detection of typically ~ 150 PN in galaxies at the distance of the Virgo cluster in one night of observations. In the same observation the radial velocities are obtained with an accuracy of ~ 20 km/s

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.

We have used the publicly available data from the 2dF Galaxy Redshift Survey and the 2dF QSO Redshift Survey to test the hypothesis that there is a periodicity in the redshift distribution of quasi-stellar objects (QSOs) found projected close to foreground galaxies. These data provide by far the largest and most homogeneous sample for such a study, yielding 1647 QSO-galaxy pairs. There is no evidence for a periodicity at the predicted frequency in log(1+z), or at any other frequency.

We report here the successful commissioning of the PN.Spectrograph, the first special-purpose instrument for the measurement of galaxy kinematics through the PN population

To quantify the distribution of bar shapes in spiral galaxies, we have analysed 113 H-band and 89 B-band galaxy images from the Ohio State University Bright Spiral Galaxy Survey. Parameters measuring bar shape and position along the Hubble sequence were obtained in each waveband. Evidence was found for a bimodality in the distribution of bar shape, implying that barred and unbarred galaxies are not just the extrema of a single distribution, and that any evolution between these two states must occur on a rapid timescale. Objective bar shapes measured in the H-band were found to be more closely related to visual classifications than B-band bar strengths, as the B-band images are somewhat compromised by localised star formation, especially in later-type systems. Galaxies were found to be more centrally concentrated in the infrared. Later type galaxies showed greater asymmetry in the optical than the infrared, presumably again due to localised star formation, but on average the bar strengths in the two bands were found to be the same.

We investigate the visibility of galactic bars and spiral structure in the distant Universe by artificially redshifting 101 B-band CCD images of local spiral galaxies from the Ohio State University Bright Spiral Galaxy Survey. Our artificially redshifted images correspond to Hubble Space Telescope I-band observations of the local galaxy sample seen at z=0.7, with integration times matching those of both the very deep Northern Hubble Deep Field data, and the much shallower Flanking Field observations. The expected visibility of galactic bars is probed in two ways: (1) using traditional visual classification, and (2) by charting the changing shape of the galaxy distribution in "Hubble space", a quantitative two-parameter description of galactic structure that maps closely on to Hubble's original tuning fork. Both analyses suggest that over 2/3 of strongly barred luminous local spirals i.e. objects classified as SB in the Third Reference Catalog) would still be classified as strongly barred at z=0.7 in the Hubble Deep Field data. Under the same conditions, most weakly barred spirals (classified SAB in the Third Reference Catalog) would be classified as regular spirals. The corresponding visibility of spiral structure is assessed visually, by comparing luminosity classifications for the artificially redshifted sample with the corresponding luminosity classifications from the Revised Shapley Ames Catalog. We find that for exposures times similar to that of the Hubble Deep Field spiral structure should be detectable in most luminous low-inclination spiral galaxies at z=0.7 in which it is present. [ABRIDGED]

This paper describes a framework for studying galaxy morphology, particularly bar strength, in a quantitative manner, and presents applications of this approach that reveal observational evidence for secular evolution in bar morphology. The distribution of bar strength in galaxies is quite strongly bimodal, suggesting that barred and unbarred systems are distinct entities, and that any evolution between these two states must occur on a relatively rapid timescale. Bars' strengths appear to be correlated with their pattern speeds, implying that these structures weaken as they start to slow, and disappear entirely before the bars have slowed significantly. There is also tantalizing evidence that bars are rare beyond a redshift of z ~ 0.7, indicating that galaxies have only recently evolved to a point where bars can readily form.

This paper investigates the detailed dynamical properties of a relatively homogeneous sample of disc-dominated S0 galaxies, with a view to understanding their formation, evolution and structure. By using high signal-to-noise ratio long-slit spectra of edge-on systems, we have been able to reconstruct the complete line-of-sight velocity distributions of stars along the galaxies' major axes. From these data, we have derived both model distribution functions (the phase density of their stars) and the approximate form of their gravitational potentials. The derived distribution functions are all consistent with these galaxies being simple disc systems, with no evidence for a complex formation history. Essentially no correlation is found between the characteristic mass scale-lengths and the photometric scale-lengths in these galaxies, suggesting that they are dark-matter dominated even in their inner parts. Similarly, no correlation is found between the mass scale-lengths and asymptotic rotation speed, implying a wide range of dark matter halo properties. By comparing their asymptotic rotation speeds with their absolute magnitudes, we find that these S0 galaxies are systematically offset from the Tully-Fisher relation for later-type galaxies. The offset in luminosity is what one would expect if star formation had been suddenly switched off a few Gyrs ago, consistent with a simple picture in which these S0s were created from ordinary later-type spirals which were stripped of their star-forming ISM when they encountered a dense cluster environment.

This paper reviews the constraints that can be placed on the shapes of disk galaxies' dark halos using the distribution and kinematics of atomic hydrogen. These data indicate that dark halos are close to axisymmetric, with their axes of symmetry co-aligned with their disk axes. They also appear to be oblate, with shortest-to-longest axis ratios displaying quite a broad range of values from ~0.2 to ~0.8. These results are consistent with the predicted shapes of halos in cold dark matter scenarios, but rule out some of the more exotic dark matter candidates. However, the total number of measurements is still depressingly small, and more data are required if halo shape is to become a powerful diagnostic for theories of galaxy formation and evolution.

(Abridged) Axisymmetric models of the Milky Way exhibit strong interrelations between the Galactic constants (R_0 and T_0), the stellar columndensity (S_*) and the shape of the dark matter (DM) halo. Here we present analytical relations that can be used to investigate the effects of the uncertain gaseous velocity dispersion on the HI flaring constraints. The contribution of cosmic rays and magnetic fields to the pressure gradients is small. A significantly flattened dark matter halo is only possible if R_0 <~ 6.8 kpc. If R_0 is larger than ~7 kpc, or T_0 >~ 170 km/s, we can rule out two DM candidates that require a highly flattened DM halo: 1) decaying massive neutrinos; and 2) a disk of cold molecular hydrogen. It is only possible to construct self-consistent models of the Galaxy based on the IAU-recommended values for the Galactic constants in the unlikely case that the the stellar columndensity is smaller than ~18 M_sun/pc^2. If we assume that the halo is oblate and S_* = 35 +/- 5 M_sun/pc^2, R_0 <~ 8 kpc and T_0 <~ 200 km/s. Combining the best kinematical and star-count estimates of S_*, we conclude that: 25 <~ S_* <~ 45 M_sun/pc^2. Kuijken & Gilmore's (1991) determination of the columndensity of matter with |z|<=1.1 kpc is robust and valid over a wide range of Galactic constants. Our mass models show that the DM density in the Galactic centre is uncertain by a factor 1000. In the Solar neighbourhood we find: rho_DM ~0.42 GeV/c^2/cm^3 or (11 +/- 5) mM_sun/pc^3 -- roughly 15% of rho_tot.

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.

These lectures are intended to provide an introduction to the rich interplay between N-body simulations and stellar-kinematic observations of galaxies. The first section describes the kinematic properties of galaxies that are accessible to observation, and gives a brief introduction to the stellar dynamics that one might wish to infer from such observations. In the next section, the use of N-body simulations in the study of galaxy dynamics is placed in its historical context. The next two sections provide case studies of the stellar-dynamical properties of elliptical and disk galaxies where N-body simulations have been used to explain kinematic data, to test new analysis techniques, and to motivate further observations. The final section contains some speculation as to how the symbiotic relation between observations and N-body simulations may develop in future studies of galaxy dynamics.

In order to establish an objective framework for studying galaxy morphology, we have developed a quantitative two-parameter description of galactic structure that maps closely on to Hubble's original tuning fork. Any galaxy can be placed in this "Hubble space", where the x-coordinate measures position along the early-to-late sequence, while the y-coordinate measures in a quantitative way the degree to which the galaxy is barred. The parameters defining Hubble space are sufficiently robust to allow the formation of Hubble's tuning fork to be mapped out to high redshifts. In the present paper, we describe a preliminary investigation of the distribution of local galaxies in Hubble space, based on the CCD imaging atlas of Frei et al. (1996). We find that barred, weakly-barred, and unbarred galaxies are remarkably well-separated on this diagnostic diagram. The spiral sequence is clearly bimodal and indeed approximates a tuning fork: strongly-barred and unbarred spirals do not simply constitute the extrema of a smooth unimodal distribution of bar strength, but rather populate two parallel sequences. Strongly barred galaxies lie on a remarkably tight sequence, strongly suggesting the presence of an underlying unifying physical process. Rather surprisingly, weakly barred systems do not seem to correspond to objects bridging the parameter space between unbarred and strongly barred galaxies, but instead form an extension of the regular spiral sequence. This relation lends support to models in which the bulges of late-type spirals originate from secular processes driven by bars.

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.

We present a method for recovering the distribution functions of edge-on thin axisymmetric disks directly from their observable kinematic properties. The most generally observable properties of such a stellar system are the line-of-sight velocity distributions of the stars at different projected radii along the galaxy. If the gravitational potential is known, then the general two-integral distribution function can be reconstructed using the shapes of the high-velocity tails of these line-of-sight distributions. If the wrong gravitational potential is adopted, then a distribution function can still be constructed using this technique, but the low-velocity parts of the observed velocity distributions will not be reproduced by the derived dynamical model. Thus, the gravitational potential is also tightly constrained by the observed kinematics.

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.

We present evidence that there is a significant correlation between the fraction of a galaxy's mass that lies in its central black hole and the age of the galaxy's stellar population. Since the absorption line indices that are used to estimate the age are luminosity weighted, they essentially measure the time since the last significant episode of star formation in the galaxy. The existence of this correlation is consistent with several theories of galaxy formation, including the currently-favoured hierarchical picture of galaxy evolution, which predicts just such a relation between black hole mass and the time since the last burst of merger-induced star formation. It is not consistent with models in which the massive black hole is primordial, and hence uncoupled from the stellar properties of the galaxy.

`Wakes' of X-ray emission have now been detected trailing behind a few (at least seven) elliptical galaxies in clusters. To quantify how widespread this phenomenon is, and what its nature might be, we have obtained a deep (70 ksec) X-ray image of the poor cluster Abell 160 using the ROSAT HRI. Combining the X-ray data with optical positions of confirmed cluster members, and applying a statistic designed to search for wake-like excesses, we confirm that this phenomenon is observed in galaxies in this cluster. The probability that the detections arise from chance is less than 0.0038. Further, the wakes are not randomly distributed in direction, but are preferentially oriented pointing away from the cluster centre. This arrangement can be explained by a simple model in which wakes arise from the stripping of their host galaxies' interstellar media due to ram pressure against the intracluster medium through which they travel.

To investigate the origins of wide-angle tailed radio sources (WATs), we have compiled a sample of these systems in Abell clusters for which X-ray data exist. Contrary to conventional wisdom, the WATs are found to be significantly displaced from the X-ray centroids of their host clusters. The bends in the WATs' radio jets are found to be oriented preferentially such that they point directly away from or toward the cluster centre, with more of the former than the latter. If this morphology is attributed to ram pressure, then the WATs are on primarily radial orbits, with more approaching the X-ray centroid than receding. There is also some evidence that the in-coming WATs are on average further from the X-ray centroid than the out-going ones. All of these observations strongly support a scenario in which WATs are created in cluster mergers.

In order to test the reliability of determinations of the shapes of galaxies' dark matter halos, we have made such measurements for the Milky Way by two independent methods, which make use of the stellar kinematics in the solar neighbourhood and the observed flaring of the Galactic HI layer to estimate the flattening of the Galactic dark halo. These techniques are found to produce a consistent estimate for the halo shape, with a shortest-to-longest axis ratio of q ~ 0.8, but only if one adopts somewhat non-standard values for the distance to the Galactic centre, R_0, and the local Galactic rotation speed, Theta_0. For consistency, one requires values of R_0 < 7.6 kpc and Theta_0 < 190 km/s. Although differing significantly from the current IAU-sanctioned values, these upper limits are consistent with all existing observational constraints. If future measurements confirm these lower values for the Galactic constants, then the validity of the gas layer flaring method will be confirmed. Further, dark matter candidates such as cold molecular gas and massive decaying neutrinos, which predict very flat dark halos with q < 0.2, will be ruled out. Conversely, if the Galactic constants were found to be close to the more conventional values, then there would have to be some systematic error in the methods for measuring dark halo shapes, so the existing modeling techniques would have to be viewed with some scepticism.

It has been suggested that the boxy and peanut-shaped bulges found in some edge-on galaxies are galactic bars viewed from the side. We investigate this hypothesis by presenting emission-line spectra for a sample of 10 edge-on galaxies that display a variety of bulge morphologies. To avoid potential biases in the classification of this morphology, we use an objective measure of bulge shape. Generally, bulges classified as more boxy show the more complicated kinematics characteristic of edge-on bars, confirming the intimate relation between the two phenomena.

The frequency of barred spiral galaxies as a function of redshift contains important information on the gravitational influence of stellar disks in their dark matter halos and also may distinguish between contemporary theories for the origin of galactic bulges. In this paper we present a new quantitative method for determining the strength of barred spiral structure, and verify its robustness to redshift-dependent effects. By combining galaxy samples from the Hubble Deep Field North with newly available data from the Hubble Deep Field South, we are able to define a statistical sample of 18 objectively-defined low-inclination barred spiral systems with $I_{814W}<23.2$ mag. Analysing the proportion of barred spiral galaxies seen as a function of redshift, we find a significant decline in the barred fraction beyond redshifts $z\simeq 0.5$. The physical significance of this effect remains unclear, but several possibilities include dynamically hotter (or increasingly dark-matter dominated) high-redshift discs, or an enhanced efficiency in bar destruction at high redshifts. By investigating the formation of the ``orthogonal'' axis of Hubble's classification tuning fork, our result complements studies of evolution in the early--late sequence, and pushes to later epochs the redshift at which the Hubble classification sequence is observed to be in place.

The disks of galaxies are primarily stellar systems, and fundamentally dynamical entities. Thus, to fully understand galactic disks, we must study their stellar kinematics as well as their morphologies. Observational techniques have now advanced to a point where quite detailed stellar-kinematic information can be extracted from spectral observations. This review presents three illustrative examples of analyses that make use of such information to study the formation and evolution of these systems: the derivation of the pattern speed of the bar in NGC 936; the calculation of the complete velocity ellipsoid of random motions in NGC 488; and the strange phenomenon of counter-rotation seen in NGC 3593.

We show that the kinematics of the shells seen around some elliptical galaxies provide a new, independent means for measuring the gravitational potentials of elliptical galaxies out to large radii. A numerical simulation of a set of shells formed in the merger between an elliptical and a smaller galaxy reveals that the shells have a characteristic observable kinematic structure, with the maximum line-of-sight velocity increasing linearly as one moves inward from a shell edge. A simple analytic calculation shows that this structure provides a direct measure of the gradient of the gravitational potential at the shell radius. In order to extract this information from attainable data, we have also derived the complete distribution of line-of-sight velocities for material within a shell; comparing the observed spectra of a shell to a stellar template convolved with this distribution will enable us to measure the gradient of the potential at this radius. Repeating the analysis for a whole series of nested shells in a galaxy allows the complete form of the gravitational potential as a function of radius to be mapped out. The requisite observations lie within reach of the up-coming generation of large telescopes.

The local stellar kinematics of the Milky Way offer a useful tool for studying the rotation curve of the Galaxy. These kinematics -- usually parameterized by the Oort constants A and B -- depend on the local gradient of the rotation curve as well as its absolute value (Theta_0), and the Sun's distance to the Galactic center (R_0). The density of interstellar gas in the Milky Way is shown to vary non-monotonically with radius, and so contributes significantly to the local gradient of the rotation curve. We have therefore calculated mass models for the Milky Way that include this component, and have derived the corresponding radial variation in the Oort constants. Between 0.9R_0 and 1.2R_0, the Oort functions, A(R) and B(R), differ significantly from the general Theta/R dependence. Various previously-inexplicable observations are shown to be consistent with these predictions. These models can explain the 40% difference between the values for 2 A R_0 derived from radial velocity data originating in the inner and outer Galaxy. They also go some way toward explaining the different shapes of the velocity ellipsoids of giant and dwarf stars in the solar neighbourhood. However, a consistent picture only emerges if one adopts small values of R_0 = 7.1 +/- 0.4 kpc and Theta_0 = 184 +/- 8 km/s. With these Galactic constants, the Milky Way's rotation curve declines slowly in the outer Galaxy; V_rot(20 kpc) = 166 kms. Our low value for R_0 agrees well with the only direct determination (7.2 +/- 0.7 kpc, Reid 1993). Using these Galactic constants, we find that the proper motion of Sgr A^* is consistent with the observational constraints. The radial velocities and proper motions of our best fit model are entirely consistent with the radial velocities of Cepheids and the Hipparcos measurements of their proper motions.

The local stellar kinematics of the Milky Way, parameterized by the Oort constants A and B, depend on the local gradient of the rotation curve, its absolute value (Theta_0), and the distance to the Galactic center (R_0). The surface density of interstellar gas in the Milky Way varies non-monotonically with radius, and so contributes significantly to the local gradient of the rotation curve, and the Oort constants. Because of this, the Oort functions A(R) and B(R) differ significantly from the dominant ~Theta_0/R dependence, in the Solar neighborhood and other locations in the Galaxy. These models may explain the ~40% difference between the values for 2 A R_0 derived from radial velocity data originating in the inner and outer Galaxy (Merrifield 1992). Incorporating these local non-linearities explains the significant differences between the Oort constants derived from nearby stars (d <= 1 kpc; Hanson 1987) and distant Cepheids (d=0.5-6 kpc; Feast & Whitelock 1997). However, a consistent picture only emerges if one adopts small values for the Galactic constants: R_0 = 7.1 +/- 0.4 kpc, and Theta_0 = 184 +/- 8 km/s. These values are consistent with most kinematical methods of determining R_0, including the proper motion of Sgr A* (Backer 1996), the direct determination of R_0 using water masers (7.2 +/- 0.7 kpc; Reid 1993), and constraints set by the shape of the Milky Way's dark halo (Olling & Merrifield 1997b).

We present the results of the deepest optically identified X-ray survey yet made. The X-ray survey was made with the ROSAT PSPC and reaches a flux limit of 1.6x10^-15 erg cm^-2 s^-1 (0.5--2.0 keV). Above a flux limit of 2x10^-15 erg cm^-2 s^-1 we define a complete sample of 70 sources of which 59 are identified. Some (5) other sources have tentative identifications and in a further 4 the X-ray error-boxes are blank to R=23 mag. At the brighter flux levels (>= 10^-14 erg cm^-2 s^-1) we confirm the results of previous less deep X-ray surveys with 84% of the sources begin QSOs. At fainter fluxes, however, the survey is dominated by a population of galaxies with narrow optical emission lines (NELGs). In addition, a number of groups and clusters of galaxies are found at intermediate fluxes. Most of these are poor systems of low X-ray luminosity and are generally found at redshifts of > 0.3. Their numbers are consistent with a zero evolutionary scenario, in contrast to the situation for high luminosity clusters at the same redshift. We discuss the significance of these results to the determination of the cosmic soft X-ray background (XRB) and show that at 2x10^-15 erg cm^-2 s^-1, we have resolved more than 50% of the background. We also briefly consider the probable importance of NELG objects to the residual background and look at some of the properties of these unusual objects.

If a galaxy resides in a cluster, then its passage through the pervasive intracluster medium will produce a detectable signature in the X-ray emission from the cluster. The simplest kinematic information that can be extracted from this signature is the galaxy's direction of motion on the plane of the sky. This paper explores the constraints on cluster dynamics that could be derived from such information. In particular, we show that it is possible to define a projected anisotropy parameter, B(R), which is directly analogous to the usual orbital anisotropy parameter. We describe an estimator for this quantity, Bhat(R), which can be derived in a robust and straightforward manner. Using a simple dynamical model, we demonstrate the ambiguity between the distribution of mass and the distribution of galaxy orbits when interpreting the traditional measures of cluster kinematics (the projected density of galaxies and their line-of-sight velocity dispersion). As an example, we show how two very different dynamical models can fit the kinematic properties of the Coma cluster. We demonstrate that the measurement of Bhat using a relatively small sample of wake directions (N_{wake} ~ 50) would provide an effective mechanism for lifting this degeneracy. Thus, by combining X-ray measurements of wake directions with number counts and line-of-sight velocities derived from optical data, it will prove possible to measure both the orbit distribution and the form of the gravitational potential in clusters of galaxies.

The distribution of shapes of galaxies' dark halos provides a basic test for models of galaxy formation. To-date, few dark halo shapes have been measured, and the results of different methods appear contradictory. Here, we add to the sample of measured shapes by calculating the flattening of the Milky Way's dark halo based on the manner in which the gas layer in the Galaxy flares with radius. We also test the validity of this technique -- which has already been applied to several other galaxies -- by comparing the inferred halo flattening to that obtained from a stellar-kinematic analysis, which can only be applied to the Milky Way. Both methods return consistent values for the shape of the Milky Way's halo, with a shortest-to-longest axis ratio for the dark matter distribution of q = 0.75 +/- 0.25. However, this consistency is only achieved if we adopt a value of R_0 = 7 +/- 1 kpc for the Sun's distance to the Galactic center. Although this value is smaller than the IAU-sanctioned R_0 = 8.5 kpc, it is quite consistent with current observations. Whatever value of R_0 is adopted, neither method returns halo parameters consistent with a disk-like mass distribution, for which q ~ 0.2. This finding rules out cold molecular gas and decaying massive neutrinos as dark matter candidates.

In this paper we report on the analysis of a deep ROSAT HRI image of the moderately-rich cluster, Abell 2634, by which we have been able to detect the X-ray emission from the galaxies in the cluster. The ICM of Abell 2634 is an order of magnitude denser than that of the Virgo cluster, and so this analysis allows us to explore the X-ray properties of individual galaxies in the richest environment yet explored. By stacking the X-ray images of the galaxies together, we have shown that their emission appears to be marginally resolved by the HRI. This extent is smaller than for galaxies in poorer environments, and is comparable to the size of the galaxies in optical light. These facts suggest that the detected X-ray emission originates from the stellar populations of the galaxies, rather than from extended hot interstellar media. Support for this hypothesis comes from placing the optical and X-ray luminosities of these galaxies in the L_B - L_X plane: the galaxies of Abell 2634 lie in the region of this plane where models indicate that all the X-ray emission can be explained by the usual population of X-ray binaries. It is therefore probable that ram pressure stripping has removed the hot gas component from these galaxies.

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.

In the deepest optically identified X-ray survey yet performed, we have identified 32 X-ray selected QSOs to a flux limit of 2x10^{-15} erg cm^{-2} s^{-1} (0.5-2 keV). The survey, performed with the ROSAT PSPC, has 89% spectroscopic completeness. The QSO log(N)-log(S) relation is found to have a break to a flat slope at faint fluxes. The surface density of QSOs at the survey limit is 230+/-40 per square degree, the largest so far of any QSO survey. We have used this survey to measure the QSO X-ray luminosity function at low luminosities (Lx<10^{44.5} erg s^{-1}) and high redshifts (1<z<2.5). The highest redshift QSO in the survey has z=3.4. Combined with the QSOs from the Einstein EMSS at bright fluxes, we find pure luminosity evolution of the form Lx proportional to (1+z)^{3.0(+0.2,-0.3)} is an adequate description of the evolution of the X-ray luminosity function at low redshifts. A redshift cutoff in the evolution is required at z=1.4 ^{+0.4}_{-0.17} (for qo=0.5). We discuss the form of this evolution, its dependence on the model assumed and the errors on the derived parameters. We show that most previous X-ray surveys, including the EMSS, are consistent with a power law luminosity evolution index of 3.0. The contribution of QSOs to the 1-2 keV cosmic X-ray background is found to be between 31% and 51%.

We have obtained high signal-to-noise spectra along the major axes of 28 S0 galaxies in order to search for the presence of disk stars on retrograde orbits. Full line-of-sight velocity distributions were extracted from the data, and the velocity distributions were modelled as arising from the superposition of populations of stars on prograde and retrograde orbits. We find no new cases in which a significant fraction of disk stars lie on retrograde orbits; an identical analysis of NGC~4550 does reveal the previously-known counter-rotating stellar disk in this system. Upper limits determined for each object indicate that no more than about 5% of the observed disk star light could arise from counter-rotating stellar components. These results suggest that previously-discovered disk galaxies with counter-rotating stars are exceptional and that (at 95% confidence) at most 10% of S0 galaxies contain significant counter-rotating populations. The most likely value for the fraction of such S0 galaxies lies closer to 1%. This result contrasts with the prevalence of counter-rotating gas in these systems; combining our new observations with existing data, we find that 24+/-8% (1-sigma error) of the gas disks in S0 galaxies counter-rotate relative to their stellar components.

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.

It has been suggested that the peanut-shaped bulges seen in some edge-on disk galaxies are produced when bars in these galaxies buckle. This paper reviews the modelling which seeks to show how bars buckle, and I present a very simple new model which captures the essential physics of this process. I then discuss the problems in establishing observationally the connection between peanut-shaped bulges and bars: confirmation of the link has proved difficult because boxy bulges are only apparent in edge-on galaxies whereas bars are only easily detectable in more face-on systems. Finally, I present a new technique which avoids this difficulty by searching for the distinctive kinematic signature of an edge-on bar; application of this method to spectra of peanut-shaped bulges reveals that they are, indeed, associated with hidden bars.

We have used the Tremaine-Weinberg method to measure the angular speed of rotation for the bar in the SB0 galaxy NGC 936. With this technique, the bar's pattern speed, Omega_p, can be derived from the luminosity and stellar-kinematic information in long-slit spectral observations taken parallel to the major axis of the galaxy. The kinematic measurement required is the mean line-of-sight velocity of all stellar light entering the slit. This quantity can only be calculated reliably if any asymmetry in the shape of the broadening function of the spectral lines is also measured, and so we present a method which allows for such asymmetry. The technique also returns a true measure of the RMS uncertainty in the estimate. Application of the analysis to a set of long-slit spectra of NGC 936 returns four separate measures of Omega_p which are mutually consistent. Combining these data produces a best estimate for the bar pattern speed of Omega_p = 60 +/- 14 km/s/kpc (assuming a distance of 16.6 Mpc). This result refines the only previous attempt to make this measurement, which yielded an estimate for Omega_p in NGC 936 of 104 +/- 37 km/s/kpc (Kent 1987). The new measurement places the co-rotation radius just beyond the end of the bar, in agreement with theoretical calculations.

It has been suggested that the peanut-shaped bulges seen in some edge-on disk galaxies are due to the presence of a central bar. Although bars cannot be detected photometrically in edge-on galaxies, we show that barred potentials produce a strong kinematic signature in the form of double-peaked line-of-sight velocity distributions with a characteristic ``figure-of-eight'' variation with radius. We have obtained spectroscopic observations of two edge-on galaxies with peanut-shaped bulges (NGC~5746 and NGC~5965), and they reveal exactly such line-of-sight velocity distributions in both their gaseous (emission line) and their stellar (absorption line) components. These observations provide strong observational evidence that peanut-shaped bulges are a by-product of bar formation.