We construct a perturbative S-matrix for interacting massive scalar fields in global de Sitter space. Our S-matrix is formulated in terms of asymptotic particle states in the far past and future, taking appropriate care for light fields whose wavefunctions decay only very slowly near the de Sitter conformal boundaries. An alternative formulation expresses this S-matrix in terms of residues of poles in analytically-continued Euclidean correlators (computed in perturbation theory), making it clear that the standard Minkowski-space result is obtained in the flat-space limit. Our S-matrix transforms properly under CPT, is invariant under the de Sitter isometries and perturbative field redefinitions, and is unitary. This unitarity implies a de Sitter version of the optical theorem. We explicitly verify these properties to second order in the coupling for a general cubic interaction, including both tree- and loop-level contributions. Contrary to other statements in the literature, we find that a particle of any positive mass may decay at tree level to any number of particles, each of arbitrary positive masses. In particular, even very light fields (in the complementary series of de Sitter representations) are not protected from tree-level decays.

We extend previous combinations of LEP and cosmological relic density constraints on the parameter space of the constrained MSSM, with universal input supersymmetry-breaking parameters, to large tan beta. We take account of the possibility that the lightest Higgs boson might weigh about 115 GeV, but also retain the possibility that it might be heavier. We include the most recent implementation of the b to s gamma constraint at large tan beta. We refine previous relic density calculations at large tan beta by combining a careful treatment of the direct-channel Higgs poles in annihilation of pairs of neutralinos chi with a complete treatment of chi - stau coannihilation, and discuss carefully uncertainties associated with the mass of the b quark. We find that coannihilation and pole annihilations allow the CMSSM to yield an acceptable relic density at large tan beta, but it is consistent with all the constraints only if m_chi > 140 (180) GeV for mu > 0 (mu < 0) for our default choices m_b(m_b) = 4.25 GeV, m_t = 175 GeV, and A_0 = 0.

Assuming that the lightest supersymmetric particle (LSP) is the lightest neutralino, we present a detailed exploration of neutralino-stau coannihilation channels, including analytical expressions and numerical results. We also include neutralino coannihilations with the selectron and smuon. We evaluate the implications of coannihilations for the cosmological relic density of the LSP, which is assumed to be stable, in the constrained minimal supersymmetric extension of the Standard Model (CMSSM), in which the soft supersymmetry-breaking parameters are universal at the supergravity GUT scale. We evaluate the changes due to coannihilations in the region of the MSSM parameter space that is consistent with the cosmological upper limit on the relic LSP density. In particular, we find that the upper limit on the mass of the neutralino is increased from about 200 GeV to about 600 GeV in the CMSSM, and estimate a qualitatively similar increase for gauginos in the general MSSM.

Matsumoto and Yoshimura [hep-ph/9910393] have argued that there are loop corrections to the number density of heavy particles (in thermal equilibrium with a gas of light particles) that are not Boltzmann suppressed by a factor of e^(-M/T) at temperatures T well below the mass M of the heavy particle. We argue, however, that their definition of the number density does not correspond to a quantity that could be measured in a realistic experiment. We consider a model where the heavy particles carry a conserved U(1) charge, and the light particles do not. The fluctuations of the net charge in a given volume then provide a measure of the total number of heavy particles in that same volume. We show that these charge fluctuations are Boltzmann suppressed (to all orders in perturbation theory). Therefore, we argue, the number density of heavy particles is also Boltzmann suppressed.

Starting from a Caldeira-Leggett model for the interaction of a system with an environment, Joichi, Matsumoto, and Yoshimura have reconsidered the derivation of the quantum Boltzmann equation. They find an extra term that accounts for the effects of virtual particles, and which drastically changes the results for relic densities of stable, weakly interacting massive particles (WIMPs), and for the decay products of unstable particles. We show, however, that this modified Boltzmann equation does not properly account for the interaction energy between the massive particles (which are decaying or annihilating) and the thermal bath of light particles. We argue that the conventional Boltzmann equation gives the correct result.

