The LIGO/Virgo collaboration has by now observed or constrained the gravitational merger rates of different classes of compact objects. We consider the possibility that the bulk of these mergers are primordial black hole (PBH) mergers of PBHs formed during the QCD epoch making up the entirety of the dark matter. Having shown in a companion paper that mergers due to the initial binary population formed in the early Universe are likely negligible, we compute current merger rates in PBH clusters in which the typical PBH resides. We consider two scenarios: (i) the PBH mass function dictated by the QCD equation of state and (ii) the PBH mass function dictated by the existence of a peak in the inflationary perturbation spectrum. In the first scenario, which is essentially parameter free, we reproduce very well the merger rates for heavy BHs, the merger rate of mass-asymmetric BHs such as GW190814, a recently discovered merger of a 23$M_{\odot}$ black hole with a $2.6M_{\odot}$ object, and can naturally explain why LIGO/Virgo has not yet observed mergers of two light PBHs from the dominant $\sim 1\,M_{\odot}$ PBH population. In the second scenario, which has some parameter freedom, we match well the observed rate of heavy PBHs, but can currently not explain the rate for mass-asymmetric events. In either case it is tantalizing that in both scenarios PBH merger rates made with a minimum of assumptions match most LIGO/Virgo observed rates very well.

The LIGO/Virgo collaboration have by now detected the mergers of ten black hole binaries via the emission of gravitational radiation. The hypothesis that these black holes have formed during the cosmic QCD epoch and make up all of the cosmic dark matter, has been rejected by many authors reasoning that, among other constraints, primordial black hole (PBH) dark matter would lead to orders of magnitude larger merger rates than observed. We revisit the calculation of the present PBH merger rate. Solar mass PBHs form clusters at fairly high redshifts, which evaporate at lower redshifts. We consider in detail the evolution of binary properties in such clusters due to three-body interactions between the two PBH binary members and a third by-passing PBH, for the first time, by full numerical integration. A Monte-Carlo analysis shows that formerly predicted merger rates are reduced by orders of magnitude due to such interactions. The natural prediction of PBH dark matter formed during the QCD epoch yields a pronounced peak around $1M_{\odot}$ with a small mass fraction of PBHs on a shoulder around $30M_{\odot}$, dictated by the well-determined equation of state during the QCD epoch. We employ this fact to make a tentative prediction of the merger rate of $\sim 30M_{\odot}$ PBH binaries, and find it very close to that determined by LIGO/Virgo. Furthermore we show that current LIGO/Virgo limits on the existence of $\sim M_{\odot}$ binaries do not exclude QCD PBHs to make up all of the cosmic dark matter. Neither do constraints on QCD PBHs from the stochastic gravitational background, pre-recombination accretion, or dwarf galaxies pose a problem. Microlensing constraints on QCD PBHs should be re-investigated. We caution, however, in this numerically challenging problem some possibly relevant effects could not be treated.

The standard cosmological model determined from the accurate cosmic microwave background measurements made by the Planck satellite implies a value of the Hubble constant $H_0$ that is $4.4$ standard deviations lower than the one determined from Type Ia supernovae. The Planck best fit model also predicts higher values of the matter density fraction $\Omega_m$ and clustering amplitude $S_8$ compared to those obtained from the Dark Energy Survey Year 1 data. Here we show that accounting for the enhanced recombination rate due to additional small-scale inhomogeneities in the baryon density may solve both the $H_0$ and the $S_8-\Omega_m$ tensions. The additional baryon inhomogeneities can be induced by primordial magnetic fields present in the plasma prior to recombination. The required field strength to solve the Hubble tension is just what is needed to explain the existence of galactic, cluster, and extragalactic magnetic fields without relying on dynamo amplification. Our results show clear evidence for this effect and motivate further detailed studies of primordial magnetic fields, setting several well-defined targets for future observations.

Primordial magnetic fields (PMFs), being present before the epoch of cosmic recombination, induce small-scale baryonic density fluctuations. These inhomogeneities lead to an inhomogeneous recombination process that alters the peaks and heights of the large-scale anisotropies of the cosmic microwave backround (CMB) radiation. Utilizing numerical compressible MHD calculations and a Monte Carlo Markov chain analysis, which compares calculated CMB anisotropies with those observed by the \textit{WMAP} and \textit{Planck} satellites, we derive limits on the magnitude of putative PMFs. We find that the {\it total} remaining present day field, integrated over all scales, cannot exceed 47 pG for scale-invariant PMFs and 8.9 pG for PMFs with a violet Batchelor spectrum at 95\% confidence level. These limits are more than one order of magnitude more stringent than any prior stated limits on PMFs from the CMB which have not accounted for this effect.

Models in which the dark matter is produced at extremely low rates from the annihilation of Standard Model particles in the early Universe allow us to explain the current dark matter relic density while easily evading the traditional experimental constraints. In scenarios where the dark matter interacts with the Standard Model via a new physics mediator, the early Universe dynamics of the dark sector can be particularly complex, as the dark matter and the mediator could be in thermal and chemical equilibrium with each other. This equilibration takes place via number-changing processes such as double Compton scattering and bremsstrahlung, whose amplitudes are cumbersome to calculate. In this paper, we show that in large regions of the parameter space, these equilibration mechanisms do not significantly affect the final dark matter relic density. In particular, for a model with a light dark photon mediator, the relic density can be reasonably estimated by considering that the dark matter is solely produced through the annihilation of Standard Model particles. This result considerably simplifies the treatment of a large class of dark matter theories, facilitating in particular the superimposition of the relic density constraints on the current and future experimental bounds.

Any neutral boson such as a dark photon or dark Higgs that is part of a non-standard sector of particles can mix with its standard model counterpart. When very weakly mixed with the Standard Model, these particles are produced in the early Universe via the freeze-in mechanism and subsequently decay back to standard model particles. In this work, we place constraints on such mediator decays by considering bounds from Big Bang nucleosynthesis and the cosmic microwave background radiation. We find both nucleosynthesis and CMB can constrain dark photons with a kinetic mixing parameter between log {\epsilon} ~ -10 to -17 for masses between 1 MeV and 100 GeV. Similarly, the dark Higgs mixing angle {\epsilon} with the Standard Model Higgs is constrained between log {\epsilon} ~ -6 to -15. Dramatic improvement on the bounds from CMB spectral distortions can be achieved with proposed experiments such as PIXIE.

We extend our recent derivation of the time evolution equations for the energy content of magnetic fields and turbulent motions for incompressible, homogeneous, and isotropic turbulence to include the case of nonvanishing helicity. These equations are subsequently numerically integrated in order to predict the present day primordial magnetic field strength and correlation length, depending on its initial helicity and magnetic energy density. We find that all prior analytic predictions for helical magnetic fields, such as the epoch when they become maximally helical and their subsequent growth of correlation length L ~ a^{1/3} and decrease of magnetic field strength B ~ a^{-1/3} with scale factor a, are well confirmed by the simulations. An initially fully helical primordial magnetic field is a factor 4 10^4 stronger at the present epoch then its nonhelical counterpart when generated during the electroweak epoch.

