Primordial black holes (PBHs) could provide the dark matter but a variety of constraints restrict the possible mass windows to $10^{16}$ - $10^{17}\,$g, $10^{20}$ - $10^{24}\,$g and $10$ - $10^{3}\,M_{\odot}$. The last possibility is contentious but of special interest in view of the recent detection of black-hole mergers by LIGO/Virgo. PBHs might have important consequences and resolve various cosmological conundra even if they have only a small fraction of the dark-matter density. In particular, those larger than $10^{3}\,M_{\odot}$ could generate cosmological structures through the `seed' or `Poisson' effect, thereby alleviating some problems associated with the standard cold dark matter scenario, and sufficiently large PBHs might provide seeds for the supermassive black holes in galactic nuclei. More exotically, the Planck-mass relics of PBH evaporations or stupendously large black holes bigger than $10^{12}\,M_{\odot}$ could provide an interesting dark component.

We update the constraints on the fraction of the Universe going into primordial black holes (PBHs) over the mass range $10^{-5}$--$10^{50}$ g. Those smaller than $\sim 10^{15}$ g would have evaporated by now due to Hawking radiation, so their abundance at formation is constrained by the effects of evaporated particles on big bang nucleosynthesis, the cosmic microwave background (CMB), the Galactic and extragalactic $\gamma$-ray and cosmic ray backgrounds and the possible generation of stable Planck mass relics. PBHs larger than $\sim 10^{15}$ g are subject to a variety of constraints associated with gravitational lensing, dynamical effects, influence on large-scale structure, accretion and gravitational waves. We discuss the constraints on both the initial collapse fraction and the current fraction of the cold dark matter in PBHs at each mass scale but stress that many of the constraints are associated with observational or theoretical uncertainties and some are no longer applicable. We also consider indirect constraints associated with the amplitude of the primordial density fluctuations, such as second-order tensor perturbations and $\mu$-distortions arising from the effect of acoustic reheating on the CMB, but these only apply if PBHs are created from the high-$\sigma$ peaks of nearly Gaussian fluctuations. Finally we discuss how the constraints are modified if the PBHs have an extended mass function, this being relevant if PBHs provide some combination of the dark matter, the LIGO/Virgo coalescences and the seeds for cosmic structure.

Stephen Hawking's contributions to the understanding of gravity, black holes and cosmology were truly immense. They began with the singularity theorems in the 1960s followed by his discovery that black holes have an entropy and consequently a finite temperature. Black holes were predicted to emit thermal radiation, what is now called Hawking radiation. He pioneered the study of primordial black holes and their potential role in cosmology. His organisation of and contributions to the Nuffield Workshop in 1982 consolidated the picture that the large-scale structure of the universe originated as quantum fluctuations during the inflationary era. Work on the interplay between quantum mechanics and general relativity resulted in his formulation of the concept of the wavefunction of the universe. The tension between quantum mechanics and general relativity led to his struggles with the information paradox concerning deep connections between these fundamental areas of physics. These achievements were all accomplished following the diagnosis during the early years of Stephen's studies as a post-graduate student in Cambridge that he had incurable motor neuron disease -- he was given two years to live. Against all the odds, he lived a further 55 years. The distinction of his work led to many honours and he became a major public figure, promoting with passion the needs of disabled people. His popular best-selling book A Brief History of Time made cosmology and his own work known to the general public worldwide. He became an icon for science and an inspiration to all.

A universal mechanism may be responsible for several unresolved cosmic conundra. The sudden drop in the pressure of relativistic matter at $W^{\pm} / Z^{0}$ decoupling, the quark-hadron transition and $e^{+}e^{-}$ annihilation enhances the probability of primordial black hole (PBH) formation in the early Universe. Assuming the amplitude of the primordial curvature fluctuations is approximately scale-invariant, this implies a multi-modal PBH mass spectrum with peaks at $10^{-6}$, $1$, $30$, and $10^{6}\,M_{\odot}$. This suggests a unified PBH scenario which naturally explains the dark matter and recent microlensing observations, the LIGO/Virgo black hole mergers, the correlations in the cosmic infrared and X-ray backgrounds, and the origin of the supermassive black holes in galactic nuclei at high redshift. A distinctive prediction of our model is that LIGO/Virgo should soon observe the merging of black holes in the mass gap between $2$ and $5\,M_{\odot}$ or above $60\,M_{\odot}$ and with low mass ratios. The newly detected events S190924h and S190930s could confirm our predictions and thereby indicate a primordial origin for the black holes. The detection of PBHs would also offer a novel way to probe the existence of new particles with mass between $1\,$MeV and $10^{10}\,$GeV.

A mechanism for generating primordial black-hole mass spectra with many spikes is proposed and investigated. This mechanism relies on the choice of non-Bunch-Davies vacua, these leading to oscillatory features in the inflationary power spectrum. This in turn generates oscillations in the primordial black-hole mass function with exponentially enhanced spikes. This ``multimodal'' effect is demonstrated for most of the well-studied models of primordial black-hole formation.

The origin of the baryon asymmetry of the Universe (BAU) and the nature of dark matter are two of the most challenging problems in cosmology. We propose a scenario in which the gravitational collapse of large inhomogeneities at the quark-hadron epoch generates both the baryon asymmetry and dark matter in the form of primordial black holes (PBHs). This is due to the sudden drop in radiation pressure during the transition from a quark-gluon plasma to non-relativistic hadrons. The collapse to a PBH is induced by fluctuations of a light spectator scalar field in rare regions and is accompanied by the violent expulsion of surrounding material, which might be regarded as a sort of "primordial supernova" . The acceleration of protons to relativistic speeds provides the ingredients for efficient baryogenesis around the collapsing regions and its subsequent propagation to the rest of the Universe. This scenario naturally explains why the observed BAU is of order the PBH collapse fraction and why the baryons and dark matter have comparable densities. The predicted PBH mass distribution ranges from sub-solar to several hundred solar masses. This is compatible with current observational constraints and could explain the rate, mass and low spin of the black hole mergers detected by LIGO-Virgo. Future observations will soon be able to test this scenario.

We revisit the cosmological and astrophysical constraints on the fraction of the dark matter in primordial black holes (PBHs) with an extended mass function. We consider a variety of mass functions, all of which are described by three parameters: a characteristic mass and width and a dark matter fraction. Various observations then impose constraints on the dark matter fraction as a function of the first two parameters. We show how these constraints relate to those for a monochromatic mass function, demonstrating that they usually become more stringent in the extended case than the monochromatic one. Considering only the well-established bounds, and neglecting the ones that depend on additional astrophysical assumptions, we find that there are three mass windows, around $4\times 10^{-17}M_\odot,$ $2\times 10^{-14}M_\odot$ and $25-100M_\odot$, where PBHs can constitute all dark matter. However, if one includes all the bounds, PBHs can only constitute of order $10\%$ of the dark matter.

