We study the phase space structure of exact quantum Wightman functions in spatially homogeneous, temporally varying systems. In addition to the usual mass shells, the Wightman functions display additional coherence shells around zero frequency $k_0=0$, which carry the information of the local quantum coherence of particle-antiparticle pairs. We find also other structures, which encode non-local correlations in time, and discuss their role and decoherence. We give a simple derivation of the cQPA formalism, a set of quantum transport equations, that can be used to study interacting systems including the local quantum coherence. We compute quantum currents created by a temporal change in a particle's mass, comparing the exact Wightman function approach, the cQPA and the semiclassical methods. We find that the semiclassical approximation, which is fully encompassed by the cQPA, works surprisingly well even for very sharp temporal features. This is encouraging for the application of semiclassical methods in electroweak baryogenesis with strong phase transitions.

It is widely believed that electroweak baryogenesis should be suppressed in strong phase transitions with fast-moving bubble walls, but this effect has never been quantitatively studied. We rederive fluid equations describing transport of particle asymmetries near the bubble wall without making the small-wall-velocity approximation. We show that the suppression of the baryon asymmetry is a smooth function of the wall speed and that there is no special behavior when crossing the sound speed barrier. Electroweak baryogenesis can thus be efficient also with strong detonations, generically associated with models with observably large gravitational waves. We also make a systematic and critical comparison of our improved transport equations to another one commonly used in the literature, based on the VEV-insertion formalism.

The dark matter annihilation channels sometimes involve sharp resonances. In such cases the usual momentum averaged approximations for computing the DM abundance may not be accurate. We develop an easily accessible momentum dependent framework for computing the DM abundance accurately and efficiently near such features. We apply the method to the case of a singlet scalar dark matter $s$ interacting with SM through higgs portal $\lambda_{\rm hs}s^2 h^2$ and compare the results with different momentum averaged methods. The accuracy of the latter depend strongly on the strength of the elastic interactions and corrections are large if WIMP has negligible interactions beyond the main annihilation channel. In the singlet scalar model however, the standard model scatterings induce an efficient kinetic equilibrium that validates the momentum averaged computation to 20 per cent accuracy. We update the current extent of the allowed region in the light singlet scalar dark matter to $m_{\rm S} \in [56,62.5]$ GeV.

We demonstrate the existence of a generic, efficient and purely gravitational channel producing a significant abundance of dark relics during reheating after the end of inflation. The mechanism is present for any inert scalar with the non-minimal curvature coupling $\xi R\chi^2$ and the relic production is efficient for natural values $\xi = {\cal O}(1)$. The observed dark matter abundance can be reached for a broad range of relic masses extending from $m \sim 1 {\rm k eV}$ to $m \sim 10^{8} {\rm GeV}$, depending on the scale of inflation and the dark sector couplings. Frustratingly, such relics escape direct, indirect and collider searches since no non-gravitational couplings to visible matter are needed.

We investigate the CMB $\mu$ distortion in models where two uncorrelated sources contribute to primordial perturbations. We parameterise each source by an amplitude, tilt, running and running of the running. We perform a detailed analysis of the distribution signal as function of the model parameters, highlighting the differences compared to single-source models. As a specific example, we also investigate the mixed inflaton-curvaton scenario. We find that the $\mu$ distortion could efficiently break degeneracies of curvaton parameters especially when combined with future sensitivity of probing the tensor-to-scalar ratio $r$. For example, assuming bounds $\mu < 0.5 \times 10^{-8}$ and $r<0.01$, the curvaton contribution should either vanish or the curvaton should dominate primordial perturbations and its slow-roll parameter $\eta_{\chi}$ is constrained to the interval $-0.007 < \eta_{\chi}< 0.045$.

Adding an extra singlet scalar $S$ to the Higgs sector can provide a barrier at tree level between a false vacuum with restored electroweak symmetry and the true one. This has been demonstrated to readily give a strong phase transition as required for electroweak baryogenesis. We show that with the addition of a fermionic dark matter particle $\chi$ coupling to $S$, a simple UV-complete model can realize successful electroweak baryogenesis. The dark matter gets a CP asymmetry that is transferred to the standard model through a $CP\ portal\ interaction$, which we take to be a coupling of $\chi$ to $\tau$ leptons and an inert Higgs doublet. The CP asymmetry induced in left-handed $\tau$ leptons biases sphalerons to produce the baryon asymmetry. The model has promising discovery potential at the LHC, while robustly providing a large enough baryon asymmetry and correct dark matter relic density with reasonable values of the couplings.

