In this paper we clarify and generalise previous work by Moser and Belbruno concerning the link between the motions in the classical Kepler problem and geodesic motion on spaces of constant curvature. Both problems can be formulated as Hamiltonian systems and the phase flow in each system is characterised by the value of the corresponding Hamiltonian and one other parameter (the mass parameter in the Kepler problem and the curvature parameter in the geodesic motion problem). Using a canonical transformation the Hamiltonian vector field for the geodesic motion problem is transformed into one which is proportional to that for the Kepler problem. Within this framework the energy of the Kepler problem is equal to (minus) the curvature parameter of the constant curvature space and the mass parameter is given by the value of the Hamiltonian for the geodesic motion problem. We work with the corresponding family of evolution spaces and present a unified treatment which is valid for all values of energy continuously. As a result, there is a correspondence between the constants of motion for both systems and the Runge-Lenz vector in the Kepler problem arises in a natural way from the isometries of a space of constant curvature. In addition, the canonical nature of the transformation guarantees that the Poisson bracket Lie algebra of constants of motion for the classical Kepler problem is identical to that associated with geodesic motion on spaces of constant curvature.

We consider the polarization arising from scattering in an envelope illuminated by a central anisotropic source. Spherical harmonics are used to describe both the light source anisotropy and the envelope density distribution functions of the scattering particles. This framework demonstrates how the net resultant polarization arises from a superposition of three basic "shape" functions: the distribution of source illumination, the distribution of envelope scatterers, and the phase function for dipole scattering. Specific expressions for the Stokes parameters and scattered flux are derived for the case of an ellipsoidal light source inside an ellipsoidal envelope, with principal axes that are generally not aligned. Two illustrative examples are considered: (a) axisymmetric mass loss from a rapidly rotating star, such as may apply to some Luminous Blue Variables, and (b) a Roche-lobe filling star in a binary system with a circumstellar envelope. As a general conclusion, the combination of source anisotropy with distorted scattering envelopes leads to more complex polarimetric behavior such that the source characteristics should be carefully considered when interpreting polarimetric data.

Scattering of starlight by dust, molecules and electrons in spiral galaxies will produce a modification of the direct intensity and a polarization in the observed light. We treat the case where the distribution of scatterers can be considered to be optically thin, and derive semi-analytic expressions for the resolved intensity and polarized intensity for Thomson, Rayleigh, and more general scattering mechanisms. These expressions are applied to a parametric model spiral galaxies. It is further shown that in the case of Thomson and Rayleigh scattering, and when scatterers and stars are distributed with rotational symmetry, the total polarized flux depends on the inclination, $i$, of the galactic axis to the line of sight according to a simple $\sin ^2 i$ law. This generalises the well known result for pointlike and spherical light sources. By using a method based on spherical harmonics, we generalise this law for more general mechanisms, and show that to good approximation, the $\sin ^2 i$ law still holds for the class of models considered.

The polarization of the CMBR scattered by galaxy clusters in the kinematic Sunyaev Zeldovich effect depends on the transverse velocity of the cluster. This polarizing effect is proportional to the transverse velocity squared, and so weaker that the change in intensity due to the radial motion in the kinematic effect. The value given by Sunyaev and Zeldovich, and which is frequently cited, underestimates the polarizing effect by a factor of ten. We show furthermore that the polarization has a strong frequency dependence. This means that the polarization should be detectable with the new generation of CMBR probes, at least for some clusters. Thus this effect offers, almost uniquely, a method of obtaining the vectorial velocity of clusters.

The measurement of the integrated optical polarization of weakly gravitationally lensed galaxies can provide considerable constraints on lens models. The method outlined depends on fact that the orientation of the direction of optical polarization is not affected by weak gravitational lensing. The angle between the semi-major axis of the imaged galaxy and the direction of integrated optical polarization thus informs one of the distortion produced by the gravitational lensing. Although the method depends on the polarimetric measurement of faint galaxies, large telescopes and improved techniques should make such measurements possible in the near future.

It has already been shown that microlensing can give rise to a non-zero variable polarisation signal. Here we use realistic simulations to demonstrate the additional information that can be gained from polarimetric observations of lensing events.

We investigate the feasibility of reconstructing the radial intensity profile of extended stellar sources by inverting their microlensed light curves. Using a simple, linear, limb darkening law as an illustration, we show that the intensity profile can be accurately determined, at least over the outer part of the stellar disc, with realistic light curve sampling and photometric errors. The principal requirement is that the impact parameter of the lens be less than or equal to the stellar radius. Thus, the analysis of microlensing events provides a powerful method for testing stellar atmosphere models.

We describe a method for determining the limb polarization and limb darkening of stars in eclipsing binary systems, by inverting photometric and polarimetric light curves. Because of the ill-conditioning of the problem, we use the Backus-Gilbert method to control the resolution and stability of the recovered solution, and to make quantitative estimates of the maximum accuracy possible. Using this method we confirm that the limb polarization can indeed be recovered, and demonstrate this with simulated data, thus determining the level of observational accuracy required to achieve a given accuracy of reconstruction. This allows us to set out an optimal observational strategy, and to critcally assess the claimed detection of limb polarization in the Algol system. The use of polarization in stars has been proposed as a diagnostic tool in microlensing surveys by Simmons et al. (1995), and we discuss the extension of this work to the case of microlensing of extended sources.