We propose a scenario in which the quantum fluctuations of a light stochastic spectator field during inflation, for instance the QCD axion, generate large curvature fluctuations in the radiation era. This inevitably leads to the existence of Hubble patches in which primordial black holes (PBHs) form at the quark-hadron epoch with a density comparable to the dark matter. The PBHs naturally have an extended mass function, with a density peaking at the Chandrasekhar mass, and extending to the mass required to explain the black hole coalescences observed by LIGO/Virgo, without violating current astrophysical constraints. At the quark-hadron epoch, the entropy production and temperature increase in PBH-forming regions provide the ingredients for a novel scenario of hot spot electroweak baryogenesis. Baryons are produced very efficiently around these collapsing pockets and then quickly diffuse to the entire Universe. This explains why the PBH collapse fraction at the QCD epoch is of order the observed baryon-to-photon ratio and why baryons and dark matter have comparable densities. If the stochastic spectator field is the QCD axion, this also explains the strong CP problem. Parameter fine-tunings are replaced by a single anthropic selection argument involving the stochasticity of the spectator field during inflation and the requirement that galaxies can form. Finally, we identify several observational predictions of our scenario that should be testable within the next few years.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

Primordial black holes (PBHs) could provide the dark matter but a variety of constraints restrict the possible mass windows to $10^{16} - 10^{17}$g, $10^{20} - 10^{24}$g and $10 - 10^3M_{\odot}$. The last possibility is of special interest in view of the recent detection of black-hole mergers by LIGO. PBHs larger than $10^3 M_{\odot}$ might have important cosmological consequences even if they have only a small fraction of the dark matter density. In particular, they could generate cosmological structures either individually through the "seed" effect or collectively through the "Poisson" effect, thereby alleviating some problems associated with the standard cold dark matter scenario.

We study the formation of primordial black holes (PBHs) in the early Universe during a period of slow reheating after inflation. We demonstrate how the PBH formation mechanism may change even before the end of the matter-dominated phase and calculate the expected PBH mass function. We find that there is a threshold for the variance of the density contrast, $\sigma_c \simeq 0.05$, below which the transition occurs even before reheating, with this having important consequences for the PBH mass function. We also show that there is a maximum cut-off for the PBH mass at around $100\,M_{\odot}$, below which the subdominant radiation bath affects PBH production, making the scenario particularly interesting for the recent LIGO observations of black hole mergers.

Primordial black holes (PBHs) could provide the dark matter in various mass windows below $10^2 M_{\odot}$ and those of $30 M_{\odot}$ might explain the LIGO events. PBHs much larger than this might have important consequences even if they provide only a small fraction of the dark matter. In particular, they could generate cosmological structure either individually through through the `seed' effect or collectively through the `Poisson' effect, thereby alleviating some problems associated with the standard CDM scenario. If the PBHs all have a similar mass and make a small contribution to the dark matter, then the seed effect dominates on small scales, in which case PBHs could seed the supermassive black holes in galactic nuclei or even galaxies themselves. If they have a similar mass and provide the dark matter, the Poisson effect dominates on all scales and the first bound clouds would form earlier than in the usual scenario, with interesting observational consequences. If the PBHs have an extended mass spectrum, which is more likely, they could fulfill all three roles - providing the dark matter, binding the first bound clouds and generating galaxies. In this case, the galactic mass function naturally has the observed form, with the galaxy mass being simply related to the black hole mass. The stochastic gravitational wave background from the PBHs in this scenario would extend continuously from the LIGO frequency to the LISA frequency, offering a potential goal for future surveys.

We completely classify Friedmann-Lema\^{i}tre-Robertson-Walker solutions with spatial curvature $K=0,\pm 1$ and equation of state $p=w\rho$, according to their conformal structure, singularities and trapping horizons. We do not assume any energy conditions and allow $\rho < 0$, thereby going beyond the usual well-known solutions. For each spatial curvature, there is an initial spacelike big-bang singularity for $w>-1/3$ and $\rho>0$, while no big-bang singularity for $w<-1$ and $\rho>0$. For $K=0$ or $-1$, $-1<w<-1/3$ and $\rho>0$, there is an initial null big-bang singularity. For each spatial curvature, there is a final spacelike future big-rip singularity for $w<-1$ and $\rho>0$, with null geodesics being future complete for $-5/3\le w<-1$ but incomplete for $w<-5/3$. For $w=-1/3$, the expansion speed is constant. For $-1<w<-1/3$ and $K=1$, the universe contracts from infinity, then bounces and expands back to infinity. For $K=0$, the past boundary consists of timelike infinity and a regular null hypersurface for $-5/3<w<-1$, while it consists of past timelike and past null infinities for $w\le -5/3$. For $w<-1$ and $K=1$, the spacetime contracts from an initial spacelike past big-rip singularity, then bounces and blows up at a final spacelike future big-rip singularity. For $w<-1$ and $K=-1$, the past boundary consists of a regular null hypersurface. The trapping horizons are timelike, null and spacelike for $w\in (-1,1/3)$, $w\in \{1/3, -1\}$ and $w\in (-\infty,-1)\cup (1/3,\infty)$, respectively. A negative energy density ($\rho <0$) is possible only for $K=-1$. In this case, for $w>-1/3$, the universe contracts from infinity, then bounces and expands to infinity; for $-1<w<-1/3$, it starts from a big-bang singularity and contracts to a big-crunch singularity; for $w<-1$, it expands from a regular null hypersurface and contracts to another regular null hypersurface.

If primordial black holes (PBHs) form directly from inhomogeneities in the early Universe, then the number in the mass range $10^5 -10^{12}M_{\odot}$ is severely constrained by upper limits to the $\mu$ distortion in the cosmic microwave background (CMB). This is because inhomogeneities on these scales will be dissipated by Silk damping in the redshift interval $5\times 10^4\lesssim z\lesssim2\times 10^6$. If the primordial fluctuations on a given mass scale have a Gaussian distribution and PBHs form on the high-$\sigma$ tail, as in the simplest scenarios, then the $\mu$ constraints exclude PBHs in this mass range from playing any interesting cosmological role. Only if the fluctuations are highly non-Gaussian, or form through some mechanism unrelated to the primordial fluctuations, can this conclusion be obviated.

We study production of primordial black holes (PBHs) during an early matter-dominated phase. As a source of perturbations, we consider either the inflaton field with a running spectral index or a spectator field that has a blue spectrum and thus provides a significant contribution to the PBH production at small scales. First, we identify the region of the parameter space where a significant fraction of the observed dark matter can be produced, taking into account all current PBH constraints. Then, we present constraints on the amplitude and spectral index of the spectator field as a function of the reheating temperature. We also derive constraints on the running of the inflaton spectral index, ${\rm d}n/{\rm d}{\rm ln}k \lesssim -0.002$, which are comparable to those from the Planck satellite for a scenario where the spectator field is absent.

