Antimatter macroscopic dark matter (macros) refers to a generic class of antimatter dark matter candidates that interact with ordinary matter primarily through annihilation with large cross-sections. A combination of terrestrial, astrophysical, and cosmological observations constrain a portion of the anti-macro parameter space. However, a large region of the parameter space remains unconstrained, most notably for nuclear-dense objects.

Macroscopic dark matter refers to a variety of dark matter candidates that would be expected to (elastically) scatter off of ordinary matter with a large geometric cross-section. A wide range of macro masses $M_X$ and cross-sections $\sigma_X$ remain unprobed. We show that over a wide region within the unexplored parameter space, collisions of a macro with a human body would result in serious injury or death. We use the absence of such unexplained impacts with a well-monitored subset of the human population to exclude a region bounded by $\sigma_X \geq 10^{-8} - 10^{-7}$ cm$^2$ and $M_X < 50$ kg. Our results open a new window on dark matter: the human body as a dark matter detector.

In the Ijjas-Steinhardt cyclic model, the universe passes through phases dominated by radiation, matter, and a dark energy scalar field, with the value of the scale factor increasing with each cycle. Since each cycle terminates in a finite time, it is straightforward to calculate the fraction of time that the universe spends in a state for which the matter and dark energy densities have comparable magnitudes; when this fraction is large, it can be taken as a solution of the coincidence problem. This solution of the coincidence problem requires a relatively short lifetime for each cycle, but unlike in the case of phantom models, there is no fixed upper bound on this lifetime. However, scalar field models satisfying the Swampland conjectures yield sufficiently short lifetimes to provide a satisfactory resolution of the coincidence problem.

We examine the thermal evolution of particle number densities in the early universe when the particles have a finite diffusion length. Assuming that annihilations are impossible when the mean separation of the particles is larger than their diffusion length, we derive a version of the Boltzmann equation for freeze out in this scenario and an approximate solution, accurate to better than 2\%. The effect of a finite diffusion length is to increase the final relic freeze-out abundance over its corresponding value when diffusion effects are ignored. When diffusion is limited only by scattering off of the thermal background, and the annihilation cross section is bounded by unitarity, a significant effect on the freeze-out abundance requires a scattering cross section much larger than the annihilation cross section. A similar effect is demonstrated when the relic particles are produced via the freeze-in mechanism, but in this case the finite diffusion length is due to the scattering of particles that annihilate into the relic particle of interest. For freeze in, the effect of a finite diffusion length is to reduce the final relic particle abundance. The effects of a finite diffusion length are most important when the scattering cross section or the relic mass are very large. While we have not found a particularly compelling example where this would affect previous results, with the current interest in new dark matter candidates it could become an important consideration.

We reexamine $k$-essence dark energy models with a scalar field $\phi$ and a factorized Lagrangian, $\mathcal L = V(\phi)F(X)$, with $X = \frac{1}{2} \nabla_\mu \phi \nabla^\mu \phi.$ A value of the equation of state parameter, $w$, near $-1$ requires either $X \approx 0$ or $dF/dX \approx 0$. Previous work showed that thawing models with $X \approx 0$ evolve along a set of unique trajectories for $w(a)$, while those with $dF/dX \approx 0$ can result in a variety of different forms for $w(a)$. We show that if $dV/d\phi$ is small and $(1/V)(dV/d\phi)$ is roughly constant, then the latter models also converge toward a single unique set of behaviors for $w(a)$, different from those with $X \approx 0$. We derive the functional form for $w(a)$ in this case, determine the conditions on $V(\phi)$ for which it applies, and present observational constraints on this new class of models. We note that $k$-essence models with $dF/dX \approx 0$ correspond to a dark energy sound speed $c_s^2 \approx 0$.

We examine the case of a random isotropic velocity field, in which one of the velocity components (the "radial" component, with magnitude $v_z$) can be measured easily, while measurement of the velocity perpendicular to this component (the "transverse" component, with magnitude $v_T$) is more difficult and requires long-time monitoring. Particularly important examples are the motion of galaxies at cosmological distances and the interpretation of Gaia data on the proper motion of stars in globular clusters and dwarf galaxies. We address two questions: what is the probability distribution of $v_T$ for a given $v_z$, and for what choice of $v_z$ is the expected value of $v_T$ maximized? We show that, for a given $v_z$, the probability that $v_T$ exceeds some value $v_0$ is $p(v_T \ge v_0 | v_z) = {p_z(\sqrt{v_0^2 + v_z^2}})/{p_z(v_z)}$, where $p_z(v_z)$ is the probability distribution of $v_z$. The expected value of $v_T$ is maximized by choosing $v_z$ as large as possible whenever $\ln p_z(\sqrt{v_z})$ has a positive second derivative, and by taking $v_z$ as small as possible when this second derivative is negative.

We reexamine big bang nucleosynthesis with large-scale baryon density inhomogeneities when the length scale of the density fluctuations exceeds the neutron diffusion length ($\sim 10^7-10^8$ cm at BBN), and the amplitude of the fluctuations is sufficiently small to prevent gravitational collapse. In this limit, the final light element abundances can be determined by simply mixing the abundances from regions with different baryon/photon ratios without interactions. We examine gaussian, lognormal, and gamma distributions for the baryon/photon ratio, $\eta $. We find that the deuterium and lithium-7 abundances increase with the RMS fluctuation in $\eta $, while the effect on helium-4 is much smaller. We show that these increases in the deuterium and lithium-7 abundances are a consequence of Jensen's inequality, and we derive analytic approximations for these abundances in the limit of small RMS fluctuations. Observational upper limits on the primordial deuterium abundance constrain the RMS fluctuation in $\eta $ to be less than $17\%$ of the mean value of $\eta $. This provides us with a new limit on the graininess of the early universe.

If the dark energy density asymptotically approaches a nonzero constant, $\rho_{DE} \rightarrow \rho_0$, then its equation of state parameter $w$ necessarily approaches $-1$. The converse is not true; dark energy with $w \rightarrow -1$ can correspond to either $\rho_{DE} \rightarrow \rho_0$ or $\rho_{DE} \rightarrow 0$. This provides a natural division of models with $w \rightarrow -1$ into two distinct classes: asymptotic $\Lambda$ ($\rho_{DE} \rightarrow \rho_0$) and pseudo-$\Lambda$ ($\rho_{DE} \rightarrow 0$). We delineate the boundary between these two classes of models in terms of the behavior of $w(a)$, $\rho_{DE}(a)$, and $a(t)$. We examine barotropic and quintessence realizations of both types of models. Barotropic models with positive squared sound speed and $w \rightarrow -1$ are always asymptotically $\Lambda$; they can never produce pseudo-$\Lambda$ behavior. Quintessence models can correspond to either asymptotic $\Lambda$ or pseudo-$\Lambda$ evolution, but the latter is impossible when the expansion is dominated by a background barotropic fluid. We show that the distinction between asymptotic $\Lambda$ and pseudo-$\Lambda$ models for $w> -1$ is mathematically dual to the distinction between pseudo-rip and big/little rip models when $w < -1$.

The asymmetry between matter and antimatter is key to the existence and nature of our Universe. A measure of the matter - antimatter asymmetry of the Universe is provided by the present value of the universal ratio of baryons (baryons minus antibaryons) to photons (or the ratio of baryons to entropy). The baryon asymmetry parameter is an important physical and cosmological parameter. But how fine tuned is it? A "natural" value for this parameter is zero, corresponding to equal amounts of matter and antimatter. Another, also possibly natural, choice for this dimensionless parameter would be of order unity, corresponding to nearly equal amounts (by number) of matter (and essentially no antimatter) and photons in every comoving volume. However, observations suggest that in the Universe we inhabit the value of this parameter is nonzero, but smaller than this natural value by some nine to ten orders of magnitude. In this contribution we review the evidence that our Universe does not contain equal amounts of matter and antimatter. Any change in the magnitude of the baryon asymmetry parameter necessarily leads to a universe with physical characteristics different from those in our own. The degree of fine tuning in the baryon asymmetry parameter is determined by the width of the range over which it can be varied and still allow for the existence of life. Our results suggest that the baryon asymmetry parameter can be varied over a very wide range without impacting the prospects for life; this result is not suggestive of fine tuning. [abridged]

We study a rapidly-oscillating scalar field with potential $V(\phi) = k|\phi|^n$ nonminally coupled to the Ricci scalar $R$ via a term of the form $(1- 8 \pi G_0 \xi \phi^2) R$ in the action. In the weak coupling limit, we calculate the effect of the nonminimal coupling on the time-averaged equation of state parameter $\gamma = (p + \rho)/\rho$. The change in $\langle \gamma \rangle$ is always negative for $n \ge 2$ and always positive for $n < 0.71$ (which includes the case where the oscillating scalar field could serve as dark energy), while it can be either positive or negative for intermediate values of $n$. Constraints on the time-variation of $G$ force this change to be infinitesimally small at the present time whenever the scalar field dominates the expansion, but constraints in the early universe are not as stringent. The rapid oscillation induced in $G$ also produces an additional contribution to the Friedman equation that behaves like an effective energy density with a stiff equation of state, but we show that, under reasonable assumptions, this effective energy density is always smaller than the density of the scalar field itself.

A change in the fundamental constants of nature or plasma effects in the early universe could stabilize $^8$Be against decay into two $^4$He nuclei. Coc et al. examined this effect on big bang nucleosynthesis as a function of $B_8$, the mass difference between two $^4$He nuclei and a single $^8$Be nucleus, and found no effects for $B_8 \le 100$ keV. Here we examine stable $^8$Be with larger $B_8$ and also allow for a variation in the rate for $^4$He + $^4$He $\longrightarrow$ $^8$Be to determine the threshold for interesting effects. We find no change to standard big bang nucleosynthesis for $B_8 < 1$ MeV. For $B_8 \gtrsim 1$ MeV and a sufficiently large reaction rate, a significant fraction of $^4$He is burned into $^8$Be, which fissions back into $^4$He when $B_8$ assumes its present-day value, leaving the primordial $^4$He abundance unchanged. However, this sequestration of $^4$He results in a decrease in the primordial $^7$Li abundance. Primordial abundances of $^7$Li consistent with observationally-inferred values can be obtained for reaction rates similar to those calculated for the present-day (unbound $^8$Be) case. Even for the largest binding energies and largest reaction rates examined here, only a small fraction of $^8$Be is burned into heavier elements, consistent with earlier studies. There is no change in the predicted deuterium abundance for any model we examined.

