Many proposals for physics beyond the Standard Model give rise to a dark sector containing many degrees of freedom. In this work, we explore the cosmological implications of the non-trivial dynamics which may arise within such dark sectors, focusing on decay processes which take place entirely among the dark constituents. First, we demonstrate that such decays can leave dramatic imprints on the resulting dark-matter phase-space distribution. In particular, this distribution need not be thermal -- it can even be multi-modal, exhibiting a non-trivial pattern of peaks and troughs as a function of momentum. We then proceed to show how these features can induce modifications to the matter power spectrum. Finally, we assess the extent to which one can approach the archaeological "inverse" problem of deciphering the properties of an underlying dark sector from the matter power spectrum. Indeed, one of the main results of this paper is a remarkably simple conjectured analytic expression which permits the reconstruction of many of the important features of the dark-matter phase-space distribution directly from the matter power spectrum. Our results therefore provide an interesting toolbox of methods for learning about, and potentially constraining, the features of non-minimal dark sectors and their dynamics in the early universe.

Dynamical Dark Matter (DDM) is an alternative framework for dark-matter physics in which the dark sector comprises a vast ensemble of particle species whose decay widths are balanced against their cosmological abundances. Previous studies of this framework have focused on a particular class of DDM ensembles --- motivated primarily by KK towers in theories with extra dimensions --- in which the density of states scales roughly as a polynomial of mass. In this paper, by contrast, we study the properties of a different class of DDM ensembles in which the density of states grows exponentially with mass. Ensembles with this Hagedorn-like property arise naturally as the "hadrons" associated with the confining phase of a strongly-coupled dark sector; they also arise naturally as the gauge-neutral bulk states of Type I string theories. We study the dynamical properties of such ensembles, and demonstrate that an appropriate DDM-like balancing between decay widths and abundances can emerge naturally --- even with an exponentially rising density of states. We also study the effective equations of state for such ensembles, and investigate some of the model-independent observational constraints on such ensembles that follow directly from these equations of state. In general, we find that such constraints tend to introduce correlations between various properties of these DDM ensembles such as their associated mass scales, lifetimes, and abundance distributions. For example, we find that these constraints allow DDM ensembles with energy scales ranging from the GeV scale all the way to the Planck scale, but the total present-day cosmological abundance of the dark sector must be spread across an increasing number of different states in the ensemble as these energy scales are dialed from the Planck scale down to the GeV scale. Numerous other correlations and constraints are also discussed.

In this paper, we apply a new statistical analysis technique, Mean Field approach to Bayesian Independent Component Analysis (MF-ICA), on galaxy spectral analysis. This algorithm can compress the stellar spectral library into a few Independent Components (ICs), and galaxy spectrum can be reconstructed by these ICs. Comparing to other algorithms which decompose a galaxy spectrum into a combination of several simple stellar populations, MF-ICA approach offers a large improvement in the efficiency. To check the reliability of this spectral analysis method, three different methods are used: (1) parameter-recover for simulated galaxies, (2) comparison with parameters estimated by other methods, and (3) consistency test of parameters from the Sloan Digital Sky Survey galaxies. We find that our MF-ICA method not only can fit the observed galaxy spectra efficiently, but also can recover the physical parameters of galaxies accurately. We also apply our spectral analysis method to the DEEP2 spectroscopic data, and find it can provide excellent fitting for those low signal-to-noise spectra.

Inflation predicts B-mode polarization with correlations that span superhorizon scales at recombination. In contrast, the correlations set up by causal sources, such as phase transitions or defects, necessarily vanish on superhorizon scales. Motivated by BICEP2's B-mode detection, we consider the prospects for measuring the inflationary superhorizon signature in future observations. We explain that the finite resolution of an experiment and the filtering of the raw data induces a transfer of spurious subhorizon power to superhorizon scales, and describe ways to correct for it. We also provide a detailed treatment of possible sources of noise in the measurement. Finally, we present forecasts for the detectability of the signal with future CMB polarization experiments.

We study the neutralino being the Lightest Supersymmetric Particle (LSP) as a cold Dark Matter (DM) candidate with a mass less than 40 GeV in the framework of the Next-to-Minimal-Supersymmetric-Standard-Model (NMSSM). We find that with the current collider constraints from LEP, the Tevatron and the LHC, there are three types of light DM solutions consistent with the direct/indirect searches as well as the relic abundance considerations: (i) A1, H1-funnels, (ii) stau coannihilation and (iii) sbottom coannihilation. Type-(i) may take place in any theory with a light scalar (or pseudo-scalar) near the LSP pair threshold; while Type-(ii) and (iii) could occur in the framework of Minimal-Supersymmetric-Standard-Model (MSSM) as well. We present a comprehensive study on the properties of these solutions and point out their immediate relevance to the experiments of the underground direct detection such as superCDMS and LUX/LZ, and the astro-physical indirect search such as Fermi-LAT. We also find that the decays of the SM-like Higgs boson may be modified appreciably and the new decay channels to the light SUSY particles may be sizable. The new light CP-even and CP-odd Higgs bosons will decay to a pair of LSPs as well as other observable final states, leading to interesting new Higgs phenomenology at colliders. For the light sfermion searches, the signals would be very challenging to observe at the LHC given the current bounds. However, a high energy and high luminosity lepton collider, such as the ILC, would be able to fully cover these scenarios by searching for events with large missing energy plus charged tracks or displaced vertices.