We present a derivation of the time evolution equations for the energy content of nonhelical magnetic fields and the accompanying turbulent flows from first principles of incompressible magnetohydrodynamics in the general framework of homogeneous and isotropic turbulence. This is then applied to the early Universe, i.e., the evolution of primordial magnetic fields. Numerically integrating the equations, we find that most of the energy is concentrated at an integral wavenumber scale k_I where the turbulence turn over time equals the Hubble time. At larger length scales L, i.e., smaller wavenumbers q = 2 \pi / L << k_I, independent of the assumed turbulent flow power spectrum, mode-mode coupling tends to develop a small q magnetic field tail with a Batchelor spectrum proportional to the fourth inverse power of L and therefore a scaling for the magnetic field of B ~ L^(-5/2).

It is shown that small-scale magnetic fields present before recombination induce baryonic density inhomogeneities of appreciable magnitude. The presence of such inhomogeneities changes the ionization history of the Universe, which in turn decreases the angular scale of the Doppler peaks and increases Silk damping by photon diffusion. This unique signature could be used to (dis)prove the existence of primordial magnetic fields of strength as small as B~10^(-11) Gauss by upcoming cosmic microwave background observations.

Cosmic magnetic fields may be generated during early cosmic phase transition, such as the QCD- or electroweak- transitions. The magnitude of the remainder of such fields at the present epoch crucially depends on the exponent $n$ of their (initially super-Hubble) large-scale tail, i.e. B(lambda) ~ lambda^(-n). It has been claimed that causality requires n=5/2, contrary to much earlier claims of n=3/2. Here we analyze this question in detail. First, we note that contrary to current belief, the large-scale magnetic field tail is not established at the phase transition itself, but rather continuoulsy evolves up to the present epoch. Neglecting turbulent flows we find n=7/2, i.e. very strongly suppressed large-scale fields. However, in the inevitable presence of turbulent flows we find that the large-scale magnetic field tail has sufficient time to evolve to that of the fluid turbulence. For white noise fluid turbulence this yields n=3/2 up to a certain scale and n=5/2 beyond for the magnetic field spectrum. This picture is also not changed when primordial viscosity and fluid flow dissipation is taken into account. Appreciable primordial magnetic fields originating from cosmic phase transitions seem thus possible.

In the pre-reheating era, following cosmic inflation and preceding radiation domination, the energy density may be dominated by an oscillating massive scalar condensate, such as is the case for quadratic chaotic inflation. We have found in a previous paper that during this period, a wide range of sub-Hubble scale perturbations are subject to a preheating instability, leading to the growth of density perturbations ultimately collapsing to form non-linear structures. We compute here the gravitational wave signal due to these structures in the linear limit and present estimates for emission in the non-linear limit due to various effects: the collapse of halos, the tidal interactions, the evaporation during the conversion of the inflaton condensate into radiation and finally the ensuing turbulent cascades. The gravitational wave signal could be rather large and potentially testable by future detectors.

After cosmic inflation and before the transition to radiation domination, the cosmic energy density may have been dominated during an extended period by an oscillating massive scalar condensate. We show that during this period, sub-Hubble scale perturbations are subject to a metric preheating instability in the narrow resonance regime. We analyze in detail both, quadratic and quartic potentials. The instability leads to the growth of density perturbations which in many cases become non-linear already before the beginning of a radiation dominated Universe. This is particularly the case when requiring a phenomenologically preferred low reheat temperature. These early structures may lead to the emission of gravitational waves and the production of primordial black holes. Furthermore, it is not clear if they could modify the prediction of linear curvature perturbations on very large scales.

We propose an alternate, calculable mechanism of dark matter genesis, "thermal freeze-in," involving a Feebly Interacting Massive Particle (FIMP) interacting so feebly with the thermal bath that it never attains thermal equilibrium. As with the conventional "thermal freeze-out" production mechanism, the relic abundance reflects a combination of initial thermal distributions together with particle masses and couplings that can be measured in the laboratory or astrophysically. The freeze-in yield is IR dominated by low temperatures near the FIMP mass and is independent of unknown UV physics, such as the reheat temperature after inflation. Moduli and modulinos of string theory compactifications that receive mass from weak-scale supersymmetry breaking provide implementations of the freeze-in mechanism, as do models that employ Dirac neutrino masses or GUT-scale-suppressed interactions. Experimental signals of freeze-in and FIMPs can be spectacular, including the production of new metastable coloured or charged particles at the LHC as well as the alteration of big bang nucleosynthesis.

We re-consider the gravitino as dark matter in the framework of the Constrained MSSM. We include several recently suggested improvements on: (i) the thermal production of gravitinos, (ii) the calculation of the hadronic spectrum from NLSP decay and (iii) the BBN calculation including stau bound-state effects. In most cases we find an upper bound on the reheating temperature $\treh\lsim {\rm a few}\times 10^7\gev$ from over-production of $^6{\rm Li}$ from bound state effects. We also find an upper limit on the stau lifetime of $3\times10^4 \sec$, which is nearly an order of magnitude larger than the simple limit $5\times 10^3 \sec$ often used to avoid the effect of bound-state catalysis. The bound on $\treh$ is relaxed to $\lsim {\rm a few}\times 10^8\gev$ when we use a more conservative bound on ${^6}{\rm Li}/{^7}{\rm Li}$, in which case a new region at small stau mass at $\sim 100\gev$ and much longer lifetimes opens up. Such a low stau mass region can be easily tested at the LHC.

We review how our current understanding of the light element synthesis during the Big Bang Nucleosynthesis era may help shed light on the identity of particle dark matter.

The Big Bang nucleosynthesis (BBN) in the presence of charged massive particles (CHAMPs) is studied in detail. All currently known effects due to the existence of bound states between CHAMPs and nuclei, including possible late-time destruction of Li6 and Li7 are included. The study sets conservative bounds on CHAMP abundances in the decay time range 3x10^2 sec - 10^12 sec. It is stressed that the production of Li6 at early times T ~ 10keV is overestimated by a factor ~ 10 when the approximation of the Saha equation for the He4 bound state fraction is utilised. To obtain conservative limits on the abundance of CHAMPs, a Monte-Carlo analysis with ~ 3x10^6 independent BBN runs, varying reaction rates of nineteen different reactions, is performed (see attached erratum, however). The analysis yields the surprising result that except for small areas in the particle parameter space conservative constraints on the abundance of decaying charged particles are currently very close to those of neutral particles. It is shown that, in case a number of heretofore unconsidered reactions may be determined reliably in future, it is conceivable that the limit on CHAMPs in the early Universe could be tightened by orders of magnitude. An ERRATUM gives limits on primordial CHAMP densities when the by Ref. Kamimura et al. recently more accurately determined CHAMP reaction rates are employed.