We consider a fermionic dark matter candidate interacting via a scalar mediator coupled with the Standard Model through a Higgs portal. We consider general setting including both scalar and pseudoscalar interactions between the scalar and fermion, and illustrate the relevant features for dark matter abundance, direct search limits and collider constraints. The case where dark matter has a self-interaction strength $\left\langle \sigma_V \right\rangle/m_\psi \sim 0.1-1 \,\mathrm{cm}^2/\mathrm{g}$ is strongly constrained, in particular by the Big Bang Nucleosynthesis. We show that these constraints can be alleviated by introducing a new light sterile neutrino $N$. The allowed region for the extended model consists of a triangle at $10 (\sin \theta)^{-2/5} \,{\rm MeV} \lesssim m_N \lesssim 1\,{\rm GeV}$.

We investigate an extension of the Standard Model containing two Higgs doublets and a singlet scalar field (2HDSM). We show that the model can have a strongly first-order phase transition and give rise to the observed baryon asymmetry of the Universe, consistent with all experimental constraints. In particular, the constraints from the electron and neutron electric dipole moments are less constraining here than in pure two-Higgs-doublet model (2HDM). The two-step, first-order transition in 2HDSM, induced by the singlet field, may lead to strong supercooling and low nucleation temperatures in comparison with the critical temperature, $T_n \ll T_c$, which can significantly alter the usual phase-transition pattern in 2HD models with $T_n \approx T_c$. Furthermore, the singlet field can be the dark matter particle. However, in models with a strong first-order transition its abundance is typically but a thousandth of the observed dark matter abundance.

We consider portal models which are ultraweakly coupled with the Standard Model, and confront them with observational constraints on dark matter abundance and isocurvature perturbations. We assume the hidden sector to contain a real singlet scalar $s$ and a sterile neutrino $\psi$ coupled to $s$ via a pseudoscalar Yukawa term. During inflation, a primordial condensate consisting of the singlet scalar $s$ is generated, and its contribution to the isocurvature perturbations is imprinted onto the dark matter abundance. We compute the total dark matter abundance including the contributions from condensate decay and nonthermal production from the Standard Model sector. We then use the Planck limit on isocurvature perturbations to derive a novel constraint connecting dark matter mass and the singlet self coupling with the scale of inflation: $m_{\rm DM}/{\rm GeV}\lesssim 0.2\lambda_{\rm s}^{\scriptscriptstyle 3/8} \left(H_*/10^{\scriptscriptstyle 11}{\rm GeV}\right)^{\scriptscriptstyle -3/2}$. This constraint is relevant in most portal models ultraweakly coupled with the Standard Model and containing light singlet scalar fields.

We consider dark matter in a minimal extension of the Standard Model (SM) which breaks electroweak symmetry dynamically and leads to a complete unification of the SM and technicolor coupling constants. The unification scale is determined to be $M_{\rm U} \approx 2.2 \times 10^{15}$ GeV and the unified coupling $\alpha_{\rm U} \approx 0.0304$. Moreover, unification strongly suggest that the technicolor sector of the model must become strong at the scale of ${\cal O}$(TeV). The model also contains a tightly constrained sector of mixing neutral fields stabilized by a discrete symmetry. We find the lightest of these states can be DM with a mass in the range $m_{\rm DM} \approx 30-800$ GeV. We find a large set of parameters that satisfy all available constraints from colliders and from dark matter search experiments. However, most of the available parameter space is within the reach of the next generation of DM search experiments. The model is also sensitive to a modest improvement in the measurement of the precision electroweak parameters.

One of the simplest models of dark matter is that where a scalar singlet field S comprises some or all of the dark matter, and interacts with the standard model through an HHSS coupling to the Higgs boson. We update the present limits on the model from LHC searches for invisible Higgs decays, the thermal relic density of S, and dark matter searches via indirect and direct detection. We point out that the currently allowed parameter space is on the verge of being significantly reduced with the next generation of experiments. We discuss the impact of such constraints on possible applications of scalar singlet dark matter, including a strong electroweak phase transition, and the question of vacuum stability of the Higgs potential at high scales.