The Planck mission has conclusively detected lensing of the Cosmic Microwave Background (CMB) radiation from foreground sources to an overall significance of greater than $25\sigma$. The high precision of this measurement motivates the development of a more complete formulation of the calculation of this effect. While most effects on the CMB anisotropies are widely studied through direct solutions of the Boltzmann equation, the non-linear effect of CMB lensing is formulated through the solutions of the geodesic equation. In this paper, we present a new formalism to the calculation of the lensing effect by \emph{directly solving the Boltzmann equation}, as we did in the calculation of the CMB anisotropies at recombination. In particular, we developed a diagrammatic approach to efficiently keep track of all the interaction terms and calculate all possible non-trivial correlations to arbitrary high orders. Using this formalism, we explicitly articulate the approximations required to recover the usual remapping approach used in current studies of the weak lensing. In addition, we point out additional unexplored corrections that are manifest in our formalism to which experiments may be sensitive. As an example, we calculate the correction to the CMB temperature power spectrum for the \emph{lens-lens} coupling effects which are neglected in standard calculations. We find that the correction is $\lesssim 0.1\%$ of the CMB temperature power spectrum for $\ell$ up to 3000 and thus is comparable to the cosmic variance.

We study cosmological perturbations by solving the governing Boltzmann and Einstein Field equations up to second order, and calculate the corresponding CMB bispectrum during recombination. We include all the second-order Liouville and collision terms, truncating the multipole hierarchy at $l=10$, consistently including all $m\neq 0$ terms when calculating the bispectrum in the flat-sky limit. At this stage, we focus on contributions at recombination and we neglect 2nd-order vector and tensor terms, lensing effects, and late time non-linear ISW. We find that the signal-to-noise for the bispectrum is 0.69 for $l_{max}=2000$, yielding an overall signal $F_{\rm NL} = 3.19$ (normalised relative to the local model). We find that the effective $f_{NL}$'s of the equilateral and local type are 5.11 and 0.88 respectively. This recombination signal will have to be taken into account in a quantitative analysis of the Planck data.

We determine the prospects for finding dark matter at the Tevatron and LHC through the production of exotic 4th generation quarks T' that decay through T' \to t X, where X is dark matter. The resulting signal of t \bar{t} + \met has not previously been considered in searches for 4th generation quarks, but there are both general and specific dark matter motivations for this signal, and with slight modifications, this analysis applies to any scenario where invisible particles are produced in association with top quarks. Current direct and indirect bounds on such exotic quarks restrict their masses to be between 300 and 600 GeV, and the dark matter's mass may be anywhere below m_T'. We simulate the signal and main backgrounds with MadGraph/MadEvent-Pythia-PGS4. For the Tevatron, we find that an integrated luminosity of 20 fb^-1 will allow 3\sigma discovery up to m_T' = 400 GeV and 95% exclusion up to m_T' = 455 GeV. For the 10 TeV LHC with 300 pb^-1, the discovery and exclusion sensitivities rise to 490 GeV and 600 GeV. These scenarios are therefore among the most promising for dark matter at colliders. Perhaps most interestingly, we find that dark matter models that can explain results from the DAMA, CDMS and CoGeNT Collaborations can be tested with high statistical significance using data already collected at the Tevatron and have extraordinarily promising implications for early runs of the LHC.