Supersymmetric extensions of the standard model of particle physics assuming the gravitino to be the lightest supersymmetric particle (LSP), and with the next-to-LSP decaying to the gravitino during Big Bang nucleosynthesis, are analyzed. Particular emphasis is laid on their potential to solve the "Li7 problem", an apparent factor 2-4 overproduction of Li7 in standard Big Bang nucleosynthesis (BBN), their production of cosmologically important amounts of Li6, as well as the resulting gravitino dark matter densities in these models. The study includes several improvements compared to prior studies. Heavy gravitinos in the constrained minimal supersymmetric standard model (CMMSM) are reanalyzed, whereas light gravitinos in gauge-mediated supersymmetry breaking scenarios (GMSB) are studied for the first time. It is confirmed that decays of NLSP staus to heavy gravitinos, while producing all the dark matter, may at the same time resolve the Li7 problem. For NLSP decay times ~ 1000 sec, such scenarios also lead to cosmologically important Li6 (and possibly Be9) abundances. However, as such scenarios require heavy > 1 TeV staus they are likely not testable at the LHC. It is found that decays of NLSP staus to light gravitinos may lead to significant Li6 (and Be9) abundances, whereas NLSP neutralinos decaying into light gravitinos may solve the Li7 problem. Though both scenarios are testable at the LHC they may not lead to the production of the bulk of the dark matter. A section of the paper outlines particle properties required to significantly reduce the Li7 abundance, and/or enhance the Li6 (and possibly Be9) abundances, by the decay of an arbitrary relic particle.

Standard Big Bang nucleosynthesis at baryonic density as inferred by WMAP implies a primordial Li7 abundance factor of two to three larger than that inferred by observations of low--metallicity halo stars. Recent observations of Li6 in halo stars suggest a pre--galactic origin of this isotope, which is exceedingly difficult to explain by putative high redshift cosmic ray populations. We investigate if one or both of these lithium problems may be solved by late decaying relic particles in the Constrained Minimal Supersymmetric Standard Model (CMSSM) coupled to gravity. Assuming that the gravitino is the lightest supersymmetric particle (LSP) we find that in large parts of the CMSSM parameter space either of these problems can be potentially solved. In particular, we also find solutions where both lithium problems may be solved simultaneously. These entail the hadronic decay of relic 1 TeV staus into 100 GeV gravitinos at 1000 sec after the Big Bang, as proposed by one of us before. Moreover, the produced gravitinos naturally contribute a large fraction, or all, to the by WMAP required dark matter density. A prediction of this model is the dark matter to be lukewarm with free--streaming lengths of the order of a 4 keV early freezing--out relic particle. Such lukewarm gravitino dark matter may have beneficial effects to the formation of galactic structures and may be constrained by the reionisation history of the Universe. The current version of the paper presents results for both cases (a) when catalytic nuclear reations are included and (b) when they are neglected.

Charged massive particles (CHAMPs), when present during the Big Bang nucleosynthesis (BBN) era, may significantly alter the synthesis of light elements when compared to a standard BBN scenario. This is due to the formation of bound states with nuclei. This paper presents a detailed numerical and analytical analysis of such CHAMP BBN. All reactions important for predicting light-element yields are calculated within the Born approximation. Three priorly neglected effects are treated in detail:(a) photodestruction of bound states due to electromagnetic cascades induced by the CHAMP decay, (b) late-time efficient destruction/production of H2, Li6, and Li7 due to reactions on charge Z=1 nuclei bound to CHAMPs, and (c) CHAMP exchange between nuclei. Each of these effects may induce orders-of-magnitude changes in the final abundance yields. The study focusses on the impact of CHAMPs on a possible simultaneous solution of the Li6 and Li7 problems. It is shown that a priorly suggested simultaneous solution of the Li6 and Li7 problems for a relic decaying at tau_x = 1000 sec is only very weakly dependent on the relic being neutral or charged, unless its hadronic branching ratio is Bh << 10^-4 very small. By use of a Monte-Carlo analysis it is shown that within CHAMP BBN the existence of further parameter space for a simultaneous solution of the Li6 and Li7 problem for long decay times tau_x > 10^6 sec seems possible but fairly unlikely.

In the framework of the Constrained MSSM we re--examine the gravitino as the lightest superpartner and a candidate for cold dark matter in the Universe. Unlike in other recent studies, we include both a thermal contribution to its relic population from scatterings in the plasma and a non--thermal one from neutralino or stau decays after freeze--out. Relative to a previous analysis [1] we update, extend and considerably improve our treatment of constraints from observed light element abundances on additional energy released during BBN in association with late gravitino production. Assuming the gravitino mass in the GeV to TeV range, and for natural ranges of other supersymmetric parameters, the neutralino region is excluded, while for smaller values of the gravitino mass it becomes allowed again. The gravitino relic abundance is consistent with observational constraints on cold dark matter from BBN and CMB in some well defined domains of the stau region but, in most cases, only due to a dominant contribution of the thermal population. This implies, depending on the gravitino mass, a large enough reheating temperature. If $\mgravitino>1$ GeV then $T_R>10^7$ GeV, if allowed by BBN and other constraints but, for light gravitinos, if $\mgravitino>100$ keV then $T_R>3\times 10^3$ GeV. On the other hand, constraints mostly from BBN imply an upper bound $T_R \lsim {a few}x 10^8\times10^9$ GeV which appears inconsistent with thermal leptogenesis. Finally, most of the preferred stau region corresponds to the physical vacuum being a false vacuum. The scenario can be partially probed at the LHC.

It is shown that gravitinos with mass m_{3/2} ~ 0.1-1 MeV may provide suitable cold dark matter candidates in scenarios of gauge mediated supersymmetry breaking (GMSB) under SO(10) grand unification coupled to supergravity, which accommodate a messenger sector of mass scale M_X ~ 10^6 GeV. This is due to the combined effects of renormalizable loop-suppressed operators and generic non-renormalizable ones governing the dilution of a pre-existing equilibrium gravitino abundance via messenger decay. The above range of gravitino and messenger masses can be accommodated in indirect GMSB scenarios. The gravitino abundance does not depend on the post-inflationary reheat temperature and it is shown that leptogenesis can generate successfully the baryon asymmetry.

Big Bang nucleosynthesis in the presence of decaying relic neutral particles is examined in detail. All non-thermal processes important for the determination of light-element abundance yields of 2H, 3H, 3He, 4He, 6Li, and 7Li are coupled to the thermonuclear fusion reactions to obtain comparatively accurate results. Predicted light-element yields are compared to observationally inferred limits on primordial light-element abundances to infer constraints on the abundances and properties of relic decaying particles with decay times in the interval 0.01 sec < tau < 10^(12) sec. Decaying particles are typically constrained at early times by 4He or 2H, at intermediate times by 6Li, and at large times by the 3He/2H ratio. Constraints are shown for a large number of hadronic branching ratios and decaying particle masses and may be applied to constrain the evolution of the early Universe.

