We propose that the state of the universe does {\it not} spontaneously violate CPT. Instead, the universe after the big bang is the CPT image of the universe before it, both classically and quantum mechanically. The pre- and post-bang epochs comprise a universe/anti-universe pair, emerging from nothing directly into a hot, radiation-dominated era. CPT symmetry selects a unique QFT vacuum state on such a spacetime, providing a new interpretation of the cosmological baryon asymmetry, as well as a remarkably economical explanation for the cosmological dark matter. Requiring only the standard three-generation model of particle physics (with right-handed neutrinos), a $\mathbb{Z}_2$ symmetry suffices to render one of the right-handed neutrinos stable. We calculate its abundance from first principles: matching the observed dark matter density requires its mass to be $4.8\times10^{8}~{\rm GeV}$. Several other testable predictions follow: (i) the three light neutrinos are Majorana and allow neutrinoless double $\beta$ decay; (ii) the lightest neutrino is massless; and (iii) there are no primordial long-wavelength gravitational waves. We mention connections to the strong CP problem and the arrow of time.

We investigate the idea that the universe before the Big Bang is the CPT reflection of the universe after the bang, so that the state of the universe does {\it not} spontaneously violate CPT. The universe before the bang and the universe after the bang may be viewed as a universe/anti-universe pair, created from nothing. The early universe is radiation dominated and inflationary energy is not required. We show how CPT selects a preferred vacuum state for quantum fields on such a cosmological spacetime. This, in turn, leads to a new view of the cosmological matter/anti-matter asymmetry, and a novel and economical explanation of the dark matter abundance. If we assume that the matter fields in the universe are described by the standard model of particle physics (including right-handed neutrinos), it is natural for one of the heavy neutrinos to be stable, and we show that in order to match the observed dark matter density, its mass must be $4.8\times10^{8}~{\rm GeV}$. We also obtain further predictions, including: (i) that the three light neutrinos are majorana; (ii) that the lightest of these is exactly massless; and (iii) that there are no primordial, long-wavelength gravitational waves.

If observations confirm BICEP2's claim of a tensor-scalar ratio $r\approx 0.2$ on CMB scales, then the inflationary consistency relation $n_{t}=-r/8$ predicts a small negative value for the tensor spectral index $n_t$. We show that future CMB polarization experiments should be able to confirm this prediction at several sigma. We also show how to properly extend the consistency relation to solar system scales, where the primordial gravitational wave density $\Omega_{gw}$ could be measured by proposed experiments such as the Big Bang Observer. This would provide a far more stringent test of the consistency relation and access much more detailed information about the early universe.

The recent BICEP2 measurement of primordial gravity waves (r = 0.2^{+0.07}_{-0.05}) appears to be in tension with the upper limit from WMAP (r<0.13 at 95% CL) and Planck (r<0.11 at 95% CL). We carefully quantify the level of tension and show that it is very significant (around 0.1% unlikely) when the observed deficit of large-scale temperature power is taken into account. We show that measurements of TE and EE power spectra in the near future will discriminate between the hypotheses that this tension is either a statistical fluke, or a sign of new physics. We also discuss extensions of the standard cosmological model that relieve the tension, and some novel ways to constrain them.

Consider the minimal renormalizable extension of the Standard Model with purely dimensionless couplings, successful electroweak symmetry breaking (via the Coleman-Weinberg mechanism) and a see-saw mechanism for neutrino mass: we will call this the Minimal Dimensionless Standard Model (MDSM). In fact, 3 closely related models fit the bill: MDSM_1, MDSM_2 and MDSM_3. We analyze the theoretical and observational constraints on these models. We argue that, when they are minimally coupled to gravity, they can accomplish several important cosmological tasks (inflation, dark matter, leptogenesis) in a way that is economical, predictive and tightly woven into the fabric of known physics. One of the models (MDSM_3), which includes an extra U(1)_{B-L} gauge symmetry, seems particularly promising.

Pulsar timing arrays (PTAs) will be sensitive to a finite number of gravitational wave (GW) "point" sources (e.g. supermassive black hole binaries). N quiet pulsars with accurately known distances d_{pulsar} can characterize up to 2N/7 distant chirping sources per frequency bin \Delta f_{gw}=1/T, and localize them with "diffraction limited" precision \delta\theta \gtrsim (1/SNR)(\lambda_{gw}/d_{pulsar}). Even if the pulsar distances are poorly known, a PTA with F frequency bins can still characterize up to (2N/7)[1-(1/2F)] sources per bin, and the quasi-singular pattern of timing residuals in the vicinity of a GW source still allows the source to be localized quasi-topologically within roughly the smallest quadrilateral of quiet pulsars that encircles it on the sky, down to a limiting resolution \delta\theta \gtrsim (1/SNR) \sqrt{\lambda_{gw}/d_{pulsar}}. PTAs may be unconfused, even at the lowest frequencies, with matched filtering always appropriate.

Radiation emitted near a black hole reaches the observer by multiple paths; and when this radiation varies in time, the time-delays between the various paths generate a "blinking" effect in the observed light curve L(t) or its auto-correlation function xi(T)= <L(t)L(t-T)>. For the particularly important "face-on" configuration (in which the hole is viewed roughly along its spin axis, while the emission comes roughly from its equatorial plane -- e.g. from the inner edge of its accretion disk, or from the violent flash of a nearby/infalling star) we calculate the blinking in detail by computing the time delay Delta t_{j}(r,a) and magnification mu_{j}(r,a) of the jth path (j=1,2,3,...), relative to the primary path (j=0), as a function of the emission radius r and black hole spin 0<a/M<1. The particular geometry and symmetry of the nearly-face-on configuration enhances and "protects" the blinking signal, making it more detectable and more independent of certain astrophysical and observational details. The effect can be surprisingly strong: e.g. for radiation from the innermost stable circular orbit ("ISCO") of a black hole of critical spin (a_{crit}/M = 0.853), the j=1,2,3 fluxes are, respectively, 27%, 2% and 0.1% of the j=0 flux.

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.