A change in the fundamental constants of nature or plasma effects in the early universe could stabilize $^8$Be against decay into two $^4$He nuclei. Coc et al. examined this effect on big bang nucleosynthesis as a function of $B_8$, the mass difference between two $^4$He nuclei and a single $^8$Be nucleus, and found no effects for $B_8 \le 100$ keV. Here we examine stable $^8$Be with larger $B_8$ and also allow for a variation in the rate for $^4$He + $^4$He $\longrightarrow$ $^8$Be to determine the threshold for interesting effects. We find no change to standard big bang nucleosynthesis for $B_8 < 1$ MeV. For $B_8 \gtrsim 1$ MeV and a sufficiently large reaction rate, a significant fraction of $^4$He is burned into $^8$Be, which fissions back into $^4$He when $B_8$ assumes its present-day value, leaving the primordial $^4$He abundance unchanged. However, this sequestration of $^4$He results in a decrease in the primordial $^7$Li abundance. Primordial abundances of $^7$Li consistent with observationally-inferred values can be obtained for reaction rates similar to those calculated for the present-day (unbound $^8$Be) case. Even for the largest binding energies and largest reaction rates examined here, only a small fraction of $^8$Be is burned into heavier elements, consistent with earlier studies. There is no change in the predicted deuterium abundance for any model we examined.

We consider a first order formulation of relativistic fluids with bulk viscosity based on a stress-energy tensor introduced by Lichnerowicz. Choosing a barotropic equation of state, we show that this theory satisfies basic physical requirements and, under the further assumption of vanishing vorticity, that the equations of motion are causal, both in the case of a fixed background and when the equations are coupled to Einstein's equations. Furthermore, Lichnerowicz's proposal does not fit into the general framework of first order theories studied by Hiscock and Lindblom, and hence their instability results do not apply. These conclusions apply to the full-fledged non-linear theory, without any equilibrium or near equilibrium assumptions. Similarities and differences between the approach here explored and other theories of relativistic viscosity, including the Mueller-Israel-Stewart formulation, are addressed. Cosmological models based on the Lichnerowicz stress-energy tensor are studied. Finally, it is also shown how the present model can be generalized to a second order formulation.

We consider Big Bang nucleosynthesis and the cosmic microwave background when both the neutrino temperature and neutrino number are allowed to vary from their standard values. The neutrino temperature is assumed to differ from its standard model value by a fixed factor from Big Bang nucleosynthesis up to the present. In this scenario, the effective number of relativistic degrees of freedom, $N_{\rm eff}^{\rm CMB}$, derived from observations of the cosmic microwave background is not equal to the true number of neutrinos, $N_\nu$. We determine the element abundances predicted by Big Bang nucleosynthesis as a function of the neutrino number and temperature, converting the latter to the equivalent value of $N_{\rm eff}^{\rm CMB}$. We find that a value of $N_{\rm eff}^{\rm CMB} \approx 3$ can be made consistent with $N_\nu = 4$ with a decrease in the neutrino temperature of $\sim 5\%$, while $N_\nu = 5$ is excluded for any value of $N_{\rm eff}^{\rm CMB}$. No observationally-allowed values for $N_{\rm eff}^{\rm CMB}$ and $N_\nu$ can solve the lithium problem.

We examine a quintessence model with a modified exponential potential given by $V(\phi) = V_0(1+e^{-\lambda \phi})$. Unlike quintessence with a standard exponential potential, our model can yield an acceptable accelerated expansion at late times, while producing a distinct "early dark energy" signature at high redshift. We determine the evolution of the equation of state parameter, $w_\phi$, and the density parameter, $\Omega_\phi$, as a function of the scale factor. The strongest constraints on the model come from cosmic microwave background observations rather than supernova data. The former give the limit $\lambda > 13$. This model predicts a value of the effective neutrino number during Big Bang nucleosynthesis larger than the standard model value. It also provides a partial solution to the coincidence problem, in that the ratio of the quintessence energy density is always within a few orders of magnitude of the background radiation or matter density from the early universe up to the present, but it does not explain why the accelerated expansion is beginning near the present day, suggesting that these two different ways of characterizing the coincidence problem are not entirely equivalent.

