We present numerical relativity simulations of cosmological scenarios in which the universe is smoothed and flattened by undergoing a phase of slow contraction and test their sensitivity to a wide range of initial conditions. Our numerical scheme enables the variation of all freely specifiable physical quantities that characterize the initial spatial hypersurface, such as the initial shear and spatial curvature contributions as well as the initial field and velocity distributions of the scalar that drives the cosmological evolution. In particular, we include initial conditions that are far outside the perturbative regime of the well-known attractor scaling solution. We complement our numerical results by analytically performing a complete dynamical systems analysis and show that the two approaches yield consistent results.

Performing a fully non-perturbative analysis using the tools of numerical general relativity, we demonstrate that a period of slow contraction is a `supersmoothing' cosmological phase that homogenizes, isotropizes and flattens the universe both classically and quantum mechanically and can do so far more robustly and rapidly than had been realized in earlier studies.

Perhaps the greatest challenge for fundamental theories based on compactification from extra dimensions is accommodating a period of accelerated cosmological expansion. Previous studies have identified constraints imposed by the existence of dark energy for two overlapping classes of compactified theories: (1) those in which the higher dimensional picture satisfies certain metric properties selected to reproduce known low energy phenomenology; or (2) those derived from string theory assuming they satisfy the Swampland conjectures. For either class, the analyses showed that dark energy is only possible if it takes the form of quintessence. In this paper, we explore the consequences for theories that belong to both classes and show that the joint constraints are highly restrictive, leaving few options.

This white paper addresses key challenges for the design of next-decade Cosmic Microwave Background (CMB) experiments, and for assessing their capability to extract cosmological information from CMB polarization. We focus here on the challenges posed by foreground emission, CMB lensing, and instrumental systematics to detect the signal that arises from gravitational waves sourced by inflation and parameterized by $r$, at the level of $r \sim 10^{-3}$ or lower, as proposed for future observational efforts. We argue that more accurate and robust analysis and simulation tools are required for these experiments to realize their promise. We are optimistic that the capability to simulate the joint impact of foregrounds, CMB lensing, and systematics can be developed to the level necessary to support the design of a space mission at $r \sim 10^{-4}$ in a few years. We make the case here for supporting such work. Although ground-based efforts present additional challenges (e.g., atmosphere, ground pickup), which are not addressed here, they would also benefit from these improved simulation capabilities.

Combining intervals of ekpyrotic (ultra-slow) contraction with a (non-singular) classical bounce naturally leads to a novel cyclic theory of the universe in which the Hubble parameter, energy density and temperature oscillate periodically, but the scale factor grows by an exponential factor from one cycle to the next. The resulting cosmology not only resolves the homogeneity, isotropy, flatness and monopole problems and generates a nearly scale invariant spectrum of density perturbations, but it also addresses a number of age-old cosmological issues that big bang inflationary cosmology does not. There may also be wider-ranging implications for fundamental physics, black holes and quantum measurement.

In this paper, we describe the first steps towards fully non-perturbative cosmology. We explain why the conventional methods used by cosmologists based on the ADM formulation are generally inadequate for this purpose and why it is advantageous instead to adapt the harmonic formulation pioneered and utilized in mathematical and numerical relativity. Here we focus on using this approach to evaluating the linear mode stability in homogeneous and nearly homogeneous backgrounds and devising a valid scheme and diagnostics for numerical computation. We also briefly touch on the relevance of these methods for extracting cosmological observables from non-perturbative simulations.

We study constraints imposed by two proposed string Swampland criteria on cosmology. These criteria involve an upper bound on the range traversed by scalar fields as well as a lower bound on $|\nabla_{\phi} V|/V$ when $V >0$. We find that inflationary models are generically in tension with these two criteria. Applying these same criteria to dark energy in the present epoch, we find that specific quintessence models can satisfy these bounds and, at the same time, satisfy current observational constraints. Assuming the two Swampland criteria are valid, we argue that the universe will undergo a phase transition within a few Hubble times. These criteria sharpen the motivation for future measurements of the tensor-to-scalar ratio $r$ and the dark energy equation of state $w$, and for tests of the equivalence principle for dark matter.

The unique CV-type meteorite Khatyrka is the only natural sample in which "quasicrystals" and associated crystalline Cu,Al-alloys, including khatyrkite and cupalite, have been found. They are suspected to have formed in the early Solar System. To better understand the origin of these exotic phases, and the relationship of Khatyrka to other CV chondrites, we have measured He and Ne in six individual, ~40-{\mu}m-sized olivine grains from Khatyrka. We find a cosmic-ray exposure age of about 2-4 Ma (if the meteoroid was <3 m in diameter, more if it was larger). The U,Th-He ages of the olivine grains suggest that Khatyrka experienced a relatively recent (<600 Ma) shock event, which created pressure and temperature conditions sufficient to form both the quasicrystals and the high-pressure phases found in the meteorite. We propose that the parent body of Khatyrka is the large K-type asteroid 89 Julia, based on its peculiar, but matching reflectance spectrum, evidence for an impact/shock event within the last few 100 Ma (which formed the Julia family), and its location close to strong orbital resonances, so that the Khatyrka meteoroid could plausibly have reached Earth within its rather short cosmic-ray exposure age.

We introduce the "wedge diagram," an intuitive way to illustrate how cosmological models with a classical (non-singular) bounce generically resolve fundamental problems in cosmology. These include the well-known horizon, flatness, and inhomogeneity problems; the small tensor-to-scalar ratio observed in the cosmic microwave background; the low entropy at the beginning of a hot, expanding phase; and the avoidance of quantum runaway. The same diagrammatic approach can be used to compare with other cosmological scenarios.

We recently showed how it is possible to use a cubic Galileon action to construct classical cosmological solutions that enter a contracting null energy condition (NEC) violating phase, bounce at finite values of the scale factor and exit into an expanding NEC-satisfying phase without encountering any singularities or pathologies. A drawback of these examples is that singular behavior is encountered at some time either just before or just after the NEC-violating phase. In this Letter, we show that it is possible to circumvent this problem by extending our method to actions that include the next order ${\cal L}_4$ Galileon interaction. Using this approach, we construct non-singular classical bouncing cosmological solutions that are non-pathological for all times.

One of the fundamental questions of theoretical cosmology is whether the universe can undergo a non-singular bounce, i.e., smoothly transit from a period of contraction to a period of expansion through violation of the null energy condition (NEC) at energies well below the Planck scale and at finite values of the scale factor such that the entire evolution remains classical. A common claim has been that a non-singular bounce either leads to ghost or gradient instabilities or a cosmological singularity. In this letter, we examine cubic Galileon theories and present a procedure for explicitly constructing examples of a non-singular cosmological bounce without encountering any pathologies and maintaining a sub-luminal sound speed for co-moving curvature modes throughout the NEC violating phase. We also discuss the relation between our procedure and earlier work.

In the context of Einstein gravity, if the null energy condition (NEC) is satisfied, the energy density in expanding space-times always decreases while in contracting space-times the energy density grows and the universe eventually collapses into a singularity. In particular, no non-singular bounce is possible. It is, though, an open question if this energy condition can be violated in a controlled way, i.e., without introducing pathologies, such as unstable negative-energy states or an imaginary speed of sound. In this paper, we will re-examine the claim that the recently proposed mimetic scenario can violate the NEC without pathologies. We show that mimetic cosmology is prone to gradient instabilities even in cases when the NEC is satisfied (except for trivial examples). Most interestingly, the source of the instability is always the Einstein-Hilbert term in the action. The matter stress-energy component does not contribute spatial gradient terms but instead makes the problematic curvature modes dynamical. We also show that mimetic cosmology can be understood as a singular limit of known, well-behaved theories involving higher-derivative kinetic terms and discuss ways of removing the instability.

The results from Planck2015, when combined with earlier observations from WMAP, ACT, SPT and other experiments, were the first observations to disfavor the "classic" inflationary paradigm. To satisfy the observational constraints, inflationary theorists have been forced to consider plateau-like inflaton potentials that introduce more parameters and more fine-tuning, problematic initial conditions, multiverse-unpredictability issues, and a new 'unlikeliness problem.' Some propose turning instead to a "postmodern" inflationary paradigm in which the cosmological properties in our observable universe are only locally valid and set randomly, with completely different properties (and perhaps even different physical laws) existing in most regions outside our horizon. By contrast, the new results are consistent with the simplest versions of ekpyrotic cyclic models in which the universe is smoothed and flattened during a period of slow contraction followed by a bounce, and another promising bouncing theory, anamorphic cosmology, has been proposed that can produce distinctive predictions.

We introduce "anamorphic" cosmology, an approach for explaining the smoothness and flatness of the universe on large scales and the generation of a nearly scale-invariant spectrum of adiabatic density perturbations. The defining feature is a smoothing phase that acts like a contracting universe based on some Weyl frame-invariant criteria and an expanding universe based on other frame-invariant criteria. An advantage of the contracting aspects is that it is possible to avoid the multiverse and measure problems that arise in inflationary models. Unlike ekpyrotic models, anamorphic models can be constructed using only a single field and can generate a nearly scale-invariant spectrum of tensor perturbations. Anamorphic models also differ from pre-big bang and matter bounce models that do not explain the smoothness. We present some examples of cosmological models that incorporate an anamorphic smoothing phase.

Ekpyrotic bouncing cosmologies have been proposed as alternatives to inflation. In these scenarios, the universe is smoothed and flattened during a period of slow contraction preceding the bounce while quantum fluctuations generate nearly scale-invariant super-horizon perturbations that seed structure in the post-bounce universe. An analysis by Tolley and Wesley (2007) showed that, for a wide range of ekpyrotic models, generating a scale-invariant spectrum of adiabatic or entropic fluctuations is only possible if the cosmological background is unstable, in which case the scenario is highly sensitive to initial conditions. In this paper, we analyze an important counterexample: a simple action that generates a Gaussian, scale-invariant spectrum of entropic perturbations during ekpyrotic contraction without requiring fine-tuned initial conditions. Based on this example, we discuss some generalizations.

We consider the cosmological consequences if a small fraction ($f\lesssim 0.1$) of the dark matter is ultra-strongly self-interacting, with an elastic self-interaction cross-section per unit mass $\sigma\gg1\ \mathrm{cm^{2}/g}$. This possibility evades all current constraints that assume that the self-interacting component makes up the majority of the dark matter. Nevertheless, even a small fraction of ultra-strongly self-interacting dark matter (uSIDM) can have observable consequences on astrophysical scales. In particular, the uSIDM subcomponent can undergo gravothermal collapse and form seed black holes in the center of a halo. These seed black holes, which form within several hundred halo interaction times, contain a few percent of the total uSIDM mass in the halo. For reasonable values of $\sigma f$, these black holes can form at high enough redshifts to grow to $\sim10^9 M_\odot$ quasars by $z \gtrsim 6$, alleviating tension within the standard $\Lambda$CDM cosmology. The ubiquitous formation of central black holes in halos could also create cores in dwarf galaxies by ejecting matter during binary black hole mergers, potentially resolving the "too big to fail" problem.

We explore a new type of entropic mechanism for generating density perturbations in a contracting phase in which there are two scalar fields, but only one has a steep negative potential. This first field dominates the energy density and is the source of the ekpyrotic equation of state. The second field has a negligible potential, but its kinetic energy density is coupled to the first field with a non-linear sigma-model type interaction. We show that for any ekpyrotic equation of state it is possible to choose the potential and the kinetic coupling such that exactly scale-invariant (or nearly scale-invariant) entropy perturbations are produced. The corresponding background solutions are stable, and the bispectrum of the entropy perturbations vanishes as no non-Gaussianity is produced during the ekpyrotic phase. Hence, the only contribution to non-Gaussianity comes from the non-linearity of the conversion process during which entropic perturbations are turned into adiabatic ones, resulting in a local non-Gaussianity parameter $f_{NL} \sim 5$.

Classic inflation, the theory described in textbooks, is based on the idea that, beginning from typical initial conditions and assuming a simple inflaton potential with a minimum of fine-tuning, inflation can create exponentially large volumes of space that are generically homogeneous, isotropic and flat, with nearly scale-invariant spectra of density and gravitational wave fluctuations that are adiabatic, Gaussian and have generic predictable properties. In a recent paper, we showed that, in addition to having certain conceptual problems known for decades, classic inflation is for the first time also disfavored by data, specifically the most recent data from WMAP, ACT and Planck2013. Guth, Kaiser and Nomura and Linde have each recently published critiques of our paper, but, as made clear here, we all agree about one thing: the problematic state of classic inflation. Instead, they describe an alternative inflationary paradigm that revises the assumptions and goals of inflation, and perhaps of science generally.

In a recent series of papers, we have shown that theories with scalar fields coupled to gravity (e.g., the standard model) can be lifted to a Weyl-invariant equivalent theory in which it is possible to unambiguously trace the classical cosmological evolution through the transition from big crunch to big bang. The key was identifying a sufficient number of finite, Weyl-invariant conserved quantities to uniquely match the fundamental cosmological degrees of freedom across the transition. In so doing we had to account for the well-known fact that many Weyl-invariant quantities diverge at the crunch and bang. Recently, some authors rediscovered a few of these divergences and concluded based on their existence alone that the theories cannot be geodesically complete. In this note, we show that this conclusion is invalid. Using conserved quantities we explicitly construct the complete set of geodesics and show that they pass continuously through the big crunch-big bang transition.

We show how to lift a generic non-scale invariant action in Einstein frame into a locally conformally-invariant (or Weyl-invariant) theory and present a new general form for Lagrangians consistent with Weyl symmetry. Advantages of such a conformally invariant formulation of particle physics and gravity include the possibility of constructing geodesically complete cosmologies. We present a conformal-invariant version of the standard model coupled to gravity, and show how Weyl symmetry may be used to obtain unprecedented analytic control over its cosmological solutions. Within this new framework, generic FRW cosmologies are geodesically complete through a series of big crunch - big bang transitions. We discuss a new scenario of cosmic evolution driven by the Higgs field in a \textquotedblleft minimal\textquotedblright% \ conformal standard model, in which there is no new physics beyond the standard model at low energies, and the current Higgs vacuum is metastable as indicated by the latest LHC data.

The large-scale structure of the universe suggests that the physics underlying its early evolution is scale-free. This was the historic motivation for the Harrison-Zel'dovich-Peebles spectrum and for inflation. Based on a hydrodynamical approach, we identify scale-free forms for the background equation-of-state for both inflationary and cyclic scenarios and use these forms to derive predictions for the spectral tilt and tensor-to-scalar ratio of primordial density perturbations. For the case of inflation, we find three classes of scale-free models with distinct predictions. Including all classes, we show that scale-free inflation predicts tensor-to-scalar ratio $r > 10^{-4}$. We show that the observationally favored class is theoretically disfavored because it suffers from an initial conditions problem and the hydrodynamical form of an unlikeliness problem similar to that identified recently for certain inflaton potentials. We contrast these results with those for scale-free cyclic models.

In bouncing cosmology, the primordial fluctuations are generated in a cosmic contraction phase before the bounce into the current expansion phase. For a nonsingular bounce, curvature and anisotropy grow rapidly during the bouncing phase, raising questions about the reliability of perturbative analysis. In this paper, we study the evolution of adiabatic perturbations in a nonsingular bounce by nonperturbative methods including numerical simulations of the nonsingular bounce and the covariant formalism for calculating nonlinear perturbations. We show that the bounce is disrupted in regions of the universe with significant inhomogeneity and anisotropy over the background energy density, but is achieved in regions that are relatively homogeneous and isotropic. Sufficiently small perturbations, consistent with observational constraints, can pass through the nonsingular bounce with negligible alteration from nonlinearity. We follow scale invariant perturbations generated in a matter-like contraction phase through the bounce. Their amplitude in the expansion phase is determined by the growing mode in the contraction phase, and the scale invariance is well preserved across the bounce.

Recent measurements at the LHC suggest that the current Higgs vacuum could be metastable with a modest barrier (height 10^{10-12}{GeV})^{4}) separating it from a ground state with negative vacuum density of order the Planck scale. We note that metastability is problematic for big bang to end one cycle, bounce, and begin the next. In this paper, motivated by the approximate scaling symmetry of the standard model of particle physics and the primordial large-scale structure of the universe, we use our recent formulation of the Weyl-invariant version of the standard model coupled to gravity to track the evolution of the Higgs in a regularly bouncing cosmology. We find a band of solutions in which the Higgs field escapes from the metastable phase during each big crunch, passes through the bang into an expanding phase, and returns to the metastable vacuum, cycle after cycle after cycle. We show that, due to the effect of the Higgs, the infinitely cycling universe is geodesically complete, in contrast to inflation.

We show that results from the Planck satellite reported in 2013 are consistent with the simplest cyclic models for natural parameter ranges i.e., order unity dimensionless coefficients, assuming the standard entropic mechanism for generating curvature perturbations. With improved precision, forthcoming results from Planck and other experiments should be able to test the parameter ranges by confirming or refuting the core predictions - i.e., no observable primordial B-mode polarization and detectable local non-gaussianity. A new prediction, given the Planck 2013 constraints on the bispectrum, is a sharp constraint on the local trispectrum parameter g_{NL}; namely, the simplest models predict it is negative, with g_{NL} < -1700.

Recent results from the Planck satellite combined with earlier observations from WMAP, ACT, SPT and other experiments eliminate a wide spectrum of more complex inflationary models and favor models with a single scalar field, as reported by the Planck Collaboration. More important, though, is that all the simplest inflaton models are disfavored statistically relative to those with plateau-like potentials. We discuss how a restriction to plateau-like models has three independent serious drawbacks: it exacerbates both the initial conditions problem and the multiverse-unpredictability problem and it creates a new difficulty that we call the inflationary "unlikeliness problem." Finally, we comment on problems reconciling inflation with a standard model Higgs, as suggested by recent LHC results. In sum, we find that recent experimental data disfavors all the best-motivated inflationary scenarios and introduces new, serious difficulties that cut to the core of the inflationary paradigm. Forthcoming searches for B-modes, non-Gaussianity and new particles should be decisive.

We study a model of a scalar field minimally coupled to gravity, with a specific potential energy for the scalar field, and include curvature and radiation as two additional parameters. Our goal is to obtain analytically the complete set of configurations of a homogeneous and isotropic universe as a function of time. This leads to a geodesically complete description of the universe, including the passage through the cosmological singularities, at the classical level. We give all the solutions analytically without any restrictions on the parameter space of the model or initial values of the fields. We find that for generic solutions the universe goes through a singular (zero-size) bounce by entering a period of antigravity at each big crunch and exiting from it at the following big bang. This happens cyclically again and again without violating the null energy condition. There is a special subset of geodesically complete non-generic solutions which perform zero-size bounces without ever entering the antigravity regime in all cycles. For these, initial values of the fields are synchronized and quantized but the parameters of the model are not restricted. There is also a subset of spatial curvature-induced solutions that have finite-size bounces in the gravity regime and never enter the antigravity phase. These exist only within a small continuous domain of parameter space without fine tuning initial conditions. To obtain these results, we identified 25 regions of a 6-parameter space in which the complete set of analytic solutions are explicitly obtained.

We consider bouncing cosmologies in which an ekpyrotic contraction phase with w >> 1 is followed by a bouncing phase with w < -1 that violates the null energy condition. The bouncing phase, induced by ghost condensation, is designed to produce a classically nonsingular bounce at a finite value of the scale factor. We show that the initial curvature and anisotropy, though diluted during the ekpyrotic phase, grow back exponentially during the bouncing phase. Moreover, curvature perturbations and anisotropy are generated by quantum fluctuations during the ekpyrotic phase. In the bouncing phase, however, an adiabatic curvature perturbation grows to dominate and contributes a blue spectrum that spoils the scale-invariance. Meanwhile, a scalar shear perturbation grows nonlinear and creates an overwhelming anisotropy that disrupts the nonsingular bounce altogether.

In inflationary models, the predicted amplitude of primordial density perturbations Q is much larger than the observed value (~10^{-5}) for natural choices of parameters. To explain the requisite exponential fine-tuning, anthropic selection is often invoked, especially in cases where microphysics is expected to produce a complex energy landscape. By contrast, we find examples of ekpyrotic models based on heterotic M-theory for which dynamical selection naturally favors the observed value of Q.

We note that the essential idea of inflation, that the universe underwent a brief period of accelerated expansion followed by a long period of decelerated expansion, can be encapsulated in a "closure condition" which relates the amount of accelerated expansion during inflation to the amount of decelerated expansion afterward. We present a protocol for systematically testing the validity of this condition observationally.

We consider non-singular bouncing cosmologies, such as the new ekpyrotic model, in which the universe undergoes a slow contraction phase with equation of state $w \gg 1$, followed by a bounce that occurs at a finite scale factor when quantum gravity corrections are still negligible. Such a non-singular bounce requires a violation of the null energy condition in which $w$ falls below -1 at some time before the bounce. In this paper, we show that a component of the adiabatic curvature perturbations, though decaying and negligible during the ekpyrotic phase, is exponentially amplified just before $w$ approaches -1, enough to spoil the scale-invariant perturbation spectrum. We discuss how the problem may be avoided, for example, in singular bounces.

We recently presented a series of dark energy theorems that place constraints on the equation of state of dark energy ($\wdark$), the ime-variation of Newton's constant ($\dot G$), and the violation of energy conditions in theories with extra dimensions. In this paper, we explore how current and future measurements of $\wdark$ and $\dot G$ can be used to place tight limits on large classes of these theories (including some of the most well-motivated examples) independent of the size of the extra dimensions. As an example, we show that models with conformally Ricci-flat metrics obeying the null energy condition (a common ansatz for Kaluza-Klein and string constructions) are highly constrained by current ata and may be ruled out entirely by future dark energy and pulsar observations.

The universe can be made flat and smooth by undergoing a phase of ultra-slow (ekpyrotic) contraction, a condition achievable with a single, canonical scalar field and conventional general relativity. It has been argued, though, that generating scale-invariant density perturbations, requires at least two scalar fields and a two-step process that first produces entropy fluctuations and then converts them to curvature perturbations. In this paper, we identify a loophole in the argument and introduce an ekpyrotic model based on a single, canonical scalar field that generates nearly scale-invariant curvature fluctuations through a purely "adiabatic mechanism" in which the background evolution is a dynamical attractor. The resulting spectrum can be slightly red with distinctive non-gaussian fluctuations.

Incorporating the QCD axion and simultaneously satisfying current constraints on the dark matter density and isocurvature fluctuations requires non-minimal fine-tuning of inflationary parameters or the axion misalignment angle (or both) for Peccei-Quinn symmetry-breaking scales $f_a > 10^{12}$ GeV. To gauge the degree of tuning in models with many axion-like fields at similar symmetry-breaking scales and masses, as may occur in string theoretic models that include a QCD axion, we introduce a figure of merit ${\cal F}$ that measures the fractional volume of allowed parameter space: the product of the slow roll parameter $\epsilon$ and each of the axion misalignment angles, $\theta_0$. For a single axion, $\mathcal{F} \lesssim 10^{-11}$ is needed to avoid conflict with observations. We show that the fine tuning of $\mathcal{F}$ becomes exponentially more extreme in the case of numerous axion-like fields. Anthropic arguments are insufficient to explain the fine tuning because the bulk of the anthropically allowed parameter space is observationally ruled out by limits on the cosmic microwave background isocurvature modes. Therefore, this tuning presents a challenge to the compatibility of string-theoretic models with light axions and inflationary cosmology.

Georges Lemaitre introduced the term "phoenix universe" to describe an oscillatory cosmology with alternating periods of gravitational collapse and expansion. This model is ruled out observationally because it requires a supercritical mass density and cannot accommodate dark energy. However, a new cyclic theory of the universe has been proposed that evades these problems. In a recent elaboration of this picture, almost the entire universe observed today is fated to become entrapped inside black holes, but a tiny region will emerge from these ashes like a phoenix to form an even larger smooth, flat universe filled with galaxies, stars, planets, and, presumably, life. Survival depends crucially on dark energy and suggests a reason why its density is small and positive today.

Non-gaussianity in the microwave background radiation is bound to play a key role in giving us clues about the physics of the very early universe. However, the associated calculations, at second and even third order in perturbation theory, tend to be complicated to the point of obscuring simple underlying physical processes. In this note, we present a simple analytic procedure for approximating the non-linearity parameters f_{NL} and g_{NL} for cyclic models in which the cosmological perturbations are generated via the entropic mechanism. Our approach is quick, physically transparent and agrees well with the results of numerical calculations.

We analyze the non-gaussian density perturbations generated in ekpyrotic/cyclic models based on heterotic M-theory. In this picture, two scalar fields produce nearly scale-invariant entropic perturbations during an ekpyrotic phase that are converted into curvature modes {\it after the ekpyrotic phase is complete} and just before the big bang. Both intrinsic non-linearity in the entropy perturbation and the conversion process contribute to non-gaussianity. The range of the non-gaussianity parameter $f_{NL}$ depends on how gradual the conversion process is and the steepness of the scalar field potential during the ekpyrotic phase. Although a wider range is possible, in principle, natural values of the ekpyrotic parameters combined with a gradual conversion process lead to values of $-60 \lesssim f_{NL} \lesssim +80$, typically much greater than slow-roll inflation but within the current observational bounds.

It has been pointed out by several groups that ekpyrotic and cyclic models generate significant non-gaussianity. In this paper, we present a physically intuitive, semi-analytic estimate of the bispectrum. We show that, in all such models, there is an intrinsic contribution to the non-gaussianity parameter f_{NL} that is determined by the geometric mean of the equation of state w_{ek} during the ekpyrotic phase and w_{c} during the phase that curvature perturbations are generated and whose value is O(100) or more times the intrinsic value predicted by simple slow-roll inflationary models, f_{NL}^{intrinsic} = O(0.1). Other contributions to f_{NL}, which we also estimate, can increase |f_{NL}| but are unlikely to decrease it significantly, making non-gaussianity a useful test of these models. Furthermore, we discuss a predicted correlation between the non-gaussianity and scalar spectral index that sharpens the test.

A new and powerful probe of the origin and evolution of structures in the Universe has emerged and been actively developed over the last decade. In the coming decade, non-Gaussianity, i.e., the study of non-Gaussian contributions to the correlations of cosmological fluctuations, will become an important probe of both the early and the late Universe. Specifically, it will play a leading role in furthering our understanding of two fundamental aspects of cosmology and astrophysics: (i) the physics of the very early universe that created the primordial seeds for large-scale structures, and (ii) the subsequent growth of structures via gravitational instability and gas physics at later times. To date, observations of fluctuations in the Cosmic Microwave Background (CMB) and the Large-Scale Structure of the Universe (LSS) have focused largely on the Gaussian contribution as measured by the two-point correlations (or the power spectrum) of density fluctuations. However, an even greater amount of information is contained in non-Gaussianity and a large discovery space therefore still remains to be explored. Many observational probes can be used to measure non-Gaussianity, including CMB, LSS, gravitational lensing, Lyman-alpha forest, 21-cm fluctuations, and the abundance of rare objects such as clusters of galaxies and high-redshift galaxies. Not only does the study of non-Gaussianity maximize the science return from a plethora of present and future cosmological experiments and observations, but it also carries great potential for important discoveries in the coming decade.

In cyclic universe models based on a single scalar field (e.g., the radion determining the distance between branes in M-theory), virtually the entire universe makes it through the ekpyrotic smoothing and flattening phase, bounces, and enters a new epoch of expansion and cooling. This stable evolution cannot occur, however, if scale-invariant curvature perturbations are produced by the entropic mechanism because it requires two scalar fields (e.g., the radion and the Calabi-Yau dilaton) evolving along an unstable classical trajectory. In fact, we show here that an overwhelming fraction of the universe fails to make it through the ekpyrotic phase; nevertheless, a sufficient volume survives and cycling continues forever provided the dark energy phase of the cycle lasts long enough, of order a trillion years. Two consequences are a new role for dark energy and a global structure of the universe radically different from that of eternal inflation.

A period of slow contraction with equation of state w > 1, known as an ekpyrotic phase, has been shown to flatten and smooth the universe if it begins the phase with small perturbations. In this paper, we explore how robust and powerful the ekpyrotic smoothing mechanism is by beginning with highly inhomogeneous and anisotropic initial conditions and numerically solving for the subsequent evolution of the universe. Our studies, based on a universe with gravity plus a scalar field with a negative exponential potential, show that some regions become homogeneous and isotropic while others exhibit inhomogeneous and anisotropic behavior in which the scalar field behaves like a fluid with w=1. We find that the ekpyrotic smoothing mechanism is robust in the sense that the ratio of the proper volume of the smooth to non-smooth region grows exponentially fast along time slices of constant mean curvature.

Although the rocket-based X-ray Quantum Calorimetry (XQC) experiment was designed for X-ray spectroscopy, the minimal shielding of its calorimeters, its low atmospheric overburden, and its low-threshold detectors make it among the most sensitive instruments for detecting or constraining strong interactions between dark matter particles and baryons. We use Monte Carlo simulations to obtain the precise limits the XQC experiment places on spin-independent interactions between dark matter and baryons, improving upon earlier analytical estimates. We find that the XQC experiment rules out a wide range of nucleon-scattering cross sections centered around one barn for dark matter particles with masses between 0.01 and 10^5 GeV. Our analysis also provides new constraints on cases where only a fraction of the dark matter strongly interacts with baryons.

We investigate whether explicit models of warped D-brane inflation are possible in string compactifications. To this end, we study the potential for D3-brane motion in a warped conifold that includes holomorphically-embedded D7-branes involved in moduli stabilization. The presence of the D7-branes significantly modifies the inflaton potential. We construct an example based on a very simple and symmetric embedding due to Kuperstein, z_1 = constant, in which it is possible to fine-tune the potential so that slow roll inflation can occur. The resulting model is rather delicate: inflation occurs in the vicinity of an inflection point, and the cosmological predictions are extremely sensitive to the precise shape of the potential.

We derive the complete spectrum of gravitational waves induced by primordial scalar perturbations ranging over all observable wavelengths. This scalar-induced contribution can be computed directly from the observed scalar perturbations and general relativity and is, in this sense, independent of the cosmological model for generating the perturbations. The spectrum is scale-invariant on small scales, but has an interesting scale-dependence on large and intermediate scales, where scalar-induced gravitational waves do not redshift and are hence enhanced relative to the background density of the Universe. This contribution to the tensor spectrum is significantly different in form from the direct model-dependent primordial tensor spectrum and, although small in magnitude, it dominates the primordial signal for some cosmological models. We confirm our analytical results by direct numerical integration of the equations of motion.

We reconsider the possibility that the observed baryon asymmetry was generated by the evaporation of primordial black holes that dominated the early universe. We present a simple derivation showing that the baryon asymmetry is insensitive to the initial black hole density and the cosmological model but is sensitive to the temperature-dependence of the CP and baryon-violating (or lepton-violating) interactions. We also consider the possibility that black holes stop evaporating and form Planck-mass remnants that act as dark matter. We show that primordial black holes cannot simultaneously account for both the observed baryon asymmetry and the (remnant) dark matter density unless the magnitude of CP violation is much greater than expected from most particle physics models. Finally, we apply these results to ekpyrotic/cyclic models, in which primordial black holes may form when branes collide. We find that obtaining the observed baryon asymmetry is compatible with the other known constraints on parameters.

We analyze a general mechanism for producing a nearly scale-invariant spectrum of cosmological curvature perturbations during a contracting phase preceding a big bang, that can be entirely described using 4d effective field theory. The mechanism, based on first producing entropic perturbations and then converting them to curvature perturbations, can be naturally incorporated in cyclic and ekpyrotic models in which the big bang is modelled as a brane collision, as well as other types of cosmological models with a pre-big bang phase. We show that the correct perturbation amplitude can be obtained and that the scalar spectral tilt n tends to range from slightly blue to red, with 0.97 < n < 1.02 for the simplest models, a range compatible with current observations but shifted by a few per cent towards the blue compared to the prediction of the simplest, large-field inflationary models.

We analyze the evolution of cosmological perturbations in the cyclic model, paying particular attention to their behavior and interplay over multiple cycles. Our key results are: (1) galaxies and large scale structure present in one cycle are generated by the quantum fluctuations in the preceding cycle without interference from perturbations or structure generated in earlier cycles and without interfering with structure generated in later cycles; (2) the ekpyrotic phase, an epoch of gentle contraction with equation of state $w\gg 1$ preceding the hot big bang, makes the universe homogeneous, isotropic and flat within any given observer's horizon; and, (3) although the universe is uniform within each observer's horizon, the global structure of the cyclic universe is more complex, owing to the effects of superhorizon length perturbations, and cannot be described in a uniform Friedmann-Robertson-Walker picture. In particular, we show that the ekpyrotic phase is so effective in smoothing, flattening and isotropizing the universe within the horizon that this phase alone suffices to solve the horizon and flatness problems even without an extended period of dark energy domination (a kind of low energy inflation). Instead, the cyclic model rests on a genuinely novel, non-inflationary mechanism (ekpyrotic contraction) for resolving the classic cosmological conundrums.

Within conventional big bang cosmology, it has proven to be very difficult to understand why today's cosmological constant is so small. In this paper, we show that a cyclic model of the universe can naturally incorporate a dynamical mechanism that automatically relaxes the value of the cosmological constant, taking account of contributions to the vacuum density at all energy scales. Because the relaxation time grows exponentially as the vacuum density decreases, nearly every volume of space spends an overwhelming majority of the time at the stage when the cosmological constant is small and positive, as observed today.

We reconsider the predictions of inflation for the spectral index of scalar (energy density) fluctuations (n_s) and the tensor/scalar ratio (r) using a discrete, model-independent measure of the degree of fine-tuning required to obtain a given combination of (n_s, r). We find that, except for cases with numerous unnecessary degrees of fine-tuning, n_s is less than 0.98, measurably different from exact Harrison-Zel'dovich. Furthermore, if n_s \gtrsim 0.95, in accord with current measurements, the tensor/scalar ratio satisfies r \gtrsim 10^{-2}, a range that should be detectable in proposed cosmic microwave background polarization experiments and direct gravitational wave searches.

Near comoving wavenumber k, the gravitational-wave background (GWB) from inflation carries information about the physical conditions near two moments in cosmic history: the moment when k ``left the horizon'' during inflation, and the moment when it ``re-entered the horizon'' after inflation. We investigate the extent to which this information can be extracted if the GWB is measured by a combination of cosmic-microwave-background (CMB) polarization experiments on large scales and space-based laser-interferometer experiments on small scales. To disentangle this information, we derive a new gravitational-wave transfer function that incorporates a number of physical effects that were treated less accurately, less generally, or were missing altogether in previous treatments. In particular, it incorporates: (i) dark energy with time-varying equation-of-state w(z); (ii) tensor anisotropic stress due to free-streaming relativistic particles in the early universe; and (iii) a variety of physical effects that cause deviations from the standard equation-of-state w=1/3 during the radiation era. Based on this transfer function, we consider the degree to which the GWB can be used to test inflation and to probe the ``primordial dark age'' between the end of inflation and the electroweak phase transition.

We consider the effect of compactification of extra dimensions on the onset of classical chaotic "Mixmaster" behavior during cosmic contraction. Assuming a universe that is well-approximated as a four-dimensional Friedmann-Robertson--Walker model (with negligible Kaluza-Klein excitations) when the contraction phase begins, we identify compactifications that allow a smooth contraction and delay the onset of chaos until arbitrarily close the big crunch. These compactifications are defined by the de Rham cohomology (Betti numbers) and Killing vectors of the compactification manifold. We find compactifications that control chaos in vacuum Einstein gravity, as well as in string theories with N = 1 supersymmetry and M-theory. In models where chaos is controlled in this way, the universe can remain homogeneous and flat until it enters the quantum gravity regime. At this point, the classical equations leading to chaotic behavior can no longer be trusted, and quantum effects may allow a smooth approach to the big crunch and transition into a subsequent expanding phase. Our results may be useful for constructing cosmological models with contracting phases, such as the ekpyrotic/cyclic and pre-big bang models.

We consider how well the dark energy equation of state $w$ as a function of red shift $z$ will be measured using current and anticipated experiments. We use a procedure which takes fair account of the uncertainties in the functional dependence of $w$ on $z$, as well as the parameter degeneracies, and avoids the use of strong prior constraints. We apply the procedure to current data from WMAP, SDSS, and the supernova searches, and obtain results that are consistent with other analyses using different combinations of data sets. The effects of systematic experimental errors and variations in the analysis technique are discussed. Next, we use the same procedure to forecast the dark energy constraints achieveable by the end of the decade, assuming 8 years of WMAP data and realistic projections for ground-based measurements of supernovae and weak lensing. We find the $2 \sigma$ constraints on the current value of $w$ to be $\Delta w_0 (2 \sigma) = 0.20$, and on $dw/dz$ (between $z=0$ and $z=1$) to be $\Delta w_1 (2 \sigma)=0.37$. Finally, we compare these limits to other projections in the literature. Most show only a modest improvement; others show a more substantial improvement, but there are serious concerns about systematics. The remaining uncertainty still allows a significant span of competing dark energy models. Most likely, new kinds of measurements, or experiments more sophisticated than those currently planned, are needed to reveal the true nature of dark energy.

We consider a picture in which the transition from a big crunch to a big bang corresponds to the collision of two empty orbifold planes approaching each other at a constant non-relativistic speed in a locally flat background space-time, a situation relevant to recently proposed cosmological models. We show that $p$-brane states which wind around the extra dimension propagate smoothly and unambiguously across the orbifold plane collision. In particular we calculate the quantum mechanical production of winding M2-branes extending from one orbifold to the other. We find that the resulting density is finite and that the resulting gravitational back-reaction is small. These winding states, which include the string theory graviton, can be propagated smoothly across the transition using a perturbative expansion in the membrane tension, an expansion which from the point of view of string theory is an expansion in {\it inverse} powers of $\alpha'$. We argue that interactions should be well-behaved because the string coupling tends to zero at the crunch. The production of massive Kaluza-Klein states should also be exponentially suppressed for small collision speeds. We contrast this good behavior with that found in previous studies of strings in Lorentzian orbifolds.

The Cyclic Model attempts to resolve the homogeneity, isotropy, and flatness problems and generate a nearly scale-invariant spectrum of fluctuations during a period of slow contraction that precedes a bounce to an expanding phase. Here we describe at a conceptual level the recent developments that have greatly simplified our understanding of the contraction phase and the Cyclic Model overall. The answers to many past questions and criticisms are now understood. In particular, we show that the contraction phase has equation of state w>1 and that contraction with w>1 has a surprisingly similar properties to inflation with w < -1/3. At one stroke, this shows how the model is different from inflation and why it may work just as well as inflation in resolving cosmological problems.

We consider how quintessence models in which the sound speed differs from the speed of light and varies with time affect the cosmic microwave background and the fluctuation power spectrum. Significant modifications occur on length scales related to the Hubble radius during epochs in which the sound speed is near zero and the quintessence contributes a non-negligible fraction of the total energy density. For the microwave background, we find that the usual enhancement of the lowest multipole moments by the integrated Sachs-Wolfe effect can be modified, resulting in suppression or bumps instead. Also, the sound speed can produce oscillations and other effects at wavenumbers $k > 10^{-2}$ h/Mpc in the fluctuation power spectrum.

We consider cosmological models with a scalar field with equation of state $w\ge 1$ that contract towards a big crunch singularity, as in recent cyclic and ekpyrotic scenarios. We show that chaotic mixmaster oscillations due to anisotropy and curvature are suppressed, and the contraction is described by a homogeneous and isotropic Friedmann equation if $w>1$. We generalize the results to theories where the scalar field couples to p-forms and show that there exists a finite value of $w$, depending on the p-forms, such that chaotic oscillations are suppressed. We show that $Z_2$ orbifold compactification also contributes to suppressing chaotic behavior. In particular, chaos is avoided in contracting heterotic M-theory models if $w>1$ at the crunch.

For a 4-dimensional spatially-flat Friedmann-Robertson-Walker universe with a scalar field $\phi(x)$, potential $V(\phi)$ and constant equation of state $w=p/\rho$, we show that an expanding solution characterized by $\epsilon=3(1+w)/2$ produces the same scalar perturbations as a contracting solution with $\hat{\epsilon}=1/\epsilon$. The same symmetry applies to both the dominant and subdominant scalar perturbation modes. This result admits a simple physical interpretation and generalizes to $d$ spacetime dimensions if we define $\epsilon \equiv [(2d-5)+(d-1)w]/(d-2)$.

We show that a gravitational interaction between the derivative of the Ricci scalar curvature and the baryon-number current dynamically breaks CPT in an expanding universe and, combined with baryon-number-violating interactions, can drive the universe towards an equilibrium baryon asymmetry that is observationally acceptable.

Recent advances in understanding the propagation of perturbations through the transition from big crunch to big bang (esp. Tolley et al. hep-th/0306109) make it possible for the first time to consider the full set of phenomenological constraints on the scalar field potential in cyclic models of the universe. We show that cyclic models require a comparable degree of tuning to that needed for inflationary models. The constraints are reduced to a set of simple design rules including "fast-roll" parameters analogous to the "slow-roll" parameters in inflation.

We present a simple, nearly model-independent estimate that yields the predictions of the simplest inflationary and ekpyrotic/cyclic models for the spectral tilt of the primordial density inhomogeneities. Remarkably, we find that the simplest models yield an identical result: n_s is approximately 0.95. For inflation, the same estimate predicts a ratio of tensor to scalar contributions to the low multipoles of the microwave background anisotropy of T/S=20%; the tensor contribution is negligible for ekpyrotic/cyclic models.

Inflation predicts a primordial gravitational wave spectrum that is slightly ``red,'' i.e., nearly scale-invariant with slowly increasing power at longer wavelengths. In this paper, we compute both the amplitude and spectral form of the primordial tensor spectrum predicted by cyclic/ekpyrotic models. The spectrum is blue and exponentially suppressed compared to inflation on long wavelengths. The strongest observational constraint emerges from the requirement that the energy density in gravitational waves should not exceed around 10 per cent of the energy density at the time of nucleosynthesis.

Dark matter, proposed decades ago as a speculative component of the universe, is now known to be the vital ingredient in the cosmos, eight times more abundant than ordinary matter, one quarter of the total energy density and the component which has controlled the growth of structure in the universe. Its nature remains a mystery, but, assuming it is comprised of weakly interacting sub-atomic particles, is consistent with large scale cosmic structure. However, recent analyses of structure on galactic and sub-galactic scales have suggested discrepancies and stimulated numerous alternative proposals. We discuss how studies of the density, demography, history and environment of smaller scale structures may distinguish among these possibilities and shed new light on the nature of dark matter.

A prescription is developed for matching general relativistic perturbations across singularities of the type encountered in the ekpyrotic and cyclic scenarios i.e. a collision between orbifold planes. We show that there exists a gauge in which the evolution of perturbations is locally identical to that in a model space-time (compactified Milne mod Z_2) where the matching of modes across the singularity can be treated using a prescription previously introduced by two of us. Using this approach, we show that long wavelength, scale-invariant, growing-mode perturbations in the incoming state pass through the collision and become scale-invariant growing-mode perturbations in the expanding hot big bang phase.

The recent announcement from the Wilkinson Microwave Anisotropy Probe (WMAP) satellite experiment combined with other recent advances in observational cosmology verifies key components of the standard cosmological model. However, we argue that there remain some significant open issues regarding the basic history and composition of the Universe and uncertainties in some of the most important parameters.

We analyze the general conditions on the equation of state $w$ required for quantum fluctuations of a scalar field to produce a scale-invariant spectrum of density perturbations, including models which (in the four dimensional effective description) bounce from a contracting to an expanding phase. We show that there are only two robust cases: $w\approx -1$ (inflation) and $w \gg 1$ (the ekpyrotic/cyclic scenario). All other cases, including the $w \approx 0$ case considered by some authors, require extreme fine-tuning of initial conditions and/or the effective potential. For the ekpyrotic/cyclic ($w \gg 1$) case, we also analyze the small deviations from scale invariance.

It has recently been argued that bubble nucleation in ekpyrotic and cyclic cosmological scenarios can lead to unacceptable inhomogeneities unless certain constraints are satisfied. In this paper we show that this is not the case. We find that bubble nucleation is completely negligible in realistic models.

We propose a cosmological model in which the universe undergoes an endless sequence of cosmic epochs each beginning with a `bang' and ending in a `crunch.' The temperature and density are finite at each transition from crunch to bang. Instead of having an inflationary epoch, each cycle includes a period of slow accelerated expansion (as recently observed) followed by slow contraction. The combination produces the homogeneity, flatness, density fluctuations and energy needed to begin the next cycle.