We study the time evolution of the quantum fluctuations of the axion field for both the QCD axion as well as axions arising in the context of supergravity and string theories. We explicitly keep track not only of the coherently oscillating zero momentum mode of the axion but also of the higher non-zero momentum modes using the full axion potential. The full axion potential makes possible two kinds of instabilities: spinodal instabilities and parametric resonance instabilities. The presence of either of these instabilities can lead to a quasi-exponential increase in the occupation of non-zero momentum modes and the build-up of the quantum fluctuations of the axions. If either of these becomes a significant effect then axions would no longer be a suitable cold dark matter candidate. Our results confirm the conventional wisdom that these effects are not significant in the setting of an expanding FRW universe and hence axions are indeed cold dark matter candidates.

We present a general statistical analysis of Gamma Ray Bursts embedded in a host population. If no host generates more than one observed burst, then we show that there is a model independent lower bound on the number of hosts, $H$, of the form $H > c B^2$, where B is the number of observed bursts, and $c$ is a constant of order one which depends on the confidence level (CL) attached to the bound. An analysis by Tegmark et al. (1996) shows that the BATSE 3B catalog of 1122 bursts is consistent with no repeaters being present, and assuming that this is indeed the case, our result implies a host population with at least H=1.2x10^6 members. Without the explicit assumption of no repeaters, a Bayesian analysis based on the results of Tegmark et al. (1996) can be performed which gives the weaker bound of $H>1.7\times 10^5$ at the 90% CL. In the light of the non-detection of identifiable hosts in the small error-boxes associated with transient counterparts to GRBs, this result gives a model independent lower bound to the number of any rare or exotic hosts. If in fact GRBs are found to be associated with a particular sub-class of galaxies, then an analysis along the lines presented here can be used to place a lower bound on the fraction of galaxies in this sub-class. Another possibility is to treat galaxy clusters (rather than individual galaxies) as the host population, provided that the angular size of each cluster considered is less than the resolution of the detector. Finally, if repeaters are ever detected in a statistically significant manner, this analysis can be readily adapted to find upper and lower limits on $H$.

We study the non-linear dynamics of quantum fields in matter and radiation dominated universes, using the non-equilibrium field theory approach combined with the non-perturbative Hartree and the large N approximations. We examine the phenomenon of explosive particle production due to spinodal instabilities and parametric amplification in expanding universes with and without symmetry breaking. For a variety of initial conditions, we compute the evolution of the inflaton, its quantum fluctuations, and the equation of state. We find explosive growth of quantum fluctuations, although particle production is somewhat sensitive to the expansion of the universe. In the large N limit for symmetry breaking scenarios, we determine generic late time solutions for any flat Friedman-Robertson-Walker cosmology. We also present a complete and numerically implementable renormalization scheme for the equation of motion and the energy momentum tensor in flat FRW cosmologies. In this scheme the renormalization constants are independent of time and of the initial conditions.

The nonlinear time evolution of the quantum fields is studied in the O(N) model for large N in a radiation dominated FRW universe, with a view towards the phenomenon of explosive particle production due to either spinodal instabilities or parametric amplification, i.e. preheating. Quantum backreaction effects due to the produced particles are included consistently within the large N approximation. We find that preheating persists when the expansion is included, although the amount of particle production is reduced compared to the values found in Minkowski space. We also see that the behavior of the evolving zero mode is very different from that in Minkowski space, though the late time behavior in all cases is determined by a sum rule that implies the existence of Goldstone bosons in the final state.