Stable particle dark matter may well originate during the decay of long-lived relic particles, as recently extensively examined in the cases of the axino, gravitino, and higher-dimensional Kaluza-Klein (KK) graviton. It is shown that in much of the viable parameter space such dark matter emerges naturally warm/hot or mixed. In particular, decay produced gravitinos (KK-gravitons) may only be considered cold for the mass of the decaying particle in the several TeV range, unless the decaying particle and the dark matter particle are almost degenerate. Such dark matter candidates are thus subject to a host of cosmological constraints on warm and mixed dark matter, such as limits from a proper reionisation of the Universe, the Lyman-alpha forest, and the abundance of clusters of galaxies.. It is shown that constraints from an early reionsation epoch, such as indicated by recent observations, may potentially limit such warm/hot components to contribute only a very small fraction to the dark matter.

The detection of non-baryonic dark matter through its gamma-ray annihilation in the centre of our galaxy has been studied. The gamma fluxes according to different models have been simulated and compared to those expected to be observed with the Alpha Magnetic Spectrometer (AMS), during a long-term mission on board of the International Space Station. Under the assumption that the dark matter halo is composed of the lightest, stable supersymmetric particle, the neutralino, the results of the simulations in the framework of mSUGRA models, show that with a cuspy dark matter halo or a clumpy halo, the annihilation gamma-ray signal would be detected by AMS. More optimistic perspectives are obtained with the Anomaly Mediated Supersymmetry Breaking (AMSB) model. The latter leads also to a cosmologically important 6Li abundance. Finally, the discovery potential for the massive Kaluza-Klein dark matter candidates has been evaluated and their detection looks feasible.

This paper investigates the parameter space of theories with gauge mediated supersymmetry breaking leading to gravitino (cold) dark matter with mass m_{3/2}= 1keV - 10MeV. We pay particular attention to the cosmological role of messenger fields. Cosmology requires that these messengers decay to the visible sector if the lightest messenger mass M_X > 30TeV. We then examine the various possible messenger number violating interactions allowed by the symmetries of the theory and by phenomenology. Late messenger decay generally results in entropy production hence in the dilution of pre-existing gravitinos. We find that in SU(5) grand unification only specific messenger-matter couplings allow to produce the required amount of gravitino dark matter particles. Gravitino dark matter with the correct abundance is however expected in larger gauge groups such as SO(10) for generic non-renormalizable messenger-matter interactions and for arbritrarily high post-inflationary reheating temperatures.

The synthesis of Li6 during the epoch of Big Bang nucleosynthesis (BBN) due to residual annihilation of dark matter particles is considered. By comparing the predicted Li6 to observations of this isotope in low-metallicity stars, generic constraints on s-wave dark matter annihilation rates into quarks, gauge bosons, and Higgs bosons are derived. It may be shown that, for example, wino dark matter in anomaly-mediated SUSY breaking scenarios with masses Mx < 250 GeV or light neutralinos with Mx < 20 GeV annihilating into light quarks are, taken face value, ruled out. These constraints may only be circumvented if significant Li6 depletion has occurred in all three low-metallicity stars in which this isotope has been observed to date. In general, scenarios invoking non-thermally generated neutralinos with enhanced annihilation rates for a putative explanation of cosmic ray positron or galactic center as well as diffuse background gamma-ray signals by present-day neutralino annihilation will have to face a stringent Li6 overproduction problem. On the other hand, it is possible that Li6 as observed in low-metallicity stars is entirely due to residual dark matter annihilation during BBN, even for neutralinos undergoing a standard thermal freeze-out.

(abridged) A detailed examination of the evolution of stochastic magnetic fields between high cosmic temperatures and the present epoch is presented. A simple analytical model matching the results of the 3D MHD simulations allows for the prediction of present day magnetic field correlation lengths and energy. Our conclusions are multi fold. (a) Initial primordial fields with only a small amount of helicity are evolving into maximally helical fields. (b) There exists a correlation between the strength of the magnetic field, B, at the peak of it's spectrum and the location of the peak, given at the present epoch by: B ~ 5x10^{-12} (L/kpc) Gauss, where L is the correlation length determined by the initial magnetic field. (c) Concerning studies of generation of cosmic microwave background (CMBR) anisotropies due to primordial magnetic fields of B~10^{-9} Gauss on ~ 10 Mpc scales, such fields are not only impossible to generate in early causal magnetogenesis scenarios but also seemingly ruled out by distortions of the CMBR spectrum due to magnetic field dissipation on smaller scales and the overproduction of cluster magnetic fields. (d) The most promising detection possibility of CMBR distortions due to primordial magnetic fields may be on much smaller scales at higher multipoles l~10^6 where the signal is predicted to be the strongest. (e) It seems possible that magnetic fields in clusters of galaxies are entirely of primordial origin, without invoking dynamo amplification. Such fields would be of (pre-collapse) strength 10^{-12} - 10^{-11} Gauss with correlation lengths in the kpc range, and would also exist in voids of galaxies.

We present results of a detailed and fully non-linear numerical and analytical investigation of magnetic field evolution from the very earliest cosmic epochs to the present. We find that, under reasonable assumptions concerning the efficiency of a putative magnetogenesis era during cosmic phase transitions, surprisingly strong magnetic fields 10^{-13} - 10^{-11} Gauss, on comparatively small scales 100 pc - 10 kpc may survive to the present. Building on prior work on the evolution of magnetic fields during the course of gravitational collapse of a cluster, which indicates that pre-collapse fields of 4\times 10^{-12} Gauss extant on small scales may suffice to produce clusters with acceptable Faraday rotation measures, we question the widely hold view that cluster magnetic fields may not be entirely of primordial origin.

Results of a detailed examination of the cascade nucleosynthesis resulting from the putative hadronic decay, evaporation, or annihilation of a primordial relic during the Big Bang nucleosynthesis (BBN) era are presented. It is found that injection of energetic nucleons around cosmic time 1000 sec may lead to an observationally favored reduction of the primordial Li7/H yield by a factor 2 - 3. Moreover, such sources also generically predict the production of the Li6 isotope with magnitude close to the as yet unexplained high Li6 abundances in low-metallicity stars. The simplest of these models operate at fractional contribution to the baryon density Omega_b h^2> 0.025, slightly larger than that inferred from standard BBN. Though further study is required, such sources, as for example due to the decay of the next-to-lightest supersymmetric particle into GeV gravitinos or the decay of an unstable gravitino in the TeV range of abundance \Omega_G h^2\sim 5\times 10^{-4} show promise to explain both the Li6 and Li7 abundances in low metallicity stars.

It is widely believed that dark matter exists within galaxies and clusters of galaxies. Under the assumption that this dark matter is composed of the lightest, stable supersymmetric particle, assumed to be the neutralino, the feasibility of its indirect detection via observations of a diffuse gamma-ray signal due to neutralino annihilation within M31 is examined.