We analyze the evolution of energy density fluctuations in cosmological scenarios with a mixture of cold dark matter and quintessence, in which the quintessence field is modeled by a constant equation of state. We obtain analytic expressions for the time evolution of the quintessence perturbations in models with light fields. The fluctuations behave analogously to a driven harmonic oscillator, where the driving term arises from the inhomogeneities in the surrounding cosmological fluid. We demonstrate that the homogeneous solution, determined by the initial conditions, is completely sub-dominant to the inhomogeneous solution for physically realistic scenarios. Thus we show that the cosmic microwave background anisotropy predicted for such models is highly insensitive to the initial conditions in the quintessence field.

The Cyclic Model is a radical, new cosmological scenario which proposes that the Universe undergoes an endless sequence of epochs which begin with a `big bang' and end in a `big crunch.' When the Universe bounces from contraction to re-expansion, the temperature and density remain finite. The model does not include a period of rapid inflation, yet it reproduces all of the successful predictions of standard big bang and inflationary cosmology. We point out numerous novel elements that have not been used previously which may open the door to further alternative cosmologies. Although the model is motivated by M-theory, branes and extra-dimensions, here we show that the scenario can be described almost entirely in terms of conventional 4d field theory and 4d cosmology.

Based on concepts drawn from the ekpyrotic scenario and M-theory, we elaborate our recent proposal of a cyclic model of the Universe. In this model, the Universe undergoes an endless sequence of cosmic epochs which begin with the Universe expanding from a `big bang' and end with the Universe contracting to a `big crunch.' Matching from `big crunch' to `big bang' is performed according to the prescription recently proposed with Khoury, Ovrut and Seiberg. The expansion part of the cycle includes a period of radiation and matter domination followed by an extended period of cosmic acceleration at low energies. The cosmic acceleration is crucial in establishing the flat and vacuous initial conditions required for ekpyrosis and for removing the entropy, black holes, and other debris produced in the preceding cycle. By restoring the Universe to the same vacuum state before each big crunch, the acceleration insures that the cycle can repeat and that the cyclic solution is an attractor.

We study the generation of density perturbations in the ekpyrotic scenario for the early universe, including gravitational backreaction. We expose interesting subtleties that apply to both inflationary and ekpyrotic models. Our analysis includes a detailed proposal of how the perturbations generated in a contracting phase may be matched across a `bounce' to those in an expanding hot big bang phase. For the physical conditions relevant to the ekpyrotic scenario, we re-obtain our earlier result of a nearly scale-invariant spectrum of energy density perturbations. We find that the perturbation amplitude is typically small, as desired to match observation.

