The cold dark matter (CDM) scenario generically predicts the existence of triaxial dark matter haloes which contain notable amounts of substructure. However, analytical halo models with smooth, spherically symmetric density profiles are routinely adopted in the modelling of light propagation effects through such objects. In this paper, we address the biases introduced by this procedure by comparing the surface mass densities of actual N-body haloes against the widely used analytical model suggested by Navarro, Frenk and White (1996) (NFW). We conduct our analysis in the redshift range of 0.0 - 1.5. In cluster sized haloes, we find that triaxiality can cause scatter in the surface mass density of the haloes up to sigma_+ = +60% and sigma_- = -70%, where the 1-sigma limits are relative to the analytical NFW model given value. Subhaloes can increase this scatter to sigma_+ = +70% and sigma_- = -80%. In galaxy sized haloes, the triaxial scatter can be as high as sigma_+ = +80% and sigma_- = -70%, and with subhaloes the values can change to sigma_+ = +40% and sigma_- = -80%. We present an analytical model for the surface mass density scatter as a function of distance to the halo centre, halo redshift and halo mass. The analytical description enables one to investigate the reliability of results obtained with simplified halo models. Additionally, it provides the means to add simulated surface density scatter to analytical density profiles. As an example, we discuss the impact of our results on the calculation of microlensing optical depths for MACHOs in CDM haloes.

If the concordance $\Lambda$CDM model is a true description of the universe, it should also properly predict the properties and structure of dark matter haloes, where galaxies are born. Using N-body simulations with a broad scale of mass and spatial resolution, we study the structure of dark matter haloes, the distribution of masses and the spatial distribution of subhaloes within the main haloes. We carry out three $\Lambda$CDM simulations with different resolutions using the AMIGA code. Dark matter haloes are identified using an algorithm that is based on the adaptive grid structure of the simulation code. The haloes we find encompass the mass scales from $10^8\mathrm{M}_{\sun}$ to $10^{15}\mathrm{M}_{\sun}$. We find that if we have to study the halo structure (search for subhaloes), the haloes have to contain at least $10^4$ particles. For such haloes, where we can resolve substructure, we determined the subhalo mass function and found that it is close to a power law with the slope -0.9 (at present time), consistent with previous studies. This slope depends slightly on the redshift and it is approximately the same for main haloes. The subhalo mass fraction ($M_{\s subH}/M_{\s MH}$) is between 0.08 and 0.2, increasing slightly with redshift and with the mass of the main halo. Its distribution is approximated using the Weibull distribution at different epochs. The mean values of subhalo mass are independent of the main halo mass. The spatial density of subhaloes, scaled to the virial radius of the main halo ($r_{\s vir}$), is independent of redshift and follows the $r^{1/3}$ rule.

Using eight dark matter haloes extracted from fully-self consistent cosmological N-body simulations, we perform microlensing experiments. A hypothetical observer is placed at a distance of 8.5 kpc from the centre of the halo measuring optical depths, event durations and event rates towards the direction of the Large Magellanic Cloud. We simulate 1600 microlensing experiments for each halo. Assuming that the whole halo consists of MACHOs, f = 1.0, and a single MACHO mass is m_M = 1.0 Msun, the simulations yield mean values of tau = 4.7^{+5.0}_{-2.2} 10^{-7} and Gamma = 1.6^{+1.3}_{-0.6} 10^{-6} events/star/year. We find that triaxiality and substructure can have major effects on the measured values so that tau and Gamma values of up to three times the mean can be found. If we fit our values of tau and Gamma to the MACHO collaboration observations (Alcock et al. 2000), we find f = 0.23^{+0.15}_{-0.13} and m_M = 0.44^{+0.24}_{-0.16}. Five out of the eight haloes under investigation produce f and m_M values mainly concentrated within these bounds.

We place observational constraints on a recently proposed Galactic population, dubbed the shroud (Gyuk & Gates 1999, Gates & Gyuk 2001). The shroud would be a very thick Galactic disk of low luminosity objects, most likely old white dwarfs, proposed to explain the optical depth seen in microlensing surveys towards the Magellanic clouds. The shroud is a simple alternative to the lenses being distributed in a classical, near-spherical dark halo; the advantage of the shroud is that it would compose only a fraction of a dark halo's total mass. In this paper, we argue that stars of the Galactic shroud would be detectable in the recent proper motion survey of Oppenheimer et al. (2001) if their absolute luminosities were brighter than M_R_59F = 19.4 or approximately M_V = 18.6. We adopt a range of simple models of the shroud's kinematics and morphology, and the colours and luminosities of its white dwarfs; via Monte-Carlo simulations, we predict the numbers expected in the Oppenheimer et al. survey, which would be clearly separated from the numbers produced by white dwarfs of the disk, thick disk and halo. The number of white dwarf detections in the proper motion survey (98) is found to be well explained by the disk, thick disk and halo. With the most conservative kinematic and density parameters for the shroud, and an absolute luminosity of the white dwarfs of M_R_59F = 17.6, we find that the proper motion survey would detect over 100 WDs, just from the shroud. For a M_R_59F = 19.4 shroud, the survey would find 5 +/- 2 peculiar objects, whereas only two white dwarfs with such characteristics are found in the original data. M_R_59F = 19.4 corresponds to M_V = 18.6 for WDs with V - I = -1.030.

Using stellar population synthesis techniques, we explore the photometric signatures of white dwarf progenitor dominated galactic halos, in order to constrain the fraction of halo mass that may be locked-up in white dwarf stellar remnants. We first construct a 10^9 M_sun stellar halo using the canonical Salpeter initial stellar mass distribution, and then allow for an additional component of low- and intermediate-mass stars, which ultimately give rise to white dwarf remnants. Microlensing observations towards the Large Magellanic Cloud, coupled with several ground-based proper motion surveys, have led to claims that in excess of 20% of the dynamical mass of the halo (10^12 M_sun) might be found in white dwarfs. Our results indicate that (1) even if only 1% of the dynamical mass of the dark halo today could be attributed to white dwarfs, their main sequence progenitors at high redshift (z ~ 3) would have resulted in halos more than 100 times more luminous than those expected from conventional initial mass functions alone, and (2) any putative halo white dwarf progenitor dominated initial mass function component, regardless of its dynamical importance, would be virtually impossible to detect at the present-day, due to its extremely faint surface brightness.

We discuss the recent discovery by Oppenheimer et al (2001) of old, cool white dwarf stars, which may be the first direct detection of Galactic halo dark matter. We argue here that the contribution of more mundane white dwarfs of the stellar halo and thick disk would contribute sufficiently to explain the new high velocity white dwarfs without invoking putative white dwarfs of the dark halo. This by no means rules out that the dark matter has been found, but it does constrain the overall contribution by white dwarfs brighter than M_V ~ 16 to significantly less than 1% of the Galactic dark matter. This work confirms a similar study by Reyle et al (2001).