The properties of universes are explored that are entirely in the interior of black holes. It is argued that these models offer a paradigm that may shed a new light on old cosmological problems. The topics that are addressed include: geometry of the universes, evolution, relation to the concordance model, dark matter, dark energy, mass ejection from black holes, galaxy models with a central black hole, Mach's principle.

We show results based on a large suite of N-Body/SPH simulations of isolated, flat dwarf galaxies, both rotating and non-rotating. The main goal is to investigate possible mechanisms to explain the observed dichotomy in radial stellar metallicity profiles of dwarf galaxies: dwarf irregulars (dIrr) and flat, rotating dwarf ellipticals (dE) generally possess flat metallicity profiles, while rounder and non-rotating dEs show strong negative metallicity gradients. These simulations show that flattening by rotation is key to reproducing the observed characteristics of flat dwarf galaxies, proving particularly efficient in erasing metallicity gradients. We propose a "centrifugal barrier mechanism" as an alternative to the previously suggested "fountain mechanism" for explaining the flat metallicity profiles of dIrrs and flat, rotating dEs. While only flattening the dark-matter halo has little influence, the addition of angular momentum slows down the infall of gas, so that star formation (SF) and the ensuing feedback are less centrally concentrated, occurring galaxy-wide. Additionally, this leads to more continuous SFHs by preventing large-scale oscillations in the SFR ("breathing"), and creates low density holes in the ISM, in agreement with observations of dIrrs. Our general conclusion is that rotation has a significant influence on the evolution and appearance of dwarf galaxies, and we suggest angular momentum as a second parameter (after galaxy mass as the dominant parameter) in dwarf galaxy evolution. Angular momentum differentiates between SF modes, making our fast rotating models qualitatively resemble dIrrs, which does not seem possible without rotation.

Recently, some intriguing results have lead to speculations whether the central density slope -- velocity dispersion anisotropy inequality (An & Evans) actually holds at all radii for spherical dynamical systems. We extend these studies by providing a complete analysis of the global slope -- anisotropy inequality for all spherical systems in which the augmented density is a separable function of radius and potential. We prove that these systems indeed satisfy the global inequality if their central anisotropy is $\beta_0\leq 1/2$. Furthermore, we present several systems with $\beta_0 > 1/2$ for which the inequality does not hold, thus demonstrating that the global density slope -- anisotropy inequality is not a universal property. This analysis is a significant step towards an understanding of the relation for general spherical systems.

The inner mass-profile of the relaxed cluster Abell 1703 is analysed by two very different strong-lensing techniques applied to deep ACS and WFC3 imaging. Our parametric method has the accuracy required to reproduce the many sets of multiple images, based on the assumption that mass approximately traces light. We test this assumption with a fully non-parametric, adaptive grid method, with no knowledge of the galaxy distribution. Differences between the methods are seen on fine scales due to member galaxies which must be included in models designed to search for lensed images, but on the larger scale the general distribution of dark matter is in good agreement, with very similar radial mass profiles. We add undiluted weak-lensing measurements from deep multi-colour Subaru imaging to obtain a fully model-independent mass profile out to the virial radius and beyond. Consistency is found in the region of overlap between the weak and strong lensing, and the full mass profile is well-described by an NFW model of a concentration parameter, $c_{\rm vir}\simeq 7.15\pm0.5$ (and $M_{vir}\simeq 1.22\pm0.15 \times 10^{15}M_{\odot}/h$). Abell 1703 lies above the standard $c$--$M$ relation predicted for the standard $\Lambda$CDM model, similar to other massive relaxed clusters with accurately determined lensing-based profiles.

In this paper we investigate whether Smoothed Particle Hydrodynamics (SPH), equipped with artificial conductivity, is able to capture the physics of density/energy discontinuities in the case of the so-called shearing layers test, a test for examining Kelvin-Helmholtz (KH) instabilities. We can trace back each failure of SPH to show KH rolls to two causes: i) shock waves travelling in the simulation box and ii) particle clumping, or more generally, particle noise. The probable cause of shock waves is the Local Mixing Instability (LMI), previously identified in the literature. Particle noise on the other hand is a problem because it introduces a large error in the SPH momentum equation. We also investigate the role of artificial conductivity (AC). Including AC is necessary for the long-term behavior of the simulation (e.g. to get $\lambda=1/2, 1$ KH rolls). In sensitive hydrodynamical simulations great care is however needed in selecting the AC signal velocity, with the default formulation leading to too much energy diffusion. We present new signal velocities that lead to less diffusion. The effects of the shock waves and of particle disorder become less important as the time-scale of the physical problem (for the shearing layers problem: lower density contrast and higher Mach numbers) decreases. At the resolution of current galaxy formation simulations mixing is probably not important. However, mixing could become crucial for next-generation simulations.

In this article we study the well-known strong lensing system SDSS J1004+4112. Not only does it host a large-separation lensed quasar with measured time-delay information, but several other lensed galaxies have been identified as well. A previously developed strong lens inversion procedure that is designed to handle a wide variety of constraints, is applied to this lensing system and compared to results reported in other works. Without the inclusion of a tentative central image of one of the galaxies as a constraint, we find that the model recovered by the other constraints indeed predicts an image at that location. An inversion which includes the central image provides tighter constraints on the shape of the central part of the mass map. The resulting model also predicts a central image of a second galaxy where indeed an object is visible in the available ACS images. We find masses of 2.5x10^13 M_O and 6.1x10^13 M_O within a radius of 60 kpc and 110 kpc respectively, confirming the results from other authors. The resulting mass map is compatible with an elliptical generalization of a projected NFW profile, with r_s = 58_{-13}^{+21} arcsec and c_vir = 3.91 +/- 0.74. The orientation of the elliptical NFW profile follows closely the orientation of the central cluster galaxy and the overall distribution of cluster members.

$N$-body simulations have unveiled several apparently universal properties of dark matter halos, including a cusped density profile, a power-law pseudo phase-space density $\rho/\sigma_r^3$, and a linear $\beta-\gamma$ relation between the density slope and the velocity anisotropy. We present a family of self-consistent phase-space distribution functions $F(E,L)$, based on the Dehnen-McLaughlin Jeans models, that incorporate these universal properties very accurately. These distribution functions, derived using a quadratic programming technique, are analytical, positive and smooth over the entire phase space and are able to generate four-parameter velocity anisotropy profiles $\beta(r)$ with arbitrary asymptotic values $\beta_0$ and $\beta_\infty$. We discuss the orbital structure of six radially anisotropic systems in detail and argue that, apart from its use for generating initial conditions for $N$-body studies, our dynamical modeling provides a valuable complementary approach to understand the processes involved in the formation of dark matter halos.

We present high spatial resolution 21cm HI observations of EA01A and EA01B, a pair of interacting post-starburst, or E+A, galaxies at z = 0.0746. Based on optical HST/WFPC2 images, both galaxies are known to display disturbed morphologies. They also appear to be linked by a bridge of stars. Previous HI observations Chang et al. (2001) had already uncovered sizable quantities of neutral gas in or near these galaxies but they lacked the spatial resolution to locate the gas with any precision within this galactic binary system. We have analysed deep, high resolution archival VLA observations of the couple. We find evidence for three gaseous tidal tails; one connected to EA01A and two emanating from EA01B. These findings confirm, independently from the optical imaging, that (i) EA01A and EA01B are actively interacting, and that, as a consequence, the starbursts that occurred in these galaxies were most likely triggered by this interaction, and that (ii) 6.6+-0.9 10^9 Msun of neutral gas are still present in the immediate vicinity of the optical bodies of both galaxies. The HI column density is lowest at the optical positions of the galaxies, suggesting that most of the neutral gas that is visible in our maps is associated with the tidal arms and not with the galaxies themselves. This might provide an explanation for the apparent lack of ongoing star formation in these galaxies.

We present new fully self-consistent models of the formation and evolution of isolated dwarf galaxies. We have used the publicly available N-body/SPH code HYDRA, to which we have added a set of star formation criteria, and prescriptions for chemical enrichment (taking into account contributions from both SNIa and SNII), supernova feedback, and gas cooling. The models follow the evolution of an initially homogeneous gas cloud collapsing in a pre-existing dark-matter halo. These simplified initial conditions are supported by the merger trees of isolated dwarf galaxies extracted from the milli-Millennium Simulation. The star-formation histories of the model galaxies exhibit burst-like behaviour. These bursts are a consequence of the blow-out and subsequent in-fall of gas. The amount of gas that leaves the galaxy for good is found to be small, in absolute numbers, ranging between 3x10^7 Msol and 6x10^7 Msol . For the least massive models, however, this is over 80 per cent of their initial gas mass. The local fluctuations in gas density are strong enough to trigger star-bursts in the massive models, or to inhibit anything more than small residual star formation for the less massive models. Between these star-bursts there can be time intervals of several Gyrs. We have compared model predictions with available data for the relations between luminosity and surface brightness profile, half-light radius, central velocity dispersion, broad band colour (B-V) and metallicity, as well as the location relative to the fundamental plane. The properties of the model dwarf galaxies agree quite well with those of observed dwarf galaxies.