We consider conditions under which a universe contracting towards a big crunch can make a transition to an expanding big bang universe. A promising example is 11-dimensional M-theory in which the eleventh dimension collapses, bounces, and re-expands. At the bounce, the model can reduce to a weakly coupled heterotic string theory and, we conjecture, it may be possible to follow the transition from contraction to expansion. The possibility opens the door to new classes of cosmological models. For example, we discuss how it suggests a major simplification and modification of the recently proposed ekpyrotic scenario.

We explore various pitfalls and challenges in determining the equation-of-state (w) of dark energy component that dominates the universe and causes the current accelerated expansion. We demonstrated in an earlier paper the existence of a degeneracy that makes it impossible to resolve well the value of w or its time-derivative with supernovae data. Here we consider standard practices, such as assuming priors that w is constant or greater than -1, and show that they also can lead to gross errors in estimating the true equation-of-state. We further consider combining measurements of the cosmic microwave background anisotropy and the Alcock-Paczynski test with supernovae data and find that the improvement in resolving the time-derivative of w is marginal, although the combination can constrain its present value perhaps to 20 percent uncertainty.

Quintessence, a time-varying energy component that may account for the accelerated expansion of the universe, can be characterized by its equation of state and sound speed. In this paper, we show that if the quintessence density is at least one percent of the critical density at the surface of last scattering the cosmic microwave background anisotropy can distinguish between models whose sound speed is near the speed of light versus near zero, which could be useful in distinguishing competing candidates for dark energy.