In this paper we explore the idea that black holes can persist in a universe that collapses to a big crunch and then bounces into a new phase of expansion. We use a scalar field to model the matter content of such a universe {near the time} of the bounce, and look for solutions that represent a network of black holes within a dynamical cosmology. We find exact solutions to Einstein's constraint equations that provide the geometry of space at the minimum of expansion and that can be used as initial data for the evolution of hyperspherical cosmologies. These solutions illustrate that there exist models in which multiple distinct black holes can persist through a bounce, and allow for concrete computations of quantities such as the black hole filling factor. We then consider solutions in flat cosmologies, as well as in higher-dimensional spaces (with up to nine spatial dimensions). We derive conditions for the black holes to remain distinct (i.e. avoid merging) and hence persist into the new expansion phase. Some potentially interesting consequences of these models are also discussed.

The possibility that the dark matter comprises primordial black holes (PBHs) is considered, with particular emphasis on the currently allowed mass windows at $10^{16}$ - $10^{17}\,$g, $10^{20}$ - $10^{24}\,$g and $1$ - $10^{3}\,M_{\odot}$. The Planck mass relics of smaller evaporating PBHs are also considered. All relevant constraints (lensing, dynamical, large-scale structure and accretion) are reviewed and various effects necessary for a precise calculation of the PBH abundance (non-Gaussianity, non-sphericity, critical collapse and merging) are accounted for. It is difficult to put all the dark matter in PBHs if their mass function is monochromatic but this is still possible if the mass function is extended, as expected in many scenarios. A novel procedure for confronting observational constraints with an extended PBH mass spectrum is therefore introduced. This applies for arbitrary constraints and a wide range of PBH formation models, and allows us to identify which model-independent conclusions can be drawn from constraints over all mass ranges. We focus particularly on PBHs generated by inflation, pointing out which effects in the formation process influence the mapping from the inflationary power spectrum to the PBH mass function. We then apply our scheme to two specific inflationary models in which PBHs provide the dark matter. The possibility that the dark matter is in intermediate-mass PBHs of $1$ - $10^{3}\,M_{\odot}$ is of special interest in view of the recent detection of black-hole mergers by LIGO. The possibility of Planck relics is also intriguing but virtually untestable.

The fraction of the Universe going into primordial black holes (PBHs) with initial mass M_* \approx 5 \times 10^{14} g, such that they are evaporating at the present epoch, is strongly constrained by observations of both the extragalactic and Galactic gamma-ray backgrounds. However, while the dominant contribution to the extragalactic background comes from the time-integrated emission of PBHs with initial mass M_*, the Galactic background is dominated by the instantaneous emission of those with initial mass slightly larger than M_* and current mass below M_*. Also, the instantaneous emission of PBHs smaller than 0.4 M_* mostly comprises secondary particles produced by the decay of directly emitted quark and gluon jets. These points were missed in the earlier analysis by Lehoucq et al. using EGRET data. For a monochromatic PBH mass function, with initial mass (1+\mu) M_* and \mu << 1, the current mass is (3\mu)^{1/3} M_* and the Galactic background constrains the fraction of the Universe going into PBHs as a function of \mu. However, the initial mass function cannot be precisely monochromatic and even a tiny spread of mass around M_* would generate a current low-mass tail of PBHs below M_*. This tail would be the main contributor to the Galactic background, so we consider its form and the associated constraints for a variety of scenarios with both extended and nearly-monochromatic initial mass functions. In particular, we consider a scenario in which the PBHs form from critical collapse and have a mass function which peaks well above M_*. In this case, the largest PBHs could provide the dark matter without the M_* ones exceeding the gamma-ray background limits.

The claim that an overdense (positive curvature) region in the early universe cannot extend beyond some maximum scale and remain part of our universe, first made 40 years ago, has recently been questioned by Kopp et al. Their analysis is elucidating and demonstrates that one cannot constrain the form of primordial density perturbations using this argument. However, the notion of a separate-universe scale still applies and it places an important upper limit on the mass of primordial black holes forming at any epoch. We calculate this scale for equations of state of the form $p = k \rho c^2$ with $-1 <k < \infty$, refining earlier calculations on account of the Kopp et al. criticisms. For $-1/3 < k < \infty$, the scale is always of order the cosmological particle horizon size, with a numerical factor depending on $k$. This confirms the earlier claim that a primordial black hole cannot be much larger than the particle horizon at formation. For $-1 < k< -1/3$, as expected for some periods in the history of the universe, the situation changes radically, in that a sufficiently large positive-curvature region produces a baby universe rather than a black hole. There is still a separate-universe scale but the interpretation of these solutions requires care.

Primordial black holes (PBHs) are of special interest because of the crucial role of quantum effects in their formation and evaporation. This means that they provide a unique probe of the early universe, high-energy physics and quantum gravity. We highlight some recent developments in the subject, including improved limits on the fraction of the Universe going into evaporating PBHs in the mass range $10^{9} - 10^{17}$ g and the possibility of using PBHs to probe a cosmological bounce.

We discuss whether black holes could persist in a universe which recollapses and then bounces into a new expansion phase. Whether the bounce is of classical or quantum gravitational origin, such cosmological models are of great current interest. In particular, we investigate the mass range in which black holes might survive a bounce and ways of differentiating observationally between black holes formed just after and just before the last bounce. We also discuss the consequences of the universe going through a sequence of dimensional changes as it passes through a bounce.

We update the constraints on the fraction of the Universe going into primordial black holes in the mass range 10^9--10^17 g associated with the effects of their evaporations on big bang nucleosynthesis and the extragalactic photon background. We include for the first time all the effects of quark and gluon emission by black holes on these constraints and account for the latest observational developments. We then discuss the other constraints in this mass range and show that these are weaker than the nucleosynthesis and photon background limits, apart from a small range 10^13--10^14 g, where the damping of cosmic microwave background anisotropies dominates. Finally we review the gravitational and astrophysical effects of nonevaporating primordial black holes, updating constraints over the broader mass range 1--10^50 g.

Several authors have claimed that the observable Hawking emission from a microscopic black hole is significantly modified by the formation of a photosphere or chromosphere around the black hole due to QED or QCD interactions between the emitted particles. Analyzing these models we identify a number of physical and geometrical effects which invalidate them. In all cases, we find that the observational signatures of a cosmic or Galactic background of black holes or an individual black hole remain essentially those of the standard Hawking model, with little change to the detection probability.