We show that a bounded, isolated quantum system of many particles in a specific initial state will approach thermal equilibrium if the energy eigenfunctions which are superposed to form that state obey {\it Berry's conjecture}. Berry's conjecture is expected to hold only if the corresponding classical system is chaotic, and essentially states that the energy eigenfunctions behave as if they were gaussian random variables. We review the existing evidence, and show that previously neglected effects substantially strengthen the case for Berry's conjecture. We study a rarefied hard-sphere gas as an explicit example of a many-body system which is known to be classically chaotic, and show that an energy eigenstate which obeys Berry's conjecture predicts a Maxwell--Boltzmann, Bose--Einstein, or Fermi--Dirac distribution for the momentum of each constituent particle, depending on whether the wave functions are taken to be nonsymmetric, completely symmetric, or completely antisymmetric functions of the positions of the particles. We call this phenomenon {\it eigenstate thermalization}. We show that a generic initial state will approach thermal equilibrium at least as fast as $O(\hbar/\Delta)t^{-1}$, where $\Delta$ is the uncertainty in the total energy of the gas. This result holds for an individual initial state; in contrast to the classical theory, no averaging over an ensemble of initial states is needed. We argue that these results constitute a new foundation for quantum statistical mechanics.

We discuss the applicability and derivation of window functions for cosmic microwave background experiments on large and intermediate angular scales. These window functions describe the response of the experiment to power in a particular mode of the fluctuation spectrum. We give general formulae, illustrated with specific examples, for the most common observing strategies.

We analyze the Gamma Ursae Minoris (GUM) and Mu Pegasii (MuP) scans of the Millimeter-wave Anisotropy eXperiment in the context of cold dark matter (CDM) models of structure formation, paying particular attention to the two-dimensional nature of the GUM scan. If all of the structure in the (foreground subtracted) data is attributed to cosmic microwave background anisotropy, then there is a detection in each scan. For a standard CDM model, the amplitudes of the signals are individually compatible with the COsmic Background Explorer measurement, but marginally inconsistent with each other.

We present results on the consistency of standard cold dark matter (CDM) models with both the COBE normalization and the data from the 4th channel of the 9-point South Pole scan. We find that CDM models are consistent with both experiments, a conclusion which is at odds with some other analyses. This is partly due to our careful treatment of the temperature autocorrelation function, but also derives from a realization that the statistical conclusions depend strongly on assumptions about the prior distribution for the normalization $Q$.

The fluctuations in the cosmic microwave background radiation may contain deviations from gaussian statistics which would be reflected in a nonzero value of three-point correlation function of $\Delta T$. However, any potential observation of the three-point function is limited by cosmic variance, even if a whole-sky map of $\Delta T$ is available. Here I derive a general formula for the cosmic variance of the three-point function in terms of integrals over the two-point function. This formula can be applied to any cosmological model and to any experimental measurement strategy. It gives a fundamental lower limit on the magnitude of a measurable three-point correlation function, and hence on the measurable amount of skewness in the distribution of $\Delta T$.

We discuss the effects of finite sky coverage and the uncertainty in extracting information about the power spectrum from experiments on small angular scales. In general the cosmic variance is enhanced by a factor of $4\pi/A$, where $A$ is the solid angle sampled by the experiment. As a rough guide, an experiment with sensitivity peaking at the $\ell$th multipole has to cover $\simgt\ell$ independent patches to have a lower ``sample variance'' than for a whole-sky measurement of the quadrupole. Our approach gives a relatively simple way of attaching an error bar to the theoretical prediction for a particular experiment, and thereby comparing theories with experimental results, without the need for computationally-intensive Maximum Likelihood or Monte Carlo calculations.

Inflationary models predict a definite, model independent, angular dependence for the three-point correlation function of $\Delta T/T$ at large angles (greater than $\sim 1^\circ$) which we calculate. The overall amplitude is model dependent and generically unobservably small, but may be large in some specific models. We compare our results with other models of nongaussian fluctuations.

In the standard inflationary scenario with inflaton potential $V(\Phi)=M^4-{1\over4}\lambda\Phi^4$, the resulting density perturbations $\delta\rho/\rho$ are proportional to $\lambda^{1/2}$. Upper bounds on $\delta\rho/\rho$ require $\lambda < 10^{-13}$. Ratra has shown that an alternative treatment of reheating results in $\delta\rho/\rho \propto \lambda^{-1}$, so that an upper bound on $\delta\rho/\rho$ does not put an obvious upper bound on $\lambda$. We verify that $\delta\rho/\rho \propto \lambda^{-1}$ is indeed a possibility, but show that $\lambda < 10^{-13}$ is still required.

A very weakly coupled scalar field with mass $m$ and initial vacuum expectation value $V$ will provide enough mass to close the universe provided $V\simeq (3\times 10^8\gev)(100\gev/m)^{1/4}$. We discuss possible models in which such a field could arise.