We propose a cosmological scenario in which the hot big bang universe is produced by the collision of a brane in the bulk space with a bounding orbifold plane, beginning from an otherwise cold, vacuous, static universe. The model addresses the cosmological horizon, flatness and monopole problems and generates a nearly scale-invariant spectrum of density perturbations without invoking superluminal expansion (inflation). The scenario relies, instead, on physical phenomena that arise naturally in theories based on extra dimensions and branes. As an example, we present our scenario predominantly within the context of heterotic M-theory. A prediction that distinguishes this scenario from standard inflationary cosmology is a strongly blue gravitational wave spectrum, which has consequences for microwave background polarization experiments and gravitational wave detectors.

We show that non-topological solitons, known as Q-balls, are promising candidates for self-interacting dark matter. They can satisfy the cross-section requirements for a broad range of masses. Unlike previously considered examples, Q-balls can stick together after collision, reducing the effective self-interaction rate to a negligible value after a few collisions per particle. This feature modifies predictions for halo formation. We also discuss the possibility that Q-balls have large interaction cross-sections with ordinary matter.

In the early 1990's, an analysis was completed by several theorists of the available mass/cross-section parameter space for unusual particle candidates to solve the dark matter problem, e.g. strongly interacting massive particles (SIMPs). This analysis found several unconstrained windows, such that for SIMP masses and cross-sections within these windows, SIMPs could still be the dominant dark matter in our Galactic halo. Since the early 1990's, some of these windows have been narrowed or closed, and some of these windows have been widened further by more careful analysis. We summarize the present state of the SIMP parameter space, and point to the cosmological salience of SIMPs as dark matter, given some of the present inadequacies of the WIMP solution to the dark matter problem.