This review addresses the issue of whether there are physically realistic self-similar solutions in which a primordial black hole is attached to an exact or asymptotically Friedmann model for an equation of state of the form $p=(\gamma-1)\rho c^2$. In the positive pressure case ($1 < \gamma < 2$), there is no such solution when the black hole is attached to an exact Friedmann background via a sonic point. However, it has been claimed that there is a one-parameter family of asymptotically Friedmann black hole solutions providing the ratio of the black hole size to the cosmological horizon size is in a narrow range above some critical value. There are also "universal" black holes in which the black hole has an apparent horizon but no event horizon. It turns out that both these types of solution are only asymptotically {\it quasi}-Friedmann, because they contain a solid angle deficit at large distances, but they are not necessarily excluded observationally. We also consider the possibility of self-similar black hole solutions in a universe dominated by a scalar field. If the field is massless, the situation resembles the stiff fluid case, so any black hole solution is again contrived, although there may still be universal black hole solutions. The situation is less clear if the scalar field is rolling down a potential and therefore massive, as in the quintessence scenario. Although no explicit asymptotically Friedmann black hole solutions of this kind are known, they are not excluded and comparison with the $0 < \gamma < 2/3$ perfect fluid case suggests that they should exist if the black hole is not too large. This implies that a black hole might grow as fast as the cosmological horizon in a quintessence-dominated universe in some circumstances, supporting the proposal that accretion onto primordial black holes may have played a role in the production of the supermassive black holes in galactic nuclei.

Searching for microlensing in M31 using automated superpixel surveys raises a number of difficulties which are not present in more conventional techniques. Here we focus on the problem that the list of microlensing candidates is sensitive to the selection criteria or "cuts" imposed and some subjectivity is involved in this. Weakening the cuts will generate a longer list of microlensing candidates but with a greater fraction of spurious ones; strengthening the cuts will produce a shorter list but may exclude some genuine events. We illustrate this by comparing three analyses of the same data-set obtained from a 3-year observing run on the INT in La Palma. The results of two of these analyses have been already reported: Belokurov et al. (2005) obtained between 3 and 22 candidates, depending on the strength of their cuts, while Calchi Novati et al. (2005) obtained 6 candidates. The third analysis is presented here for the first time and reports 10 microlensing candidates, 7 of which are new. Only two of the candidates are common to all three analyses. In order to understand why these analyses produce different candidate lists, a comparison is made of the cuts used by the three groups...

Motivated by the cosmological wormhole solutions obtained from our recent numerical investigations, we provide a definition of a wormhole which applies to dynamical situations. Our numerical solutions do not have timelike trapping horizons but they are wormholes in the sense that they connect two or more asymptotic regions. Although the null energy condition must be violated for static wormholes, we find that it can still be satisfied in the dynamical context. Two analytic solutions for a cosmological wormhole connecting two Friedmann universes without trapping horizons are presented.

Several authors, most notably Heckler, have claimed that the observable Hawking emission from a microscopic black hole is significantly modified by the formation of a photosphere around the black hole due to QED or QCD interactions between the emitted particles. In this paper we analyze these claims and identify a number of physical and geometrical effects which invalidate these scenarios. We point out two key problems. First, the interacting particles must be causally connected to interact, and this condition is satisfied by only a small fraction of the emitted particles close to the black hole. Second, a scattered particle requires a distance ~ E/m_e^2 for completing each bremsstrahlung interaction, with the consequence that it is improbable for there to be more than one complete bremsstrahlung interaction per particle near the black hole. These two effects have not been included in previous analyses. We conclude that the emitted particles do not interact sufficiently to form a QED photosphere. Similar arguments apply in the QCD case and prevent a QCD photosphere (chromosphere) from developing when the black hole temperature is much greater than Lambda_QCD, the threshold for QCD particle emission. Additional QCD phenomenological arguments rule out the development of a chromosphere around black hole temperatures of order Lambda_QCD. In all cases, the observational signatures of a cosmic or Galactic halo background of primordial black holes or an individual black hole remain essentially those of the standard Hawking model, with little change to the detection probability. We also consider the possibility, as proposed by Belyanin et al. and D. Cline et al., that plasma interactions between the emitted particles form a photosphere, and we conclude that this scenario too is not supported.

We discuss a variety of bremsstrahlung processes associated with charged particles emitted by evaporating black holes. We show that such particles produce a negligible number of bremsstrahlung photons from their scattering off each other, though at low frequencies inner bremsstrahlung photons dominate over the direct Hawking emission of photons. This analysis and the further analysis of the accompanying paper invalidate Heckler's claim that sufficiently hot evaporating black holes form QED photospheres.

Based on the asymptotic analysis of ordinary differential equations, we classify all spherically symmetric self-similar solutions to the Einstein equations which are asymptotically Friedmann at large distances and contain a perfect fluid with equation of state $p=(\gamma -1)\mu$ with $0<\gamma<2/3$. This corresponds to a ``dark energy'' fluid and the Friedmann solution is accelerated in this case due to anti-gravity. This extends the previous analysis of spherically symmetric self-similar solutions for fluids with positive pressure ($\gamma>1$). However, in the latter case there is an additional parameter associated with the weak discontinuity at the sonic point and the solutions are only asymptotically ``quasi-Friedmann'', in the sense that they exhibit an angle deficit at large distances. In the $0<\gamma<2/3$ case, there is no sonic point and there exists a one-parameter family of solutions which are {\it genuinely} asymptotically Friedmann at large distances. We find eight classes of asymptotic behavior: Friedmann or quasi-Friedmann or quasi-static or constant-velocity at large distances, quasi-Friedmann or positive-mass singular or negative-mass singular at small distances, and quasi-Kantowski-Sachs at intermediate distances. The self-similar asymptotically quasi-static and quasi-Kantowski-Sachs solutions are analytically extendible and of great cosmological interest. We also investigate their conformal diagrams. The results of the present analysis are utilized in an accompanying paper to obtain and physically interpret numerical solutions.

We use a combination of numerical and analytical methods, exploiting the equations derived in a preceding paper, to classify all spherically symmetric self-similar solutions which are asymptotically Friedmann at large distances and contain a perfect fluid with equation of state $p=(\gamma -1)\mu$ with $0<\gamma<2/3$. The expansion of the Friedmann universe is accelerated in this case. We find a one-parameter family of self-similar solutions representing a black hole embedded in a Friedmann background. This suggests that, in contrast to the positive pressure case, black holes in a universe with dark energy can grow as fast as the Hubble horizon if they are not too large. There are also self-similar solutions which contain a central naked singularity with negative mass and solutions which represent a Friedmann universe connected to either another Friedmann universe or some other cosmological model. The latter are interpreted as self-similar cosmological white hole or wormhole solutions. The throats of these wormholes are defined as two-dimensional spheres with minimal area on a spacelike hypersurface and they are all non-traversable because of the absence of a past null infinity.