Models of a rotating universe have been studied widely since G{\"o}del \cite{1}, who showed an example that is consistent with General Relativity (GR). By now, the possibility of a rotating universe has been discussed comprehensively in the framework of some types of Bianchi's models, such as Type V, VII and IX \cite{2,3}, and different approaches have been proposed to constrain the rotation. Recent discoveries of some non-Gaussian properties of the Cosmic Microwave Background Anisotropies (CMBA) \cite{nG1,nG2,nG3,nG4,nG5,nG6,nG7}, such as the suppression of the quadrupole and the alignment of some multipoles draw attention to some Bianchi models with rotation \cite{bi1,bi2}. However, cosmological data, such as those of the CMBA, strongly prefer a homogeneous and isotropic model. Therefore, it is of interest to discuss the rotation of the universe as a perturbation of the Robertson-Walker metric, to constrain the rotating speed by cosmological data and to discuss whether it could be the origin of the non-Gaussian properties of the CMBA mentioned above. Here, we derive the general form of the metric (up to 2nd-order perturbations) which is compatible with the rotation perturbation in a flat $\Lambda$-CDM universe. By comparing the 2nd-order Sachs-Wolfe effect \cite{4,5,6,7,8} due to rotation with the CMBA data, we constrain the angular speed of the rotation to be less than $10^{-9}$ rad yr$^{-1}$ at the last scattering surface. This provides the first constraint on the shear-free rotation of a $\Lambda$CDM universe.

In the standard framework of cosmology, primordial density fluctuations are assumed to have an isotropic Gaussian distribution. We search for deviations from this assumption in the WMAP data for the low $l$ modes of Cosmic Microwave Background Anisotropies (CMBA), by studying the directions of the z-axis that maximize the $l=m$ modes and the resulting amplitudes of these modes. We find a general alignment of the directions for $l=2$ to 10 modes to within 1/4 of the northern hemisphere. This alignment can be regarded as a generalization of the recently discovered alignment of the $l=2$ and 3 modes - the so-called `Axis of Evil'. Furthermore, we find abnormally high (low) powers in the $l=m=6$, 12 - 17 ($l=m=5$) modes; the probabilities for having the anomalous amplitudes of the $l=m=5$, 6, 17 modes are about 0.1%, 1% and 1% respectively according to the Gaussian conjecture. The alignment and anomalous amplitudes for these low $l$ modes are very robust against foreground contamination or different cleaning strategies, suggesting a cosmological origin and possibly supporting a non-standard inflation.

The potential for neutrino telescopes to discover charged stau production in neutrino-nucleon interactions in Earth depends in part on the stau lifetime and range. In some supersymmetric scenarios, the next lightest supersymmetric particle is a stau with a decay length on the scale of 10 km. We evaluate the electromagnetic energy loss as a function of energy and stau mass. The energy loss parameter $\beta$ scales as the inverse stau mass for the dominating electromagnetic processes, photonuclear and $e^+e^-$ pair production. The range can be parameterized as a function of stau mass, initial energy and minimum final energy. In comparison to earlier estimates of the stau range, our results are as much as a factor of two larger, improving the potential for stau discovery in neutrino telescopes.

We investigate supergravity models in which the lightest supersymmetric particle (LSP) is a stable gravitino. We assume that the next-lightest supersymmetric particle (NLSP) freezes out with its thermal relic density before decaying to the gravitino at time t ~ 10^4 s - 10^8 s. In contrast to studies that assume a fixed gravitino relic density, the thermal relic density assumption implies upper, not lower, bounds on superpartner masses, with important implications for particle colliders. We consider slepton, sneutrino, and neutralino NLSPs, and determine what superpartner masses are viable in all of these cases, applying CMB and electromagnetic and hadronic BBN constraints to the leading two- and three-body NLSP decays. Hadronic constraints have been neglected previously, but we find that they provide the most stringent constraints in much of the natural parameter space. We then discuss the collider phenomenology of supergravity with a gravitino LSP. We find that colliders may provide important insights to clarify BBN and the thermal history of the Universe below temperatures around 10 GeV and may even provide precise measurements of the gravitino's mass and couplings.

Dark matter may be composed of superWIMPs, superweakly-interacting massive particles produced in the late decays of other particles. We focus here on the well-motivated supersymmetric example of gravitino LSPs. Gravitino superWIMPs share several virtues with the well-known case of neutralino dark matter: they are present in the same supersymmetric frameworks (supergravity with R-parity conservation) and naturally have the desired relic density. In contrast to neutralinos, however, gravitino superWIMPs are impossible to detect by conventional dark matter searches, may explain an existing discrepancy in Big Bang nucleosynthesis, predict observable distortions in the cosmic microwave background, and imply spectacular signals at future particle colliders.

Dark matter may be composed of superWIMPs, superweakly-interacting massive particles produced in the late decays of other particles. We focus on the case of gravitinos produced in the late decays of sleptons or sneutrinos and assume they are produced in sufficient numbers to constitute all of non-baryonic dark matter. At leading order, these late decays are two-body and the accompanying energy is electromagnetic. For natural weak-scale parameters, these decays have been shown to satisfy bounds from Big Bang nucleosynthesis and the cosmic microwave background. However, sleptons and sneutrinos may also decay to three-body final states, producing hadronic energy, which is subject to even more stringent nucleosynthesis bounds. We determine the three-body branching fractions and the resulting hadronic energy release. We find that superWIMP gravitino dark matter is viable and determine the gravitino and slepton/sneutrino masses preferred by this solution to the dark matter problem. In passing, we note that hadronic constraints disfavor the possibility of superWIMPs produced by neutralino decays unless the neutralino is photino-like.