It is shown that there exist large classes of BPS vacua in heterotic M-theory which have negative tension on the visible orbifold plane, positive tension on the hidden plane and positive tension, physical five-branes in the bulk space. Explicit examples of such vacua are presented. Furthermore, it is demonstrated that the ratio, beta/|alpha|, of the bulk five-brane tension to the visible plane tension can, for several large classes of such vacua, be made arbitrarily small. Hence, it is straightforward to find vacua with the properties required in the examples of the Ekpyrotic theory of cosmology - a visible brane with negative tension and beta/|alpha| small. This contradicts recent claims in the literature.

We respond to the criticisms by Kallosh, Kofman and Linde concerning our proposal of the ekpyrotic universe scenario. We point out a number of errors in their considerations and argue that, at this stage, the ekpyrotic model is a possible alternative to inflationary cosmology as a description of the very early universe.

Supernova searches have been been suggested as a method for determining precisely the current value and time variation of the equation of state, $w$, of the dark energy component responsible for the accelerated expansion of the Universe. We show that the method is fundamentally limited by the fact that luminosity dista nce depends on $w$ through a multiple integral relation that smears out information about $w$ and its time variatio n. The effect degrades the resolution of $w$ that can be obtained from current d ata.

Spergel and Steinhardt have recently proposed the concept of dark matter with strong self-interactions as a means to address numerous discrepancies between observations of dark matter halos on subgalactic scales and the predictions of the standard collisionless dark matter picture. We review the motivations for this scenario and discuss some recent, successful numerical tests. We also discuss the possibility that the dark matter interacts strongly with ordinary baryonic matter, as well as with itself. We present a new analysis of the experimental constraints and re-evaluate the allowed range of cross-section and mass.