The inert doublet dark matter model has recently gained attention as a possible means of facilitating a strongly first order electroweak phase transition (EWPT), as needed for baryogenesis. We extend previous results by considering the regime where the DM is heavier than half the Higgs mass, and its relic density is determined by annihilation into W, Z and Higgs bosons. We find a large natural region of parameter space where the EWPT is strongly first order, while the lightest inert doublet state typically contributes only 0.1-3% of the total dark matter. Despite this small density, its interactions with nucleons are strong enough to be directly detectable given a factor of 5 improvement over the current sensitivity of XENON100. A 10% decrease in the branching ratio for Higgs decays to two photons is predicted.

If the Higgs boson H couples to a singlet scalar S via lambda_m |H|^2 S^2, a strong electroweak phase transition can be induced through a large potential barrier that exists already at zero temperature. In this case properties of the phase transition can be computed analytically. We show that electroweak baryogenesis can be achieved using CP violation from a dimension-6 operator that couples S to the top-quark mass, suppressed by a new physics scale that can be well above 1 TeV. Moreover the singlet is a dark matter candidate whose relic density is < 3% of the total dark matter density, but which nevertheless interacts strongly enough with nuclei (through Higgs exchange) to be just below the current XENON100 limits. The DM mass is predicted to be in the range 80-160 GeV.

The halo model (HM) describes the inhomogeneous universe as a collection of halos. The full nonlinear power spectrum of the universe is well approximated by the HM, whose prediction can be easily computed without lengthy numerical simulations. This makes the HM a useful tool in cosmology. Here we explore the lensing properties of the HM by use of the stochastic gravitational lensing (sGL) method. We obtain for the case of point sources exact and simple integral expressions for the expected value and variance of the lensing convergence, which encode detailed information about the internal halo properties. In particular a wide array of observational biases can be easily incorporated and the dependence of lensing on cosmology is properly taken into account. This simple setup should be useful for a quick calculation of the power spectrum and the related lensing observables, which can play an important role in the extraction of cosmological parameters from, e.g., SNe observations. Finally, we discuss the probability distribution function of the HM which encodes more information than the first two moments and can more strongly constrain the large-scale structures of the universe. To check the accuracy of our modelling we compare our predictions to the results from the Millennium Simulation.

We revise and extend the stochastic approach to cumulative weak lensing (hereafter the sGL method) first introduced in Ref. [1]. Here we include a realistic halo mass function and density profiles to model the distribution of mass between and within galaxies, galaxy groups and galaxy clusters. We also introduce a modeling of the filamentary large-scale structures and a method to embed halos into these structures. We show that the sGL method naturally reproduces the weak lensing results for the Millennium Simulation. The strength of the sGL method is that a numerical code based on it can compute the lensing probability distribution function for a given inhomogeneous model universe in a few seconds. This makes it a useful tool to study how lensing depends on cosmological parameters and its impact on observations. The method can also be used to simulate the effect of a wide array of systematic biases on the observable PDF. As an example we show how simple selection effects may reduce the variance of observed PDF, which could possibly mask opposite effects from very large scale structures. We also show how a JDEM-like survey could constrain the lensing PDF relative to a given cosmological model. The updated turboGL code is available at turboGL.org.

We reanalyze the supernovae data from the Union Compilation including the weak lensing effects caused by inhomogeneities. We compute the lensing probability distribution function for each background solution described by the parameters Omega_M, Omega_L and w in the presence of inhomogeneities, approximately modeled with a single-mass population of halos. We then perform a likelihood analysis in the space of FLRW-parameters and compare our results with the standard approach. We find that the inclusion of lensing can move the best-fit model significantly towards the cosmic concordance of the flat LCDM model, improving the agreement with the constraints coming from the cosmic microwave background and baryon acoustic oscillations.