Grand unified theories based on large groups (with rank greater or equal to 6 are a natural context for dark matter models. They contain Standard-Model-singlet fermions that could be dark matter candidates, and can contain new non-abelian interactions whose sphalerons convert baryons, leptons, and dark matter into each other, "cogenerating" a dark matter asymmetry comparable to the baryon asymmetry. In this paper it is shown that the same non-abelian interactions can "pre-annihilate" the symmetric component of heavy dark matter particles, which then decay late into light stable dark matter particles that inherit their asymmetry. It is shown that such decays can come from d=5 operators that are Planck or GUT suppressed. We derive constraints on such models and present a simple realization based on the group SU(7).

We apply the Alcock-Paczynski (AP) test to the stacked voids identified using the large-scale structure galaxy catalog from the Baryon Oscillation Spectroscopic Survey (BOSS). This galaxy catalog is part of the Sloan Digital Sky Survey (SDSS) Data Release 12 and is the final catalog of SDSS-III. We also use 1000 mock galaxy catalogs that match the geometry, density, and clustering properties of the BOSS sample in order to characterize the statistical uncertainties of our measurements and take into account systematic errors such as redshift space distortions. For both BOSS data and mock catalogs, we use the ZOBOV algorithm to identify voids, we stack together all voids with effective radii of 30-100Mpc/h in the redshift range 0.43-0.7, and we accurately measure the shape of the stacked voids. Our tests with the mock catalogs show that we measure the stacked void ellipticity with a statistical precision of 2.6%. We find that the stacked voids in redshift space are slightly squashed along the line of sight, which is consistent with previous studies. We repeat this measurement of stacked void shape in the BOSS data assuming several values of Omega_m within the flat LCDM model, and we compare to the mock catalogs in redshift space in order to perform the AP test. We obtain a constraint of $\Omega_m = 0.38^{+0.18}_{-0.15}$ at the 68% confidence level from the AP test. We discuss the various sources of statistical and systematic noise that affect the constraining power of this method. In particular, we find that the measured ellipticity of stacked voids scales more weakly with cosmology than the standard AP prediction, leading to significantly weaker constraints. We discuss how AP constraints will improve in future surveys with larger volumes and densities.

We derive general conditions for the existence of stable scaling solutions for the evolution of noncanonical quintessence, with a Lagrangian of the form $\mathcal{L}(X,\phi)=X^{\alpha}-V(\phi)$, for power-law and exponential potentials when the expansion is dominated by a background barotropic fluid. Our results suggest that in most cases, noncanonical quintessence with such potentials does not yield interesting models for the observed dark energy. When the scaling solution is not an attractor, there is a wide range of model parameters for which the evolution asymptotically resembles a zero-potential solution with equation of state parameter $w = 1/(2\alpha -1)$, and oscillatory solutions are also possible for positive power-law potentials; we derive the conditions on the model parameters which produce both types of behavior. We investigate thawing noncanonical models with a nearly-flat potential and derive approximate expressions for the evolution of $w(a)$. These forms for $w(a)$ differ in a characteristic way from the corresponding expressions for canonical quintessence.

We present a cosmic void catalog using the large-scale structure galaxy catalog from the Baryon Oscillation Spectroscopic Survey (BOSS). This galaxy catalog is part of the Sloan Digital Sky Survey (SDSS) Data Release 12 and is the final catalog of SDSS-III. We take into account the survey boundaries, masks, and angular and radial selection functions, and apply the ZOBOV void finding algorithm to the galaxy catalog. We identify a total of 10,643 voids. After making quality cuts to ensure that the voids represent real underdense regions, we obtain 1,228 voids with effective radii spanning the range 20-100Mpc/h and with central densities that are, on average, 30% of the mean sample density. We release versions of the catalogs both with and without quality cuts. We discuss the basic statistics of voids, such as their size and redshift distributions, and measure the radial density profile of the voids via a stacking technique. In addition, we construct mock void catalogs from 1000 mock galaxy catalogs, and find that the properties of BOSS voids are in good agreement with those in the mock catalogs. We compare the stellar mass distribution of galaxies living inside and outside of the voids, and find no significant difference. These BOSS and mock void catalogs are useful for a number of cosmological and galaxy environment studies.

We consider cosmological solutions to general relativity with a single barotropic fluid, where the pressure is a general function of the density, $p = f(\rho)$. We derive conditions for static and oscillating solutions and provide examples, extending earlier work to these simpler and more general single-fluid cosmologies. Generically we expect such solutions to suffer from instabilities, through effects such as quantum fluctuations or tunneling to zero size. We also find a classical instability ("no-go" theorem) for oscillating solutions of a single barotropic perfect fluid due to a necessarily negative squared sound speed.

We examine the Chevallier-Polarski-Linder (CPL) parametrization, in the context of quintessence and barotropic dark energy models, to determine the subset of such models to which it can provide a good fit. The CPL parametrization gives the equation of state parameter $w$ for the dark energy as a linear function of the scale factor $a$, namely $w = w_0 + w_a(1-a)$. In the case of quintessence models, we find that over most of the $w_0$, $w_a$ parameter space the CPL parametrization maps onto a fairly narrow form of behavior for the potential $V(\phi)$, while a one-dimensional subset of parameter space, for which $w_a = \kappa (1+w_0)$, with $\kappa$ constant, corresponds to a wide range of functional forms for $V(\phi)$. For barotropic models, we show that the functional dependence of the pressure on the density, up to a multiplicative constant, depends only on $w_i = w_a + w_0$ and not on $w_0$ and $w_a$ separately. Our results suggest that the CPL parametrization may not be optimal for testing either type of model.

We calculate finite-temperature corrections to the decay rate of a generic neutral (pseudo)scalar particle that decays into (pseudo)scalars or fermion-antifermion pairs. The ratio of the finite-temperature decay rate to the zero-temperature decay rate is presented. Thermal effects are largest in the limit where the decaying particle is nonrelativistic but with a mass well below the background temperature, but significant effects are possible even when we relax the former assumption. Thermal effects are reduced for the case of nonzero momentum of the decaying particle. We discuss cosmological scenarios under which significant finite-temperature corrections to the decay rate can be achieved.

We examine quintessence models for dark energy in which the scalar field, $\phi$, evolves near the vicinity of a local maximum or minimum in the potential $V(\phi)$, so that $V(\phi)$ be approximated by a quadratic function of $\phi$ with no linear term. We generalize previous studies of this type by allowing the initial value of $d \phi/dt$ to be nonzero. We derive an analytic approximation for $w(a)$ and show that it is in excellent agreement with numerical simulations for a variety of scalar field potentials having local minima or maxima. We derive an upper bound on the present-day value of $w$ as a function of the other model parameters and present representative limits on these models from observational data. This work represents a final generalization of previous studies using linear or quadratic approximations for $V(\phi)$.

We examine the cosmological consequences of an alternative to the standard expression for bulk viscosity, one which was proposed to avoid the propagation of superluminal signals without the necessity of extending the space of variables of the theory. The Friedmann equation is derived for this case, along with an expression for the effective pressure. We find solutions for the evolution of the density of a viscous component, which differs markedly from the case of conventional Eckart theory; our model evolves toward late-time phantom-like behavior with a future singularity. Entropy production is addressed, and some similarities and differences to approaches based on the Mueller-Israel-Stewart theory are discussed.

We use cosmological N-body simulations to investigate whether measurements of the moments of large-scale structure can yield constraints on primordial non-Gaussianity. We measure the variance, skewness, and kurtosis of the evolved density field from simulations with Gaussian and three different non-Gaussian initial conditions: a local model with f_NL=100, an equilateral model with f_NL=-400, and an orthogonal model with f_NL=-400. We show that the moments of the dark matter density field differ significantly between Gaussian and non-Gaussian models. We also make the measurements on mock galaxy catalogs that contain galaxies with clustering properties similar to those of luminous red galaxies (LRGs). We find that, in the case of skewness and kurtosis, galaxy bias reduces the detectability of non-Gaussianity, though we can still clearly discriminate between different models in our simulation volume. However, in the case of the variance, galaxy bias greatly amplifies the detectability of non-Gaussianity. In all cases we find that redshift distortions do not significantly affect the detectability. When we restrict our measurements to volumes equivalent to the Sloan Digital Sky Survey II (SDSS-II) or Baryon Oscillation Spectroscopic Survey (BOSS) samples, the probability of detecting a departure from the Gaussian model is high by using measurements of the variance, but very low by using only skewness and kurtosis measurements. We find that skewness and kurtosis measurements are never likely to yield useful constraints on primordial non-Gaussianity, but future surveys should be large enough to place meaningful constraints using measurements of the galaxy variance.(Abridged)

The anapole moment is the only allowed electromagnetic moment for Majorana fermions. Fermionic dark matter acquiring an anapole can have a standard thermal history and be consistent with current direct detection experiments. In this paper, we calculate the collider monojet signatures of anapole dark matter and show that the current LHC results exclude anapole dark matter with mass less than 100 GeV, for an anapole coupling that leads to the correct thermal relic abundance.

We examine models in which the accelerated expansion of the universe is driven by a scalar field rolling near an inflection point in the potential. For the simplest such models, in which the potential is of the form V(\phi) = V_0 + V_3 (\phi-\phi_0)^3, the scalar field can either evolve toward \phi = \phi_0 at late times, yielding an asymptotic de Sitter expansion, or it can transition through the inflection point, producing a transient period of acceleration. We determine the parameter ranges which produce each of these two possibilities and also map out the region in parameter space for which the equation of state of the scalar field is close to -1 at all times up to the present, mimicking \LambdaCDM. We show that the latter can be consistent with either eternal or transient acceleration. More complicated inflection point models are also investigated.

We consider dark matter (DM) that interacts with ordinary matter exclusively through an electromagnetic anapole, which is the only allowed electromagnetic form factor for Majorana fermions. We show that unlike DM particles with an electric or magnetic dipole moment, anapole dark matter particles annihilate exclusively into fermions via purely p-wave interactions, while tree-level annihilations into photons are forbidden. We calculate the anapole moment needed to produce a thermal relic abundance in agreement with cosmological observations, and show that it is consistent with current XENON100 detection limits on the DM-nucleus cross-section for all masses, while lying just below the detection threshold for a mass ~ 30-40 GeV.

Short baseline neutrino experiments such as LSND and MiniBooNE seem to suggest the existence of light sterile neutrinos. Meanwhile, current cosmic microwave background (CMB) and big bang nucleosynthesis (BBN) measurements place an upper bound on the effective number of light neutrinos, $N_{eff}$ and the PLANCK satellite will measure $N_{eff}$ to a much higher accuracy and further constrain the number of sterile neutrinos allowed. We demonstrate that if an MeV dark matter particle couples more strongly to electrons and/or photons than to neutrinos, then p-wave annihilation after neutrino decoupling can reduce the value of $N_{eff}$ inferred from BBN and PLANCK. This mechanism can accommodate two eV sterile neutrinos even if PLANCK observes $N_{eff}$ as low as the standard model theoretical value of 3.046, and a large neutrino asymmetry is not needed to obtain the correct primordial element abundances. The dark matter annihilation also weakens the cosmological upper bounds on the neutrino masses, and we derive a relationship between the change in these bounds and the corresponding change in $N_{eff}$. Dark matter with an electric dipole moment or anapole moment is a natural candidate that exhibits the desired properties for this mechanism. Coincidentally, a dark matter particle with these properties and lighter than 3 MeV is precisely one that can explain the 511 keV gamma-ray line observed by INTEGRAL. We show that the addition of two eV sterile neutrinos allows this kind of dark matter to be lighter than 3 MeV, which is otherwise ruled out by the CMB bound on $N_{eff}$ if only active neutrinos are considered.

Thermal dark matter that couples more strongly to electrons and photons than to neutrinos will heat the electron-photon plasma relative to the neutrino background if it becomes nonrelativistic after the neutrinos decouple from the thermal background. This results in a reduction in N_eff below the standard-model value, a result strongly disfavored by current CMB observations. Taking conservative lower bounds on N_eff and on the decoupling temperature of the neutrinos, we derive a bound on the dark matter particle mass of m_\chi > 3-9 MeV, depending on the spin and statistics of the particle. For p-wave annihilation, our limit on the dark matter particle mass is stronger than the limit derived from distortions to the CMB fluctuation spectrum produced by annihilations near the epoch of recombination.

We examine cyclic phantom models for the universe, in which the universe is dominated sequentially by radiation, matter, and a phantom dark energy field, followed by a standard inflationary phase. Since this cycle repeats endlessly, the Universe spends a substantial portion of its lifetime in a state for which the matter and dark energy densities have comparable magnitudes, thus ameliorating the coincidence problem. We calculate the fraction of time that the universe spends in such a coincidental state and find that it is nearly the same as in the case of a phantom model with a future big rip. In the limit where the dark energy equation of state parameter, w, is close to -1, we show that the fraction of time, f, for which the ratio of the dark energy density to the matter density lies between r_1 and r_2, is f = -(1+w) ln [(\sqrt{r_2} + \sqrt{1+r_2})/(\sqrt{r_1} + \sqrt{1+r_1})].

Dark energy models with various scenarios of evolution are considered from the viewpoint of the formalism for the equation of state. It is shown that these models are compatible with current astronomical data. Some of the models presented here evolve arbitrarily close to $\Lambda$CDM up to the present, but diverge in the future into a number of different possible asymptotic states, including asymptotic de-Sitter (pseudo-rip) evolution, little rips with disintegration of bound structures, and various forms of finite-time future singularities. Therefore it is impossible from observational data to determine whether the universe will end in a future singularity or not. We demonstrate that the models under consideration are stable for a long period of time (billions of years) before entering a Little Rip/Pseudo-Rip induced dissolution of bound structures or before entering a soft finite-time future singularity. Finally, the physical consequences of Little Rip, Type II, III and Big Crush singularities are briefly compared.

If we assume that the cosmic energy density will remain constant or strictly increase in the future, then the possible fates for the universe can be divided into four categories based on the time asymptotics of the Hubble parameter H(t): the cosmological constant, for which H(t) = constant, the big rip, for which H(t) goes to infinity at finite time, the little rip, for which H(t) goes to infinity as time goes to infinity, and the pseudo-rip, for which H(t) goes to a constant as time goes to infinity. Here we examine the last of these possibilities in more detail. We provide models that exemplify the pseudo-rip, which is an intermediate case between the cosmological constant and the little rip. Structure disintegration in the pseudo-rip depends on the model parameters. We show that pseudo-rip models for which the density and Hubble parameter increase monotonically can produce an inertial force which does not increase monotonically, but instead peaks at a particular future time and then decreases.

Cosmic microwave background (CMB) observations suggest the possibility of an extra dark radiation component, while the current evidence from big bang nucleosynthesis (BBN) is more ambiguous. Dark radiation from a decaying particle can affect these two processes differently. Early decays add an additional radiation component to both the CMB and BBN, while late decays can alter the radiation content seen in the CMB while having a negligible effect on BBN. Here we quantify this difference and explore the intermediate regime by examining particles decaying during BBN, i.e., particle lifetimes \tau_X satisfying 0.1 sec < \tau_X < 1000 sec. We calculate the change in the effective number of neutrino species, N_{eff}, as measured by the CMB, \Delta N_{CMB}, and the change in the effective number of neutrino species as measured by BBN, \Delta N_{BBN}, as a function of the decaying particle initial energy density and lifetime, where \Delta N_{BBN} is defined in terms of the number of additional two-component neutrinos needed to produce the same change in the primordial helium-4 abundance as our decaying particle. As expected, for short lifetimes (\tau_X < 0.1 sec), the particles decay before the onset of BBN, and \Delta N_{CMB} = \Delta N_{BBN}, while for long lifetimes (\tau_X > 1000 sec), \Delta N_{BBN} is dominated by the energy density of the nonrelativistic particles before they decay, so that \Delta N_{BBN} remains nonzero and becomes independent of the particle lifetime. By varying both the particle energy density and lifetime, one can obtain any desired combination of \Delta N_{BBN} and \Delta N_{CMB}, subject to the constraint that \Delta N_{CMB} >= \Delta N_{BBN}. We present limits on the decaying particle parameters derived from observational constraints on \Delta N_{CMB} and \Delta N_{BBN}.

We examine in more detail specific models which yield a little rip cosmology, i.e., a universe in which the dark energy density increases without bound but the universe never reaches a finite-time singularity. We derive the conditions for the little rip in terms of the inertial force in the expanding universe and present two representative models to illustrate in more detail the difference between little rip models and those which are asymptotically de Sitter. We derive conditions on the equation of state parameter of the dark energy to distinguish between the two types of models. We show that coupling between dark matter and dark energy with a little rip equation of state can alter the evolution, changing the little rip into an asymptotic de Sitter expansion. We give conditions on minimally-coupled phantom scalar field models and on scalar-tensor models that indicate whether or not they correspond to a little rip expansion. We show that, counterintuitively, despite local instability, a little-rip cosmology has an infinite lifetime.

We examine the evolution of quintessence models with potentials satisfying (V'/V)^2<<1 and V"/V<<1, in the case where the initial field velocity is nonzero. We derive an analytic approximation for the evolution of the equation of state parameter, w, for the quintessence field. We show that such models are characterized by an initial rapid freezing phase, in which the equation of state parameter w decreases with time, followed by slow thawing evolution, for which w increases with time. These models resemble constant-V models at early times but diverge at late times. Our analytic approximation gives results in excellent agreement with exact numerical evolution.

We examine models in which the dark energy density increases with time (so that the equation-of-state parameter w satisfies w < -1), but w approaches -1 asymptotically, such that there is no future singularity. We refine previous calculations to determine the conditions necessary to produce this evolution. Such models can display arbitrarily rapid expansion in the near future, leading to the destruction of all bound structures (a "little rip"). We determine observational constraints on these models and calculate the point at which the disintegration of bound structures occurs. For the same present-day value of w, a big rip with constant w disintegrates bound structures earlier than a little rip.

Models that lead to a cosmological stiff fluid component, with a density $\rho_S$ that scales as $a^{-6}$, where $a$ is the scale factor, have been proposed recently in a variety of contexts. We calculate numerically the effect of such a stiff fluid on the primordial element abundances. Because the stiff fluid energy density decreases with the scale factor more rapidly than radiation, it produces a relatively larger change in the primordial helium-4 abundance than in the other element abundances, relative to the changes produced by an additional radiation component. We show that the helium-4 abundance varies linearly with the density of the stiff fluid at a fixed fiducial temperature. Taking $\rho_{S10}$ and $\rho_{R10}$ to be the stiff fluid energy density and the standard density in relativistic particles, respectively, at $T = 10$ MeV, we find that the change in the primordial helium abundance is well-fit by $\Delta Y_p = 0.00024(\rho_{S10}/\rho_{R10})$. The changes in the helium-4 abundance produced by additional radiation or by a stiff fluid are identical when these two components have equal density at a "pivot temperature", $T_*$, where we find $T_* = 0.55$ MeV. Current estimates of the primordial $^4$He abundance give the constraint on a stiff fluid energy density of $\rho_{S10}/\rho_{R10} < 30$.

A $\Lambda$CDM model with dark matter that decays into inert relativistic energy on a timescale longer than the Hubble time will produce an expansion history that can be misinterpreted as stable dark matter with time-varying dark energy. We calculate the corresponding spurious equation of state parameter, $\widetilde w_\phi$, as a function of redshift, and show that the evolution of $\widetilde w_\phi$ depends strongly on the assumed value of the dark matter density, erroneously taken to scale as $a^{-3}$. Depending on the latter, one can obtain models that mimic quintessence ($\widetilde w_\phi > -1$), phantom models ($\widetilde w_\phi < -1$) or models in which the equation of state parameter crosses the phantom divide, evolving from $\widetilde w_\phi > -1$ at high redshift to $\widetilde w_\phi < -1$ at low redshift. All of these models generically converge toward $\widetilde w_\phi \approx -1$ at the present. The degeneracy between the $\Lambda$CDM model with decaying dark matter and the corresponding spurious quintessence model is broken by the growth of density perturbations.

We propose a class of simple dark energy models which predict a late-time dark radiation component and a distinctive time-dependent equation of state $w(z)$ for redshift $z < 3$. The dark energy field can be coupled strongly enough to Standard Model particles to be detected in colliders, and the model requires only modest additional particle content and little or no fine-tuning other than a new energy scale of order milli-electron volts.

We reexamine the evolution of thermal relic particle abundances for the case where the interaction rate depends on the particle velocities. For the case of Sommerfeld enhancement, we show that the standard analytic approximation, modified in a straightforward way, provides an estimate of the relic particle abundance that is accurate to within 10% (in comparison to less than 1% error for the non-Sommerfeld-enhanced case). We examine the effect of kinetic decoupling on relic particle abundances when the interaction rate depends on the velocity. For the case of pure p-wave annihilation, the effect of kinetic decoupling is an increase in the relic abundance, but the effect is negligible when the kinetic decoupling temperature is much less than the chemical decoupling temperature. For the case of Sommerfeld-enhanced s-wave annihilations, after kinetic decoupling occurs, annihilations continue to change the particle abundance down to arbitrarily low temperatures, until either matter domination begins or the Sommerfeld effect cuts off. We derive analytic approximations to give the final relic particle abundances for both of these cases.

Any multivariate distribution can be uniquely decomposed into marginal (1-point) distributions, and a function called the copula, which contains all of the information on correlations between the distributions. The copula provides an important new methodology for analyzing the density field in large-scale structure. We derive the empirical 2-point copula for the evolved dark matter density field. We find that this empirical copula is well-approximated by a Gaussian copula. We consider the possibility that the full n-point copula is also Gaussian and describe some of the consequences of this hypothesis. Future directions for investigation are discussed.

We derive slow-roll conditions for thawing k-essence with a separable Lagrangian $p(X,\phi)=F(X)V(\phi)$. We examine the evolution of the equation of state parameter, $w$, as a function of the scale factor $a$, for the case where $w$ is close to -1. We find two distinct cases, corresponding to $X \approx 0$ and $F_X \approx 0$, respectively. For the case where $X\approx0$ the evolution of $\phi$ and hence $w$ is described by only two parameters, and $w(a)$ is model-independent and coincides with similar behavior seen in thawing quintessence models. This result also extends to non-separable Lagrangians where $X\approx0$. For the case $F_X \approx 0$, an expression is derived for $w(a)$, but this expression depends on the potential $V(\phi)$, so there is no model-independent limiting behavior. For the $X \approx 0$ case, we derive observational constraints on the two parameters of the model, $w_0$ (the present-day value of $w$), and the $K$, which parametrizes the curvature of the potential. We find that the observations sharply constrain $w_0$ to be close to -1, but provide very poor constraints on $K$.

We examine the class of barotropic fluid models of dark energy, in which the pressure is an explicit function of the density, p = f(\rho). Through general physical considerations we constrain the asymptotic past and future behaviors and show that this class is equivalent to the sum of a cosmological constant and a decelerating perfect fluid, or "aether", with w_{AE}\ge0. Barotropic models give substantially disjoint predictions from quintessence, except in the limit of \LambdaCDM. They are also interesting in that they simultaneously can ameliorate the coincidence problem and yet "predict" a value of w\approx-1.

We examine dark energy models in which a phantom field $\phi$ is rolling near a local minimum of its potential $V(\phi)$.We require that $(1/V)(dV/d\phi) \ll 1$, but $(1/V)(d^2 V/d\phi^2)$ can be large. Using techniques developed in the context of hilltop quintessence, we derive a general expression for $w$ as a function of the scale factor, and as in the hilltop case, we find that the dynamics of the field depend on the value of $(1/V)(d^2 V/d\phi^2)$ near the mimimum. Our general result gives a value for $w$ that is within 1% of the true (numerically-derived) value for all of the particular cases examined. Our expression for $w(a)$ reduces to the previously-derived phantom slow-roll result of Sen and Scherrer in the limit where the potential is flat, $(1/V)(dV/d\phi) \ll 1$.

We propose a simple dark energy model with the following properties: the model predicts a late-time dark radiation component that is not ruled out by current observational data, but which produces a distinctive time-dependent equation of state w(z) for z < 3. The dark energy field can be coupled strongly enough to Standard Model particles to be detected in colliders, and the model requires only modest additional particle content and little or no fine-tuning other than a new energy scale of order milli-electron volts.

We examine the time variation of a previously-uninvestigated fundamental dimensionless constant. Constraints are placed on this time variation using historical measurements. A model is presented for the time variation, and it is shown to lead to an accelerated expansion for the universe. Directions for future research are discussed.

We examine dark energy models in which a quintessence or a phantom field, $\phi$, rolls near the vicinity of a local minimum or maximum, respectively, of its potential $V(\phi)$. Under the approximation that $(1/V)(dV/d\phi) \ll 1$, [although $(1/V)(d^2 V/d\phi^2)$ can be large], we derive a general expression for the equation of state parameter $w$ as a function of the scale factor for these models. The dynamics of the field depends on the value of $(1/V)(d^2 V/d\phi^2)$ near the extremum, which describes the potential curvature. For quintessence models, when $(1/V)(d^2 V/d\phi^2)<3/4$ at the potential minimum, the equation of state parameter $w(a)$ evolves monotonically, while for $(1/V)(d^2 V/d\phi^2)>3/4$, $w(a)$ has oscillatory behavior. For phantom fields, the dividing line between these two types of behavior is at $(1/V)(d^2 V/d\phi^2) = -3/4$. Our analytical expressions agree within 1% with the exact (numerically-derived) behavior, for all of the particular cases examined, for both quintessence and phantom fields. We present observational constraints on these models.

We examine hilltop quintessence models, in which the scalar field is rolling near a local maximum in the potential, and w is close to -1. We first derive a general equation for the evolution of the scalar field in the limit where w is close to -1. We solve this equation for the case of hilltop quintessence to derive w as a function of the scale factor; these solutions depend on the curvature of the potential near its maximum. Our general result is in excellent agreement (delta w < 0.5%) with all of the particular cases examined. It works particularly well (delta w < 0.1%) for the pseudo-Nambu-Goldstone Boson potential. Our expression for w(a) reduces to the previously-derived slow-roll result of Sen and Scherrer in the limit where the curvature goes to zero. Except for this limiting case, w(a) is poorly fit by linear evolution in a.

Oscillating scalar fields, with an oscillation frequency much greater than the expansion rate, have been proposed as models for dark energy. We examine these models, with particular emphasis on the evolution of the ratio of the oscillation frequency to the expansion rate. We show that this ratio always increases with time if the dark energy density declines less rapidly than the background energy density. This allows us to classify oscillating dark energy models in terms of the epoch at which the oscillation frequency exceeds the expansion rate, which is effectively the time at which rapid oscillations begin. There are three basic types of behavior: early oscillation models, in which oscillations begin during the matter-dominated era, late oscillation models, in which oscillations begin after scalar-field domination, and non-oscillating models. We examine a representative set of models (those with power-law potentials) and determine the parameter range giving acceptable agreement with the supernova observations. We show that a subset of all three classes of models can be consistent with the observational data.

Observational limits on the high-energy neutrino background have been used to place general constraints on dark matter that annihilates only into standard model particles. Dark matter particles that annihilate into neutrinos will also inevitably branch into electromagnetic final states through higher-order tree and loop diagrams that give rise to charged leptons, and these charged particles can transfer their energy into photons via synchrotron radiation or inverse Compton scattering. In the context of effective field theory, we calculate the loop-induced branching ratio to charged leptons and show that it is generally quite large, typically >1%, when the scale of the dark matter mass exceeds the electroweak scale, M_W. For a branching fraction >3%, the synchrotron radiation bounds on dark matter annihilation are currently stronger than the corresponding neutrino bounds in the interesting mass range from 100 GeV to 1 TeV. For dark matter masses below M_W, our work provides a plausible framework for the construction of a model for "neutrinos only" dark matter annihilations.

We examine phantom dark energy models produced by a field with a negative kinetic term and a potential that satisfies the slow roll conditions: [(1/V)(dV/dphi)]^2 << 1 and (1/V)(d^2 V/dphi^2) << 1. Such models provide a natural mechanism to produce an equation of state parameter, w, slightly less than -1 at present. Using techniques previously applied to quintessence, we show that in this limit, all such phantom models converge to a single expression for w(a), which is a function only of the present-day values of Omega_phi and w. This expression is identical to the corresponding behavior of w(a) for quintessence models in the same limit. At redshifts z < 1, this limiting behavior is well fit by the linear parametrization, w=w_0 + w_a(1-a), with w_a \approx -1.5(1+w_0).

The thawing quintessence model with a nearly flat potential provides a natural mechanism to produce an equation of state parameter, w, close to -1 today. We examine the behavior of such models for the case in which the potential satisfies the slow roll conditions: [(1/V)(dV/dphi)]^2 << 1 and (1/V)(d^2 V/dphi^2) << 1, and we derive the analog of the slow-roll approximation for the case in which both matter and a scalar field contribute to the density. We show that in this limit, all such models converge to a unique relation between 1+w, Omega_phi, and the initial value of (1/V)(dV/dphi). We derive this relation, and use it to determine the corresponding expression for w(a), which depends only on the present-day values for w and Omega_phi. For a variety of potentials, our limiting expression for w(a) is typically accurate to within delta w < 0.005 for w<-0.9. For redshift z < 1, w(a) is well-fit by the Chevallier-Polarski-Linder parametrization, in which w(a) is a linear function of a.

Big Bang nucleosynthesis requires a fine balance between equations of state for photons and relativistic fermions. Several corrections to equation of state parameters arise from classical and quantum physics, which are derived here from a canonical perspective. In particular, loop quantum gravity allows one to compute quantum gravity corrections for Maxwell and Dirac fields. Although the classical actions are very different, quantum corrections to the equation of state are remarkably similar. To lowest order, these corrections take the form of an overall expansion-dependent multiplicative factor in the total density. We use these results, along with the predictions of Big Bang nucleosynthesis, to place bounds on these corrections.

We examine flat models containing a dark matter component and an arbitrary dark energy component, subject only to the constraint that the dark energy satisfies the weak energy condition. We determine the constraints that these conditions place on the evolution of the Hubble parameter with redshift, H(z), and on the scaling of the coordinate distance with redshift, r(z). Observational constraints on H(z) are used to derive an upper bound on the current matter density. We demonstrate how the weak energy condition constrains fitting functions for r(z).

We demonstrate that as we extrapolate the current $\Lambda$CDM universe forward in time, all evidence of the Hubble expansion will disappear, so that observers in our "island universe" will be fundamentally incapable of determining the true nature of the universe, including the existence of the highly dominant vacuum energy, the existence of the CMB, and the primordial origin of light elements. With these pillars of the modern Big Bang gone, this epoch will mark the end of cosmology and the return of a static universe. In this sense, the coordinate system appropriate for future observers will perhaps fittingly resemble the static coordinate system in which the de Sitter universe was first presented.

We demonstrate that in a vacuum-energy-dominated expansion phase, surprisingly neither the decay of matter nor matter-antimatter annihilation into relativistic particles can ever cause radiation to once again dominate over matter in the future history of the universe.

We examine phantom dark energy models derived from a scalar field with a negative kinetic term for which V(phi) approaches infinity asymptotically. All such models can be divided into three classes, corresponding to an equation of state parameter w with asymptotic behavior w -> -1, w -> w_0 < -1, and w -> infinity. We derive the conditions on the potential V(phi) which lead to each of these three types of behavior. For models with w -> -1, we derive the conditions on V(phi) which determine whether or not such models produce a future big rip. Observational constraints are derived on two classes of these models: power-law potentials with V(phi) = lambda phi^alpha (with alpha positive or negative) and exponential potentials of the form V(phi) = beta e^{lambda phi^alpha}. It is shown that these models spend more time in a state with Omega_m ~ Omega_phi than do corresponding models with a constant value of w, thus providing a more satisfactory solution to the coincidence problem.

We derive conditions for stable tracker solutions for both quintessence and k-essence in a general cosmological background, H^2 \propto f(\rho). We find that tracker solutions are possible only when \eta = d ln f /d ln \rho is constant, aside from a few special cases, which are enumerated. Expressions for the quintessence or k-essence equation of state are derived as a function of \eta and the equation of state of the dominant background component.

We examine the behavior of dark energy models in the plane defined by w (the equation of state parameter for the dark energy) and w' (the derivative of w with respect to the logarithm of the scale factor). For non-phantom barotropic fluids with positive squared sound speed, we find that w' < 3w(w+1), the opposite of the bound on quintessence models previously derived by Caldwell and Linder. Thus, these barotropic models and quintessence models for the dark energy occupy disjoint regions in the w - w' plane. We also derive two new bounds for quintessence models in the w - w' plane: the first is a general bound for any scalar field with a monotonic potential, while the second improves on the Caldwell-Linder bound for tracker quintessence models. Observationally distinguishing barotropic models from quintessence models requires \sigma(w') < 1+w.

We revisit a scenario in which the cosmological constant is cancelled by the potential energy of a slowly evolving scalar field, or "cosmon". The cosmon's evolution is tied to the cosmological constant by a feedback mechanism. This feedback is achieved by an unconventional coupling of the cosmon field to the Ricci curvature scalar. The solutions show that the effective cosmological constant evolves approximately as $t^{-2}$ and remains always of the same order as the density of ordinary matter and radiation. Newton's constant varies on cosmological time scales, with $\dot{G}_N/G_N \ll 1/t$. $G_N$ could have been somewhat different, and possibly smaller, at the time of Big Bang nucleosynthesis.

The generalized Chaplygin gas is characterized by the equation of state p = - A/rho^alpha, with alpha > -1 and w > -1. We generalize this model to allow for the cases where alpha < -1 or w < -1. This generalization leads to three new versions of the generalized Chaplygin gas: an early phantom model in which w << -1 at early times and asymptotically approaches w = -1 at late times, a late phantom model with w \approx -1 at early times and w -> - \infty at late times, and a transient model with w \approx -1 at early times and w -> 0 at late times. We consider these three cases as models for dark energy alone and examine constraints from type Ia supernovae and from the subhorizon growth of density perturbations. The transient Chaplygin gas model provides a possible mechanism to allow for a currently accelerating universe without a future horizon, while some of the early phantom models produce w < -1 without either past or future singularities.

We study here the effect of a varying G on the evolution of the early Universe and, in particular, on primordial nucleosynthesis. This variation of G is modelled using the Brans-Dicke theory as well as a more general class of scalar-tensor theories. Modified nucleosynthesis codes are used to investigate this effect and the results obtained are used to constrain the parameters of the theories. We extend previous studies of primordial nucleosynthesis in scalar-tensor theories by including effects which can cause a slow variation of G during radiation domination and by including a late-time accelerating phase to the Universe's history. We include a brief discussion on the epoch of matter-radiation equality in Brans-Dicke theory, which is also of interest for determining the positions of the cosmic microwave background power-spectrum peaks.

Phantom dark energy models, with w < -1, are characterized by a future singularity and therefore a finite lifetime for the universe. Because the future singularity is triggered by the onset of dark-energy domination, the universe spends a significant fraction of its total lifetime in a state for which the dark energy and matter densities are roughly comparable. We calculate, as a function of w, the fraction of the total lifetime of the universe for which the dark energy and matter densities differ by less than the ratio r_0 in either direction. For r_0 = 10, this fraction varies from 1/3 to 1/8 as w varies from -1.5 to -1.1; the fraction is smaller for smaller values of r_0. This result indicates that the coincidence problem is significantly ameliorated in phantom-dominated cosmologies with a future singularity.

We examine some cosmological consequences of gravity coupling with different strength to fermions and bosons. We show that this leads to a different perturbation of the standard picture of primordial nucleosynthesis than the addition of extra neutrino types or overall scaling of the value of G. Observed abundances of deuterium and helium-4 place bounds on the ratio of the bosonic gravitational constant (G_B) to the fermionic gravitational constant (G_F) of 0.45 < G_B/G_F < 0.92 at 1-sigma, and 0.33 < G_B/G_F < 1.10 at 2-sigma. A value of G_B < G_F can reconcile the current "tension" between the abundances of deuterium and helium-4 predicted by primordial nucleosynthesis. We comment briefly on other cosmological effects.

We examine k-essence models in which the Lagrangian p is a function only of the derivatives of a scalar field phi and does not depend explicity on phi. The evolution of phi for an arbitrary functional form for p can be given in terms of an exact analytic solution. For quite general conditions on the functional form of p, such models can evolve to a state characterized by a density scaling with the scale factor as rho = rho_0 + rho_1(a/a_0)^{-3}, but with a sound speed c_s^2 << 1 at all times. Such models can serve as a unified model for dark matter and dark energy, while avoiding the problems of the generalized Chaplygin gas models, which are due to a non-negligible sound speed in these models. A dark energy component with c_s << 1 serves to suppress cosmic microwave background fluctuations on large angular scales.

We examine cosmic microwave background constraints on variations in the recombination history of the universe. We use a very general extension to the standard model of recombination characterized by two parameters, a and b, which multiply the overall rates of recombination and ionization and the binding energies of hydrogen respectively. We utilize WMAP temperature and TE cross-correlation data, along with temperature data at smaller scales from ACBAR and CBI, to place constraints on these parameters. The range of recombination histories which gives the best fit to the data is -0.6 < log(a) < 0.5 and 0.85 < b < 1.15, forming a diagonal region from low-a/high-b to high-a/low-b. We find z_rec = 1055 +- 25 and \delta z = 198 +- 5 at 68% confidence for the range of recombination models which are best able to reproduce the observed CMB power spectra. Standard recombination provides an acceptable fit to the data; while there is still room for non-standard recombination, the allowed variation is small.

The current uncertainty in Newton's constant, G_N, is of the order of 0.15%. For values of the baryon to photon ratio consistent with both cosmic microwave background observations and the primordial deuterium abundance, this uncertainty in G_N corresponds to an uncertainty in the primordial 4He mass fraction, Y_P, of +-1.3 x 10^{-4}. This uncertainty in Y_P is comparable to the effect from the current uncertainty in the neutron lifetime, which is often treated as the dominant uncertainty in calculations of Y_P. Recent measurements of G_N seem to be converging within a smaller range; a reduction in the estimated error on G_N by a factor of 10 would essentially eliminate it as a source of uncertainty in the calculation of the primordial 4He abundance.

Quintessence models based on a scalar field, phi, with an inverse power law potential display simple tracking behavior at early times, when the quintessence energy density, rho_phi, is sub-dominant. At late times, when rho_phi becomes comparable to the matter density, the evolution of phi diverges from its scaling behavior. We calculate the first order departure of phi from its tracker solution at low redshift. Our results for the evolution of phi, rho_phi, Omega_phi, and w are suprisingly accurate even down to z=0. We find that w and Omega_phi are related linearly to first order. We also derive a semi-analytic expression for w(z) which is accurate to within a few percent. Our analytic techniques are potentially applicable to any quintessence model in which the quintessence component comes to dominate at late times.

We derive constraints on the time variation of the Higgs vacuum expectation value $<\phi>$ through the effects on Big Bang nucleosynthesis (BBN) and the cosmic microwave background (CMB). In the former case, we include the (previously-neglected) effect of the change in the deuteron binding energy, which alters both the $^4$He and deuterium abundances significantly. We find that the current BBN limits on the relative change in \higgs are $-(0.6 - 0.7) \times 10^{-2} < Delta<\phi>/<\phi> < (1.5 - 2.0) \times 10^{-2}$, where the exact limits depend on the model we choose for the dependence of the deuteron binding energy on \higgs.The limits from the current CMB data are much weaker.

We examine constraints from Big Bang nucleosynthesis on type II Randall-Sundrum brane cosmologies with both a dark radiation component and a quadratic term that depends on the 5-dimensional Planck mass, M_5. Using limits on the abundances of deuterium and helium-4, we calculate the allowed region in the M_5-dark radiation plane and derive the precise BBN bound on M_5 alone with no dark radiation: M_5 > 13 TeV.

We consider the effect of inhomogeneous neutrino degeneracy on Big Bang nucleosynthesis for the case where the distribution of neutrino chemical potentials is given by a Gaussian. The chemical potential fluctuations are taken to be isocurvature, so that only inhomogeneities in the electron chemical potential are relevant. Then the final element abundances are a function only of the baryon-photon ratio $\eta$, the effective number of additional neutrinos $\Delta N_\nu$, the mean electron neutrino degeneracy parameter $\bar \xi$, and the rms fluctuation of the degeneracy parameter, $\sigma_\xi$. We find that for fixed $\eta$, $\Delta N_\nu$, and $\bar \xi$, the abundances of helium-4, deuterium, and lithium-7 are, in general, increasing functions of $\sigma_\xi$. Hence, the effect of adding a Gaussian distribution for the electron neutrino degeneracy parameter is to decrease the allowed range for $\eta$. We show that this result can be generalized to a wide variety of distributions for $\xi$.

We study the effect on the cosmic microwave background (CMB) anisotropy and large scale structure (LSS) power spectrum of a scattering interaction between cold dark matter and baryons. This scattering alters the CMB anisotropy and LSS spectrum through momentum transfer between the cold dark matter particles and the baryons. We find that current CMB observations can put an upper limit on the scattering cross section which is comparable with or slightly stronger than previous disk heating constraints at masses greater than 1 GeV, and much stronger at smaller masses. When large-scale structure constraints are added to the CMB limits, our constraint is more stringent than this previous limit at all masses. In particular, a dark matter-baryon scattering cross section comparable to the ``Spergel-Steinhardt'' cross section is ruled out for dark matter mass greater than 1 GeV.

The dark energy that appears to produce the accelerating expansion of the universe can be characterized by an equation of state p=w\rho with w<-1/3. A number of observational tests have been proposed to study the value or redshift dependence of w, including SN Ia distances, the Sunyaev-Zel'dovich effect, cluster abundances, strong and weak gravitational lensing, galaxy and quasar clustering, galaxy ages, the \lya forest, and CMB anisotropies. The proposed observational tests based on these phenomena measure either the distance-redshift relation d(z), the Hubble parameter H(z), the age of the universe t(z), the linear growth factor D_1(z), or some combination of these quantities. We compute the evolution of these four observables, and of the combination H(z)d(z) that enters the Alcock-Paczyznski anisotropy test, in models with constant w, in quintessence models with some simple forms of the potential V(\phi), and in toy models that allow more radical time variations of w. Measurement of any of these quantities to precision of a few percent is generally sufficient to discriminate between w=-1 and w=-2/3. However, the time-dependence predicted in quintessence models is extremely difficult to discern because the quintessence component is dynamically unimportant at the redshifts where w departs substantially from its low-z value. Even for the toy models that allow substantial changes in w at low redshift, there is always a constant-w model that produces very similar evolution of all of the observables simultaneously. We conclude that measurement of the effective equation of state of the dark energy may be achieved by several independent routes in the next few years, but that detecting time-variation in this equation of state will prove very difficult except in specialized cases.

Recent cosmic microwave background (CMB) results from BOOMERANG, MAXIMA, and DASI provide cosmological constraints on new physics that can be competitive with those derived from Big Bang Nucleosynthesis (BBN). In particular, both CMB and BBN can be used to place limits on models involving neutrino degeneracy and additional relativistic degrees of freedom. However, for the case of the CMB, these constraints are, in general, sensitive to the assumed priors. We examine the CMB and BBN constraints on such models and study the sensitivity of ``new physics" to the assumed priors. If we add a constraint on the age of the universe ($t_0 \ga 11$ Gyr), then for models with a cosmological constant, the range of baryon densities and neutrino degeneracy parameters allowed by the CMB and BBN is fairly robust: $\eta_{10} = 6.0 \pm 0.6$, \deln $\la 6$, $\xi_e \la 0.3$. In the absence of new physics, models without a cosmological constant are only marginally compatible with recent CMB observations (excluded at the 93% confidence level).

We use the spherical collapse (SC) approximation to derive expressions for the smoothed redshift-space probability distribution function (PDF), as well as the $p$-order hierarchical amplitudes $S_p$, in both real and redshift space. We compare our results with numerical simulations, focusing on the $\Omega=1$ standard CDM model, where redshift distortions are strongest. We find good agreement between the SC predictions and the numerical PDF in real space even for $\sigma_L \simgt 1$, where $\sigma_L$ is the linearly-evolved rms fluctuation on the smoothing scale. In redshift space, reasonable agreement is possible only for $\sigma_L \simlt 0.4$. Numerical simulations also yield a simple empirical relation between the real-space PDF and redshift-space PDF: we find that for $\sigma \simlt 1$, the redshift space PDF, P[\delta_z], is, to a good approximation, a simple rescaling of the real space PDF, P[\delta], i.e., P[\delta/\sigma] d[\delta/\sigma] = P[\delta_z/\sigma_z] d[\delta_z/\sigma_z], where $\sigma$ and \sigma_z are the real-space and redshift-space rms fluctuations, respectively. This result applies well beyond the validity of linear perturbation theory, and it is a good fit for both the standard CDM model and the Lambda-CDM model. It breaks down for SCDM at $\sigma \approx 1$, but provides a good fit to the \Lambda-CDM models for $\sigma$ as large as 0.8.

In extended quintessence models, a scalar field which couples to the curvature scalar R provides most of the energy density of the universe. We point out that such models can also lead naturally to a decrease in the primordial abundance of helium-4, relieving the tension which currently exists between the primordial helium-4 abundance inferred from observations and the amount predicted by Standard Big Bang Nucleosynthesis (SBBN) corresponding to the observed deuterium abundance. Using negative power-law potentials for the quintessence field, we determine the range of model parameters which can lead to an interesting reduction in the helium-4 abundance, and we show that it overlaps with the region allowed by other constraints on extended quintessence models.

Interest in non-standard recombination scenarios has been spurred by recent cosmic microwave background (CMB) results from BOOMERANG and MAXIMA, which show an unexpectedly low second acoustic peak, resulting in a best-fit baryon density that is 50% larger than the prediction of big-bang nucleosynthesis (BBN). This apparent discrepancy can be avoided if the universe has a non-standard ionization history in which the recombination of hydrogen is significantly delayed relative to the standard model. While future CMB observations may eliminate this discrepancy, it is useful to develop a general framework for analyzing non-standard ionization histories. We develop such a framework, examining non-standard models in which the hydrogen binding energy E_b and the overall expression for the time rate of change of the ionized fraction of electrons are multiplied by arbitrary factors. This set of models includes a number of previously-proposed models as special cases. We find a wide range of models with delayed recombination that are able to fit the CMB data with a baryon density in accordance with BBN, but there are even allowed models with earlier recombination than in the standard model. A generic prediction of these models is that the third acoustic CMB peak should be very low relative to what is found in the standard model. This is the case even for the models with earlier recombination than in the standard model, because here the third peak is lowered by an increased diffusion damping at recombination relative to the standard model. Interestingly, the specific height of the third peak depends sensitively on the model parameters, so that future CMB measurements will be able to distinguish between different non-standard recombination scenarios.

Energy injected into the Cosmic Microwave Background at redshifts z<10^6 will distort its spectrum permanently. In this paper we discuss the distortion caused by annihilations of relic particles. We use the observational bounds on deviations from a Planck spectrum to constrain a combination of annihilation cross section, mass, and abundance. For particles with (s-wave) annihilation cross section, <\sigma|v|>=\sigma_0, the bound is f[(\sigma_0/6e-27cm^3/s)(\Omega_{X\bar{X}}h^2)^2]/(m_X/MeV)<0.2, where m_X is the particle mass, \Omega_{X\bar{X}} is the fraction of the critical density the particle and its antiparticle contribute if they survive to the present time, h=H_0/(100km/s/Mpc), H_0 is the Hubble constant, and f is the fraction of the annihilation energy that interacts electromagnetically. We also compute the less stringent limits for p-wave annihilation. We update other bounds on residual annihilations and compare them to our CMB bound.

It has been shown by many independent studies that the cold dark matter scenario produces singular galactic dark halos, in strong contrast with observations. Possible remedies are that either the dark matter is warm so that it has significant thermal motion or that the dark matter has strong self interactions. We have combined these ideas to calculate the linear mass power spectrum and the spectrum of cosmic microwave background (CMB) fluctuations for self-interacting warm dark matter. Our results indicate that such models have more power on small scales than is the case for the standard warm dark matter model, with a CMB fluctuation spectrum which is nearly indistinguishable from standard cold dark matter. This enhanced small-scale power may provide better agreement with the observations than does standard warm dark matter.

A time variation in the Higgs vacuum expectation value alters the electron mass and thereby changes the ionization history of the universe. This change produces a measurable imprint on the pattern of cosmic microwave background (CMB) fluctuations. The nuclear masses and nuclear binding energies, as well as the Fermi coupling constant, are also altered, with negligible impact on the CMB. We calculate the changes in the spectrum of the CMB fluctuations as a function of the change in the electron mass. We find that future CMB experiments could be sensitive to |\Delta m_e/m_e| \sim |\Delta G_F/G_F| \sim 10^{-2} - 10^{-3}. However, we also show that a change in the electron mass is nearly, but not exactly, degenerate with a change in the fine-structure constant. If both the electron mass and the fine-structure constant are time-varying, the corresponding CMB limits are much weaker, particularly for l < 1000.

We examine Big Bang nucleosynthesis (BBN) in the case of inhomogenous neutrino degeneracy, in the limit where the fluctuations are sufficiently small on large length scales that the present-day element abundances are homogeneous. We consider two representive cases: degeneracy of the electron neutrino alone, and equal chemical potentials for all three neutrinos. We use a linear programming method to constrain an arbitrary distribution of the chemical potentials. For the current set of (highly-restrictive) limits on the primordial element abundances, homogeneous neutrino degeneracy barely changes the allowed range of the baryon-to-photon ratio. Inhomogeneous degeneracy allows for little change in the lower bound on the baryon-to-photon ratio, but the upper bound in this case can be as large as 1.1 \times 10^{-8} (only electron neutrino degeneracy) or 1.0 \times 10^{-9} (equal degeneracies for all three neutrinos). For the case of inhomogeneous neutrino degeneracy, we show that there is no BBN upper bound on the neutrino energy density, which is bounded in this case only by limits from structure formation and the cosmic microwave background.

A massive neutrino which decays after recombination (t > 10^{13} sec) into relativistic decay products produces an enhanced integrated Sachs-Wolfe effect, allowing constraints to be placed on such neutrinos from present cosmic microwave background anisotropy data. Previous treatments of this problem have approximated the decay products as an additional component of the neutrino background. This approach violates energy-momentum conservation, and we show that it leads to serious errors for some neutrino masses and lifetimes. We redo this calculation more accurately, by correctly incorporating the spatial distribution of the decay products. For low neutrino masses and long lifetimes, we obtain a much smaller distortion in the CMB fluctuation spectrum than have previous treatments. We combine these new results with a recent set of CMB data to exclude the mass and lifetime range m_h > 100 eV, \tau > 10^{12} sec. Masses as low as 30 eV are excluded for a narrower range in lifetime.

Any time variation in the fine-structure constant alters the ionization history of the universe and therefore changes the pattern of cosmic microwave background fluctuations. We calculate the changes in the spectrum of these fluctuations as a function of the change in the fine-structure constant, and we find that these changes are dominated by the change in the redshift of recombination due to the shift in the binding energy of hydrogen. We estimate the accuracy with which the next generation of cosmic microwave background experiments might constrain any variation in the fine-structure constant at z of roughly 1000. We find that such experiments could potentially be sensitive to $|\Delta \alpha/\alpha| \sim 10^{-2} - 10^{-3}$.

An attractive method of obtaining an effective cosmological constant at the present epoch is through the potential energy of a scalar field. Considering models with a perfect fluid and a scalar field, we classify all potentials for which the scalar field energy density scales as a power-law of the scale factor when the perfect fluid density dominates. There are three possibilities. The first two are well known; the much-investigated exponential potentials have the scalar field mimicking the evolution of the perfect fluid, while for negative power-laws, introduced by Ratra and Peebles, the scalar field density grows relative to that of the fluid. The third possibility is a new one, where the potential is a positive power-law and the scalar field energy density decays relative to the perfect fluid. We provide a complete analysis of exact solutions and their stability properties, and investigate a range of possible cosmological applications.

We examine simulations of the formation of domain walls, cosmic strings, and monopoles on a cubic lattice, in which the topological defects are assumed to lie at the zeros of a piecewise constant 1, 2, or 3 component Gaussian random field, respectively. We derive analytic expressions for the corresponding topological defect densities in the continuum limit and show that they fail to agree with simulation results, even when the fields are smoothed on small scales to eliminate lattice effects. We demonstrate that this discrepancy, which is related to a classic geometric fallacy, is due to the anisotropy of the cubic lattice, which cannot be eliminated by smoothing. This problem can be resolved by linearly interpolating the field values on the lattice, which gives results in good agreement with the continuum predictions. We use this procedure to obtain a lattice-free estimate (for Gaussian smoothing) of the fraction of the total length of string in the form of infinite strings: $f_\infty = 0.716 \pm 0.015$.

The pattern of anisotropies in the Cosmic Microwave Background depends upon the masses and lifetimes of the three neutrino species. A neutrino species of mass greater than 10 eV with lifetime between 10^{13} sec and 10^{17} sec leaves a very distinct signature (due to the integrated Sachs-Wolfe effect): the anisotropies at large angles are predicted to be comparable to those on degree scales. Present data exclude such a possibility and hence this region of parameter space. For $m_\nu \simeq 30$ eV, $\tau \simeq 10^{13}$ sec, we find an interesting possibility: the Integrated Sachs Wolfe peak produced by the decaying neutrino in low-$\Omega$ models mimics the acoustic peak expected in an $\Omega = 1$ model.

While it is well-known that "biased galaxy formation" can increase the strength of galaxy clustering, it is less clear whether straightforward biasing schemes can change the shape of the galaxy correlation function on large scales. Here we consider "local" biasing models, in which the galaxy density field $\delta_g$ at a point is a function of the matter density field $\delta$ at that point: $\delta_g = f(\delta)$. We consider both deterministic biasing, where $f$ is simply a function, and stochastic biasing, in which the galaxy density is a random variable whose distribution depends on the matter density: $\delta_g=X(\delta)$. We show that even when this mapping is performed on a highly nonlinear density field with a hierarchical correlation structure, the correlation function $\xi$ is simply scaled up by a constant, as long as $\xi << 1$. In stochastic biasing models, the galaxy autocorrelation function behaves exactly as in deterministic models, with the mean value of $X$ for a given value of $\delta$ taking the role of the deterministic bias function. We extend our results to the power spectrum P(k), showing that for sufficiently small k, the effect of local biasing is equivalent to the multiplication of P(k) by a constant, with the addition of a constant term. If a cosmological model predicts a large-scale mass correlation function in conflict with the shape of the observed galaxy correlation function, then the model cannot be rescued by appealing to a complicated but local relation between galaxies and mass.

We simulate the formation of cosmic strings at the zeros of a complex Gaussian field with a power spectrum $P(k) \propto k^n$, specifically addressing the issue of the fraction of length in infinite strings. We make two improvements over previous simulations: we include a non-zero random background field in our box to simulate the effect of long-wavelength modes, and we examine the effects of smoothing the field on small scales. The inclusion of the background field significantly reduces the fraction of length in infinite strings for $n < -2$. Our results are consistent with the possibility that infinite strings disappear at some $n = n_c$ in the range $-3 \le n_c < -2.2$, although we cannot rule out $n_c = -3$, in which case infinite strings would disappear only at the point where the mean string density goes to zero. We present an analytic argument which suggests the latter case. Smoothing on small scales eliminates closed loops on the order of the lattice cell size and leads to a ``lattice-free" estimate of the infinite string fraction. As expected, this fraction depends on the type of window function used for smoothing.

We test tree-level perturbation theory for Gaussian initial conditions with power spectra $P(k)\propto k^n$ by comparing the probability distribution function (PDF) for the density predicted by the Local Lagrangian Approximation (LLA) with the results of numerical gravitational clustering simulations. Our results indicate that our approximation correctly reproduces the evolved density PDF for $-3 \leq n \leq-1$ power spectra up to the weakly nonlinear regime, while it shows marginal agreement for power indices n=0 and +1 in the linear regime and poor agreement beyond this point. This suggests that tree-level perturbation theory (as realized in the Local Lagrangian Approximation) can accurately predict the density distribution function for $-3 \leq n \leq -1$ but fails for $n \ge 0$.

We compute the skewness $t_3$ and the corresponding hierarchical amplitude $T_3$ of the divergence of the velocity field for arbitrary non-Gaussian initial conditions. We find that $T_3$ qualitatively resembles the corresponding hierarchical amplitude for the density field, $S_3$, in that it contains a term proportional to the initial skewness, which decays inversely as the linear growth factor, plus a constant term which differs from the corresponding Gaussian term by a complex function of the initial three- and four- point functions. We extend the results for $S_3$ and $T_3$ with non-Gaussian initial conditions to evolved fields smoothed with a spherical tophat window function. We show that certain linear combinations, namely $S_3 + {1 \over 2} T_3$, $S_3 + T_3$, and $s_3 + t_3$, lead to expressions which are much simpler, for non-Gaussian initial conditions, than $S_3$ and $T_3$ (or $s_3$ and $t_3$) considered separately.

We examine local Lagrangian approximations for the gravitational evolution of the density distribution function. In these approximations, the final density at a Lagrangian point q at a time t is taken to be a function only of t and of the initial density at the same Lagrangian point. A general expression is given for the evolved density distribution function for such approximations, and we show that the vertex generating function for a local Lagrangian mapping applied to an initially Gaussian density field bears a simple relation to the mapping itself. Using this result, we design a local Lagrangian mapping which reproduces nearly exactly the hierarchical amplitudes given by perturbation theory for gravitational evolution. When extended to smoothed density fields and applied to Gaussian initial conditions, this mapping produces a final density distribution function in excellent agreement with full numerical simulations of gravitational clustering. We also examine the application of these local Lagrangian approximations to non-Gaussian initial conditions.

We present new upper and lower bounds to the primordial abundances of deuterium and helium-3 based on observational data from the solar system and the interstellar medium. Independent of any model for the primordial production of the elements we find (at the 95\% C.L.): $1.5 \times 10^{-5} \le (D/H)_P \le 10.0 \times 10^{-5}$ and $(^3He/H)_P \le 2.6\times 10^{-5}$. When combined with the predictions of standard big bang nucleosynthesis, these constraints lead to a 95\% C.L. bound on the primordial abundance of deuterium: $(D/H)_{best} = (3.5^{+2.7}_{-1.8})\times 10^{-5}$. Measurements of deuterium absorption in the spectra of high redshift QSOs will directly test this prediction. The implications of this prediction for the primordial abundances of helium-4 and lithium-7 are discussed, as well as those for the universal density of baryons.

A new evaluation of the constraint on the number of light neutrino species (N_nu) from big bang nucleosynthesis suggests a discrepancy between the predicted light element abundances and those inferred from observations, unless the inferred primordial 4He abundance has been underestimated by 0.014 +/- 0.004 (1 sigma) or less than 10% (95%C.L.) of 3He survives stellar processing. With the quoted systematic errors in the observed abundances and a conservative chemical evolution parameterization, the best fit to the combined data is N_nu = 2.1 +/- 0.3 (1 sigma) and the upper limit is N_nu < 2.6 (95% C.L.). The data are inconsistent with the Standard Model (N_nu = 3) at the 98.6% C.L.

We examine inhomogeneous primordial nucleosynthesis for {\it arbitrary} distributions $f$ of the baryon-to-photon ratio $\eta$, in the limit where neither particle diffusion nor gravitational collapse is important. By discretizing $f(\eta)$ and using linear programming, we show that for a set of $m$ constraints on the primordial element abundances, the maximum and minimum possible values of $\bar \eta$ (the final mean value of $\eta$) are given when $f(\eta)$ is a sum of at most $m+1$ delta functions. Our linear programming results indicate that when $f$ is taken to be an arbitrary function, there is no lower bound on $\bar \eta$, while the upper bound is essentially the homogeneous upper bound.

A decaying particle that dominates the energy density of the universe and then decays with a lifetime $\sim 10^{-1}$ sec into electromagnetically-interacting particles can heat the background electron neutrinos more than the $\mu$ and $\tau$ neutrinos, resulting in a decrease in the primordial $^4$He abundance. Previous calculations of this effect assumed sudden decoupling of the neutrinos from the electromagnetic plasma. Here we repeat this calculation using an improved (although still approximate) treatment of the energy flow from the blackbody radiation into the various neutrino components. We find that with our improved calculation, the reduction in the primordial $^4$He is not as large as had been previously thought, but it occurs over a larger range in decaying particle lifetime $\tau$. The reduction in $^4$He can be as large as $\Delta Y = -0.01$ for $\tau = 0.1 - 0.7$ sec.

We examine Gaussian isocurvature hot dark matter (massive neutrino) models for large-scale structure in which the initial density perturbations are produced in the baryons with a power--law spectrum $P_B(k) = Ak^{n_B}$. We calculate the linearly-evolved power spectrum and cosmic microwave fluctuations. We find that models with only isocurvature perturbations are inconsistent with observations of damped Ly$\alpha$ systems and COBE constraints on the power index. However, models which contain a mixture of adiabatic and isocurvature perturbations can be made consistent with COBE, galaxy surveys and damped Ly$\alpha$ systems.Isocurvature hot dark matter models also produce a bias between baryons and neutrinos even in the linear regime. We find that this ``natural bias'' can increase the baryon fraction in small scale objects like damped Ly$\alpha$ systems, but it has no effect on cluster scales.

We consider the growth of density perturbations in the presence of self--interacting dark matter, SIDM, proposed by Carlson, Machacek and Hall (1992). We determine the range of values for the coupling constant $\lambda$ and the particle mass $m^\prime$, for which the power spectrum lies in the ``allowed" range based on constraints from the IRAS galaxy survey and damped Lyman--$\alpha $ systems. Our results show that no combination of parameters can meet both limits. We consider constraints on the $\2-2$ scatterings which keep the SIDM particles in pressure equilibrium, and we show that if such interactions maintain pressure equilibrium down to the present, they will be strong enough to disrupt galaxy mergers and may lead to stripping of galaxy halos as galaxies move through the dark matter background of these particles. Hence, we also investigate the evolution of large-scale structure in the SIDM model when the particles drop out of pressure equilibrium at some higher redshift. The resulting free-streaming leads to an additional suppression of small-scale perturbations, but it does not significantly affect our results.

The Sachs-Wolfe temperature fluctuations produced by primordial density perturbations are proportional to the potential field \phi, which is a weighted integral over the density field \delta. Because of the central limit theorem, \phi can be approximately Gaussian even when \delta is non-Gaussian. Using the Wold representation for non-Gaussian density fields, \delta(\rvec) = \int f(|\rvec - \rvec^\prime|) \Delta(\rvec^\prime) d^3 \rvec^\prime, we find conditions on \Delta and f for which \phi must have a Gaussian one-point distribution, while \delta can be non-Gaussian. Sufficient (but not necessary) conditions are that the density field have a power spectrum (which determines f) of P(k) \propto k^n, with -2 < n \le +1, and that \Delta(\rvec) be non-Gaussian with no long-range correlations. Thus, there is an infinite set of non-Gaussian density fields which produce a nearly Gaussian one-point distribution for the Sachs-Wolfe effect.

We examine the multi--point velocity field for non--Gaussian models as a probe of non--Gaussian behavior. The two--point velocity correlation is not a useful indicator of a non--Gaussian density field, since it depends only on the power spectrum, even for non--Gaussian models. However, we show that the distribution of velocity differences $\bfv_1 - \bfv_2$, where $\bfv_1$ and $\bfv_2$ are measured at the points $\bfr_1$ and $\bfr_2$, respectively, is a good probe of non--Gaussian behavior, in that $p(\bfv_1 - \bfv_2)$ tends to be non--Gaussian whenever the density field is non--Gaussian. As an example, we examine the behavior of $p(\bfv_1 - \bfv_2)$ for non--Gaussian seed models, in which the density field is the convolution of a distribution of points with a set of density profiles. We apply these results to the global texture model.