We recently introduced the concept of "k-essence" as a dynamical solution for explaining naturally why the universe has entered an epoch of accelerated expansion at a late stage of its evolution. The solution avoids fine-tuning of parameters and anthropic arguments. Instead, k-essence is based on the idea of a dynamical attractor solution which causes it to act as a cosmological constant only at the onset of matter-domination. Consequently, k-essence overtakes the matter density and induces cosmic acceleration at about the present epoch. In this paper, we present the basic theory of k-essence and dynamical attractors based on evolving scalar fields with non-linear kinetic energy terms in the action. We present guidelines for constructing concrete examples and show that there are two classes of solutions, one in which cosmic acceleration continues forever and one in which the acceleration has finite duration.

We present a comparison of halo properties in cosmological simulations of collisionless cold dark matter (CDM) and self-interacting dark matter (SIDM) for a range of dark matter cross sections. We find, in agreement with various authors, that CDM yields cuspy halos that are too centrally concentrated as compared to observations. Conversely, SIDM simulations using a Monte Carlo N-body technique produce halos with significantly reduced central densities and flatter cores with increasing cross section. We introduce a concentration parameter based on enclosed mass that we expect will be straightforward to determine observationally, unlike that of Navarro, Frenk & White, and provide predictions for SIDM and CDM. SIDM also produces more spherical halos than CDM, providing possibly the strongest observational test of SIDM. We discuss our findings in relation to various relevant observations as well as SIDM simulations of other groups. Taking proper account of simulation limitations, we find that a dark matter cross section per unit mass of sigma_DM ~= 10^{-23}-10^{-24} cm^2/GeV is consistent with all current observational constraints.