We consider a model where electroweak symmetry breaking is driven by Technicolor dynamics with minimal particle content required for walking coupling and saturation of global anomalies. Furthermore, the model features three additional Weyl fermions singlet under Technicolor interactions, which provide for a one-loop unification of the Standard Model gauge couplings. Among these extra matter fields exists a possible candidate for weakly interacting dark matter. We evaluate the relic densities and find that they are sufficient to explain the cosmological observations and avoid the experimental limits from earth-based searches. Hence, we establish a non-supersymmetric framework where hierarchy and naturality problems are solved, coupling constant unification is achieved and a plausible dark matter candidate exists.

We study a superweakly interacting dark matter particle motivated by minimal walking technicolor theories. Our WIMP is a mixture of a sterile state and a state with the charges of a standard model fourth family neutrino. We show that the model can give the right amount of dark matter over a range of the WIMP mass and mixing angle. We compute bounds on the model parameters from the current accelerator data including the oblique corrections to the precision electroweak parameters, as well as from cryogenic experiments, Super-Kamiokande and from the IceCube experiment. We show that consistent dark matter solutions exist which satisfy all current constraints. However, almost the entire parameter range of the model lies within the the combined reach of the next generation experiments.

We study the weak gravitational lensing effects caused by a stochastic distribution of dark matter halos. We develop a simple approach to calculate the magnification probability distribution function which allows us to easily compute the magnitude bias and dispersion for an arbitrary data sample and a given universe model. As an application we consider the effects of single-mass large-scale cosmic inhomogeneities to the SNe magnitude-redshift relation, and conclude that such structures could bias the PDF enough to affect the extraction of cosmological parameters from the limited size of present-day SNe data samples. We also release turboGL, a simple and very fast (<= 1s) Mathematica code based on the method here presented.

We study the impact of cosmic inhomogeneities on the interpretation of SNe observations. We build an inhomogeneous universe model that can confront supernova data and yet is reasonably well compatible with the Copernican Principle. Our model combines a relatively small local void, that gives apparent acceleration at low redshifts, with a meatball model that gives sizeable lensing (dimming) at high redshifts. Together these two elements, which focus on different effects of voids on the data, allow the model to mimic the concordance model.

We study the impact of cosmic inhomogeneities on the interpretation of observations. We build an inhomogeneous universe model without dark energy that can confront supernova data and yet is reasonably well compatible with the Copernican Principle. Our model combines a relatively small local void, that gives apparent acceleration at low redshifts, with a meatball model that gives sizeable lensing (dimming) at high redshifts. Together these two elements, which focus on different effects of voids on the data, allow the model to mimic the concordance model.

We derive quantum kinetic equations for scalar fields undergoing coherent evolution either in time (coherent particle production) or in space (quantum reflection). Our central finding is that in systems with certain space-time symmetries, quantum coherence manifests itself in the form of new spectral solutions for the dynamical 2-point correlation function. This spectral structure leads to a consistent approximation for dynamical equations that describe coherent evolution in presence of decohering collisions. We illustrate the method by solving the bosonic Klein problem and the bound states for the nonrelativistic square well potential. We then compare our spectral phase space definition of particle number to other definitions in the nonequilibrium field theory. Finally we will explicitly compute the effects of interactions to coherent particle production in the case of an unstable field coupled to an oscillating background.

We consider stability properties of spherically symmetric spacetimes of stars in metric f(R) gravity. We stress that these not only depend on the particular model, but also on the specific physical configuration. Typically configurations giving the desired $\gamma_{\rm PPN} \approx 1$ are strongly constrained, while those corresponding to $\gamma_{\rm PPN} \approx 1/2$ are less affected. Furthermore, even when the former are found strictly stable in time, the domain of acceptable static spherical solutions typically shrinks to a point in the phase space. Unless a physical reason to prefer such a particular configuration can be found, this poses a naturalness problem for the currently known metric f(R) models for accelerating expansion of the Universe.

We derive quantum kinetic equations for fermions in a homogeneous time-dependent background in presence of decohering collisions, by use of the Schwinger-Keldysh CTP-formalism. The quantum coherence (between particles and antiparticles) is found to arise from new spectral solutions for the dynamical 2-point correlation function in the mean field limit. The physical density matrix $\rho$ and its dynamics is shown to be necessarily dependent on the extrenous information on the system, and expressions that relate $\rho$ to fundamental coherence functions and fermionic particle and antiparticle numbers are derived. For an interacting system we demonstrate how smooth decoherence effects are induced by collisions. As special applications we study the production of unstable particles during the preheating stage of the inflation and an evolution of an initially quantum $\rho$ towards a statistical limit including decoherence and thermalisation.