It is widely believed that dark matter exists within galaxies and clusters of galaxies. Under the assumption that this dark matter is composed of the lightest, stable supersymmetric particle, assumed to be the neutralino, the feasibility of its indirect detection via observations of a diffuse gamma-ray signal due to neutralino annihilations within M31 is examined. To this end, first the dark matter halo of the close spiral galaxy M31 is modeled from observations, then the resultant gamma-ray flux is estimated within supersymmetric model configurations. We conclude that under favorable conditions such as the rapid accretion of neutralinos on the central black hole in M31 and/or the presence of many clumps inside its halo with $r^{-3/2}$ inner profiles, a neutralino annihilation gamma-ray signal is marginally detectable by the ongoing collaboration CELESTE.

The distant globular cluster Palomar 13 has been found to have a very high mass-to-light ratio and its profile can be well fitted either by a King model with a tail, or with a NFW model. This cluster may be the first case of the many clumps predicted by CDM simulations that would not be disrupted by the galactic halo potential. We make the hypothesis that Pal 13 is made of neutralinos and run the DarkSuspect code to estimate the high-energy photon flux due to the annihilation of neutralinos through various channels in some benchmark scenarios. These low fluxes may be used as targets to be reached in proposals for future ground-based high altitude Cerenkov telescopes.

Microlensing searches aim to detect compact halo dark matter via its gravitational lensing effect on stars within the Large Magellanic Cloud. The most recent results have led to the claim that roughly one fifth of the galactic halo dark matter may be in the form of compact, solar-mass objects. We analyze this hypothesis by considering the goodness-of-fit of the best-fit halo dark matter solutions to the observational data. We show that the distribution of the durations of the observed microlensing events is narrower than that expected to result from a standard halo lens population at 90 to 95% confidence. Only when the fraction of expected background (non-halo) events is significantly increased does the discrepancy between the observed and theoretical event duration distributions disappear. This indicates that it is possible that most of the lenses responsible for the observed microlensing events are not located within the Milky Way halo.

It is widely believed that the cosmic baryon density may be obtained by inferring deuterium abundances in low-metallicity quasar absorption line systems. The implicit assumptions which enter this argument are critically assessed. In particular, the production of deuterium in non-standard Big Bang nucleosynthesis scenarios, the putative production of deuterium in astrophysical environments, and the possible destruction of deuterium via non-standard chemical evolution are discussed.

This article gives a brief qualitative description of the possible evolution of the early Universe between the end of an inflationary epoch and the end of Big Bang nucleosynthesis. After a general introduction, establishing the minimum requirements cosmologists impose on this cosmic evolutionary phase, namely, successful baryogenesis, the production of cosmic dark matter, and successful light-element nucleosynthesis, a more detailed discussion on some recent developments follows. This latter includes the physics of preheating, the putative production of (alternative) dark matter, and the current status of Big Bang nucleosynthesis.

The primordial helium abundance Y_p is commonly inferred from abundance determinations in low-metallicity extragalactic HII-regions. Such determinations may be subject to systematic uncertainties that are investigated here. Particular attention is paid to two effects: icf-corrections for ``imperfect'' ionization structure leading to significant amounts of (unobservable) neutral helium or hydrogen and ``tcf''-corrections due to non-uniform temperature. Model HII-regions with a large number of parameters are constructed and it is shown that required corrections are almost exclusively functions of two physical parameters: the number of helium- to hydrogen- ionizing photons in the illuminating continuum Q(He0)/Q(H0), and the ratio of width to radius delta r_S/r_S of the Stromgren sphere. For clouds of sufficient He-ionizing photons Q(He0)/Q(H0) >~ 0.15 and non-negligible width of the Stromgren sphere a significant overestimate of helium abundances may result. Such clouds show radiation softness parameters in the range -0.4 <~ log(eta) <~ 0.3 coincident with the range of eta in observed HII-regions. Existing data of HII-regions indeed seem to display a correlation which is consistent with a typical ~ 2-4% overestimate of helium abundances due to these effects. In case such an interpretation prevails, and in the absence of other compensating effects, a significant downward revision of Y_p may result. It is argued that caution should be exercised regarding the validity of commonly quoted error bars on Y_p.

We examine the Big Bang nucleosynthesis (BBN) process in the presence of small-scale baryon inhomogeneities. Primordial abundance yields for D, He4, Li6, Li7, Be9, and B11 are computed for wide ranges of parameters characterizing the inhomogeneities taking account of all relevant diffusive and hydrodynamic processes. These calculations may be of interest due to (a) recent observations of the anisotropies in the cosmic microwave background radiation favoring slightly larger baryonic contribution to the critical density, Omega_b, than allowed by a standard BBN scenario and (b) new observational determinations of Li6 and Be9 in metal-poor halo stars. We find considerable parameter space in which production of D and He4 is in agreement with observational constraints even for Omega_b h^2 a factor 2-3 larger than the Omega_b inferred from standard BBN. Nevertheless, in this parameter space synthesis of Li7 in excess of the inferred Li7 abundance on the Spite plateau results. Production of Li6, Be9, and B11 in inhomogeneous BBN scenarios is still typically well below the abundance of these isotopes observed in the most metal-poor stars to date thus neither confirming nor rejecting inhomogeneous BBN. In an appendix we summarize results of a reevaluation of baryon diffusion constants entering inhomogeneous BBN calculations.

It has been recently suggested that stable, supersymmetric B-balls formed in the early universe could not only be the dark matter at the present epoch, but also be responsible for baryogenesis by their partial evaporation at high temperatures. We reinvestigate the efficiency of B-ball baryogenesis and find it to be limited by the diffusion of baryon number away from the B-balls. Successful baryogenesis may only occur for B-balls with charges Q\simle 10^{20} - 5\times 10^{23}, which is close to the observational lower limits on the Q of a significant B-ball dark matter component. We also present some cosmological constraints on the abundances of larger B-balls in the early universe.

Spatially varying primordial magnetic fields may be efficiently dissipated prior to the epoch of recombination due to the large viscosity of the baryon-photon fluid. We show that this dissipation may result in observable chemical potential mu and Compton y distortions in the cosmic microwave background radiation (CMB) spectrum. Current upper limits on mu and y from FIRAS constrain magnetic fields to have strength B_0 < 3\times 10^{-8}Gauss (scaled to the present) between comoving coherence length \approx 400 pc and \approx 0.6 Mpc. This represents the strongest upper limit on small-scale primordial magnetic fields to date.