We consider the possible existence of self-similar solutions containing black holes in a Friedmann background with a stiff fluid or a scalar field. We carefully study the relationship between the self-similar equations in these two cases and emphasize the crucial role of the similarity horizon. We show that there is no self-similar black hole solution surrounded by an exact or asymptotically flat Friedmann background containing a massless scalar field. This result also applies for a scalar field with a potential, providing the universe is decelerating. However, if there is a potential and the universe is accelerating (as in the quintessence scenario), the result only applies for an exact Friedmann background. This extends the result previously found in the stiff fluid case and strongly suggests that accretion onto primordial black holes is ineffective even during scalar field domination. It also contradicts recent claims that such black holes can grow appreciably by accreting quintessence. Appreciable growth might be possible with very special matter fields but this requires {\it ad hoc} and probably unphysical conditions.

Self-similar models are important in general relativity and other fundamental theories. In this paper we shall discuss the ``similarity hypothesis'', which asserts that under a variety of physical circumstances solutions of these theories will naturally evolve to a self-similar form. We will find there is good evidence for this in the context of both spatially homogenous and inhomogeneous cosmological models, although in some cases the self-similar model is only an intermediate attractor. There are also a wide variety of situations, including critical pheneomena, in which spherically symmetric models tend towards self-similarity. However, this does not happen in all cases and it is it is important to understand the prerequisites for the conjecture.

Recent developments in the study of primordial black holes (PBHs) are reviewed, with particular emphasis on their formation and evaporation. It is still not clear whether PBHs formed but, if they did, they could provide a unique probe of the early Universe, gravitational collapse, high energy physics and quantum gravity. Indeed their study may place interesting constraints on the physics relevant to these areas even if they never existed.

The POINT-AGAPE collaboration is carrying out a search for gravitational microlensing toward M31 to reveal galactic dark matter in the form of MACHOs (Massive Astrophysical Compact Halo Objects) in the halos of the Milky Way and M31. A high-threshold analysis of 3 years of data yields 6 bright, short--duration microlensing events, which are confronted to a simulation of the observations and the analysis. The observed signal is much larger than expected from self lensing alone and we conclude, at the 95% confidence level, that at least 20% of the halo mass in the direction of M31 must be in the form of MACHOs if their average mass lies in the range 0.5-1 M$_\odot$. This lower bound drops to 8% for MACHOs with masses $\sim 0.01$ M$_\odot$. In addition, we discuss a likely binary microlensing candidate with caustic crossing. Its location, some 32' away from the centre of M31, supports our conclusion that we are detecting a MACHO signal in the direction of M31.

We report the results of our statistical study of the POINT-AGAPE CNe catalogue. The first Sloan $r'$ and $i'$ analysis of the MMRD relationship in M31 is performed using the resulting POINT-AGAPE CN catalogue. Good fits are produced to the MMRD in the two filters. However, we are unable to verify the $t_{15}$ relationship for either Sloan filter. The subsequent analysis of our automated CN detection pipeline has provided us with the most thorough knowledge of the completeness of a CN survey to-date. In addition, the large field of view of the survey has permitted us to probe the outburst rate well into the galactic disk, unlike previous CCD imaging surveys. Using this analysis we are able to probe the CN distribution of M31 and evaluate the global nova rate. Using models of the galactic surface brightness of M31, we show that the observed CN distribution consists of a separate bulge and disk population. We also show that the M31 bulge CN eruption rate per unit $r'$ flux is more than five times greater than that of the disk. We deduce a global M31 CN rate of $65^{+16}_{-15}$ year$^{-1}$, a value much higher than found by previous surveys. Using the global rate, we derive a M31 bulge rate of $38^{+15}_{-12}$ year$^{-1}$ and a disk rate of $27^{+19}_{-15}$ year$^{-1}$. Given our understanding of the completeness and an analysis of other sources of error, we conclude that the true global nova rate of M31 is at least 50% higher than was previously thought and this has consequent implications for the presumed CN rate in the Milky Way. We deduce a Galactic bulge rate of $14^{+6}_{-5}$ year$^{-1}$, a disk rate of $20^{+14}_{-11}$ year$^{-1}$ and a global Galactic rate of $34^{+15}_{-12}$ year$^{-1}$.

We discuss a new class of solutions to the Einstein equations which describe a primordial black hole (PBH) in a flat Friedmann background. Such solutions arise if a Schwarzschild black hole is patched onto a Friedmann background via a transition region. They are possible providing the black hole event horizon is larger than the cosmological apparent horizon. Such solutions have a number of strange features. In particular, one has to define the black hole and cosmological horizons carefully and one then finds that the mass contained within the black hole event horizon decreases when it is larger than the Friedmann cosmological apparent horizon, although its area always increases. These solutions involve two distinct future null infinities and are interpreted as the conversion of a white hole into a black hole. Although such solutions may not form from gravitational collapse in the same way as standard PBHs, there is nothing unphysical about them, since all energy and causality conditions are satisfied. Their conformal diagram is a natural amalgamation of the Kruskal diagram for the extended Schwarzschild solution and the conformal diagram for a black hole in a flat Friedmann background. In this paper, such solutions are obtained numerically for a spherically symmetric universe containing a massless scalar field, but it is likely that they exist for more general matter fields and less symmetric systems.

We provide precise constraints on the size of any black holes forming in the early Universe for a variety of formation scenarios. In particular, we prove that the size of the apparent horizon of a primordial black hole formed by causal processes in a flat Friedmann universe is considerably smaller than the cosmological apparent horizon size for an equation of state $p=k\rho$ ($1/3<k<1$). This also applies for a stiff equation of state ($k=1$) or for a massless scalar field. The apparent horizon of a primordial black hole formed through hydrodynamical processes is also considerably smaller than the cosmological apparent horizon for $0<k\le 1$. We derive an expression for the maximum size which an overdense region can have without being a separate closed universe rather than part of our own. Newtonian argument shows that a black hole smaller than the cosmological horizon can never accrete much.

The evolution of primordial black holes in a flat Friedmann universe with a massless scalar field is investigated in fully general relativistic numerical relativity. A primordial black hole is expected to form with a scale comparable to the cosmological apparent horizon, in which case it may go through an initial phase with significant accretion. However, if it is very close to the cosmological apparent horizon size, the accretion is suppressed due to general relativistic effects. In any case, it soon gets smaller than the cosmological horizon and thereafter it can be approximated as an isolated vacuum solution with decaying mass accretion. In this situation the dynamical and inhomogeneous scalar field is typically equivalent to a perfect fluid with a stiff equation of state $p=\rho$. The black hole mass never increases by more than a factor of two, despite recent claims that primordial black holes might grow substantially through accreting quintessence. It is found that the gravitational memory scenario, proposed for primordial black holes in Brans-Dicke and scalar-tensor theories of gravity, is highly unphysical.