A new density matrix and corresponding quantum kinetic equations are introduced for fermions undergoing coherent evolution either in time (coherent particle production) or in space (quantum reflection). A central element in our derivation is finding new spectral solutions for the 2-point Green's functions written in the Wigner representation, that are carrying the information of the quantum coherence. Physically observable density matrix is then defined from the bare singular 2-point function by convoluting it with the extrenous information about the state of the system. The formalism is shown to reproduce familiar results from the Dirac equation approach, like Klein problem and nonlocal reflection from a mass wall. The notion of the particle number in the presence of quantum coherence is shown to be particularily transparent in the current picture. We extend the formalism to the case of mixing fields and show how the usual flavour mixing and oscillation of neutrinos emerges again from a singular shell structure. Finally, we show how the formalism can be extended to include decohering interactions.

We study the interior spacetimes of stars in the Palatini formalism of f(R) gravity and derive a generalized Tolman-Oppenheimer-Volkoff and mass equation for a static, spherically symmetric star. We show that matching the interior solution with the exterior Schwarzschild-De Sitter solution in general gives a relation between the gravitational mass and the density profile of a star, which is different from the one in General Relativity. These modifications become neglible in models for which $\delta F(R) \equiv \partial f/\partial R - 1$ is a decreasing function of R however. As a result, both Solar System constraints and stellar dynamics are perfectly consistent with $f(R) = R - \mu^4/R$.

We study both analytically and numerically the gravitational fields of stars in f(R) gravity theories. We derive the generalized Tolman-Oppenheimer-Volkov equations for these theories and show that in metric f(R) models the Parameterized Post-Newtonian parameter $\gamma_{\rm PPN} = 1/2$ is a robust outcome for a large class of boundary conditions set at the center of the star. This result is also unchanged by introduction of dark matter in the Solar System. We find also a class of solutions with $\gamma_{\rm PPN} \approx 1$ in the metric $f(R)=R-\mu^4/R$ model, but these solutions turn out to be unstable and decay in time. On the other hand, the Palatini version of the theory is found to satisfy the Solar System constraints. We also consider compact stars in the Palatini formalism, and show that these models are not inconsistent with polytropic equations of state. Finally, we comment on the equivalence between f(R) gravity and scalar-tensor theories and show that many interesting Palatini f(R) gravity models can not be understood as a limiting case of a Jordan-Brans-Dicke theory with $\omega \to -3/2$.

We explore in detail a dilatonic scalar-tensor theory of gravity inspired by large extra dimensions, where a radion field from compact extra dimensions gives rise to quintessence in our 4-dimensional world. We show that the model can give rise to other types of cosmologies as well, some more akin to $k$-essence and possibly variants of phantom dark energy. In our model the field (or radius) stabilization arises from quantum corrections to the effective 4D Ricci scalar. We then show that various constraints nearly determine the model parameters, and give an example of a quintessence-type cosmology consistent with observations. We show that the upcoming SNAP-experiment would easily distinguish the present model from a constant $\Lambda$ model with an \emph{equal} amount of dark energy, but that the SNAP-data alone will not be able distinguish it from a $\Lambda$ model with about 5% less dark energy.

We review some astrophysical and cosmological properties and implications of neutrino masses and mixing angles. These include: constraints based on the relic density of neutrinos, limits on their masses and lifetimes, BBN limits on mass parameters, neutrinos and supernovae, and neutrinos and high energy cosmic rays.

We consider the dynamics of fermions with a spatially varying mass which couple to bosons through a Yukawa interaction term and perform a consistent weak coupling truncation of the relevant kinetic equations. We then use a gradient expansion and derive the CP-violating source in the collision term for fermions which appears at first order in gradients. The collisional sources together with the semiclassical force constitute the CP-violating sources relevant for baryogenesis at the electroweak scale. We discuss also the absence of sources at first order in gradients in the scalar equation, and the limitations of the relaxation time approximation.

