We constrain the most likely range of masses for the Milky Way and M31 using an application of the Numerical Action Method (NAM) that optimizes the fit to observed parameters over a large ensemble of NAM-generated solutions. Our 95% confidence level mass ranges, 1.5-4.5 x 10^12 m_sun for MW and 1.5-5.5 x 10^12 m_sun for M31, are consistent with the upper range of estimates from other methods and suggests that a larger proportion of the total mass becomes detectable when the peculiar motions of many nearby satellites are taken into account in the dynamical analysis. We test the method against simulated Local Group catalogs extracted from the Millennium Run to confirm that mass predictions are consistent with actual galaxy halo masses.

We test the ability of the numerical action method (NAM) to recover the individual orbit histories of mass tracers in an expanding universe in a region of radius 26Mpc/h, given the masses and redshift-space coordinates at the present epoch. The mass tracers are represented by dark matter haloes identified in a high resolution N-body simulation of the standard LCDM cosmology. Since previous tests of NAM at this scale have traced the underlying distribution of dark matter particles rather than extended haloes, our study offers an assessment of the accuracy of NAM in a scenario which more closely approximates the complex dynamics of actual galaxy haloes. We show that NAM can recover present-day halo distances with typical errors of less than 3 per cent, compared to 5 per cent errors assuming Hubble flow distances. The total halo mass and the linear bias were both found to be constained at the 50 per cent level. The accuracy of individual orbit reconstructions was limited by the inability of NAM, in some instances, to correctly model the positions of haloes at early times solely on the basis of the redshifts, angular positions, and masses of the haloes at the present epoch. Improvements in the quality of NAM reconstructions may be possible using the present-day three-dimensional halo velocities and distances to further constrain the dynamics. This velocity data is expected to become available for nearby galaxies in the coming generations of observations by SIM and GAIA.

Analysis of the peculiar velocities of galaxies should take account of the uncertainties in both redshifts and distances. We show how this can be done by a numerical application of the action principle. The method is applied to an improved catalog of the galaxies and tight systems of galaxies within 4$h_{75}^{-1}$ Mpc, supplemented with a coarser sample of the major concentrations at 4$h_{75}^{-1}$ Mpc to 20$h_{75}^{-1}$ Mpc distance. Inclusion of this outer zone improves the fit of the mass tracers in the inner zone to their measured redshifts and distances, yielding best fits with reduced $\chi ^2$ in redshift and distance in the range 1.5 to 2. These solutions are based on the assumption that the galaxies in and near the Local Group trace the mass, and a powerful test would be provided by observations of proper motions of the nearby galaxies. Predicted transverse galactocentric velocities of some of the nearby galaxies are confined to rather narrow ranges of values, and are on the order of 100 km~s$^{-1}$, large enough to be detected and tested by the proposed SIM and GAIA satellite missions.

We examine the correlational property of $\lya$ clouds in detail and compare it to that of mass and galaxies. We show that the correlation strength of $\lya$ clouds is somewhat weaker than that of the underlying matter, which in turn is weaker than that of galaxies (biased galaxy formation). On the scales probed, $10-300$km/s, for $\lya$ clouds, we find that higher density, higher optical depth, higher column density regions are more strongly clustered than lower density, lower optical depth, lower column density regions. Thus, a consistent picture seems to emerge: the correlation strength for a given set of objects is positively correlated with their characteristic global density and the differences among the correlations of galaxies, $\lya$ clouds and mass reflect the differences in density that each trace. Significant positive correlations with a strength of $0.1-1.0$ are found for $\lya$ clouds in the velocity range $50-300$km/s. Among the correlational measures examined, an optical depth correlation function (Equation 5) proposed here may serve as the best correlational measure. It reasonably faithfully represents the true correlation of the underlying matter, enabling a better indication of both matter correlation and the relationship between galaxies and $\lya$ clouds. Furthermore, it appears to be an alternative to the conventional line-line correlation function with the virtue that it does not require ambiguous post-observation fitting procedures such as those commonly employed in the conventional line-finding methods. Neither does it depend sensitively on the observational resolution (e.g., FWHM), insofaras the clouds are resolved Conveniently, it can be easily measured with the current observational sensitivity without being contaminated significantly by the presence of noise.