Porcupines are networks of gravitational wave detectors in which the detectors and the distances between them are short relative to the gravitational wavelengths of interest. Perfect porcupines are special configurations whose sensitivity to a gravitational plane wave is independent of the propagation direction or polarization of the wave. I develop the theory of porcupines, including the optimal estimator \hat{h}^{ij} for the gravitational wave field; useful formulae for the spin-averaged and rotationally-averaged SNR^{2}; and a simple derivation of the properties of perfect porcupines. I apply these results to the interesting class of ``simple'' porcupines, and mention some open problems.

Perfect porcupines are specially-configured networks of gravitational wave detectors, in the limit that the individual detectors and the distances between them are short relative to the gravitational wavelengths of interest. They have beautiful properties which make them ideal gravitational wave telescopes. I present the most important cases explicitly. For a network of one-arm detectors (like "AGIS" \cite{Dimopoulos_et_al}), the minimal perfect porcupine has 6 detectors, oriented along the 6 diameters of a regular icosahedron. For a network of two-arm detectors (like the equal-arm Michelson interferometers LIGO/VIRGO) the minimal perfect porcupine is a certain 5 detector configuration.

In a recent paper arXiv:0709.0299, we introduced a spin expansion that provides a simple yet powerful way to understand aspects of binary black hole (BBH) merger. This approach relies on the symmetry properties of initial and final quantities like the black hole mass m, kick velocity {\bf k}, and spin vector {\bf s}, rather than a detailed understanding of the merger dynamics. In this paper, we expand on this proposal, examine how well its predictions agree with current simulations, and discuss several future directions that would make it an even more valuable tool. The spin expansion yields many new predictions, including several exact results that may be useful for testing numerical codes. Some of these predictions have already been confirmed, while others await future simulations. We explain how a relatively small number of simulations -- 10 equal-mass simulations, and 16 unequal-mass simulations -- may be used to calibrate all of the coefficients in the spin expansion up to second order at the minimum computational cost. For a more general set of simulations of given covariance, we derive the minimum-variance unbiased estimators for the spin expansion coefficients. We discuss how this calibration would be interesting and fruitful for general relativity and astrophysics. Finally, we sketch the extension to eccentric orbits.

We regard binary black hole (BBH) merger as a map from a simple initial state (two Kerr black holes, with dimensionless spins {\bf a} and {\bf b}) to a simple final state (a Kerr black hole with mass m, dimensionless spin {\bf s}, and kick velocity {\bf k}). By expanding this map around {\bf a} = {\bf b} = 0 and applying symmetry constraints, we obtain a simple formalism that is remarkably successful at explaining existing BBH simulations. It also makes detailed predictions and suggests a more efficient way of mapping the parameter space of binary black hole merger. Since we rely on symmetry rather than dynamics, our expansion complements previous analytical techniques.

We derive a general master equation relating the gravitational-wave observables r and Omega_gw(f). Here r is the tensor-to-scalar ratio, constrained by cosmic-microwave-background (CMB) experiments; and Omega_gw(f) is the energy spectrum of primordial gravitational-waves, constrained e.g. by pulsar-timing measurements, laser-interferometer experiments, and Big Bang Nucleosynthesis (BBN). Differentiating the master equation yields a new expression for the tilt d(ln Omega_gw(f))/d(ln f). The relationship between r and Omega_gw(f) depends sensitively on the uncertain physics of the early universe, and we show that this uncertainty may be encapsulated (in a model-independent way) by two quantities: w_hat(f) and nt_hat(f), where nt_hat(f) is a certain logarithmic average over nt(k) (the primordial tensor spectral index); and w_hat(f) is a certain logarithmic average over w_tilde(a) (the effective equation-of-state in the early universe, after horizon re-entry). Here the effective equation-of-state parameter w_tilde(a) is a combination of the ordinary equation-of-state parameter w(a) and the bulk viscosity zeta(a). Thus, by comparing constraints on r and Omega_gw(f), one can obtain (remarkably tight) constraints in the [w_hat(f), nt_hat(f)] plane. In particular, this is the best way to constrain (or detect) the presence of a ``stiff'' energy component (with w > 1/3) in the early universe, prior to BBN. Finally, although most of our analysis does not assume inflation, we point out that if CMB experiments detect a non-zero value for r, then we will immediately obtain (as a free by-product) a new upper bound w_hat < 0.55 on the logarithmically averaged effective equation-of-state parameter during the ``primordial dark age'' between the end of inflation and the start of BBN.

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.

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)$.

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.

We investigate a class of models for dark matter and/or negative-pressure, dynamical dark energy consisting of ``spintessence,'' a complex scalar field $\phi$ spinning in a U(1)-symmetric potential $V(\phi)=V(|\phi|)$. As the Universe expands, the field spirals slowly toward the origin. The internal angular momentum plays an important role in the cosmic evolution and fluctuation dynamics. We outline the constraints on a cosmic spintessence field, describing the properties of the potential necessary to sustain a viable dark energy model, making connections with quintessence and self-interacting and fuzzy cold dark matter. Possible implications for the coincidence problem, baryogenesis, and cosmological birefringence, and generalizations of spintessence to models with higher global symmetry and models in which the symmetry is not exact are also discussed.