We study the dynamics of active-sterile neutrino oscillations in the early universe using full momentum-dependent quantum-kinetic equations. These equations are too complicated to allow for an analytical treatment, and numerical solution is greatly complicated due to very pronounced and narrow structures in the momentum variable introduced by resonances. Here we introduce a novel dynamical discretization of the momentum variable which overcomes this problem. As a result we can follow the evolution of neutrino ensemble accurately well into the stable growing phase. Our results confirm the existence of a "chaotic region" of mixing parameters, for which the final sign of the asymmetry, and hence the SBBN prediction of He(4)-abundance cannot be accurately determined.

We perform a systematic gradient expansion on kinetic equations and derive the CP-violating semiclassical force for fermions propagating in presence of a CP-violating wall at a first order electroweak phase transition. The force appears at order h-bar in the flow term of the kinetic equation and agrees with the semiclassical force used for baryogenesis computations. In particular we consider the force for charginos in both the MSSM and NMSSM. We then study the continuity equations for the vector and axial vector currents and stress the role of the latter as the one containing baryogenesis sources. We also show that there is no CP-violating force for bosons to order h-bar in gradient expansion.

We consider active-sterile neutrino oscillations in the early universe in an inhomogeneous isocurvature background. We show that very small initial baryonic seed-inhomogeneities can trigger a growth of very large amplitude spatial fluctuations in lepton asymmetry. Domains of varying asymmetry are observed to persist for a long time despite dissipative effects. The space dependent asymmetry profiles give rise to MSW-resonances within the domain boundaries, enhancing dramatically the equilibration of the sterile neutrino species. According to our one-dimensional toy-model, the effect is so strong that almost the entire parameter space where exponential growth of asymmetry occurs would be ruled out by nucleosynthesis.

We study the fermion propagator in a spatially varying classical background field, and show that, contrary to common wisdom, it may get nontrivial gradient corrections already at the first order in derivative expansion. This occurs whenever the fermion self-energy acquires a spatially (or temporally) varying pseudoscalar term, a simple example of which is given by a complex mass term m(x)= m_R + i gamma_5 m_I. Such effective mass terms arise for example in extensions of the Standard Model during the electroweak transition, and they are crucial in providing the CP-violation necessary for electroweak baryogenesis.

We consider numerically the growth of neutrino asymmetry in active-sterile neutrino oscillations in the early universe. It is shown that the final sign of the asymmetry can be highly sensitive to small variations of the oscillation parameters. We find regions which are completely or partially chaotic, but also regions where the sign remains very robust. The consequences for atmospheric neutrino oscillations and primordial nucleosynthesis are then discussed. In the completely chaotic region the predicted He(4)-abundance has an inherent arbitrariness \Delta Y \simeq 0.01.

We consider the effect of the presence of a hypermagnetic field at the electroweak phase transition. Screening of the Z-component inside a bubble of the broken phase delays the phase transition and makes it stronger first order. We show that the sphaleron constraint can be evaded for m_H up to 100 GeV if a B_Y\gsim 0.3 T^2 exists at the time of the EW phase transition, thus resurrecting the possibility for baryogenesis within the minimal standard model (provided enough CP violation can be obtained). We estimate that for m_H\gsim 100 GeV the Higgs condendsate behaves like a type II superconductor with Z-vortices penetrating the bubble. Also, for such high Higgs masses the minimum B_Y field required for a strong first order phase transition is large enough to render the W-field unstable towards forming a condensate which changes the simple picture of the symmetry breaking.

We consider inhomogeneous big bang nucleosynthesis in light of the present observational situation. Different observations of He-4 and D disagree with each other, and depending on which set of observations one uses, the estimated primordial He-4 corresponds to a lower baryon density in standard big bang nucleosynthesis than what one gets from deuterium. Recent Kamiokande results rule out a favorite particle physics solution to this tension between He-4 and D. Inhomogeneous nucleosynthesis can alleviate this tension, but the more likely solution is systematics in the observations. The upper limit to Omega_b from inhomogeneous nucleosynthesis is higher than in standard nucleosynthesis, given that the distance scale of the inhomogeneity is near the optimal value, which maximizes effects of neutron diffusion. Possible sources of baryon inhomogeneity include the QCD and electroweak phase transitions. The distance scale of the inhomogeneities arising from the electroweak transition is too small for them to have a large effect on nucleosynthesis, but the effect may still be larger than some of the other small corrections recently incorporated to SBBN codes.