The POINT-AGAPE collaboration surveyed M31 with the primary goal of optical detection of microlensing events, yet its data catalogue is also a prime source of lightcurves of variable and transient objects, including classical novae (CNe). A reliable means of identification, combined with a thorough survey of the variable objects in M31, provides an excellent opportunity to locate and study an entire galactic population of CNe. This paper presents a set of 440 neural networks, working in 44 committees, designed specifically to identify fast CNe. The networks are developed using training sets consisting of simulated novae and POINT-AGAPE lightcurves, in a novel variation on K-fold cross-validation. They use the binned, normalised power spectra of the lightcurves as input units. The networks successfully identify 9 of the 13 previously identified M31 CNe within their optimal working range (and 11 out of 13 if the network error bars are taken into account). They provide a catalogue of 19 new candidate fast CNe, of which 4 are strongly favoured.

Recent developments in the study of primordial black holes (PBHs) will be reviewed, with particular emphasis on their formation and evaporation. PBHs could provide a unique probe of the early Universe, gravitational collapse, high energy physics and quantum gravity. Indeed their study may place interesting constraints on the physics relevant to these areas even if they never formed. In the "early Universe" context, particularly useful constraints can be placed on inflationary scenarios, especially if evaporating PBHs leave stable Planck-mass relicts. In the "gravitational collapse" context, the existence of PBHs could provide a unique test of the sort of critical phenomena discovered in recent numerical calculations. In the "high energy physics" context, information may come from gamma-ray bursts (if a subset of these are generated by PBH explosions) or from cosmic rays (if some of these derive from evaporating PBHs). In the "quantum gravity" context, the formation and evaporation of small black holes could lead to observable signatures in cosmic ray events and accelerator experiments, providing there are extra dimensions and providing the quantum gravity scale is around a TeV.

An automated search is carried out for microlensing events using a catalogue of 44554 variable superpixel lightcurves derived from our three-year monitoring program of M31. Each step of our candidate selection is objective and reproducible by a computer. Our search is unrestricted, in the sense that it has no explicit timescale cut. So, it must overcome the awkward problem of distinguishing long-timescale microlensing events from long-period stellar variables. The basis of the selection algorithm is the fitting of the superpixel lightcurves to two different theoretical models, using variable star and blended microlensing templates. Only if microlensing is preferred is an event retained as a possible candidate. Further cuts are made with regard to (i) sampling, (ii) goodness of fit of the peak to a Paczynski curve, (iii) consistency of the microlensing hypothesis with the absence of a resolved source, (iv) achromaticity, (v) position in the colour-magnitude diagram and (vi) signal-to-noise ratio. Our results are reported in terms of first-level candidates, which are the most trustworthy, and second-level candidates, which are possible microlensing but have lower signal-to-noise and are more questionable. The pipeline leaves just 3 first-level candidates, all of which have very short full-width half-maximum timescale (<5 days) and 3 second-level candidates, which have timescales of 31, 36 and 51 days respectively. We also show 16 third-level lightcurves, as an illustration of the events that just fail the threshold for designation as microlensing candidates. They are almost certainly mainly variable stars. Two of the 3 first-level candidates correspond to known events (PA 00-S3 and PA 00-S4) already reported by the POINT-AGAPE project. The remaining first-level candidate is new.

The POINT-AGAPE survey is an optical search for gravitational microlensing events towards the Andromeda Galaxy (M31). As well as microlensing, the survey is sensitive to many different classes of variable stars and transients. Here we describe the automated detection and selection pipeline used to identify M31 classical novae (CNe) and we present the resulting catalogue of 20 CN candidates observed over three seasons. CNe are observed both in the bulge region as well as over a wide area of the M31 disk. Nine of the CNe are caught during the final rise phase and all are well sampled in at least two colours. The excellent light-curve coverage has allowed us to detect and classify CNe over a wide range of speed class, from very fast to very slow. Among the light-curves is a moderately fast CN exhibiting entry into a deep transition minimum, followed by its final decline. We have also observed in detail a very slow CN which faded by only 0.01 mag day$^{-1}$ over a 150 day period. We detect other interesting variable objects, including one of the longest period and most luminous Mira variables. The CN catalogue constitutes a uniquely well-sampled and objectively-selected data set with which to study the statistical properties of classical novae in M31, such as the global nova rate, the reliability of novae as standard-candle distance indicators and the dependence of the nova population on stellar environment. The findings of this statistical study will be reported in a follow-up paper.

The POINT-AGAPE collaboration has been monitoring M31 for three seasons with the Wide Field Camera on the Isaac Newton Telescope. In each season, data are taken for one hour per night for roughly sixty nights during the six months that M31 is visible. The two fields of view straddle the central bulge, northwards and southwards. We have calculated the locations, periods and amplitudes of 35414 variable stars in M31 as a by-product of our microlensing search. The variables are classified according to their period and amplitude of variation. They are classified into population I and II Cepheids, Miras and semi-regular long-period variables. The population I Cepheids are associated with the spiral arms, while the central concentration of the Miras and long-period variables varies noticeably, the stars with brighter (and shorter) variations being much more centrally concentrated. A crucial role in the microlensing experiment is played by the asymmetry signal. It was initially assumed that the variable stars would not be a serious problem as their distributions would be symmetric. We demonstrate that this assumption is not correct. We find that differential extinction associated with the dust lanes causes the variable star distributions to be asymmetric. The size and direction of the asymmetry of the variable stars is measured as a function of period and amplitude of variation. The implications of this discovery for the successful completion of the microlensing experiments towards M31 are discussed. (Abridged)

Recent developments in the study of primordial black holes (PBHs) will be reviewed, with particular emphasis on their formation and evaporation. PBHs could provide a unique probe of the early Universe, gravitational collapse, high energy physics and quantum gravity. Indeed their study may place interesting constraints on the physics relevant to these areas even if they never formed.

The lightcurve of PA-99-N2, one of the recently announced microlensing candidates towards M31, shows small deviations from the standard Paczynski form. We explore a number of possible explanations, including correlations with the seeing, the parallax effect and a binary lens. We find that the observations are consistent with an unresolved RGB or AGB star in M31 being microlensed by a binary lens. We find that the best fit binary lens mass ratio is about one hundredth, which is one of most extreme values found for a binary lens so far. If both the source and lens lie in the M31 disk, then the standard M31 model predicts the probable mass range of the system to be 0.02-3.6 solar masses (95 % confidence limit). In this scenario, the mass of the secondary component is therefore likely to be below the hydrogen-burning limit. On the other hand, if a compact halo object in M31 is lensing a disk or spheroid source, then the total lens mass is likely to lie between 0.09-32 solar masses, which is consistent with the primary being a stellar remnant and the secondary a low mass star or brown dwarf. The optical depth (or alternatively the differential rate) along the line of sight toward the event indicates that a halo lens is more likely than a stellar lens provided that dark compact objects comprise no less than 15 per cent (or 5 per cent) of haloes.