The cluster lens Cl 0024+1654 is undoubtedly one of the most beautiful examples of strong gravitational lensing, providing five large images of a single source with well-resolved substructure. Using the information contained in the positions and the shapes of the images, combined with the null space information, a non-parametric technique is used to infer the strong lensing mass map of the central region of this cluster. This yields a strong lensing mass of 1.60x10^14 M_O within a 0.5' radius around the cluster center. This mass distribution is then used as a case study of the monopole degeneracy, which may be one of the most important degeneracies in gravitational lensing studies and which is extremely hard to break. We illustrate the monopole degeneracy by adding circularly symmetric density distributions with zero total mass to the original mass map of Cl 0024+1654. These redistribute mass in certain areas of the mass map without affecting the observed images in any way. We show that the monopole degeneracy and the mass-sheet degeneracy together lie at the heart of the discrepancies between different gravitational lens reconstructions that can be found in the literature for a given object, and that many images/sources, with an overall high image density in the lens plane, are required to construct an accurate, high-resolution mass map based on strong-lensing data.

Lopsidedness is common in spiral galaxies. Often, there is no obvious external cause, such as an interaction with a nearby galaxy, for such features. Alternatively, the lopsidedness may have an internal cause, such as a dynamical instability. In order to explore this idea, we have developed a computer code that searches for self-consistent perturbations in razor-thin disc galaxies and performed a thorough mode-analysis of a suite of dynamical models for disc galaxies embedded in an inert dark-matter halo with varying amounts of rotation and radial anisotropy. Models with two equal-mass counter-rotating discs and fully rotating models both show growing lopsided modes. For the counter-rotating models, this is the well-known counter-rotating instability, becoming weaker as the net rotation increases. The m=1 mode of the maximally rotating models, on the other hand, becomes stronger with increasing net rotation. This rotating m=1 mode is reminiscent of the eccentricity instability in near-Keplerian discs. To unravel the physical origin of these two different m=1 instabilities, we studied the individual stellar orbits in the perturbed potential and found that the presence of the perturbation gives rise to a very rich orbital behaviour. In the linear regime, both instabilities are supported by aligned loop orbits. In the non-linear regime, other orbit families exist that can help support the modes. In terms of density waves, the counter-rotating m=1 mode is due to a purely growing Jeans-type instability. The rotating m=1 mode, on the other hand, grows as a result of the swing amplifier working inside the resonance cavity that extends from the disc center out to the radius where non-rotating waves are stabilized by the model's outwardly rising Q-profile.

The inversion of a gravitational lens system is, as is well known, plagued by the so-called mass-sheet degeneracy: one can always rescale the density distribution of the lens and add a constant-density mass-sheet such that the, also properly rescaled, source plane is projected onto the same observed images. For strong lensing systems, it is often claimed that this degeneracy is broken as soon as two or more sources at different redshifts are available. This is definitely true in the strict sense that it is then impossible to add a constant-density mass-sheet to the rescaled density of the lens without affecting the resulting images. However, often one can easily construct a more general mass distribution -- instead of a constant-density sheet of mass -- which gives rise to the same effect: a uniform scaling of the sources involved without affecting the observed images. We show that this can be achieved by adding one or more circularly symmetric mass distributions, each with its own center of symmetry, to the rescaled mass distribution of the original lens. As it uses circularly symmetric distributions, this procedure can lead to the introduction of ring shaped features in the mass distribution of the lens. In this paper, we show explicitly how degenerate inversions for a given strong lensing system can be constructed. It then becomes clear that many constraints are needed to effectively break this degeneracy.