We consider the cosmology of a pair of domain walls bounding a five-dimensional bulk space-time with negative cosmological constant, in which the distance between the branes is not fixed in time. Although there are strong arguments to suggest that this distance should be stabilized in the present epoch, no such constraints exist for the early universe and thus non-static solutions might provide relevant inflationary scenarios. We find the general solution for the standard ansatz where the bulk is foliated by planar-symmetric hypersurfaces. We show that in all cases the bulk geometry is that of anti-de Sitter (AdS_5). We then present a geometrical interpretation for the solutions as embeddings of two de Sitter (dS_4) surfaces in AdS_5, which provide a simple interpretation of the physical properties of the solutions. A notable feature explained in the analysis is that two-way communication between branes expanding away from one another is possible for a finite amount of time, after which communication can proceed in one direction only. The geometrical picture also shows that our class of solutions (and related solutions in the literature) are not completely general, contrary to some claims. We then derive the most general solution for two walls in AdS_5. This includes novel cosmologies where the brane tensions are not constrained to have opposite signs. The construction naturally generalizes to arbitrary FRW cosmologies on the branes.

Increasing evidence suggests that most of the energy density of the universe consists of a dark energy component with negative pressure, a ``cosmological constant" that causes the cosmic expansion to accelerate. In this paper, we address the puzzle of why this component comes to dominate the universe only recently rather than at some much earlier epoch. We present a class of theories based on an evolving scalar field where the explanation is based entirely on internal dynamical properties of the solutions. In the theories we consider, the dynamics causes the scalar field to lock automatically into a negative pressure state at the onset of matter-domination such that the present epoch is the earliest possible time, consistent with nucleosynthesis restrictions, when it can start to dominate.

Cosmological models with cold dark matter composed of weakly interacting particles predict overly dense cores in the centers of galaxies and clusters and an overly large number of halos within the Local Group compared to actual observations. We propose that the conflict can be resolved if the cold dark matter particles are self-interacting with a large scattering cross-section but negligible annihilation or dissipation. In this scenario, astronomical observations may enable us to study dark matter properties that are inaccessible in the laboratory

In this paper, we show how the generic coupling of moduli to the kinetic energy of ordinary matter fields results in a cosmological mechanism that influences the evolution and stability of moduli. As an example, we reconsider the problem of stabilizing the dilaton in a non-perturbative potential induced by gaugino condensates. A well-known difficulty is that the potential is so steep that the dilaton field tends to overrun the correct minimum and to evolve to an observationally unacceptable vacuum. We show that the dilaton coupling to the thermal energy of matter fields produces a natural mechanism for gently relaxing the dilaton field into the correct minimum of the potential without fine-tuning of initial conditions. The same mechanism is potentially relevant for stabilizing other moduli fields.

We present a comprehensive study of the observational constraints on spatially flat cosmological models containing a mixture of matter and quintessence --- a time varying, spatially inhomogeneous component of the energy density of the universe with negative pressure. Our study also includes the limiting case of a cosmological constant. Low red shift constraints include the Hubble parameter, baryon fraction, cluster abundance, age of the universe, bulk velocity and shape of the mass power spectrum; intermediate red shift constraints are due to type 1a supernovae, gravitational lensing, the Ly-a forest, and the evolution of large scale structure; high red shift constraints are based on cosmic microwave background temperature anisotropy. Mindful of systematic errors, we adopt a conservative approach in applying these constraints. We determine that quintessence models in which the matter density parameter is $0.2 \ls \Omega_m \ls 0.5$ and the effective, density-averaged equation of state is $-1 \le w \ls -0.2$, are consistent with the most reliable, current low red shift and CMB observations at the $2\sigma$ level. Factoring in the constraint due to type 1a SNe, the range for the equation of state is reduced to $-1 \le w \ls -0.4$, where this range represents models consistent with each observational constraint at the 2$\sigma$ level or better (concordance analysis). A combined maximum likelihood analysis suggests a smaller range, $-1 \le w \ls -0.6$. We find that the best-fit and best-motivated quintessence models lie near $\Omega_m \approx 0.33$, $h \approx 0.65$, and spectral index $n_s=1$, with an effective equation of state $w \approx -0.65$ for ``tracker'' quintessence and $w=-1$ for ``creeper'' quintessence. (abstract shortened)

We propose a simple and accurate method for computing analytically the mass correlation function for cold dark matter and scale-free models that fits N-body simulations over a range that extends from the linear to the strongly non-linear regime. The method, based on the dynamical evolution of the pair conservation equation, relies on a universal relation between the pair-wise velocity and the smoothed correlation function valid for high and low density models, as derived empirically from N-body simulations. An intriguing alternative relation, based on the stable-clustering hypothesis, predicts a power-law behavior of the mass correlation function that disagrees with N-body simulations but conforms well to the observed galaxy correlation function if negligible bias is assumed. The method is a useful tool for rapidly exploring a wide span of models and, at the same time, raises new questions about large scale structure formation.

The "cosmic triangle" is introduced as a way of representing the past, present, and future status of the universe. Our current location within the cosmic triangle is determined by the answers to three questions: How much matter is in the universe? Is the expansion rate slowing down or speeding up? And, is the universe flat? A review of recent observations suggests a universe that is lightweight (matter density about one-third the critical value), is accelerating, and is flat. The acceleration implies the existence of cosmic dark energy that overcomes the gravitational self-attraction of matter and causes the expansion to speed up.

We argue that the brane-world picture with matter-fields confined to 4-d domain walls and with gravitational interactions across the bulk disallows adding an arbitrary constant to the low-energy, 4-d effective theory -- which finesses the usual cosmological constant problem. The analysis also points to difficulties in stabilizing moduli fields; as an alternative, we suggest scenarios in which the moduli motion is heavily damped by various cosmological mechanisms and varying ultra-slowly with time.

Recently, we introduced the notion of "tracker fields," a form of quintessence which has an attractor-like solution. Using this concept, we showed how to construct models in which the ratio of quintessence to matter densities today is independent of initial conditions. Here we apply the same idea to the standard cold dark matter component in cases where it is composed of oscillating fields. Combining these ideas, we can construct a model in which quintessence, cold dark matter, and ordinary matter all contribute comparable amounts to the total energy density today irrespective of initial conditions.

A substantial fraction of the energy density of the universe may consist of quintessence in the form of a slowly-rolling scalar field. Since the energy density of the scalar field generally decreases more slowly than the matter energy density, it appears that the ratio of the two densities must be set to a special, infinitesimal value in the early universe in order to have the two densities nearly coincide today. Recently, we introduced the notion of tracker fields to avoid this initial conditions problem. In the paper, we address the following questions: What is the general condition to have tracker fields? What is the relation between the matter energy density and the equation-of-state of the universe imposed by tracker solutions? And, can tracker solutions explain why quintessence is becoming important today rather than during the early universe?

Recent observations suggest that a large fraction of the energy density of the universe has negative pressure. One explanation is vacuum energy density; another is quintessence in the form of a scalar field slowly evolving down a potential. In either case, a key problem is to explain why the energy density nearly coincides with the matter density today. The densities decrease at different rates as the universe expands, so coincidence today appears to require that their ratio be set to a specific, infinitessimal value in the early universe. In this paper, we introduce the notion of a "tracker field", a form of quintessence, and show how it may explain the coincidence, adding new motivation for the quintessence scenario.

Some form of missing energy may account for the difference between the observed cosmic matter density and the critical density. Two leading candidates are a cosmological constant and quintessence (a time-varying, inhomogenous component with negative pressure). We show that an ideal, full-sky cosmic background anisotropy experiment may not be able to distinguish the two, even when non-linear effects due to gravitational lensing are included. Due to this ambiguity, microwave background experiments alone may not determine the matter density or Hubble constant very precisely. We further show that degeneracy may remain even after considering classical cosmological tests and measurements of large scale structure.

The abundance of rich clusters is a strong constraint on the mass power spectrum. The current constraint can be expressed in the form $\sigma_8 \Omega_m^{\gamma} = 0.5 \pm 0.1$ where $\sigma_8$ is the $rms$ mass fluctuation on 8 $h^{-1}$ Mpc scales, $\Omega_m$ is the ratio of matter density to the critical density, and $\gamma$ is model-dependent. In this paper, we determine a general expression for $\gamma$ that applies to any models with a mixture of cold dark matter plus cosmological constant or quintessence (a time-evolving, spatially-inhomogeneous component with negative pressure) including dependence on the spectral index $n$, the Hubble constant $h$, and the equation-of-state of the quintessence component $w$. The cluster constraint is combined with COBE measurements to identify a spectrum of best-fitting models. The constraint from the evolution of rich clusters is also discussed.

We examine the possibility that a significant component of the energy density of the universe has an equation-of-state different from that of matter, radiation or cosmological constant ($\Lambda$). An example is a cosmic scalar field evolving in a potential, but our treatment is more general. Including this component alters cosmic evolution in a way that fits current observations well. Unlike $\Lambda$, it evolves dynamically and develops fluctuations, leaving a distinctive imprint on the microwave background anisotropy and mass power spectrum.

An alternative to the standard cold dark matter model has been recently proposed in which a significant fraction of the energy density of the universe is due to a dynamical scalar field ($Q$) whose effective equation-of-state differs from that of matter, radiation or cosmological constant ($\Lambda$). In this paper, we determine how the Q-component modifies the primordial inflation gravitational wave (tensor metric) contribution to the cosmic microwave background anisotropy and, thereby, one of the key tests of inflation.