We present detailed numerical calculations of the light element abundances synthesized in a Universe consisting of matter- and antimatter- domains, as predicted to arise in some electroweak baryogenesis scenarios. In our simulations all relevant physical effects, such as baryon-antibaryon annihilations, production of secondary particles during annihilations, baryon diffusion, and hydrodynamic processes are coupled to the nuclear reaction network. We identify two dominant effects, according to the typical spatial dimensions of the domains. Small antimatter domains are dissipated via neutron diffusion prior to He4 synthesis at T_He4 \approx 80 keV, leading to a suppression of the primordial He4 mass fraction. Larger domains are dissipated below T_He4 via a combination of proton diffusion and hydrodynamic expansion. In this case the strongest effects on the elemental abundances are due to anti-p He4 annihilations, leading to an overproduction of He3 relative to H2 and to overproduction of Li6 via non-thermal nuclear reactions. Both effects may result in light element abundances deviating substantially from the standard Big Bang Nucleosynthesis yields and from the observationally inferred values. This allows us to derive stringent constraints on the antimatter parameters. For some combinations of the parameters, one may obtain both, low H2 and low He4, at a common value of the cosmic baryon density, a result seemingly favored by current observational data.

I consider the synthesis of 6Li due to the decay of relic particles, such as gravitinos or moduli, after the epoch of Big Bang Nucleosynthesis. The synthesized 6Li/H ratio may be compared to 6Li/H in metal-poor stars which, in the absence of stellar depletion of 6Li, yields significantly stronger constraints on relic particle densities than the usual consideration of overproduction of 3He. Production of 6Li during such an era of non-thermal nucleosynthesis may also be regarded as a possible explanation for the relatively high 6Li/H ratios observed in metal-poor halo stars.

Standard Big Bang nucleosynthesis predicts an essentially zero primordial metallicity. I speculate on possible metal (i.e. nucleon number A\geq 12) production in scenarios of inhomogeneous Big Bang nucleosynthesis. It is conceivable, though not necessarily probable, that some primordial metallicity is synthesized if a small fraction of all cosmic baryons reside in very high-density regions. Such conditions could possibly result from the evaporation of some baryon-number carrying soliton prior to the epoch of Big Bang nucleosynthesis.

We consider the generation of super-horizon metric fluctuations during an epoch of preheating in the presence of a scalar field \chi quadratically coupled to the inflaton. We find that the requirement of efficient broad resonance is concomitant with a severe damping of super-horizon \delta\chi quantum fluctuations during inflation. Employing perturbation theory with backreaction included as spatial averages to second order in the scalar fields and in the metric, we argue that the usual inflationary prediction for metric perturbations on scales relevant for structure formation is not strongly modified.

We present a numerical investigation of the gravitational collapse of horizon-size density fluctuations to primordial black holes (PBHs) during the radiation-dominated phase of the Early Universe. The collapse dynamics of three different families of initial perturbation shapes, imposed at the time of horizon crossing, is computed. The perturbation threshold for black hole formation, needed for estimations of the cosmological PBH mass function, is found to be $\delta_{\rm c} \approx 0.7$ rather than the generally employed $\delta_{\rm c} \approx 1/3$, if $\delta$ is defined as $\Delta M/\mh$, the relative excess mass within the initial horizon volume. In order to study the accretion onto the newly formed black holes, we use a numerical scheme that allows us to follow the evolution for long times after formation of the event horizon. In general, small black holes (compared to the horizon mass at the onset of the collapse) give rise to a fluid bounce that effectively shuts off accretion onto the black hole, while large ones do not. In both cases, the growth of the black hole mass owing to accretion is insignificant. Furthermore, the scaling of black hole mass with distance from the formation threshold, known to occur in near-critical gravitational collapse, is demonstrated to apply to primordial black hole formation.

Primordial black holes (PBHs) may form in the early universe when pre-existing adiabatic density fluctuations enter into the cosmological horizon and recollapse. It has been suggested that PBH formation may be facilitated when fluctuations enter into the horizon during a strongly first-order phase transition which proceeds in approximate equilibrium. We employ general-relativistic hydrodynamics numerical simulations in order to follow the collapse of density fluctuations during first-order phase transitions. We find that during late stages of the collapse fluctuations separate into two regimes, an inner part existing exclusively in the high-energy density phase with energy density $\epsilon_{\rm h}$, surrounded by an outer part which exists exclusively in the low-energy density phase with energy density $\epsilon_{\rm h}-L$, where $L$ is the latent heat of the transition. We confirm that the fluctuation density threshold $\delta\epsilon /\epsilon$ required for the formation of PBHs during first-order transitions decreases with increasing $L$ and falls below that for PBH formation during ordinary radiation dominated epochs. Our results imply that, in case PBHs form at all in the early universe, their mass spectrum is likely dominated by the approximate horizon masses during epochs when the universe undergoes phase transitions.

We investigate Big Bang nucleosynthesis (hereafter, BBN) in a cosmic environment characterized by a distribution of small-scale matter/antimatter domains. Production of antimatter domains in a baryo-asymmetric universe is predicted in some electroweak baryogenesis scenarios. We find that cosmic antimatter domains of size exceeding the neutron-diffusion length at temperature T approx. 1 MeV significantly affect the light-element production. Annihilation of antimatter preferentially occurs on neutrons such that antimatter domains may yield a reduction of the He-4 abundance relative to a standard BBN scenario. In the limiting case, all neutrons will be removed before the onset of light-element production, and a universe with net baryon number but without production of light elements results. In general, antimatter domains spoil agreement between BBN abundance yields and observationally inferred primordial abundances limits which allows us to derive limits on their presence in the early universe. However, if only small amounts of antimatter are present, BBN with low deuterium and low He-4, as seemingly favored by current observational data, is possible.

The recent discovery of critical phenomena arising in gravitational collapse near the threshold of black hole formation is used to estimate the initial mass function of primordial black holes (PBHs). It is argued that the universal scaling relation between black hole mass and initial perturbation found for a variety of collapsing space-times also applies to PBH formation, indicating the possibility of the formation of PBHs with masses much smaller than one horizon mass. Owing to the natural fine-tuning of initial conditions by the exponential decline of the probability distribution for primordial density fluctuations, sub-horizon mass PBHs are expected to form at all epochs. This result suggests that the constraints on the primordial fluctuation spectrum based on the abundance of PBHs at different mass scales may have to be revisited.

A brief overview of non-standard big bang nucleosynthesis (BBN) scenarios is presented. Trends and results of the light-element nucleosynthesis in BBN scenarios with small-scale or large-scale inhomogeneity, the presence of antimatter domains, stable or unstable massive neutrinos, neutrino oscillations, neutrino degeneracy, or massive decaying particles are summarized.

Observations by the MACHO collaboration indicate that a significant fraction of the galactic halo dark matter may be in form of compact objects with masses $M\sim 0.5M_{\odot}$. Identification of these objects as red or white dwarfs is problematic due to stringent observational upper limits on such dwarf populations. Primordial black hole (PBH) formation from pre-existing density fluctuations is facilitated during the cosmic QCD transition due to a significant decrease in pressure forces. For generic initial density perturbation spectra this implies that essentially all PBHs may form with masses close to the QCD-horizon scale, $M_h^{QCD}\sim 1M_{\odot}$. It is possible that such QCD PBHs contribute significantly to the closure density today. I discuss the status of theoretical predictions for the properties of QCD PBH dark matter. Observational signatures of and constraints on a cosmic solar mass PBH population are also discussed.