The POINT-AGAPE collaboration is currently searching for massive compact halo objects (MACHOs) towards the Andromeda galaxy (M31). The survey aims to exploit the high inclination of the M31 disk, which causes an asymmetry in the spatial distribution of M31 MACHOs. Here, we investigate the effects of halo velocity anisotropy and flattening on the asymmetry signal using simple halo models. For a spherically symmetric and isotropic halo, we find that the underlying pixel-lensing rate in far-disk M31 MACHOs is more than 5 times the rate of near-disk events. We find that the asymmetry is increased further by about 30% if the MACHOs occupy radial orbits rather than tangential orbits, but is substantially reduced if the MACHOs lie in a flattened halo. However, even for haloes with a minor-to-major axis ratio q = 0.3, the numbers of M31 MACHOs in the far-side outnumber those in the near-side by a factor of ~2. We show that, if positional information is exploited in addition to number counts, then the number of candidate events required to confirm asymmetry for a range of flattened and anisotropic halo models is achievable, even with significant contamination by variable stars and foreground microlensing events. For pixel-lensing surveys which probe a representative portion of the M31 disk, a sample of around 50 candidates is likely to be sufficient to detect asymmetry within spherical haloes, even if half the sample is contaminated, or to detect asymmetry in haloes as flat as q = 0.3 provided less than a third of the sample comprises contaminants. We also argue that, provided its mass-to-light ratio is less than 100, the recently observed stellar stream around M31 is not problematic for the detection of asymmetry. (Abstract slightly abridged.)

We have carried out a survey of the Andromeda galaxy for unresolved microlensing (pixel lensing). We present a subset of four short timescale, high signal-to-noise microlensing candidates found by imposing severe selection criteria: the source flux variation exceeds the flux of an R=21 magnitude star and the full width at half maximum timescale is less than 25 days. Remarkably, in three out of four cases, we have been able to measure or strongly constrain the Einstein crossing time of the event. One event, which lies projected on the M31 bulge, is almost certainly due to a stellar lens in the bulge of M31. The other three candidates can be explained either by stars in M31 and M32 or by MACHOs.

In cosmological models with a varying gravitational constant, it is not clear whether primordial black holes preserve the value of $G$ at their formation epoch. We investigate this question by using the Tolman-Bondi model to study the evolution of a background scalar field when a black hole forms from the collapse of dust in a flat Friedmann universe. Providing the back reaction of the scalar field on the metric can be neglected, we find that the value of the scalar field at the event horizon very quickly assumes the background cosmological value. This suggests that there is very little gravitational memory.

We report the discovery of a microlensing candidate projected 2'54" from the center of M32, on the side closest to M31. The blue color (R-I= 0.00 +/- 0.14) of the source argues strongly that it lies in the disk of M31, while the proximity of the line of sight to M32 implies that this galaxy is the most likely host of the lens. If this interpretation is correct, it would confirm previous arguments that M32 lies in front of M31. We estimate that of order one such event or less should be present in the POINT-AGAPE data base. If more events are discovered in this direction in a dedicated experiment, they could be used to measure the mass function of M32 up to an unknown scale factor. By combining microlensing observations of a binary-lens event with a measurement of the M31-M32 relative proper motion using the astrometric satellites SIM or GAIA, it will be possible to measure the physical separation of M31 and M32, the last of the six phase-space coordinates needed to assign M32 an orbit.

We discuss recent developments in the study of primordial black holes, focussing particularly on their formation and quantum evaporation. Such studies can place important constraints on models of the early Universe. An especially interesting development has been the realization that such constraints may be severely modified if the value of the gravitational "constant" G varies with cosmological epoch, a possibility which arises in many scenarios for the early Universe. The nature of the modification depends upon whether the value of $G$ near a black hole maintains the value it had at its formation epoch (corresponding to gravitational memory) or whether it tracks the background cosmological value. This is still uncertain but we discuss various approaches which might help to resolve the issue.

Cosmological nucleosynthesis calculations imply that many of the baryons in the Universe must be dark. We discuss the likelihood that some of these dark baryons may reside in galaxies as Massive Compact Halo Objects (MACHOs), the remnants of a first generation of pregalactic or protogalactic stars. Various candidates have been proposed for such remnants and we review the many types of observations which can be used to detect or exclude them. Claims to have found positive evidence for some of the candidates have generally turned out to be spurious or questionable, so the status of the MACHO scenario remains controversial. However, it would be premature to reject MACHOs altogether and further observations are likely to resolve the issue soon.

It has been pointed out that the value of the gravitational ``constant'' in the early universe may be different from that at present. In that case, it was conjectured that primordial black holes may ``remember'' the value of the gravitational constant in the early universe. The present analysis shows that this is not the case, at least in certain contexts.

We report the discovery of a short-duration microlensing candidate in the northern field of the POINT-AGAPE pixel lensing survey towards M31. The full-width half-maximum timescale is very short, just 1.8 days. Almost certainly, the source star has been identified on Hubble Space Telescope archival images, allowing us to infer an Einstein crossing time of 10.4 days, a maximum magnification of about 18, and a lens-source proper motion greater than 0.3 microarcsec/day. The event lies projected at 8' from the center of M31, which is beyond the bulk of the stellar lens population. The lens is likely to reside in one of three locations. It may be a star in the M31 disk, or a massive compact halo object (Macho) in either M31 or the Milky Way. The most probable mass is 0.06 solar masses for an M31 Macho, 0.02 solar masses for a Milky Way Macho and 0.2 solar masses for an M31 stellar lens. Whilst the stellar interpretation is plausible, the Macho interpretation is the most probable for halo fractions above 20%.

Cosmological nucleosynthesis calculations imply that many of the baryons in the Universe must be dark. We discuss the likelihood that some of these dark baryons may reside in the discs or halos of galaxies. If they were in the form of compact objects, they would then be natural MACHO candidates, in which case they are likely to be the remnants of a first generation of pregalactic or protogalactic Population III stars. Various candidates have been proposed for such remnants - brown dwarfs, red dwarfs, white dwarfs, neutron stars or black holes - and we review the many types of observations (including microlensing searches) which can be used to constrain or exclude them.

