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.

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.

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.

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.