We examine the evolution of magnetic fields in an expanding fluid composed of matter and radiation with particular interest in the evolution of cosmic magnetic fields. We derive the propagation velocities and damping rates for relativistic and non-relativistic fast and slow magnetosonic, and Alfv\'en waves in the presence of viscous and heat conducting processes. The analysis covers all MHD modes in the radiation diffusion and the free-streaming regimes. When our results are applied to the evolution of magnetic fields in the early universe, we find that cosmic magnetic fields are damped from prior to the epoch of neutrino decoupling up to recombination. Our findings have multifold implications for cosmology. The dissipation of magnetic field energy into heat during the epoch of neutrino decoupling ensures that most magnetic field configurations generated in the very early universe satisfy big bang nucleosynthesis constraints. Further dissipation before recombination constrains models in which primordial magnetic fields give rise to galactic magnetic fields or density perturbations. Finally, the survival of Alfv\'en and slow magnetosonic modes on scales well below the Silk mass may be of significance for the formation of structure on small scales (abridged).

We employ observationally determined intrinsic velocity widths and column densities of damped Lyman-alpha systems at high redshift to investigate the distribution of baryons in protogalaxies within the context of a standard cold dark matter model. We proceed under the assumption that damped Lyman alpha systems represent a population of cold, rotationally supported, protogalactic disks and that the abundance of protogalactic halos is well approximated by a cold dark matter model with critical density and vanishing cosmological constant. Using conditional cross sections to observe a damped system with a given velocity width and column density, we compare observationally inferred velocity width and column density distributions to the corresponding theoretically determined distributions for a variety of disk parameters and CDM normalizations. In general, we find that the observations can not be reproduced by the models for most disk parameters and CDM normalizations. Whereas the column density distribution favors small disks with large neutral gas fraction, the velocity width distribution favors large and thick disks with small neutral gas fraction. The possible resolutions of this problem in the context of this CDM model may be: (1) an increased contribution of rapidly rotating disks within massive dark matter halos to damped Lyman-alpha absorption or (2) the abandoning of simple disk models within this CDM model for damped Lyman-alpha systems at high redshift. Here the first possibility may be achieved by supposing that damped Lya system formation only occurs in halos with fairly large circular velocities and the second possibility may result from a large contribution of mergers and double-disks to damped Lya absorption at high redshift.