Cosmological nucleosynthesis calculations imply that there should be both non-baryonic and baryonic dark matter. Recent data suggest that some of the non-baryonic dark matter must be "hot" (i.e. massive neutrinos) and there may also be evidence for "cold" dark matter (i.e. WIMPs). If the baryonic dark matter resides in galactic halos, it is likely to be in the form of compact objects (i.e. MACHOs) and these would probably be the remnants of a first generation of pregalactic or protogalactic Population III stars. Many candidates have been proposed - brown dwarfs, red dwarfs, white dwarfs or black holes - and at various times each of these has been in vogue. We review the many types of observations which can be used to constrain or exclude both baryonic and non-baryonic dark matter candidates.

We discuss the various ways in which primordial black holes may have formed in the early Universe and how the effects of such black holes can be used to place constraints on cosmological models. We show that such constraints may be severely modified if the value of the gravitational "constant" G varies with cosmological epoch, a possibility which arises in many scenarios for the early Universe. The nature of the modification depends upon whether the value of G near a black hole maintains the value it had at its formation epoch (corresponding to gravitational memory) or whether it tracks the background cosmological value. This is still uncertain but we discuss various approaches which might help to resolve the issue.

POINT-AGAPE is an Anglo-French collaboration which is employing the Isaac Newton Telescope (INT) to conduct a pixel-lensing survey towards M31. In this paper we investigate what we can learn from pixel-lensing observables about the MACHO mass and fractional contribution in M31 and the Galaxy for the case of spherically-symmetric near-isothermal haloes. We employ detailed pixel-lensing simulations which include many of the factors which affect the observables. For a maximum MACHO halo we predict an event rate in V of up to 100 per season for M31 and 40 per season for the Galaxy. However, the Einstein radius crossing time is generally not measurable and the observed full-width half-maximum duration provides only a weak tracer of lens mass. Nonetheless, we find that the near-far asymmetry in the spatial distribution of M31 MACHOs provides significant information on their mass and density contribution. We present a likelihood estimator for measuring the fractional contribution and mass of both M31 and Galaxy MACHOs which permits an unbiased determination to be made of MACHO parameters, even from data-sets strongly contaminated by variable stars. If M31 does not have a significant population of MACHOs in the mass range 0.001-1 Solar masses strong limits will result from the first season of INT observations. Simulations based on currently favoured density and mass values indicate that, after three seasons, the M31 MACHO parameters should be constrained to within a factor four uncertainty in halo fraction and an order of magnitude uncertainty in mass (90% confidence). Interesting constraints on Galaxy MACHOs may also be possible. For a campaign lasting ten years, comparable to the lifetime of current LMC surveys, reliable estimates of MACHO parameters in both galaxies should be possible. (Abridged)

We classify all spherically symmetric perfect fluid solutions of Einstein's equations with equation of state p/mu=a which are self-similar in the sense that all dimensionless variables depend only upon z=r/t. For a given value of a, such solutions are described by two parameters and they can be classified in terms of their behaviour at large and small distances from the origin; this usually corresponds to large and small values of z but (due to a coordinate anomaly) it may also correspond to finite z. We base our analysis on the demonstration that all similarity solutions must be asymptotic to solutions which depend on either powers of z or powers of lnz. We show that there are only three similarity solutions which have an exact power-law dependence on z: the flat Friedmann solution, a static solution and a Kantowski-Sachs solution (although the latter is probably only physical for a<0). For a>1/5, there are also two families of solutions which are asymptotically (but not exactly) Minkowski: the first is asymptotically Minkowski as z tends to infinity and is described by one parameter; the second is asymptotically Minkowski at a finite value of z and is described by two parameters. A complete analysis of the dust solutions is given, since these can be written down explicitly and elucidate the link between the z>0 and z<0 solutions. Solutions with pressure are then discussed in detail; these share many of the characteristics of the dust solutions but they also exhibit new features.

If compact baryonic objects contribute significantly to the dark matter in our Galaxy, their mass function will present vital clues for galaxy formation theories and star formation processes in the early Universe. Here we discuss what one might expect to learn about the mass function of Galactic dark matter from microlensing and from direct searches in the infrared and optical wavebands. Current microlensing results from the \eros\/ collaboration already constrain halo mass functions which extend below $10^{-4}~\sm$, whilst recent \hst\/ observations place strong constraints on disc and halo dark matter mass functions extending above $0.1~\sm$. Infrared observations should either detect or constrain objects larger than $0.01~\sm$ in the near future. Objects below $0.01~\sm$ should be detectable through microlensing, although the prospects of determining their mass function depend critically on a number of factors.

Blue primordial power spectra with a spectral index $n>1$ can lead to a significant production of primordial black holes in the very early Universe. The evaporation of these objects leads to a number of observational consequences and a model independent upper limit of $n \approx 1.4$. In some cases this limit is strengthened to $n=1.3$. Such limits may be employed to define the boundary to the region of parameter space consistent with generalized inflationary predictions. [To appear in Proceedings of the CASE WESTERN CMB WORKSHOP, April 22-24 1994. Figures available on request from J.H.Gilbert@qmw.ac.uk]

Blue primordial power spectra have spectral index $n>1$ and arise naturally in the recently proposed hybrid inflationary scenario. An observational upper limit on {\em n} is derived by normalizing the spectrum at the quadrupole scale and considering the possible overproduction of Planck mass relics formed in the final stage of primordial black hole evaporation. In the inflationary Universe with the maximum reheating temperature compatible with the observed quadrupole anisotropy, the upper limit is $n=1.4$, but it is slightly weaker for lower reheat temperatures. This limit applies over 57 decades of mass and is therefore insensitive to cosmic variance and any gravitational wave contribution to the quadrupole anisotropy. It is also independent of the dark matter content of the Universe and therefore the bias parameter. In some circumstances, there may be an extended dust-like phase between the end of inflation and reheating. In this case, primordial black holes form more abundantly and the upper limit is $n=1.3$.

Brown dwarfs, stars with insufficient mass to burn hydrogen, could contribute to the dark matter in the Galactic disk, galactic halos or even a background critical density. We consider the detectability of such brown dwarfs in various scenarios, extending previous work by allowing for the possibility that they may have an extended mass spectrum or be clumped into dark clusters. We investigate the constraints placed on such scenarios by the \iras survey. Whilst an extrapolation of the mass function of visible disk stars makes it unlikely that brown dwarfs comprise all of the proposed disk dark matter, \iras does not exclude brown dwarfs providing the dark matter in our own halo or a cosmological background. Neither does it improve on existing dynamical constraints on the mass and radius of brown dwarf clusters in our halo. Future satellites such as \iso and \sirtf will either detect brown dwarfs or brown dwarf clusters or else severely constrain their contribution to the dark matter.