We derive a closed-form, analytical expression for the spectrum of long-wavelength density perturbations in inflationary models with two (or more) inflaton degrees of freedom that is valid in the slow-roll approximation. We illustrate several classes of potentials for which this expression reduces to a simple, algebraic expression.

We examine the theoretical foundations of standard methods for computing density perturbations in inflationary models. We find that: (1) the time-delay formalism (introduced by Guth and Pi, 1982) is only valid when inflation is well-described by the de Sitter solution and the equation-of-state is nearly unchanging; and, (2) the horizon-crossing/Bessel approximation extends to non-exponential inflation, but only if the equation-of-state is changing slowly. Integration of the gauge-invariant perturbation equations mode-by-mode is the only method reliable for general models. For models with rapidly varying equation-of-state, the correction leads to significantly different predictions for the microwave background anisotropy. An important corollary is that methods proposed for "reconstruction" of the inflaton potential from anisotropy data are unreliable for general models.

Parametric resonance has been discussed as a mechanism for copious particle production following inflation. Here we present a simple and intuitive calculational method for estimating the efficiency of parametric amplification as a function of parameters. This is important for determining whether resonant amplification plays an important role in the reheating process. We find that significant amplification occurs only for a limited range of couplings and interactions.

A proposed strategy for determining the deceleration parameter entails measuring the deviation from a linear (Hubble) distance-red shift relation. However, even at moderate red shifts, z > 0.2, the deviation does not depend on q_0 alone, but also on the equation-of-state of the universe.

Observational constraints guide one forcefully to examine models in which the matter density is substantially less than critical density. Particularly noteworthy are those which are consistent with inflation. For these models, microwave background anisotropy, large-scale structure measurements, direct measurements of the Hubble constant, $H_0$, and the closure parameter, $\Omega_{\rm Matter}$, ages of stars and a host of more minor facts are all consistent with a spatially flat model having significant cosmological constant $\Omega_{\Lambda} = 0.65 \pm 0.1$, $\Omega_{\rm Matter} = 1 - \Omega_{\Lambda}$ (in the form of ``cold dark matter") and a small tilt: $0.8< n < 1.2$. (Read comments at top of file to find out how to obtain figures.)

An exploratory study of the cosmology of moduli in string theory. Moduli are argued to be natural inflaton fields and lead to a robust inflationary cosmology in which inflation takes place at the top of domain walls. The amplitude of microwave background fluctuations constrains the dynamics responsible for inflation to take place at a higher scale than supersymmetry breaking. Models explaining this difference in scales and also preventing the dilaton from running to infinity are proposed. The problem of dilaton domination of the energy density of the universe is not resolved.

Observational tests during the next decade may determine if the evolution of the Universe can be understood from fundamental physical principles, or if special initial conditions, coincidences, and new, untestable physical laws must be invoked. The inflationary model of the Universe is an important example of a predictive cosmological theory based on physical principles. In this talk, we discuss the distinctive fingerprint that inflation leaves on the cosmic microwave background anisotropy. We then suggest a series of five milestone experimental tests of the microwave background which could determine the validity of the inflationary hypothesis within the next decade. The paper is a Review based on a Plenary talk given at the Snowmass Workshop on Particle Astrophysics and Cosmology, 1995 It will appear in the Proceedings edited by E. Kolb and R.Peccei. Software package for computing filter functions and band power estimates available thru world-wide-web at http://dept.physics.upenn.edu/~www/astro-cosmo/ .

A survey of solutions to the cosmological moduli problem in string theory. The only extant proposal which may work is Intermediate Scale Inflation as proposed by Randall and Thomas. Supersymmetry preserving dynamics which could give large masses to the moduli is strongly constrained by cosmology and requires the existence of string vacuum states possessing properties different from those of any known vacuuum. Such a mechanism cannot give mass to the dilaton unless there are cancellations between different exponentially small contributions to the superpotential. Our investigation also shows that stationary points of the effective potential with negative vacuum energy do not correspond to stationary solutions of the equations of postinflationary cosmology. This suggests that supersymmetry breaking is a requirement for a successful inflationary cosmology.

We examine the possibility that an epoch of inflationary expansion induces high-frequency oscillations of Newton's constant, $G$. The effect occurs because inflation can shift the expectation value of a non-minimally coupled, Brans-Dicke-like field away from the minimum of its effective potential. At some time after inflation ends, the field begins to oscillate, resulting in periodic variations in $G$. We find conditions for which the oscillation energy would be sufficient to close the universe, consistent with all known constraints from cosmology and local tests of general relativity.

The CMB anisotropy depends sensitively upon the slope and amplitude of primordial density and gravitational wave fluctuations, the baryon density, the Hubble constant, the cosmological constant, the ionization history, {\it etc.} We report on recent work showing how well multi-scale measurements of the anisotropy power spectrum can resolve these factors. We identify a hypersurface in cosmic parameter space that can be accurately localized by observations, but along which the likelihood will hardly vary. Other cosmic observations will be needed to break this degeneracy.

The cosmic microwave background anisotropy is sensitive to the slope and amplitude of primordial energy density and gravitational wave fluctuations, the baryon density, the Hubble constant, the cosmological constant, the ionization history, {\it etc.} In this Letter, we examine the degree to which these factors can be separately resolved from combined small- and large-angular scale anisotropy observations. We isolate directions of degeneracy in this cosmic parameter space, but note that other cosmic observations can break the degeneracy. (Correction to Eq. 4)

The contribution of gravitational wave (tensor metric) and energy density (scalar metric) fluctuations to the cosmic microwave background polarization is computed by numerically solving the relativistic radiation transfer equations. We find that the tensor contribution is significant only at large angular scales (multipoles $\ell \lta 40$). For standard recombination, the tensor contribution can dominate at $\ell \lta 40$; however, the effect would be difficult to detect since the total (scalar plus tensor) polarization is $< 1$\%. For models with late reionization, the total large angular scale polarization is large ($\sim 7-9$\%), but the tensor fraction is negligibly small. Hence, polarization may be useful for discriminating ionization history, but is much less promising as a means for detecting tensor fluctuations.

Long-wavelength gravitational waves can induce significant temperature anisotropy in the cosmic microwave background. Distinguishing this from anisotropy induced by energy density fluctuations is critical for testing inflationary cosmology and theories of large-scale structure formation. We describe full radiative transport calculations of the two contributions and show that they differ dramatically at angular scales below a few degrees. We show how anisotropy experiments probing large- and small-angular scales can combine to distinguish the imprint due to gravitational waves.

Extended inflation is a promising new approach to implementing the inflationary universe scenario. This paper reviews recent advances including a new, more robust mechanism for ending extended inflation, a new prediction for the density fluctuation spectrum generated by extended inflation, and the discovery that extended inflation can produce gravitational waves that can significantly add to the cosmic background anisotropy. Invited Talk at the Journees Relativistes, Amsterdam, May 14-16, 1992

Could COBE DMR be detecting the imprint from a spectrum of gravitational waves generated during inflation? The conventional inflationary prediction had been that the cosmic microwave anisotropy is dominated by energy density fluctuations generated during inflation and that the gravitational waves contribute negligibly. In this paper, we report on recent work (in collaboration with R. Davis, H. Hodges, and M. Turner) that has shown that the conventional wisdom may be wrong; specifically, gravitational waves may dominate the anisotropy in inflationary models where the spectrum of perturbations deviates significantly from scale invariance (e.g., extended and power-law inflation models and extreme versions of chaotic inflation). If gravitational waves do dominate at the large-angular scales measured by COBE DMR, the expectation and interpretation of anisotropies on small-angular scales is profoundly altered. Invited Paper for Proceedings of the Journees Relativistes, Amsterdam, May 14-16, 1992

We consider whether current notions about superstring theory below the Planck scale are compatible with cosmology. We find that the anticipated form for the dilaton interaction creates a serious roadblock for inflation and makes it unlikely that the universe ever reaches a state with zero cosmological constant and time-independent gravitational constant.

Inflation creates both scalar (density) and tensor (gravity wave) metric perturbations. We find that the tensor mode contribution to the CMB anisotropy on large-angular scales can only exceed that of the scalar mode in models where the spectrum of perturbations deviates significantly from scale invariance (e.g., extended and power-law inflation models and extreme versions of chaotic inflation). If the tensor mode dominates at large-angular scales, then the value of $\Delta T/T$ predicted on $1^\circ$ is less than if the scalar mode dominates, and, for cold dark matter models, $b>1$ can be made consistent with the COBE DMR results.

Two qualitatively different modes of ending superluminal expansion are possible in extended inflation. One mode, different from the one envoked in most extended models to date, easily avoids making big bubbles that distort the cosmic microwave background radiation (CMBR). In this mode, the spectrum of density fluctuations is found to be scale-free, $P(k) \propto k^n$, where $n$ might lie anywhere between 0.5 and 1.0 (whereas, previously, it appeared that the range $1.0> n \gtsim 0.84$ was disallowed).