We present new estimates of ages and metallicities, based on FORS/VLT optical (4400-5500A) spectroscopy, of 16 dwarf elliptical galaxies (dE's) in the Fornax Cluster and in Southern Groups. These dE's are more metal-rich and younger than previous estimates based on narrow-band photometry and low-resolution spectroscopy. For our sample we find a mean metallicity [Z/H] = -0.33 dex and mean age 3.5 Gyr, consistent with similar samples of dE's in other environments (Local Group, Virgo). Three dE's in our sample show emission lines and very young ages. This suggests that some dE's formed stars until a very recent epoch and were self-enriched by a long star formation history. Previous observations of large near-infrared (~8500A) Ca II absorption strengths in these dE's are in good agreement with the new metallicity estimates, solving part of the so-called Calcium puzzle.

Galaxies acting as gravitational lenses are surrounded by, at most, a handful of images. This apparent paucity of information forces one to make the best possible use of what information is available to invert the lens system. In this paper, we explore the use of a genetic algorithm to invert in a non-parametric way strong lensing systems containing only a small number of images. Perhaps the most important conclusion of this paper is that it is possible to infer the mass distribution of such gravitational lens systems using a non-parametric technique. We show that including information about the null space (i.e. the region where no images are found) is prerequisite to avoid the prediction of a large number of spurious images, and to reliably reconstruct the lens mass density. While the total mass of the lens is usually constrained within a few percent, the fidelity of the reconstruction of the lens mass distribution depends on the number and position of the images. The technique employed to include null space information can be extended in a straightforward way to add additional constraints, such as weak lensing data or time delay information.

We study the locus of dwarf and giant early and late-type galaxies on the Tully-Fisher relation (TFR), the stellar mass Tully-Fisher relation (sTFR) and the so-called baryonic or HI gas+stellar mass Tully-Fisher relation (gsTFR). We show that early-type and late-type galaxies, from dwarfs to giants, trace different yet approximately parallel TFRs. Surprisingly, early-type and late-type galaxies trace a single yet curved sTFR over a range of 3.5 orders of magnitude in stellar mass. Moreover, all galaxies trace a single, linear gsTFR, over 3.5 orders of magnitude in HI gas+stellar mass. Dwarf ellipticals, however, lie slightly below the gsTFR. This may indicate that early-type dwarfs, contrary to the late-types, have lost their gas, e.g. by galactic winds or ram-pressure stripping. Overall, environment only plays a secondary role in shaping these relations, making them a rather ``clean'' cosmological tool. LCDM simulations predict roughly the correct slopes for these relations.

We present a method to investigate the radial stability of a spherical anisotropic system that hosts a central supermassive black hole (SBH). Such systems have never been tested before for stability, although high anisotropies have been considered in the dynamical models that were used to estimate the masses of the central putative supermassive black holes. A family of analytical anisotropic spherical Hernquist models with and without a black hole were investigated by means of N-body simulations. A clear trend emerges that the supermassive black hole has a significant effect on the overall stability of the system, i.e. an SBH with a mass of a few percent of the total mass of the galaxy can prevent or reduce the bar instabilities in anisotropic systems. Its mass not only determines the strength of the instability reduction, but also the time in which this occurs. These effects are most significant for models with strong radial anisotropies. Furthermore, our analysis shows that unstable systems with similar SBH but with different anisotropy radii evolve differently: highly radial systems become oblate, while more isotropic models tend to form into prolate structures. In addition to this study, we also present a Monte-Carlo algorithm to generate particles in spherical anisotropic systems.

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

We present a search for CO(1->0) emission in three Local Group dwarf irregular galaxies: IC5152, the Phoenix dwarf, and UGCA438, using the ATNF Mopra radio telescope. Our scans largely cover the optical extent of the galaxies and the stripped HI cloud West of the Phoenix dwarf. Apart from a tentative but non-significant emission peak at one position in the Phoenix dwarf, no significant emission was detected in the CO spectra of these galaxies. For a velocity width of 6 km/s, we derive 4sigma upper limits of 0.03 K km/s, 0.04 K km/s and 0.06 K km/s for IC5152, the Phoenix dwarf and UGCA438, respectively. This is an improvement of over a factor of 10 compared with previous observations of IC5152; the other two galaxies had not yet been observed at millimeter wavelengths. Assuming a Galactic CO-to-H_2 conversion factor, we derive upper limits on the molecular gas mass of 6.2 x 10^4 M_sun, 3.7 x 10^3 M_sun and 1.4 x 10^5 M_sun for IC5152, the Phoenix dwarf and UGCA438, respectively. We investigate two possible causes for the lack of CO emission in these galaxies. On the one hand, there may be a genuine lack of molecular gas in these systems, in spite of the presence of large amounts of neutral gas. However, in the case of IC5152 which is actively forming stars, molecular gas is at least expected to be present in the star forming regions. On the other hand, there may be a large increase in the CO-to-H_2 conversion factor in very low-metallicity dwarfs (-2 <= [Fe/H] <= -1), making CO a poor tracer of the molecular gas content in dwarf galaxies.

We present an updated investigation of the relation between large scale disk circular velocity, v_c, and bulge velocity dispersion, sigma_c. New bulge velocity dispersions are measured for a sample of 11 low surface brightness (LSB) and 7 high surface brightness (HSB) spiral galaxies for which v_c is known from published optical or HI rotation curves. We find that, while LSB galaxies appear to define the upper envelope of the region occupied by HSB galaxies (having relatively larger v_c for any given sigma_c), the distinction between LSB and HSB galaxies in the v_c-sigma_c plane becomes less pronounced for sigma_c <= 80 km/s. We conclude that either the scatter of the v_c-sigma_c relation is a function of v_c (and hence galaxy mass) or that the character of the v_c-sigma_c relation changes at v_c ~ 80 km/s. Some inplications of our findings are discussed.

We observed the CO(3-2) emission of the emission-line regions HubbleI, HubbleV, HubbleX, Holmberg 18, and the stellar emission-line object S28 in NGC6822 with the ESO Atacama Pathfinder Experiment (APEX) 12m telescope as part of its science verification. The very low system temperature of 130-180K enabled us to achieve detections in 4 single pointings and in a high spatial resolution 70''x70'' map of HubbleV. We compare the spectra with HI observations, obtained with the Australia Telescope Compact Array, of the same regions. In combination with previous multi-line CO observations, we perform a preliminary investigation of the physical conditions in HubbleV using a simple LTE model. We estimate the mass of the HubbleV region and the H_2/I_CO(3-2) conversion factor. Also, we show that HubbleV is located very near the line-width versus size relation traced by the Milky Way and LMC molecular clouds.

We present deep single-dish HI observations of a sample of six nearby E+A galaxies (0.05<z<0.1). A non-negligible fraction of a local sample of E+As are detected in HI. In four galaxies, we have detected up to a few times 10^9 Msun of neutral gas. These E+A galaxies are almost as gas-rich as spiral galaxies with comparable luminosities. There appears to exist no direct correlation between the amount of HI present in an E+A galaxy and its star-formation rate as traced by radio continuum emission. Moreover, the end of the starburst does not necessarily require the complete exhaustion of the neutral gas reservoir. Most likely, an intense burst of star formation consumed the dense molecular clouds, which are the sites of massive star formation. This effectively stops star formation, even though copious amounts of diffuse neutral gas remain. The remaining HI reservoir may eventually lead to further episodes of star formation. This may indicate that some E+As are observed in the inactive phase of the star-formation duty cycle.

We present three-integral dynamical models for the three Local Group dwarf elliptical galaxies: NGC147, NGC185, and NGC205. It is the first time models are constructed for the Local Group dEs that allow for the presence of dark matter at large radii and that are constrained by kinematics out to at least one half-light radius. The B-band mass-to-light ratios of all three galaxies are rather similar, (M/L)_B = 3-4. Within the inner two half-light radii, about 40-50% of the mass is in the form of dark matter, so dEs contain about as much dark matter as bright ellipticals. NGC205 and NGC147 have a relatively isotropic velocity dispersion tensor within the region where the internal dynamics are strongly constrained by the data. Our estimated inclination for NGC185 is i ~ 50 degrees because in that case the model has an intrinsic flattening close to the peak of the intrinsic shape distribution of dEs and it, like the best fitting models for NGC147 and NGC205, is nearly isotropic. We also show that the dynamical properties of the bright nucleus of NGC205 are not unlike those of a massive globular cluster.

Aims: We aim at showing that the broad-band SED characteristics of our sample of post-AGB stars are best interpreted, assuming the circumstellar dust is stored in Keplerian rotating passive discs. Methods: We present a homogeneous and systematic study of the Spectral Energy Distributions (SEDs) of a sample of 51 post-AGB objects. The selection criteria to define the whole sample were tuned to cover the broad-band characteristics of known binary post-AGB stars. The whole sample includes 20 dusty RV Tauri stars from the General Catalogue of Variable Stars (GCVS). We supplemented our own Geneva optical photometry with literature data to cover a broad range of fluxes from the UV to the far-IR. Results: All the SEDs display very similar characteristics: a large IR excess with a dust excess starting near the sublimation temperature, irrespective of the effective temperature of the central star. Moreover, when available, the long wavelength fluxes show a black-body slope indicative of the presence of a component of large mm sized grains. Conclusions: We argue that in all systems, gravitationally bound dusty discs are present. The discs must be puffed-up to cover a large opening angle for the central star and we argue that the discs have some similarity with the passive discs detected around young stellar objects. We interpret the presence of a disc to be a signature for binarity of the central object, but this will need confirmation by long-term monitoring of the radial velocities. We argue that dusty RV Tauri stars are those binaries which happen to be in the Population II instability strip.

We present upper-limits on the masses of the putative central intermediate-mass black holes in two nearby Galactic globular clusters: 47Tuc (NGC104), the second brightest Galactic globular cluster, and NGC6397, a core-collapse globular cluster and, with a distance of 2.7 kpc, quite possibly the nearest globular cluster, using a technique suggested by T. Maccarone. These mass estimates have been derived from 3sigma upper limits on the radio continuum flux at 1.4 GHz, assuming that the putative central black hole accretes the surrounding matter at a rate between 0.1% and 1% of the Bondi accretion rate. For 47Tuc, we find a 3sigma upper limit of 2060 - 670 solar masses, depending on the actual accretion rate of the black hole and the distance to 47Tuc. For NGC6397, which is closer to us, we derive a 3sigma upper limit of 1290 - 390 solar masses. While estimating mass upper-limits based on radio continuum observations requires making assumptions about the gas density and the accretion rate of the black hole, their derivation does not require complex and time consuming dynamical modeling. Thus, this method offers an independent way of estimating black hole masses in nearby globular clusters. If, generally, central black holes in stellar systems accrete matter faster than 0.1% of the Bondi accretion rate, then these results indicate the absence of black holes in these globular clusters with masses as predicted by the extrapolated M-sigma relation.

We present a non-parametric technique to infer the projected-mass distribution of a gravitational lens system with multiple strong-lensed images. The technique involves a dynamic grid in the lens plane on which the mass distribution of the lens is approximated by a sum of basis functions, one per grid cell. We used the projected mass densities of Plummer spheres as basis functions. A genetic algorithm then determines the mass distribution of the lens by forcing images of a single source, projected back onto the source plane, to coincide as well as possible. Averaging several tens of solutions removes the random fluctuations that are introduced by the reproduction process of genomes in the genetic algorithm and highlights those features common to all solutions. Given the positions of the images and the redshifts of the sources and the lens, we show that the mass of a gravitational lens can be retrieved with an accuracy of a few percent and that, if the sources sufficiently cover the caustics, the mass distribution of the gravitational lens can also be reliably retrieved. A major advantage of the algorithm is that it makes full use of the information contained in the radial images, unlike methods that minimise the residuals of the lens equation, and is thus able to accurately reconstruct also the inner parts of the lens.

We present HI 21cm line observations, obtained with the Australia Telescope Compact Array, of two dwarf elliptical galaxies (dEs) in the Fornax cluster: FCC032 and FCC336. The optical positions and velocities of these galaxies place them well within the Fornax cluster. FCC032 was detected at the 3sigma significance level with a total HI flux density of 0.66+-0.22 Jy km s^-1 or an HI mass of 5.0+-1.7 10^7 h_75^-2 M_sun. Based on our deep Halpha+[NII] narrow-band images, obtained with FORS2 mounted on the VLT, this dE was already known to contain 600-1800 h_75^-2 M_sun of ionised Hydrogen (depending on the relative strengths of the Halpha and [NII] emission lines). Hence, this is the first study of the complex, multi-phase interstellar medium of a dE outside the Local Group. FCC336 was detected at the same significance level: 0.37+-0.10 Jy km s^-1 or a total HI mass of 2.8+-0.7 10^7 h_75^-2,M_sun. Using a compilation of HI data of dwarf galaxies, we find that the observed high HI-mass boundary of the distribution of dIrrs, BCDs, and dEs in a log(L_B) versus log (M_HI) diagram is in good agreement with a simple chemical evolution model with continuous star formation. The existence of many gas-poor dEs (undetected at 21cm) suggest that the environment (or more particularly, a galaxy's orbit within a cluster) also plays a crucial role in determining the amount of gas in present-day dEs. E.g., FCC032 and FCC336 are located in the sparsely populated outskirts of the Fornax cluster. This is in agreement with HI surveys of dEs in the Virgo Cluster and an Halpha survey of the Fornax Cluster, which also tend to place gas-rich dwarf galaxies in the cluster periphery.

We present SKIRT, a three-dimensional Monte Carlo radiative transfer code developed to study dusty galaxies. We discuss SKIRT's most important characteristics and present a number of applications. In particular, we focus on the kinematical aspect of SKIRT. We demonstrate that dust attenuation mimics the presence of dark matter around elliptical galaxies and that it severely affects the rotation curves of edge-on galaxies.

We present extended Spectral Energy Distributions (SEDs) of seven classical RV Tauri stars, using newly obtained submillimetre continuum measurements and Geneva optical photometry supplemented with literature data. The broad-band SEDs show a large IR excess with a black-body slope at long wavelengths in six of the seven stars, R Sct being the noticeable exception. This long wavelength slope is best explained assuming the presence of a dust component of large grains in the circumstellar material. We show that the most likely distribution of the circumstellar dust around the six systems is that the dust resides in a disc. Moreover, very small outflow velocities are needed to explain the presence of dust near the sublimation temperature and we speculate that the discs are Keplerian. The structure and evolution of these compact discs are as yet not understood but a likely prerequisite for their formation is that the dusty RV Tauri stars are binaries.

This paper is the first in a series in which we present the results of an ESO Large Program on the kinematics and internal dynamics of dwarf elliptical galaxies (dEs). We investigate the relations between the parameters that quantify the structure and internal dynamics of dEs such as the Faber_Jackson relation and the Fundamental Plane (FP). We show that the dE sequences in the various diagrams are disjunct from those traced by bright and intermediate-luminosity elliptical galaxies and bulges of spirals. It appears that semi-analytical models (SAMs) are able to reproduce the position of the dEs in those diagrams. While these findings are clearly a success for the hierarchical-merging picture of galaxy formation, they do not necessarily invalidate the alternative ``harassment'' scenario, which posits that dEs stem from perturbed and stripped late-type disk galaxies that entered clusters and groups of galaxies about 5 Gyr ago.

It is generally agreed that the rotational velocities of spiral galaxies indicate that a lot of matter is unseen. Hitherto, the majority of these mass models are all based on spherical symmetry, although this assumption is obviously not the one that is closest to reality. We present a new method to derive dynamical mass models from HI-data based on an axisymmetric geometry. To validate our method, we started a project with the aim to create such models from a combined sample of late-type Low Surface Brightness galaxies and early-type High Surface Brightness spiral galaxies. The whole sample spans a wide range of bulge-to-disk ratios and will ultimately allow us to investigate the orbital distribution as a function of morphology.

We investigate the relation between the dark matter distribution of galaxies and their central supermassive black holes which is suggested by the v_c-sigma relation. Since early-type galaxies appear to have larger black holes than late-type ones, we look for an equivalent pattern in the dark matter distribution as a function of Hubble type. To achieve our goal we use a state-of-the-art modelling code that allows a variety of geometries to be fitted to a combination of radio and optical observations of galaxies with different morphology.

Using an N-body/SPH code, we explore the scenario in which a dwarf elliptical galaxy (dE) is subjected to ram-pressure stripping due to the intracluster medium (ICM). Our simulations show that while (i) smaller dEs lose their ISM almost immediately after entering the cluster, (ii) more massive dEs are able to retain their gas for considerable timespans.

The true nature of the Hyades and Sirius superclusters is still an open question. In this contribution, we confront Eggen's hypothesis that they are cluster remnants with the results of a kinematic analysis of more than 6000 K and M giants in the solar neighbourhood. This analysis includes new radial velocity data from a large survey performed with the Coravel spectrometer, complemented by Hipparcos parallaxes and Tycho-2 proper motions (Famaey et al. 2004). A maximum-likelihood method, based on a bayesian approach, has been applied to the data, in order to make full use of all the available data (including less precise parallaxes) and to derive the properties of the different kinematic subgroups. Two such subgroups can be identified with the Hyades and Sirius superclusters. Stars belonging to them span a very wide range of age, which is difficult to account for in Eggen's scenario. These groups are thus most probably "dynamical streams" related to the dynamical perturbation by spiral waves rather than to cluster remnants. In this scenario, the Hyades and Ursa Major clusters just happen to be in the Hyades and Sirius streams, which are purely dynamical features that have nothing to do with the remnants of more massive primordial clusters. This mechanism could be the key to understanding the presence of an old metal-rich population, and of many exoplanetary systems in our neighbourhood. Moreover, a strong spiral pattern seems to be needed in order to yield such prominent streams. Since spiral structure is usually baryonic, this would leave very little room for dark matter. This may be an indication that the era of the dark-matter paradigm explaining the dynamics of the Galaxy may come to an end, and is being superseded by modified gravity.

We discuss the kinematical structure of a two-parameter family of isotropic models with a central black hole. The family contains the slope of the central density cusp and the relative black hole mass as parameters. Most of the basic kinematical quantities of these models can be expressed analytically. This family contains three distinct models where also the distribution function, differential energy distribution and spatial LOSVDs can be expressed completely analytically. Each of these models show a drastically different behaviour of the distribution function. Although the effect of a black hole on the distribution function is very strong, in particular for models with a shallow density cusp, the differential energy distribution is only marginally affected. We discuss in detail the effects of a central black hole on the LOSVDs. The projected velocity dispersion increases with black hole mass at small projected radii, but the effect of a black hole on the shape of the LOSVDs (characterized by the h4 parameter) is less straightforward to interpret. Too much reliance on the wings of the LOSVDs and the value of the h4 parameter to determine black hole masses might hence be dangerous.

The availability of the Hipparcos Catalogue triggered many kinematic and dynamical studies of the solar neighbourhood. Nevertheless, those studies generally lacked the third component of the space velocities, i.e., the radial velocities. This work presents the kinematic analysis of 5952 K and 739 M giants in the solar neighbourhood which includes for the first time radial velocity data from a large survey performed with the CORAVEL spectrovelocimeter. It also uses proper motions from the Tycho-2 catalogue, which are expected to be more accurate than the Hipparcos ones. The UV-plane constructed from these data for the stars with precise parallaxes reveals a rich small-scale structure, with several clumps corresponding to the Hercules stream, the Sirius moving group, and the Hyades and Pleiades superclusters. A maximum-likelihood method, based on a bayesian approach, has been applied to the data, in order to make full use of all the available stars and to derive the kinematic properties of these subgroups. Isochrones in the Hertzsprung-Russell diagram reveal a very wide range of ages for stars belonging to these groups. These groups are most probably related to the dynamical perturbation by transient spiral waves rather than to cluster remnants. A possible explanation for the presence of young clusters in the same area of the UV-plane is that they have been put there by the spiral wave associated with their formation, while the kinematics of the older stars of our sample has also been disturbed by the same wave. The term "dynamical stream" for the kinematic groups is thus more appropriate than the traditional term "supercluster" since it involves stars of different ages, not born at the same place nor at the same time.

We identified 21 new Planetary Nebula (PN) candidates in the Sculptor Group galaxy NGC 55. We determined a most likely distance of 2.00 +/- 0.2 Mpc using the Planetary Nebulae Luminosity Function (PNLF) method. The distance to NGC 55 is larger than previously determined distances, which means that the Sculptor Group is a bit further away from the Local Group than previously thought. The distance to NGC 55 is again similar to the distance of NGC 300, adding support to the suggestion that these galaxies form a bound pair.

We present R-band and H-alpha+[N II] narrow-band imaging of FCC032, FCC206 and FCCB729, three dwarf elliptical galaxies (dEs) in the Fornax Cluster. These dEs contain significant amounts of ionized gas. FCC032 harbours a large ionized gas complex, consisting of several individual clouds, a superbubble and a filament that extends away from the galaxy centre. The ionized gas structures observed in FCC032 bear a strong resemblance to those observed in more gas-rich and more fiercely star-forming dwarf galaxies. FCC206, a very low surface brightness dE, contains one faint extended emission region, and two compact clouds. In FCCB729, the only nucleated galaxy in this sample, one of the ionized gas clouds coincides with the stellar nucleus. We derive ionized gas masses of a few 100 to 1000 solar masses for these galaxies. This brings our sample of dEs with ionized gas with H-alpha+[N II] narrow-band imaging to five objects. The ionized gas morphologies in these galaxies range from pure nuclear emission peaks to extended emission complexes. This morphological diversity could also indicate a diversity in ionizing processes in dEs with a warm interstellar medium (active galactic nuclei, starbursts, post-AGB stars,...). Inside FCC206, four star clusters appear to be on the verge of merging to form a nucleus in this presently non-nucleated dE. Understanding the formation of nuclei in dEs could give us important clues to the formation of super-massive black holes.

{We present, for the first time, photometric and kinematical evidence, obtained with FORS2 on the VLT, for the existence of kinematically decoupled cores (KDCs) in two dwarf elliptical galaxies; FS76 in the NGC5044 group and FS373 in the NGC3258 group. Both kinematically peculiar subcomponents rotate in the same sense as the main body of their host galaxy but betray their presence by a pronounced bump in the rotation velocity profiles at a radius of about 1". The KDC in FS76 rotates at 10+/-3km/s, with the host galaxy rotating at 15+/-6km/s; the KDC in FS373 has a rotation velocity of 6+/-2km/s while the galaxy itself rotates at 20+/-5km/s. FS373 has a very complex rotation velocity profile with the velocity changing sign at 1.5 R_e. The velocity and velocity dispersion profiles of FS76 are asymmetric at larger radii. This could be caused by a past gravitational interaction with the giant elliptical NGC5044, which is at a projected distance of 50kpc. We argue that these decoupled cores are most likely not produced by mergers in a group or cluster environment because of the prohibitively large relative velocities. A plausible alternative is offered by flyby interactions between a dwarf elliptical or its disky progenitor and a massive galaxy. The tidal forces during an interaction at the relative velocities and impact parameters typical for a group environment exert a torque on the dwarf galaxy that, according to analytical estimates, transfers enough angular momentum to its stellar envelope to explain the observed peculiar kinematics.

We investigate the relation between circular velocity vc and bulge velocity dispersion sigma in spiral galaxies, based on literature data and new spectroscopic observations. We find a strong, nearly linear vc-sigma correlation with a negligible intrinsic scatter, and a striking agreement with the corresponding relation for elliptical galaxies. The least massive galaxies (sigma < 80 km/s) significantly deviate from this relation. We combine this vc-sigma correlation with the well-known MBH-sigma relation to obtain a tight correlation between circular velocity and supermassive black hole mass, and interpret this as observational evidence for a close link between supermassive black holes and the dark matter haloes in which they presumably formed. Apart from being an important ingredient for theoretical models of galaxy formation and evolution, the relation between MBH and circular velocity has the potential to become an important practical tool in estimating supermassive black hole masses in spiral galaxies.

We present a family of spherical models for elliptical galaxies and bulges consisting of a stellar component and a central black hole. All models in this family share the same stellar density profile, which has a steep central cusp. The gravitational potential of each models is a linear combination of the potential generated selfconsistently by the stars and the potential of a central black hole. The relative importance of these two contributions is a free parameter in the models. Assuming an isotropic dynamical structure, almost all kinematical properties of these models can be calculated analytically. In particular, they form the first simple dynamical models for galaxies with a central black hole where the distribution function and differential energy distribution can be written completely in terms of elementary functions only. We also present various extensions of this family to models with anisotropic orbital structures. Also for these dynamical models, the distribution function and its moments can be expressed completely in terms of elementary functions. This family is useful for a large range of applications, in particular to generate initial conditions for N-body and hydrodynamical simulations to model galactic nuclei with a central black hole.

In this paper we present a method to analyse absorption line spectra of a galaxy designed to determine the stellar dynamics and the stellar populations by a direct fit to the spectra. This paper is the first one to report on the application of the method to data. The modelling results in the knowledge of distribution functions that are sums of basis functions. The practical implementation of the method is discussed and a new type of basis functions is introduced. With this method, a dynamical model for NGC 3258 is constructed. This galaxy can be successfully modelled with a potential containing 30% dark matter within 1r_e with a mass of 1.6x10^11 M_o. The total mass within 2r_e is estimated as 5x10^11 M_o, containing 63% dark matter. The model is isotropic in the centre, is radially anisotropic between 0.2 and 2 kpc (0.88 r_e) and becomes tangentially anisotropic further on. The photometry reveals the presence of a dust disk near the centre.

In this paper we adopt a method to analyse absorption line spectra from elliptical galaxies that determines the dynamics of the galaxy and at the same time offers a way to study the stellar populations in that galaxy by a direct fit to the spectra. The result of the modelling is a distribution function for each stellar population that is considered in the fit. The method is described in detail in an accompanying paper. This paper reports on a dynamical stellar population study in NGC3258, based on spectra in two different wavelength regions, the near-IR Ca II triplet around 8600A and the Ca H and K lines around 3900A. These absorption lines have discriminating power toward various stellar types. The dynamical model shows an increase in dwarfs (represented by G2V stars) toward the centre. Most of the rotation in the model is delivered by the giants (represented by M1III stars). Moreover, the different models that were considered indicate that establishing a potential for a galaxy is dependent on the wavelength range used for the modelling.

We present new velocity dispersion measurements for a set of 12 spiral galaxies and use them to derive a more accurate V_c - sigma relation which holds for a wide morphological range of galaxies. Combined with the M_BH - sigma relation, this relation can be used as a tool to estimate supermassive black hole (SMBH) masses by means of the asymptotic circular velocity. Together with the Tully-Fisher relation, it serves as a constraint for galaxy formation and evolution models.

We present central CaT, PaT, and CaT* indices for a sample of fifteen dwarf elliptical galaxies (dEs). Twelve of these have CaT* ~ 7 A and extend the negative correlation between the CaT* index and central velocity dispersion sigma, which was derived for bright ellipticals (Es), down to 20 < sigma < 55 km/s. For five dEs we have independent age and metallicity estimates. Four of these have CaT* ~ 7 A, much higher than expected from their low metallicities (-1.5 < [Z/H] < -0.5). The observed anti-correlation of CaT* as a function of sigma or Z is in flagrant disagreement with theory. We discuss some of the amendments that have been proposed to bring the theoretical predictions into agreement with the observed CaT*-values of bright Es and how they can be extended to incorporate also the observed CaT*-values of dEs. Moreover, 3 dEs in our sample have CaT* ~ 5 A, as would be expected for metal-poor stellar systems. Any theory for dE evolution will have to be able to explain the co-existence of low-CaT* and high-CaT* dEs at a given mean metallicity. This could be the first direct evidence that the dE population is not homogeneous, and that different evolutionary paths led to morphologically and kinematically similar but chemically distinct objects.

We present SKIRT (Stellar Kinematics Including Radiative Transfer), a new Monte Carlo radiative transfer code that allows the calculation of the observed stellar kinematics of a dusty galaxy. The code incorporates the effects of both absorption and scattering by interstellar dust grains, and calculates the Doppler shift of the emerging radiation exactly by taking into account the velocities of the emitting stars and the individual scattering dust grains. The code supports arbitrary distributions of dust through a cellular approach, whereby the integration through the dust is optimized by means of a novel efficient trilinear interpolation technique. We apply our modelling technique to calculate the observed kinematics of realistic models for dusty disc galaxies. We find that the effects of dust on the mean projected velocity and projected velocity dispersion are severe for edge-on galaxies. For galaxies which deviate more than a few degrees from exactly edge-on, the effects are already strongly reduced. As a consequence, dust attenuation cannot serve as a possible way to reconcile the discrepancy between the observed shallow slopes of the inner rotation curves of LSB galaxies and the predictions of CDM cosmological models. For face-on galaxies, the velocity dispersion increases with increasing dust mass due to scattering, but the effects are limited, even for extended dust distributions. Finally, we show that serious errors can be made when the individual velocities of the dust grains are neglected in the calculations.

We present new velocity dispersion measurements of sample of 12 spiral galaxies for which extended rotation curves are available. These data are used to refine a recently discovered correlation between the circular velocity and the central velocity dispersion of spiral galaxies. We find a slightly steeper slope for our larger sample, we confirm the negligible intrinsic scatter on this correlation, and we find a striking agreement with a corresponding relation for elliptical galaxies. We combine this correlation with the well-known MBH-sigma relation to obtain a tight correlation between the circular velocities of galaxies and the masses of the supermassive black holes they host. This correlation is the observational evidence for an intimate link between dark matter haloes and supermassive black holes. Apart from being an important ingredient for theoretical models of galaxy formation and evolution, the relation between MBH and circular velocity can serve as a practical tool to estimate black hole masses in spiral galaxies.

Current error estimates on kinematic parameters are based on the assumption that the data points in the spectra follow a Poisson distribution. For realistic data that have undergone several steps in a reduction process, this is generally not the case. Neither is the noise distribution independent in adjacent pixels. Hence, the error estimates on the derived kinematic parameters will (in most cases) be smaller than the real errors. In this paper we propose a method that makes a diagnosis of the characteristics of the observed noise The method also offers the possibility to calculate more realistic error estimates on kinematic parameters. The method was tested on spectroscopic observations of NGC3258. In this particular case, the realistic errors are almost a factor of 2 larger than the errors based on least squares statistics.

We present photometric and kinematic evidence for the presence of stellar disks, seen practically edge-on, in two Fornax dwarf galaxies, FCC204 (dS0(6)) and FCC288 (dS0(7)). This is the first time such structures have been identified in Fornax dwarfs. FCC2088 has only a small bulge and a bright flaring and slightly warped disk that can be traced out to 23" from the center (2.05 kpc for H_0=75 km/s/Mpc). FCC204's disk can be traced out to 20" (1.78 kpc). This galaxy possesses a large bulge. These results can be compared to the findings of Jerjen et al. (2000) and Barazza et al. (2002) who discovered nucleated dEs with spiral and bar features in the Virgo Cluster.

We present a set of three-component Stackel potentials defined by five parameters and designed to model the Milky Way. We review the fundamental constraints that any model of the Milky Way must satisfy, including the most recent ones derived from Hipparcos data, and we study how the parameters of the presented potentials can vary in order to match these constraints. Five different valid potentials are presented and analyzed in detail: they are designed to be confronted with kinematical surveys in the future, by the construction of three-integral analytic distribution functions.

We present Halpha+[NII] narrow-band imaging of FCC046 and FCC207, two dES in the Fornax Cluster. B-R and B-I color-maps clearly show signs of dust-absorption in FCC207. FCC046 has a very bright blue nucleus, offset by about 1.1'' with respect to the outer isophotes. Moreover, FCC046 shows a pronounced lopsided shape. The emitting regions differ considerably between the two galaxies. Whereas FCC207 has only one central emission region, FCC046 also contains fainter emission regions. Based on broad-band colours, its disturbed shape and its very bright nucleus, FCC046 is akin to the class of amorphous dwarfs. The central emission regions of both galaxies are resolved and we estimate their diameters at about 60 pc. Their Halpha luminosities can be explained as photo-ionisation by post-AGB stars in an old population. Some of the emission regions in FCC046 are resolved and have diameters of the order of 50-150 pc and H$\alpha$ luminosities of the order of 10^30 W, comparable to supernova remnants or nebulae around Wolf-Rayet stars. Hence, FCC046 is clearly undergoing star-formation while for FCC207 the case is not as clearcut. We estimate the mass of the HII gas in FCC046 at M(HII) = 40-150 Mo. The ionised-gas content of FCC207 is somewhat higher : M(HII) = 60-190 Mo.

Simple analytical models, such as the Hernquist model, are very useful tools to investigate the dynamical structure of galaxies. Unfortunately, most of the analytical distribution functions are either isotropic or of the Osipkov-Merritt type, and hence basically one-dimensional. We present three different families of anisotropic distribution functions that self-consistently generate the Hernquist potential-density pair. These families have constant, increasing and decreasing anisotropy profiles respectively, and can hence represent a wide variety of orbital structures. For all of the models presented, the distribution function and the velocity dispersions can be written in terms of elementary functions. These models are ideal tools for a wide range of applications, in particular to generate the initial conditions for N-body or Monte Carlo simulations.

We present new equilibrium component distribution functions that depend on three analytic integrals in a Stackel potential, and that can be used to model stellar discs of galaxies. These components are generalizations of two-integral ones and can thus provide thin discs in the two-integral approximation. Their most important properties are the partly analytical expression for their moments, the disc-like features of their configuration space densities (exponential decline in the galactic plane and finite extent in the vertical direction) and the anisotropy of their velocity dispersions. We further show that a linear combination of such components can fit a van der Kruit disc.

This paper is the third one in a series, intended to investigate how the observed kinematics of elliptical galaxies are affected by dust attenuation. In Paper I and Paper II, we investigated the effects of dust absorption; here we extend our modelling in order to include the effects of scattering. We describe how kinematical information can be combined with the radiative transfer equation, and present a Monte Carlo code that can handle kinematical information in an elegant way. Compared to the case where only absorption is taken into account, we find that dust attenuation considerably affects the observed kinematics when scattering is included. For the central lines of sight, dust can either decrease or increase the central observed velocity dispersion. The most important effect of dust attenuation, however, is found at large projected radii. The kinematics at these lines of sight are strongly affected by photons scattered into these lines of sight, which were emitted by high-velocity stars in the central regions of the galaxy. These photons bias the LOSVDs towards high line-of-sight velocities, and significantly increase the observed velocity dispersion and LOSVD shape parameters. These effects are similar to the expected kinematical signature of a dark matter halo, such that dust attenuation may form an alternative explanation for the usual stellar kinematical evidence for dark matter halos around elliptical galaxies. We apply our results to discuss several other topics in galactic dynamics, where we feel dust attenuation should be taken into account. In particular, we argue that the kinematics observed at various wavelengths can help to constrain the spatial distribution of dust in elliptical galaxies.

In this paper, the influence of binary stars on the measured kinematics of dwarf galaxies is investigated. Using realistic distributions of the orbital parameters (semi-major axis, eccentricity, ...), analytical expressions are derived for the changes induced by the presence of binary stars in the measured velocity moments of low-mass galaxies (such as the projected velocity dispersion and the 4th order Gauss-Hermite coefficient h4). It is shown that there is a noticeable change in the observed velocity dispersion if the intrinsic velocity dispersion of a galaxy is of the same order as the binary velocity dispersion. The kurtosis of the line-of-sight velocity distribution (LOSVD) is affected even at higher values of the intrinsic velocity dispersion.

We present new tools to establish axisymmetric equilibrium models of the Milky Way. The models we wish to establish are pairs (V,F) where V is the gravitational potential generated by the whole mass distribution including the dark matter, and F is the distribution function in phase space for late-type tracer stars. We present a set of Stackel potentials that fit some fundamental parameters of the Milky Way (mass density in the solar neighbourhood and Oort constants). Then we define new component distribution functions that can be combined with these potentials in order to reproduce kinematical data like those that will be provided by GAIA.

Hierarchical models of galaxy formation and various observational evidence suggest that elliptical galaxies are, like disk galaxies, embedded in massive dark matter halos. Stellar kinematics are considered the most important tracer for this dark halo at a few effective radii. Using detailed modeling techniques, several authors have recently presented stellar kinematical evidence of a dark halo for a number of elliptical galaxies. In these modeling techniques, dust attenuation (absorption and scattering of starlight by dust grains) has not been taken into account. Nevertheless, elliptical galaxies contain a significant amount of interstellar dust, which affects all observable quantities, including the observed kinematics. We constructed a set of dynamical models for elliptical galaxies, in which dust attenuation is included through a Monte Carlo technique. We find that a dust component, shallower than the stellar distribution and with an optical depth of order unity, affects the observed kinematics significantly, in the way that it mimics the presence of a dark halo. If such dust distributions are realistic in elliptical galaxies, we are faced with a new mass-dust degeneracy. Taking dust into account in dynamical modeling procedures will hence reduce or may even eliminate the need for a dark matter halo at a few effective radii.

This is the first of two papers where we present the analysis of an HST NICMOS-Cam2 near-infared (NIR) snapshot survey in the F160W (H) filter for a sample of 78 spiral galaxies selected from the UGC and ESOLV catalogs. For 69 of these objects we provide nuclear color information derived by combining the H data either with additional NICMOS F110W (J) images or with V WFPC2/HST data. Here we present the NIR images and the optical-NIR color maps. We focus our attention on the properties of the photometrically-distinct `nuclei' which are found embedded in most of the galaxies, and provide measurements of their half-light radii and magnitudes in the H (and when available, in the J) band. We find that: (i) In the NIR, the nuclei embedded in the bright early- to intermediate-type galaxies span a much larger range in brightness than the nuclei which are typically found embedded in bulgeless late-type disks; (ii) Nuclei are found in both non-barred and barred hosts, in large-scale ($\gta$ 1 kpc) as well as in nuclear (up to a few 100pc) bars; (iii) There is a significant increase in half-light radius with increasing luminosity of the nucleus in the early/intermediate types; (iv) The nuclei of early/intermediate-type spirals cover a large range of optical-NIR colors, from $V-H\approx$ -0.5 to 3. Some nuclei are bluer and others redder than the surrounding galaxy, indicating the presence of activity or reddening by dust in many of these systems; (v) Some early/intermediate nuclei are elongated and/or slightly offset from the isophotal center of the host galaxy. On average, however, these nuclei appear as centered, star-cluster-like structures similar to those which are found in the late-type disks. [abridged]

We present surface brightness profiles for 56 of the 78 spiral galaxies observed in the HST/{\tt NICMOS/Camera-2} F160W snapshot survey introduced in Carollo et al. (2001; paper I), as well as surface brightness profiles for 23 objects out of the 41 that were also observed in the F110W filter. We fit these surface brightness profiles with the Nuker law of Lauer et al. (1995), and use the smooth analytical descriptions of the data to compute the average nuclear stellar cusp slopes $<\gamma>$ in the 0.1$''-0.5''$ radial range. Our main result is the startling similarity between the nuclear stellar cusp slopes $<\gamma>$ in the near-infared compared to those derived in the visual passband. This similarity has several implications: (1.) There are typically little or no optical-NIR {\it global} color gradients, and thus no global stellar population variations, inside $\sim 50$-100 pc from the nucleus in nearby spirals. (2.) The large observed range of the strength of the nuclear stellar cusps seen in the HST optical study of spiral galaxies reflects a physical difference between galaxies, and is not an artifact caused by nuclear dust and/or recent star formation. (3.) The dichotomy between $R^{1/4}$ bulges, with steep nuclear stellar cusps $<\gamma > \sim 1$, and exponential bulges, with shallow nuclear stellar cusps $<\gamma > < 0.3$, is also not an artifact of the effects of dust or recent star formation. (4) A surrounding massive disk appears to have no effect on the rise of the stellar density distribution within the innermost hundred pc of the

We present major and minor axis kinematics out to 2 half-light radii for the bright (M_B=-16.7) dwarf elliptical FS76, a member of the NGC5044 group. Its velocity dispersion is 46 \pm 2 km/s in the center and rises to 70 \pm 10 km/s at half-light radius. Beyond 1R_e the dispersion starts to fall again. The maximum rotation velocity is 15 \pm 6 km/s, about the value expected for an oblate isotropic rotator with the same flattening as FS76 (dE1). Hence, FS76 is the first dE discovered so far that is not flattened predominantly by anisotropy. Using dynamical models, we estimate the mass-to-light ratio to be between 3.2 and 9.1, consistent with predictions based on CDM cosmological scenarios for galaxy formation.

We investigate the influence of scattering and geometry on the attenuation curve in disc galaxies. We investigate both qualitatively and quantitatively which errors are made by either neglecting or approximating scattering, and which uncertainties are introduced due to a simplification of the star-dust geometry. We find that the magnitude of these errors depends on the inclination of the galaxy, and in particular that for face-on galaxies, the errors due to an improper treatment of scattering dominate those due to an imprecise star-dust geometry. Therefore we argue that in all methods which aim at determining the opacity of disc galaxies, scattering should be taken into account in a proper way.

Accurate photometric and kinematic modelling of disc galaxies requires the inclusion of radiative transfer models. Due to the complexity of the radiative transfer equation (RTE), sophisticated techniques are required. Various techniques have been employed for the attenuation in disc galaxies, but a quantitative comparison of them is difficult, because of the differing assumptions, approximations and accuracy requirements which are adopted in the literature. In this paper, we present an unbiased comparison of four methods to solve the RTE, in terms of accuracy, efficiency and flexibility. We apply them all on one problem that can serve as a first approximation of large portions of disc galaxies: a one-dimensional plane-parallel geometry, with both absorption and multiple scattering taken into account, with an arbitrary vertical distributions of stars and dust and an arbitrary angular redistribution of the scattering. We find that the spherical harmonics method is by far the most efficient way to solve the RTE, whereas both Monte Carlo simulations and the iteration method, which are straightforward to extend to more complex geometries, have a cost which is about 170 times larger.

Disk galaxies contain a large amount of interstellar dust, which affects the projection of kinematic quantities. We investigate in detail the effects of dust extinction on the mean projected velocity and the projected velocity dispersion. We use our results to construct a general strategy to determine the dynamical structure of disk galaxies, with the aim to constrain their mass distribution and dynamical history.

We have analysed $V$, $H$ and $J$ HST images for a sample of early- to late-type spiral galaxies, and reported elsewhere the statistical frequency of $R^{1/4}$-law and exponential bulges in our sample as a function of Hubble type, and the frequency of occurrence and structural properties of the resolved central nuclei hosted by intermediate- to late-type bulges and disks (see references in the text). Here we use these data to show that: {\it (i)} The $V-H$ color distribution of the $R^{1/4}$ bulges peaks around $<V-H> \sim 1.3$, with a sigma $\Delta (V-H) \sim 0.1$ magnitudes. Assuming a solar metallicity, these values correspond to stellar ages of $\approx 6 \pm 3$ Gyrs. In contrast, the $V-H$ color distribution of the exponential bulges peaks at $<V-H> \sim 0.9$ and has a sigma $\Delta (V-H) \sim 0.4$ mags. This likely implies significantly smaller ages and/or lower metallicities for (a significant fraction of the stars in) the exponential bulges compared to the $R^{1/4}$-law spheroids. {\it (ii)} Most of the central nuclei hosted by the exponential bulges have $V-H$ and $J-H$ colors which are compatible with relatively unobscured stellar populations. Assuming no or little dust effects, ages $\gta 1$ Gyrs are suggested for these nuclei, which in turn imply masses of about a few 10$^6$ to a few 10$^7$ M$_\odot$, i.e., sufficient to dissolve progenitor bars with masses consistent with those inferred for the exponential bulges by their luminosities. {\it (iii)} While different bulge-nucleus pairs cover a large range of $V-H$ colors, each bulge-nucleus pair has quite similar $V-H$ colors, and thus possibly similar stellar populations. (abridged)

Observations have demonstrated the presence of substantial amounts of interstellar dust in elliptical galaxies, most of which is believed to be distributed diffusely over the galaxy. Absorption by dust grains has a major impact on the transfer of stellar radiation, and affects the projection of each physical (and kinematic) quantity. In a previous paper, we have investigated the effects of a diffusely distributed dust component on the observed kinematics of spherical galaxies. In this paper we investigate the effect of not taking dust into account in dynamical modelling procedures. We use a set of semi-analytical dusty galaxy models to create synthetic LOSVD data sets, which we then model as if no dust were present. We find some considerable differences between the best-fitting models and the input models, and we find that these differences are dependent on the orbital structure of the input galaxy. For radial and isotropic models on the one hand, we find that the dynamical mass of the models decreases nearly linearly with optical depth, with an amount of 5 per cent per optical depth unit, whereas their orbital structure is hardly affected. For tangential models on the other hand, the dynamical mass decrease is smaller, but their orbital structure is affected: their distribution functions appears less tangentially anisotropic. For all models the mass-to-light ratio will be underestimated, by a factor of around 20 per cent per optical depth unit. We discuss these results in the light of the limited effects of dust extinction on the LOSVDs, as obtained in paper I, and conclude that the determination of the dynamical mass and the kinematic structure of galaxies... (for full abstract, see paper)

We present dynamical distribution functions for evolved stars in the inner galactic plane. We use an axisymmetric, two-component Stackel potential that satisfies recent constraints on the galactic potential, amongst others a slightly declining local rotation curve. We show that this potential is adequate to model stellar-kinematic samples with radial extent ranging from \~100 pc to ~5 kpc in the Galaxy. The two-integral model that gives the best fit to the first three projected moments provides a very good global representation of the data but fails to reproduce the central dispersion, the central apparent scaleheight and the cylindrical rotation at intermediate longitudes. All these features are fitted well by a three--integral model. We discuss various properties of the 2I- and 3I models and the implications for galactic structure. A somewhat thicker disk component is needed to explain the distribution of older AGB stars in the plane; this component also fits the kinematics at higher latitudes better. We find that the Disk and the Bulge, as traced by AGB stars, are very similar dynamically and could well be one and the same component. There is a dynamically distinct component in the inner 100 pc of the Bulge, however.

Observations show that early-type galaxies contain a considerable amount of interstellar dust, most of which is believed to exist as a diffusely distributed component. We construct a four-parameter elliptical galaxy model in order to investigate the effects of such a smooth absorbing component on the projection of kinematic quantities, such as the line profiles and their moments. We investigate the dependence on the optical depth and on the dust geometry. Our calculations show that both the amplitude and the morphology of these quantities can be significantly affected. Dust effects should therefore be taken in consideration when interpreting photometric and kinematic properties, and correlations that utilize these quantities.

We use a set of Stackel potentials to obtain a local approximation for an effective third integral in axisymmetric systems. We present a study on the feasibility and effectiveness of this approach. We have applied it to three trial potentials of various flattenings, corresponding to nearly ellipsoidal, disky and boxy density isophotes. In all three cases, a good fit to the potential requires only a small set of Stackel potentials, and the associated Stackel third integral provides a very satisfactory, yet analytically simple, approximation to the trial potentials effective third integral.

We present new dynamical models of the S0 galaxy N3115, making use of the available published photometry and kinematics as well as of two-dimensional TIGER spectrography. We first examined the kinematics in the central 40 arcsec in the light of two integral f(E,J) models. Jeans equations were used to constrain the mass to light ratio, and the central dark mass whose existence was suggested by previous studies. The even part of the distribution function was then retrieved via the Hunter & Qian formalism. We thus confirmed that the velocity and dispersion profiles in the central region could be well fit with a two-integral model, given the presence of a central dark mass of ~10^9 Msun. However, no two integral model could fit the h_3 profile around a radius of 25 arcsec where the outer disc dominates the surface brightness distribution. Three integral analytical models were therefore built using a Quadratic Programming technique. These models showed that three integral components do indeed provide a reasonable fit to the kinematics, including the higher Gauss-Hermite moments. Again, models without a central dark mass failed to reproduce the observed kinematics in the central arcseconds. This clearly supports the presence of a nuclear black hole of at least 6.5 10^8 Msun in the centre of NGC 3115. These models were finally used to estimate the importance of the dark matter in the outer part of NGC 3115, suggested by the flat stellar rotation curve observed by Capaccioli et al. (1993).

A variety of modeling techniques is used with surface photometry from the literature and recently acquired high-accuracy stellar kinematic data to constrain the three-dimensional mass distribution in the luminous cuspy elliptical galaxy NGC 1700. First, we model the radial velocity field and photometry, and, using a Bayesian technique, estimate the triaxiality T and short-to-long axis ratio c in five concentric annuli between approximately 1 and 3 effective radii. The results are completely consistent with T being constant inside about 2.5 r_e (36 arcsec; 6.7/h kpc). Adding an assumption of constant T as prior information gives an upper limit of T < 0.16 (95% confidence); this relaxes to T < 0.22 if it is also assumed that there is perfect alignment between the angular momentum and the galaxy's intrinsic short axis. Near axisymmetry permits us then to use axisymmetric models to constrain the radial mass profile. Using the Jeans (moment) equations, we demonstrate that 2-integral, constant-M/L models cannot fit the data; but a 2-integral model in which the cumulative enclosed M/L increases by a factor of roughly 2 from the center out to 12/h kpc can. Three-integral models constructed by quadratic programming show that, in fact, no constant-M/L model is consistent with the kinematics. Anisotropic 3-integral models with variable M/L, while not uniquely establishing a minimum acceptable halo mass, imply, as do the moment models, a cumulative M/L_B approximately 10 h at 12/h kpc. We conclude that NGC 1700 represents the best stellar dynamical evidence to date for dark matter in elliptical galaxies.

Galaxy spectra are a rich source of kinematical information since the shapes of the absorption lines reflect the movement of stars along the line-of-sight. We present a technique to directly build a dynamical model for a galaxy by fitting model spectra, calculated from a dynamical model, to the observed galaxy spectra. Using synthetic spectra from a known galaxy model we demonstrate that this technique indeed recovers the essential dynamical characteristics of the galaxy model. Moreover, the method allows a statistically meaningful error analysis on the resulting dynamical quantities.

A quasi self-consistent model for a barred structure in the central regions of our Galaxy is used to calculate the signature of such a triaxial structure on the kinematical properties of star samples. We argue that, due to the presence of a velocity dispersion, such effects are much harder to detect in the stellar component than in the gas. It might be almost impossible to detect stellar kinematical evidence for a bar using only l-v diagrams, if there is no a priori knowledge of the potential. Therefore, we propose some test parameters that can easily be applied to observed star samples, and that also incorporate distances or proper motions. We discus the diagnostic power of these tests as a function of the sample size and the bar strength. We conclude that about 1000 stars would be necessary to diagnose triaxiality with some statistical confidence.

We present a method for solving the two-dimensional linearized collisionless Boltzmann equation using Fourier expansion along the orbits. It resembles very much solutions present in the literature, but it differs by the fact that everything is performed in coordinate space instead of using action-angle variables. We show that this approach, though less elegant, is both feasible and straightforward. This approach is then incorporated in a matrix method in order to calculate self-consistent modes, using a set of potential-density pairs which is obtained numerically. We investigated the stability of some unperturbed disks having an almost flat rotation curve, an exponential disk and a non-zero velocity dispersion. The influence of the velocity dispersion, halo mass and anisotropy on the stability is further discussed.

A family of triaxial St\"ackel potential-density pairs is introduced. With the help of a Quadratic Programming method, a linear combination of potential-density pairs of this family which fits a given projected density distribution can be built. This deprojection strategy can be used to model the potentials of triaxial elliptical galaxies with or without dark halos. Besides, we show that the expressions for the St\"ackel triaxial density and potential are considerably simplified when expressed in terms of divided differences, which are convenient numerically. We present an example of triaxial deprojection for the galaxy NGC~5128 whose photometry follows the de Vaucouleurs law.

We present a dynamical analysis of 673 galactic Planetary Nebulae, using a two-integral axisymmetric model with a Kuzmin-Kutuzov St\"{a}ckel potential. The method fits the kinematics to the projected moments of a distribution function, by means of Quadratic Programming. The 2.2 $\mu$m COBE brightness map has been used after correction for the interstellar extinction as a projected star counts map in the modeling, because it constitutes a galactic distribution view of evolved red populations which are considered to be the progenitors of PNe. The model we have obtained provides a 2-integral distribution function for the COBE 2.2 $\mu$m map, and thus {\it a fortiori} a deprojection of it, which allows moreover the identification of all the major Galactic components. We derive the density laws for them. The projected velocity dispersions are not well fitted though, especially in the disk, which points at the likely presence of a third integral. If this result can be confirmed by additional data, this would mean that for the first time the presence and importance of a third integral on a global scale is demonstrated.

Photometry and kinematics of the giant elliptical galaxy NGC~5128 (Centaurus~A) based on planetary nebulae observations (Hui~\etal 1995) are used to build dynamical models which allow us to infer the presence of a dark matter halo. To this end, we apply a Quadratic Programming method. Constant mass-to-light ratio models fail to reproduce the major axis velocity dispersion measurements at large radii: the profile of this kind of models falls off too steeply when compared to the observations, clearly suggesting the necessity of including a dark component in the halo. By assuming a mass-to-light ratio which is increasing with radius, the model satisfactorily matches the observations. The total mass for the best fit model is $\sim4\times10^{11}M_\odot$ of which about 50\% is dark matter. However, models with different total masses and dark halos are also consistent with the data; we estimate that the total mass of Cen~A within 50~kpc may vary between $3\times10^{11}M_\odot$ and $5\times10^{11}M_\odot$. The best fit model consists of 75\% of stars rotating around the short axis $z$ and 25\% of stars rotating around the long axis $x$. Finally, the morphology of the projected velocity field is analyzed using Statler's classification criteria (Statler 1991). We find that the appearance of our velocity field is compatible with a type 'Nn' or 'Nd'.

Using the photometry and the kinematic data, part of which acquired with the NTT, we constructed a 2-integral and a 3-integral distribution function for NGC 4697. We detected a nuclear dust lane at $3.4''$ or 0.4 kpc from the centre. A comparison of the Lucy-deprojection method and the multi-gaussian expansion method showed that the latter offers a number of advantages, among which a qualitatively better fit and better controllable convergence. Futhermore, we show that a 2-integral model cannot be excluded for as far as the Satoh parameter $k$ is allowed to vary. The 3-integral model is based on a St\"ackel third integral. The original potential is retained where appropriate, in order to minimise errors due to the St\"ackel approximation. On the whole, the 3-integral model produces more satisfying results than the 2-integral model, and is better in handling some peculiarities in the kinematic data. For both models we show detailed contour plots for all deprojected moments up to order 2. The mass-to-light ratio appears to be a fairly well determined $(M/L)_B=4.8h_{50}$, and is much better constrained than would be the case for spherical anisotropic models with comparable data.

We explore a series expansion method to calculate the modes of oscillations for a variety of uniformly rotating finite disks, either with or without a dark halo. Since all models have the same potential, this survey focuses on the role of the distribution function in stability analyses. We show that the stability behaviour is greatly influenced by the structure of the unperturbed distribution, particularly by its energy dependence. In addition we find that uniformly rotating disks with a halo in general can feature spiral-like instabilities.

Along with data on radial velocities, more and more data on proper motions of individual stars in globular clusters are becoming available. Their usage was until now rather limited. It was mostly restricted to determining cluster membership of individual stars and to determining the spatial velocity of the cluster. We study the two dimensional distribution of the proper motions, in order to clarify the relation between the projection of the velocities and the dynamical structure. Obviously some assumptions on the dynamics of the system have to be made. In this study we chose a Plummer model. This complements an earlier study about line profiles. The three velocity dispersions in the main directions are determined and compared with the three spherical velocity dispersions. Two parameters analogous to Binney's anisotropy parameter $\beta$ are defined. Moreover, since the proper motion components are distance dependent and the line of sight velocity dispersion is distance independent, we show explicitly how to determine distances.

We construct analytical phase-space solutions for perturbations of flat disks by performing a power series expansion for the radius and the velocity coordinates. We show that this approach translates into an elegant mathematical formulation which is easy to use for a wide variety of distribution functions, for as far as resonances do not play a role, such as is the case for potentials which are close to quadratic. As a testcase, the method is applied on the Kalnajs disks. The results obtained are in full agreement with the analytical solutions of the mode analysis. The strongest advantages of this method are its independence of the mathematical complexity of the unperturbed distribution, the degree of detail with which the solutions can be calculated and its computational straightforwardness. On the contrary, power series solutions are not suitable for describing regions where resonant orbits occur, which we therefore exclude in this paper. We used the technique to analyse perturbations in the central regions of a galaxy, tracking the dynamical consequences of a Galactic bar on the kinematics of the solar neighbourhood (Hipparcos). We showed how the orientation and strength of the bar is related to the properties of the velocity ellipsoid in our model.