In this paper we outline the framework of mathematical statistics with which one may study the properties of galaxy distance estimators. We describe, within this framework, how one may formulate the problem of distance estimation as a Bayesian inference problem, and highlight the crucial question of how one incorporates prior information in this approach. We contrast the Bayesian approach with the classical `frequentist' treatment of parameter estimation, and illustrate -- with the simple example of estimating the distance to a single galaxy in a redshift survey -- how one can obtain a significantly different result in the two cases. We also examine some examples of a Bayesian treatment of distance estimation -- involving the definition of Malmquist corrections -- which have been applied in recent literature, and discuss the validity of the assumptions on which such treatments have been based.

Microlensing of extended stellar sources in the LMC and Galaxy by low mass lenses can produce variable polarisation. The characteristics of the polarisation and flux profiles can provide considerable information about the lens geometry, and help determine the importance of such low mass objects as a dark matter component.

Much interest has been generated recently by the ongoing MACHO, EROS and OGLE projects to identify gravitationally lensed stars from the Large Magellanic Cloud and Galactic bulge, and the positive identification of several events (Alcock et al, (1993), Aubourg et al, (1993) and Udalski et al, (1993)). The rate at which such events are found should provide considerable information about the distribution of the low mass compact objects, brown dwarfs etc. responsible for the lensing, and their relative importance as a component of the cosmological `dark matter'. We show measurement of the polarisation of starlight during these events can yield considerably more information about the lensing objects than was previously considered possible. Furthermore, the consideration of extended sources is shown to have a significant effect on the interpretation of the profiles and statistics of the events.

Much current interest in cosmology has revolved around the peculiar velocity field of galaxies. Many methods have been proposed to derive this, probably the most prominent being POTENT (Dekel et al. 1990). However, this is prone to many forms of systematic error or bias. Here, we examine some of the methods currently used to minimise these biases, particularly Malmquist corrections, and highlight some of their pitfalls. As an alternative we propose an iterative scheme based on Monte Carlo simulations of POTENT that shows noticeable improvement over the other methods considered.

Although Potent purports to use only radial velocities in reconstructing the potential velocity field of galaxies, the derivation of transverse components is implicit in the smoothing procedures adopted. Thus the possibility arises of using nonradial line integrals to derive a smoothed velocity field. For an inhomogeneous galaxy distribution the optimal path for integration need not be radial, and can be obtained using max-flow algorithms. In this paper we describe how one may use Dijkstra's algorithm to obtain this optimal path and velocity field, and present the results of applying the algorithm to a realistic spatial distribution of galaxies. These results show that the method has limited effect due to the large smoothing scales employed in Potent. However, the viability of the technique is demonstrated and, finally, we discuss other possible methods involving averaging over an ensemble of non-radial paths for improving a potential velocity field derived from redshifts.

The statistical properties of galaxy distance estimators are studied and a rigorous framework is developed for identifying and removing the effects of Malmquist bias due to obsevational selection. The prescription of Schechter (1980) for defining unbiased distance estimators is extended to more general -- and more realistic -- cases. The derivation of `optimal' unbiased distance estimators of minimum dispersion, by utilising information from additional -- suitably correlated -- observables, is discussed and the results applied to a calibrating sample from the Fornax cluster, as used in the Mathewson spiral galaxy redshift survey. The optimal distance estimator derived from I-band magnitude, diameter and 21cm line width has an intrinsic scatter which is 25 \% smaller than that of the Tully-Fisher relation quoted for this calibrating sample. (Figures are available on request).

In this paper we review the two main approaches to the problem of Malmquist bias which have been adopted in the cosmology literature, and show how these two formulations of the problem represent fundamentally different views of the nature of probability. We discuss the assumptions upon which both approaches are based and indicate some of their limitations. In particular we identify a basic flaw in the definition of homogeneous and inhomogeneous Malmquist corrections as they have frequently been applied in the literature, and indicate how this flaw may be corrected. (Figures available on request).

We investigate the effect of using different distance estimators on the recovery of the peculiar velocity field of galaxies using Potent. An inappropriate choice of distance estimator will give rise to spurious flows. We discuss methods of minimising these biases and the levels of accuracy required of distance estimators to retrieve velocity fields to a given standard.

Although Potent purports to use only radial velocities in retrieving the potential velocity field of galaxies, the derivation of transverse components is implicit in the smoothing procedures. Thus the possibility of using nonradial line integrals to derive the velocity field arises. In the case of inhomogeneous distributions of galaxies, the optimal path for integration need not be radial, and can be obtained by using max-flow algorithms. In this paper we present the results of using Dijkstra's algorithm to obtain this optimal path and velocity field.