We calculate the baryon asymmetry generated at the electroweak phase transition in the minimal supersymmetric standard model, treating the particles in a WKB approximation in the bubble wall background. A set of diffusion equations for the particle species relevant to baryon generation, including source terms arising from the CP violation associated with the complex phase $\delta$ of the $\mu$ parameter, are derived from Boltzmann equations, and solved. The conclusion is that $\delta$ must be $\gsim 0.1$ to generate a baryon asymmetry consistent with nucleosynthesis. We compare our results to several other recent computations of the effect, arguing that some are overestimates.

A computation of nucleosynthesis bounds on the masses of long-lived Dirac and Majorana neutrinos is reviewed. In particular an explicit treatment of the ``differential heating'' of the $\nue$ and $\bnue$ ensembles due to the residual out-of-equilibrium annihilations of decoupled heavy neutrinos is included. The effect is found to be considerably weaker than recently reported by Dolgov et al. For example, the bounds for a Dirac tau neutrino are $\mnt < 0.37$ MeV or $\mnt > 25$ MeV (for $\dNu > 1$), whereas the present laboratory bound is $\mnt < 23.1$ MeV.

We examine in detail how BBN theory is constrained, and what predictions it can make, when using only the most model-independent observational constraints. We avoid the uncertainties and model-dependencies that necessarily arise when solar neighborhood D and \he3 abundances are used to infer primordial D and \he3 via chemical and stellar evolution models. Instead, we use \he4 and \li7, thoroughly examining the effects of possible systematic errors in each. Via a likelihood analysis, we find near perfect agreement between BBN theory and the most model-independent data. Given this agreement, we then {\it assume} the correctness of BBN to set limits on the single parameter of standard BBN, the baryon-to-photon ratio, and to predict the primordial D and \he3 abundances. We also repeat our analysis including recent measurements of D/H from quasar absorption systems and find that the near perfect agreement between theory and observation of the three isotopes, D, \he4 and \li7 is maintained. These results have strong implications for the chemical and stellar evolution of the light elements, in particular for \he3. In addition, our results (especially if the D/H measurements are confirmed) have implications for the stellar depletion of \li7. Finally, we set limits on the number \nnu\ of neutrino flavors, using an analysis which carefully and systematically includes all available experimental constraints. The value \nnu = 3.0 fits best with BBN and a 95\% CL upper limit of \nnu $\la 4$ is established.

We compute the nucleosynthesis bounds on the masses of stable Dirac and Majorana neutrinos by solving an evolution equation network comprising of all neutrino species which in the Dirac case includes different helicity states as separate species. We will not commit ourselves to any particular value of the nucleosynthesis bound on effective number of light neutrino degrees of freedom $N_\nu$, but present all our mass bounds as {\em functions} of $\dNu$. For example, we find that the excluded region in the mass of a Majorana $\mu$- or $\tau$- neutrino, $0.31$ MeV $< m_{\nu}^M < 52$ MeV corresponding to a bound $\dNu < 0.3$ gets relaxed to $0.93$ MeV $< m_{\nu}^M < 31$ MeV if $\dNu < 1.0$ is used instead. For the Dirac neutrinos this latter constraint gives the upper limits (for $T_{\rm QCD} = 100$ MeV): $\mnm < 0.31$ MeV and $\mnt < 0.37$ MeV. Also, the constraint $\dNu <1$ {\it allows} a stable Dirac neutrino with $m^D_{\nu_\tau} > 22$ MeV.

We show that there is no one-loop enhancement of the rate for a light neutrino to decay into a lighter neutrino plus a majoron, contrary to a recent claim. Thus the light neutrinos must satisfy the cosmological bound of having masses less than 35 eV in the singlet majoron model, or else violate the constraint imposed by galaxy formation. In the latter case, $\nu_\tau$ could have a mass between 40 and 500 keV, while satisfying all other cosmological constraints.