We consider the formation of horizon-size primordial black holes (PBH's) from pre-existing density fluctuations during cosmic phase transitions. It is pointed out that the formation of PBH's should be particularly efficient during the QCD epoch due to a substantial reduction of pressure forces during adiabatic collapse, or equivalently, a significant decrease in the effective speed of sound during the color-confinement transition. Our considerations imply that for generic initial density perturbation spectra PBH mass functions are expected to exhibit a pronounced peak on the QCD-horizon mass scale $\sim 1 M_{\odot}$. This mass scale is roughly coincident with the estimated masses for compact objects recently observed in our galactic halo by the MACHO collaboration. Black holes formed during the QCD epoch may offer an attractive explanation for the origin of halo dark matter evading possibly problematic nucleosynthesis and luminosity bounds on baryonic halo dark matter.

We give an improved estimate of primordial magnetic fields generated during cosmological first order phase transitions. We examine the charge distribution at the nucleated bubble wall and its dynamics. We consider instabilities on the bubble walls developing during the phase transition. It is found that damping of these instabilities due to viscosity and heat conductivity caused by particle diffusion can be important in the QCD phase transition, but is probably negligible in the electroweak transition. We show how such instabilities together with the surface charge densities on bubble walls excite magnetic fields within a certain range of wavelengths. We discuss how these magnetic seed fields may be amplified by MHD effects in the turbulent fluid. The strength and spectrum of the primordial magnetic field at the present time for the cases where this mechanism was operative during the electroweak or the QCD phase transition are estimated. On a 10 Mpc comoving scale, field strengths of the order 10**(-29) G for electroweak and 10**(-20) G for QCD, could be attained for reasonable phase transition parameters.

Detections of deuterium in high redshift Lyman limit absorption systems along the line of sight to QSOs promise to reveal the primordial deuterium abundance. At present, the deuterium abundances (D/H) derived from the very few systems observed are significantly discordant. Assuming the validity of all the data, if this discordance does not reflect intrinsic primordial inhomogeneity, then it must arise from processes operating after the primordial nucleosynthesis epoch. We consider processes which might lead to significant deuterium production/destruction, yet allow the cloud to mimick a chemically unevolved system. These processes include, for example, anomalous/stochastic chemical evolution and D/{$^4$He} photo-destruction. In general, we find it unlikely that these processes could have altered significantly (D/H) in Lyman limit clouds. We argue that chemical evolution scenarios, unless very finely tuned, cannot account for significant local deuterium depletion since they tend to overproduce $^{12}$C, even when allowance is made for possible outflow. Similarly, D/{$^4$He} photo-destruction schemes engineered to locally produce or destroy deuterium founder on the necessity of requiring an improbably large $\gamma$-ray source density. Future observations of (D/H) in Lyman limit systems may provide important insight into the initial conditions for the primordial nucleosynthesis process, early chemical evolution, and the galaxy formation process.

We examine effects on primordial nucleosynthesis from a truly random spatial distribution in the baryon-to-photon ratio ($\eta$). We generate stochastic fluctuation spectra characterized by different spectral indices and root-mean-square fluctuation amplitudes. For the first time we explicitly calculate the effects of baryon diffusion on the nucleosynthesis yields of such stochastic fluctuations. We also consider the collapse instability of large-mass-scale inhomogeneities. Our results are generally applicable to any primordial mechanism producing fluctuations in $\eta$ which can be characterized by a spectral index. In particular, these results apply to primordial isocurvature baryon fluctuation (PIB) models. The amplitudes of scale-invariant baryon fluctuations are found to be severely constrained by primordial nucleosynthesis. However, when the $\eta$ distribution is characterized by decreasing fluctuation amplitudes with increasing length scale, surprisingly large fluctuation amplitudes on the baryon diffusion scale are allowed.

We consider the production of $^3$He and $^2$H by $^4$He photodisintegration initiated by non-thermal energy releases during early cosmic epochs. We find that this process cannot be the predominant source of primordial $^2$H since it would result in anomalously high $^3$He/D ratios in conflict with standard chemical evolution assumptions. We apply this fact to constrain topological defect models of highest energy cosmic ray (HECR) production. Such models have been proposed as possible sources of ultrahigh energy particles and gamma-rays with energies above $10^{20}$eV. The constraints on these models derived from $^4$He-photodisintegration are compared to corresponding limits from spectral distortions of the cosmic microwave background radiation (CMBR) and from the observed diffuse gamma-ray background. It is shown that for reasonable primary particle injection spectra superconducting cosmic strings, unlike ordinary strings or annihilating monopoles, cannot produce the HECR flux at the present epoch without violating at least the $^4$He-photodisintegration bound. The constraint from the diffuse gamma-ray background rules out the dominant production of HECR by the decay of Grand Unification particles in models with cosmological evolution assuming standard fragmentation functions. Constraints on massive black hole induced photodisintegration are also discussed.

We study big bang nucleosynthesis in the presence of large mass-scale, non-linear entropy fluctuations. Overdense regions, with masses above the local baryon-Jeans mass, are expected to collapse and form condensed objects. Surviving nucleosynthesis products therefore tend to originate from underdense regions. We compute expected surviving light element ($^2$H, $^3$He, $^4$He, $^7$Li) abundance yields for a variety of stochastic fluctuation spectra. In general, we find that spectra with significant power in fluctuations on length scales below that of the local baryon Jeans mass produce nucleosynthesis yields which are in conflict with observationally inferred primordial abundances. However, when this small scale structure is absent or suppressed, and the collapse efficiency of overdense regions is high, there exists a range of fluctuation spectral characteristics which meet all primordial abundance constraints. In such models abundance constraints can be met even when the pre-collapse baryonic fraction of the closure density is $\Omega_b\approx 0.2h^{-2}$($h$ is the Hubble parameter in units of 100 km\ s$^{-1}$Mpc$^{-1}$). Nucleosynthesis in these models is characterized by high $^2$H/H and low $^4$He mass fraction relative to a homogeneous big bang at a given value of $\Omega_bh^2$. A potentially observable signature of these models is the production of intrinsic primordial abundance variations on baryon mass-scales up to $10^{10}M_{\odot}-10^{12}_{\odot}$.

The formalism by Press and Schechter (PS) is often used to infer number densities of virialized objects of mass M (e.g. quasars, galaxies, clusters of galaxies, etc.) from a count of initially overdense regions in a Gaussian density perturbation field. We reanalyze the PS-formalism by explicitly counting underdense regions which are embedded within overdense regions, so called cloud-in-clouds. In contrast to the original PS-formalism, our revised analysis automatically accounts for all the cosmic material. We find that mass distribution functions for virialized objects are altered by the proper solution of the cloud-in-cloud problem. These altered distribution functions agree much better with distribution functions inferred form N-body simulations than the original PS-distribution functions.

We investigate the transport of baryon number across phase boundaries in a putative first order QCD-phase transition. Two independent phenomenological models are employed to estimate the baryon penetrability at the phase boundary: chromoelectric flux tube models; and an analogy to baryon-baryon coalescence in nuclear physics. Our analysis indicates that baryon transport across phase boundaries may be order of magnitude more efficient than other work has suggested. We discuss the substantial uncertainties involved in estimating baryon penetrability at phase boundaries.

Primordial nucleosynthesis calculations are shown to be able to provide constraints on electroweak baryogenesis which produce a highly inhomogeneous distribution of the baryon-to-photon ratio. Such baryogenesis scenarios overproduce 4He and/or 7Li and can be ruled out whenever a fraction f<3*10e-6(100 GeV/T)^3 of nucleated bubbles of broken-symmetry phase contributes > 10% of the baryon number within the horizon volume.

We show that a class of inhomogeneous big bang nucleosynthesis models exist which yield light-element abundances in agreement with observational constraints for baryon-to-photon ratios significantly smaller than those inferred from standard homogeneous big bang nucleosynthesis (HBBN). These inhomogeneous nucleosynthesis models are characterized by a bimodal distribution of baryons in which some regions have a local baryon-to-photon ratio eta=3*10e-10, while the remaining regions are baryon-depleted. HBBN scenarios with primordial (2H+3He)/H<9*10e-5 necessarily require that most baryons be in a dark or non-luminous form, although new observations of a possible high deuterium abundance in Lyman-alpha clouds may relax this requirement somewhat. The models described here present another way to relax this requirement and can even eliminate any lower bound on the baryon-to-photon ratio.

We present a detailed study of inhomogeneous Big Bang nucleosynthesis where, for the first time, nuclear reactions are coupled to all significant fluctuation dissipation processes. Theses processes include neutrino heat transport, baryon diffusion, photon diffusive heat transport, and hydrodynamic expansion with photon-electron Thomson drag. Light element abundance yields are presented for broad ranges of initial amplitudes and length scales for spherically condensed fluctuations. The $^2$H, $^3$He, $^4$He, and $^7$Li nucleosynthesis yields are found to be inconsistent with observationally inferred primordial abundances for all but very narrow ranges of fluctuation characteristics. Rapid hydrodynamic expansion of fluctuations late in the nucleosynthesis epoch results in significant destruction of $^7$Li ($^7$Be) only if the baryonic conytribution to the closure density ($\Omega_b$) is less than or comparable to the upper limit on this quantity from homogeneous Big Bang nucleosynthesis. This implies that $^7$Li overproduction will peclude an increase on the upper limit for $\Omega_b$ for any inhomogeneous nucleosynthesis scenarios employing spherically condensed fluctuations.

We show that primordial nucleosynthesis in baryon inhomogeneous big-bang models can lead to significant heavy-element production while still satisfying all the light-element abundance constraints including the low lithium abundance observed in population II stars. The parameters which admit this solution arise naturally from the process of neutrino induced inflation of baryon inhomogeneities prior to the epoch of nucleosynthesis. These solutions entail a small fraction of baryons (\le 2\%) in very high density regions with local baryon-to-photon ratio $\eta^h\approx 10^{-4}$, while most baryons are at a baryon-to-photon ratio which optimizes the agreement with light-element abundances. The model would imply a unique signature of baryon inhomogeneities in the early universe, evidenced by the existence of primordial material containing heavy-element products of proton and alpha- burning reactions with an abundance of $[Z]\sim -6 to -4$.

We examine the damping of non-linear sub-horizon scale entropy fluctuations in early epochs of the universe ($T\approx 100$ GeV to $T\approx 1$keV) by neutrino, baryon, and photon induced dissipative processes. Results of numerical evolution calculations are presented for broad ranges of initial fluctuation amplitudes and length scales. These calculations include a detailed treatment of neutrino inflation, neutron and proton diffusion, photon diffusive heat transport, and hydrodynamic expansion with photon-electron Thomson drag. Neutrino inflation is treated both in the diffusive heat transport regime where neutrinos are optically thick on the length scales of the fluctuations, and in the homogeneous heating regime where neutrinos are optically thin on these scales. We find considerable convergence in amplitude evolution for appreciable ranges in initial fluctuation length scales and amplitudes. Fluctuations produced with the right characteristics at very early times ($T$ $^>_{\sim}$ 100 GeV) are found to survive through the nucleosynthesis epoch.